kkd.rar_matlab 数据_变形_变形预测_预测 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