内容简介作者简介 - 科学与工程计算系

16
14
12
10
8
6
4
2
y=3*exp(-0.4*x)+12*exp(-3.2*x)
y est =2.8967*exp(-0.38045*x) + 11.983*exp(-3.0985*x)
0
0 0.5 1 1.5 2 2.5 3 3.5 4
图 4.10-3
24
chi2=17.7066
freedom=17
Q=0.40757
−a1x
−a2
x
【例 4.10.2.2-2】以 y = b1e
+ b2e
为模型,从受污染的数据中,使用伪线性法估
计 b = [ b1
b2
] ,使用 fminsearch 估计 a = [ a1
a2
] 。x, y 原始数据同上例。
(1)
[twoexps2.m]
function E=twoexps2(a,x,y,b)
%twoexps2.m
x=x(:);y=y(:);Y=b(1)*exp(-a(1)*x)+b(2)*exp(-a(2)*x);
E=sum((y-Y).^2);
(2)
[exm041022_2.m]
%exm041022_2.m
k_noise=0.3;
[x,y,STDY]=xydata(k_noise);
a0=[1 2]';
options=optimset('fminsearch');
options.TolX=0.001;
options.Display='off';
%
while 1
Mb=exp(-x*a0');
b=Mb\y;
a=fminsearch(@twoexps2,a0,options,x,y,b);
r=norm(a-a0)/norm(a);
if r