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 ...

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

4Â¥: Originally posted by xiong1021 at 2017-02-23 22:40:31
лл½â´ð£¬¿ÉÄܽá¹ûÎÞ·¨¾«È·°É£¬ÎÒÓÃfsolveÊÔÁËÒ²Êdzõʼֵ²»Ò»Ñù×îºó½á¹ûÆ«²î½Ï´ó¡£
¶ÔÁË£¬ÇëÎÊÏÂÀïÃæµÄ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