7楼: Originally posted by dingd at 2011-12-05 10:57:31:
这种动力学微分方程拟合，1stOpt很简单好使啊！
[code]
Variable T, yA, yB, yC, yD, yE;
ODEFunction yA'=-(k1+k2+k3+k4)*yA;
            yB'=k1*yA-(k5+k6+k7)*yB;
            yC'=k2*yA+k5*yB-(k8+k9)* ...

Parameter x,y,z,k2,k3,k4,k5,k6,k7,k8,k9,k10,k11;
Variable t,y(1:4);
ODEFunction   Y1'=x*k4*y5*y7+k6*y5+k10*y2*y7-0.1942*k10*y1*y4;
Y2'=x*k3*y5*y7+k5*y5-k10*y2*y7+0.1942*k10*y1*y4-k11*y2*y4+2.7624E-6*k11*y3*y7;
Y3'=k7*y5+k11*y2*y4-2.7624E-6*k11*y3*y7;
Y4'=(x-z+y/2)*k3*y5*y7+(2x-z+y/2)*k4*y5*y7+k8*y5+k10*y2*y7-0.1942*k10*y1*y4-3*k11*y2*y4+8.2872E-6*k11*y3*y7;
Y5'=k2*y7-k3*y5*y7-k4*y5*y7-k9*y5;
Y6'=-k2*y6;
Y7'=-(2x-z)*k4*y5*y7-k10*y2*y7+0.1942*k10*y1*y4+k11*y2*y4-2.7624E-6*k11*y3*y7;
Data;
//t, y1, y2, y3, y4
2.5                3.8                2.5                0.5                0.5
5                5.8                3.3                1.0                0.6
7.5                5.8                2.7                1.2                0.7
10                5.4                2.3                1.2                0.8
15                5.7                2.4                1.2                1.3
30                8.5                2.6                2.3                2.4

12楼: Originally posted by dingd at 2011-12-12 21:21:52:
下面是正确代码，待求参数数目比数据行数还多，解估计不是唯一的，你自己用1stOpt慢慢跑下吧：
[code]
Parameter x,y,z,k2,k3,k4,k5,k6,k7,k8,k9,k10,k11;
Variable t,y(1:4);
ODEFunction   Y1'=x*k4*y5*y7+ ...

2楼: Originally posted by leedobb at 2011-03-21 15:36:02
% 哦，你的题目是这个啊，
% 这只能依靠数值解，matlab肯定没有现成的程序吧。
%
% 你这样做吧。构造一个能量函数（其实就是最小二乘法啦）
%
% F =  (yA(t2)- YYA(t2,k1,k2,k3...,k9) ) ^ 2
% ..嗯，算 ...