five.m

给定数的最小二乘拟合以及在给定两组数据下用矩法估计以及用拟和估计得对比

clc;clear;
load m5_1.dat;
% 第一问
for i=1:200
    x=rand(1,1);
    for j=1:100
        y(j)=(m5_1(j)-3*exp(x))^2;
        if j==1
            continue;
        else
            y(j)=y(j)+y(j-1);
        end
    end
    figure(1);
    plot(x,y(j),'*'),hold on;
    title('J(x)的图象')
end

% 第二问
for i=0:0.001:1
    for j=1:100
        y(j)=(m5_1(j)*exp(i)-3*exp(2*i))*(-6);
        if j==1
            continue;
        else
            y(j)=y(j)+y(j-1);
        end
    end
    figure(2);
    plot(i,y(j)),hold on;
    title('J(x)求导后的图象')
end  



% 第三问
for i=0:0.001:1
    for j=1:100
        y(j)=(m5_1(j)*exp(i)-3*exp(2*i))*(-6);
        if j==1
            continue;
        else
            y(j)=y(j)+y(j-1);
        end
    end
    if y(j)>=0
       break;
    end
    z=i;
    figure(3);
    plot(i,y(j)),hold on;
    title('求极值')
end  
z=i


% 第四问
for i=1:0.001:2.71
    for j=1:100
        y(j)=(m5_1(j)-3*i)*(-6);
        if j==1
            continue;
        else
            y(j)=y(j)+y(j-1);
        end
    end
    if y(j)>=0
       break;
    end
    zer=i;
    figure(4);
    plot(i,y(j)),hold on;
end  





% 第5问
x=log(zer)



%第六问


% 第七问
x=(z+log(zer))/2