funs1.m

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

function m=Funs1(randa)
syms x 
a=randa(1);
b=randa(2);
c=randa(3);
d=randa(4);
e=randa(5);
 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(a+b.*x+c.*x.^2+d.*x.^3+e.*x.^4);
%Fx=taylor(Fx,5,x)
%digits(4)
%aa=sym(1000)

mm1=taylor(int(Fx,x))
m(1)=subs(mm1,'x','1000')-subs(mm1,'x','0')-1

mm2=taylor(int(Fx.*x,x))
m(2)=subs(mm2,'x','1000')-subs(mm2,'x','0')-k1

mm3=taylor(int(Fx.*x.^2,x))
m(3)=subs(mm3,'x','1000')-subs(mm3,'x','0')-k2

mm4=taylor(int(Fx.*x.^3,x))
m(4)=subs(mm4,'x','1000')-subs(mm4,'x','0')-k3


mm5=taylor(int(Fx.*x.^4,x))
m(5)=subs(mm5,'x','1000')-subs(mm5,'x','0')-k4