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

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