请帮我看看这个程序那个地方出问题了，多谢多谢！

maoking01 · 发表于 2006-3-24 01:33:36

<

>请帮我看看这个程序那个地方出问题了，多谢多谢！
<

>这个程序（matlab)是用来进行非线性参数求解的，我使用了另一个程序改编的，可是就是老出错。请哪位高手帮忙看看。
<

>function KineticsEst6 % 动力学ODE方程模型的参数估计 % 
clear all clc tspan = [0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100]; x0 = [7.7339; 0; 0;0]; k0 = [0.02 0.005 1.25 0.123 0.14 0.13 0.21 0.234 0.536]; lb = [0 0 0 0 0 0 0 0 0]; ub = [10 10 10 10 10 10 10 10 10];
KineticsData2; yexp = Kinetics(:,2:5)
% 使用函数lsqnonlin()进行参数估计 [k,resnorm,residual,exitflag,output,lambda,jacobian]= ... lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,x0,yexp); ci = nlparci(k,residual,jacobian); fprintf('\n\n使用函数lsqnonlin()估计得到的参数值为:\n') fprintf('\tk1 = %.4f ± %.4f\n',k(1),ci(1,2)-k(1)) fprintf('\tk2 = %.4f ± %.4f\n',k(2),ci(2,2)-k(2)) fprintf('\tk3 = %.4f ± %.4f\n',k(3),ci(3,2)-k(3)) fprintf('\tk3 = %.4f ± %.4f\n',k(4),ci(4,2)-k(4)) fprintf('\tk3 = %.4f ± %.4f\n',k(5),ci(5,2)-k(5)) fprintf('\tk3 = %.4f ± %.4f\n',k(6),ci(6,2)-k(6)) fprintf('\tk3 = %.4f ± %.4f\n',k(7),ci(7,2)-k(7)) fprintf('\tk3 = %.4f ± %.4f\n',k(8),ci(8,2)-k(8)) fprintf('\tk3 = %.4f ± %.4f\n',k(9),ci(9,2)-k(9)) fprintf(' The sum of the squares is: %.1e\n\n',resnorm)
% ------------------------------------------------------------------ function f=ObjFunc(k,tspan,x0,yexp) % 目标函数 [t Xsim]=ode45(@KineticsEqs1,tspan,x0,[],k); ysim(:,1)=Xsim(2:end,1); ysim(:,2)=Xsim(2:end,2); ysim(:,3)=Xsim(2:end,3); ysim(:,4)=Xsim(2:end,4); f=[ysim(:,1)-yexp(:,1);ysim(:,2)-yexp(:,2);ysim(:,3)-yexp(:,3);ysim(:,4)-yexp(:,4)];
% ------------------------------------------------------------------ function dCdt = KineticsEqs1(t,c,k) % ODE模型方程
r1 = k(1)*c(1)/(k(4)*c(1)+k(5)*c(2)+k(6)*c(3)+k(7)); r2 = k(2)*c(2)/(k(4)*c(1)+k(5)*c(2)+k(6)*c(3)+k(7)); r3 = k(3)*c(2)/((k(4)*c(1)+k(5)*c(2)+k(6)*c(3)+k(8))^0.574); r4=k(3)*c(3)/(k(4)*c(1)+k(5)*c(2)+k(6)*c(3)+k(9));
dCAdt=-r1-r2; dCBdt=r1+r4-r3; dCCdt=r2-r4; dCDdt=r3; %dCAdt = -r1; %dCBdt = 0.5*r1-r2 + r3; %dCCdt = 0.5*r1-r3;
dCdt = [dCAdt; dCBdt; dCCdt;dCDdt];

		自动登录	找回密码
密码			注-册-帐-号