funs2.m

基于极大嫡谱估计准则对变形数据进行预测 基于MATLAB平台编制了相应程序

function m=Funs2(randa)
syms x
 y=[18.85 22.44 20.09 21.36 18.2 20.9 22.04 22.02 22.01 19.72 12.05 18.13 22.74 17.75 19.83 21.57 20.62 21.82 19.86 21.11 22 20.74 18.7 20.28 20.26 21.17 21.37 21.08 19.15];
for i=1:29
    y1(i)=y(i);
    y2(i)=y(i)^2;
    y3(i)=y(i)^3;
    y4(i)=y(i)^4;
end    
k1=mean(y1);
k2=mean(y2);
k3=mean(y3);
k4=mean(y4);
Fx=exp(randa(1)+randa(2).*x+randa(3).*x.^2+randa(4).*x.^3+randa(5).*x.^4);
%Fx=taylor(exp(randa(1)+randa(2).*x+randa(3).*x.^2+randa(4).*x.^3+randa(5).*x.^4),8,x);

%fx1=inline(Fx); 
fx1=Fx
mm1=comsimpon(Fx,0,1000,10000);
m(1)=mm1-1;

fx2=inline(Fx*x);
mm2=comsimpon(fx2,0,1000,10000);
m(2)=mm2-k1;

fx3=inline(Fx*x^2);
mm3=comsimpon(fx3,0,1000,10000);
m(3)=mm3-k2;

fx4=inline(Fx*x^3);
mm4=comsimpon(fx4,0,1000,10000);
m(4)=mm4-k3;

fx5=inline(Fx*x^4);
mm5=comsimpon(fx5,0,1000,10000);
m(5)=mm5-k4;

m=double(m)