1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28
| clc; clear;
x=0:0.2:10; y=0.25*x+20*sin(x);
p5=polyfit(x,y,5); y5=polyval(p5,x);
p8=polyfit(x,y,8); y8=polyval(p8,x);
p60=polyfit(x,y,60); y60=polyval(p60,x);
figure plot(x,y,'ro'); hold on; plot(x,y5,'b--'); plot(x,y8,'b:'); plot(x,y60,'r-.'); xlabel('x'); ylabel('y'); legend('原始数据','5阶多项式','8阶多项式','60阶多项式');
|