www.gusucode.com > matlab通信工程仿真源码（张德丰等编著）程序书籍 > matlab_code/matlab通信工程仿真源码（张德丰等编著）/第9章/li9_26.m

    clear;
x0=-1+2*[0:10]/10;
y0=1./(1+25*x0.^2);
x=-1:0.01:1;
y1=1./(1+25*x.^2);
p2=polyfit(x0,y0,2);
y2=polyval(p2,x);
p5=polyfit(x0,y0,5);
y5=polyval(p5,x);
p8=polyfit(x0,y0,8);
y8=polyval(p8,x);
p10=polyfit(x0,y0,10);
y10=polyval(p10,x);
plot(x,y1,x,y2,'r:',x,y5,'p',x,y8,'-',x,y10,'-.');
legend('原函数','二次拟合','五次拟合','八次拟合','十次拟合');