ch8example20prog1.m

来自「清华大学出版社 邵玉斌编写的《通信系统建模与仿真实例分析》一书的所有MATLAB」· M 代码 · 共 15 行

% ch8example20prog1.m
function ch8example19prog12()
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_lsq,Jmin] = lsqcurvefit(@fun2,[1.1,0.9,8],xs,ys)% 最小二乘法拟合
[a_lsq,r,J] =nlinfit(xs,ys,@fun2,[1.1,0.9,8])% 最小二乘法拟合
yfit=fun2(a_lsq,x);
plot(x,yfit,'.');
legend('原始函数y=f(x)','样本点','最小二乘法拟合');
function yhat=fun2(a,x)
yhat=a(1)./(a(2)+a(3).*x.^2);