2014全国大学生数学建模竞赛A题论文示范(15)

仅供参考

end if k==5 break end

fprintf('**********************************\n'); end

3．问题四中模型：非线性回归法（源代码：p5.m）

xy0=[120 1200 0.6 10 120 930 0.5 300 120 840 0.4 172 120 760 0.5 65 120 1100 0.3 60 120 900 0.6 90 120 890 0.7 560 ];

x=xy0(:,[2:4]); y=xy0(:,1);

huaxue=@(beta,x)

(beta(1)-x(:,1).*beta(1).*beta(2)./x(:,2)).*(beta(3).*x(:,1)-x(:,3)); beta0=[0.05,0.005,0.02]';

[beta,r,j]=nlinfit(x,y,huaxue,beta0) betaci=nlparci(beta,r,'jacobian',j)

[yhat,delta]=nlpredci(huaxue,x,beta,r,'jacobian',j)

运行结果： beta =

0.0001 1.0243 -0.3253 r =

28.4765 -8.0898 13.1246 65.7630 -55.1525 54.3847 -3.4178 j =

1.0e+006 *

0.8056 0.0001 -0.0003 1.1274 0.0001 -0.0002 0.9407 0.0001 -0.0002 0.4774 0.0001 -0.0001 1.5416 0.0002 -0.0005