《参数拟合.pptVIP

超定方程写成矩阵表达就是确定 作图比较 1:描出原始点散点图 plot(tdata,cdata,’*’, tdata,f) legend(‘原始散点图’,’拟合曲线’) 借助fminsearch指令进行非线性最小二乘估计 （1）编写两个M函数文件：仿真数据产生文件和二乘残差计算文件。 function [x,y]=xydata(k_noise) %xydata.m 产生仿真数据 x=[0:0.2:4]; %自变量采样向量 yo=3*exp(-0.4*x)+12*exp(-3.2*x); %产生未受干扰的函数量 rand(seed,234) %均布随机发生器初始化y_noise=k_noise*(rand(size(x))-0.5); %产生均布随机数 y=yo+y_noise; %生成受干扰函数量 构建脚本M函数作为目标控制函数 function E=twoexps(a,x,y) %控制残差 x=x(:);y=y(:); Y=a(1)*exp(-a(3)*x)+a(2)*exp(-a(4)*x); %计算估计函数值 E=sum((y-Y).^2); %二乘残差 运行主要命令计算未