mypolyfit.m

来自「自己编写的matlab程序 数值分析里的程序」· M 代码 · 共 23 行

for i=1:50
   X(i)=i;
end
Y=fun(X);
A=polyfit(X,Y,1);
syms x
ya=vpa(A(1)+A(2)*x,6)
c=polyval(A,X);
plot(X,Y,'bs',X,c,'r')

B=polyfit(X,Y,2);
yb=vpa(B(1)+B(2)*x+B(3)*x^2,6)
hold on
plot(X,polyval(B,X),'y')
C=polyfit(X,Y,5);
yc=0;
for i=1:length(C)
    yc=yc+C(i)*x^(i-1);
end
vpa(yc,6)
hold on
plot(X,polyval(C,X),'g')