function quesion_27
clear;clc;
%2017-2-23
data=[13.449 15.750 20.68 24.83;
      258 270 285 300];
x0=[10 20 1 1;
   300 200 1 10];
for i=1:2
    res(i,:) = fsolve(@fun,x0(i,:))
end
    function y=fun(x)
        % u1=x(1) u2=x(2) p1=x(3)  ss=x(4)
        for j=1:4
            y(j)=integral(@(t) x(3)*normpdf(t,x(1),x(4))+(1-x(3))*normpdf(t,x(2),x(4)),-inf,data(i,j))-j/5;
        end
    end
end

3楼: Originally posted by pdl9527 at 2017-02-23 21:55:43
唯一的问题是对初值的要求很高你可以尝试改变一下x0的值，可以发现如果初值和结果偏差很大将不能得到正确结果。


function quesion_27
clear;clc;
%2017-2-23
data=;
x0=;
for i=1:2
    res(i, = fsol ...

4楼: Originally posted by xiong1021 at 2017-02-23 22:40:31
谢谢解答，可能结果无法精确吧，我用fsolve试了也是初始值不一样最后结果偏差较大。
对了，请问下里面的p要求得是大于0小于1的约束怎么加进去啊
...

function quesion_27
clear;clc;
%2017-2-23
data=[13.449 15.750 20.68 24.83;
      258 270 285 300];
x0=[10 20 1 1;
   300 200 0.5 10];
for i=1:2
    res(i,:) = lsqnonlin(@fun,x0(i,:),[-inf -inf 0 -inf],[+inf +inf 1 +inf])
end
    function y=fun(x)
        % u1=x(1) u2=x(2) p1=x(3)  ss=x(4)
        for j=1:4
            y(j)=integral(@(t) x(3)*normpdf(t,x(1),x(4))+(1-x(3))*normpdf(t,x(2),x(4)),-inf,data(i,j))-j/5;
        end
    end
end

5楼: Originally posted by pdl9527 at 2017-02-23 23:16:08
要添加边界条件可以用最小二乘法。

function quesion_27
clear;clc;
%2017-2-23
data=;
x0=;
for i=1:2
    res(i, = lsqnonlin(@fun,x0(i,,,)
end
    function y=fun(x)
        % u1=x(1) u2=x ...