ch8example18.m

通信程序源代码和模型通信程序源代码和模型通信程序源代码和模型

% ch8example18.m
x=-1:0.01:1;
y=1./(1+9*x.^2);
plot(x,y,'k');hold on;                 % 原始函数曲线
xs=[-1, -0.4, -0.1,0,0.3,0.7];         % 样本位置
ys=1./(1+9*xs.^2);
plot(xs,ys,'o');hold on;               % 样本点
p=polyfit(xs,ys,5);
yfit=polyval(p,x);
plot(x,yfit,'--');
legend('原始函数y=f(x)','样本点','拟合曲线');