ch8example19prog1.m

清华大学出版社 邵玉斌编写的《通信系统建模与仿真实例分析》一书的所有MATLAB和SIMULINK代码

% ch8example19prog1.m
function ch8example19prog1()
x=-1:0.01:1;
y=1./(1+9*x.^2);
plot(x,y,'k');hold on;                 % 原始函数曲线
xs=-1:0.35:1;                          % 样本位置
ys=1./(1+9*xs.^2);
plot(xs,ys,'o');hold on;               % 样本点
a_poly=polyfit(xs,ys,5)
yfit=polyval(a_poly,x);                % 多项式拟合曲线
plot(x,yfit,'-.');
[a_lsq,Jmin] = lsqcurvefit(@fun1,[1,1,1,1,1,1],xs,ys)% 最小二乘法拟合
yfit=fun1(a_lsq,x);
plot(x,yfit,'.');
legend('原始函数y=f(x)','样本点','多项式拟合曲线','最小二乘法拟合');
function yhat=fun1(a,x)
yhat=a(1)*x.^5+a(2)*x.^4+a(3)*x.^3+a(4)*x.^2+a(5)*x+a(6);