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

chi2(n)=sum((y-ye).^2)/dy/dy;
end
Q=1-chi2cdf(chi2,freedom);
%
subplot(1,2,1),plot(1:6,abs(chi2-freedom),'b')
xlabel('阶次'),title('chi 2 与自由度')
subplot(1,2,2),plot(1:6,Q,'r',1:6,ones(1,6)*0.5)
xlabel('阶次'),title('Q 与 0.5 线')
16
14
12
10
8
6
4
2
chi 2与 自由度
0
0 2 4 6
阶次
图 4.10-1
22
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
Q 与 0.5 线
0
0 2 4 6
阶次
%
clf,plot(x,y,'b+');axis([0,1,1,6]);hold on
errorbar(x,YE{3},D{3},'r');hold off
title('较适当的三阶拟合')
text(0.1,5.5,['chi2=' num2str(chi2(3)) '~' int2str(freedom(3))])
text(0.1,5,['freedom=' int2str(freedom(3))])
text(0.6,1.7,['Q=' num2str(Q(3)) '~0.5'])
6
5.5
5
4.5
4
3.5
3
2.5
2
1.5
chi2=4.9707~7
freedom=7
较适 当的三阶 拟合
Q=0.66353~0.5
1
0 0.2 0.4 0.6 0.8 1
图 4.10-2
A{3},DA{3}
ans =
0.6993 1.2005 1.8869 2.1077
ans =
1.9085 2.9081 1.2142 0.1333
4.10.2 非线性最小二乘估计
4.10.2.1 伪线性最小二乘
4.10.2.2 借助 fminsearch 指令进行非线性最小二乘估计