funs2.m

来自「基于极大嫡谱估计准则对变形数据进行预测 基于MATLAB平台编制了相应程序」· M 代码 · 共 31 行

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;
end    
k1=mean(y1);
k2=mean(y2);
k3=mean(y3);
Fx=exp(randa(1)+randa(2).*x+randa(3).*x.^2+randa(4).*x.^3);
%Fx=taylor(exp(randa(1)+randa(2).*x+randa(3).*x.^2+randa(4).*x.^3),8,x);

fx1=inline(Fx); 
mm1=quad(fx1,-1000,1000);
m(1)=mm1-1;

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

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

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

m=double(m)