24小时热门版块排行榜    

查看: 231  |  回复: 0

whqs8426212

铜虫 (正式写手)

[求助] 求助版主 1stopt 方程的拟合(含若干未知参数),麻烦了 谢谢啦

****************实验数据*************
t         0,                  5,                10,                    20,                30,                      60,           90,                   120,             150,                  180

C(1)   0.00446429 ,0.004108578,0.00390681,0.003674009,0.003399052,0.002833664,0.002330345,0.001822586,0.001385776,0.001002198




初值就是t=0时,C(1)= 0.00446429, C(2)=0, C(3)=0, C(4)=0,C(5)=0,C(6)=0,C(7)=0, C(8)=0, C(9)=0,C(10)=0, C(11)=0,  C(12)=0,  C(13)=0, C(14)=0,   C(15)=0, C(16)=0, C(17)=10^(-7), C(18)=0,C(19)=0
***********函数*************************

fuction:
dC(1)= -k(5)*C(1)*C(16)-k(6)*C(1)*C(15)-k(15)*C(1)*C(15)-k(16)*C(1)*C(15);
dC(2)=-k(7)*C(2)*C(15)-k(8)*C(2)*C(15)-k(9)*C(2)*C(15);
dC(3)=k(7)*C(2)*C(15)+k(15)*C(1)*C(15)-k(17)*C(3)*C(15);
dC(4)=k(8)*C(2)*C(15)-k(12)*C(4)*C(15);
dC(5)=k(10)*C(2)*C(15)-k(11)*C(6)*C(15);
dC(6)=k(11)*C(6)*C(15)-k(12)*C(4)*C(15);
dC(7)=k(10)*C(2)*C(15)+k(11)*C(6)*C(15)-k(13)*C(7)*C(15)+k(21)*C(10)*C(15)-k(22)*C(7)*C(15);
dC(8)=k(13)*C(7)*C(15)-k(14)*C(8)*C(15);
dC(9)=k(16)*C(1)*C(15)-k(18)*C(9)*C(15);
dC(10)=k(17)*C(3)*C(15)-k(21)*C(10)*C(15);
dC(11)=k(18)*C(9)*C(15)-k(19)*C(11)*C(15);
dC(12)=k(19)*C(11)*C(15)-k(20)*C(12)*C(15);
dC(13)=k(20)*C(12)*C(15)-k(23)*C(13)*C(15);
dC(14)=k(7)*C(2)*C(15)+k(12)*C(4)*C(15)+k(11)*C(6)*C(15) +k(15)*C(1)*C(15)+k(16)*C(1)*C(15);
dC(15)=k(1)*C(17)-k(2)*C(15)-k(3)*C(15)-k(6)*C(1)*C(15)-k(7)*C(2)*C(15)-k(8)*C(2)*C(15)-k(9)*C(2)*C(15)-k(10)*C(2)*C(15)-k(11)*C(6)*C(15)-k(12)*C(4)*C(15)-k(13)*C(7)*C(15)-k(14)*C(8)*C(15)-k(15)*C(1)*C(15)-k(16)*C(1)*C(15)-k(17)*C(3)*C(15)-k(18)*C(9)*C(15)-k(19)*C(11)*C(15)-k(20)*C(12)*C(15)-k(21)*C(10)*C(15)-k(22)*C(7)*C(15)-k(23)*C(13)*C(15);
dC(16)=k(2)*C(15)-k(5)*C(1)*C(16);
dC(17)=-k(1)*C(17)-k(2)*C(15)-k(3)*C(15);
dC(18)=0.5*k(3)*C(15)+k(4)*C(16);
dC(19)=k(1)*C(17)+k(2)*C(15)+k(3)*C(15);
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 whqs8426212 的主题更新
信息提示
请填处理意见