lanxiang 发表于 2011-10-19 09:29

帮忙运行下1stopt 程序?(MATLAB没有解出结果)

问题描述:4个待拟合未知参数 n,v,c,s(实际有8个,这里先能拟合好4个就很不错了),五个已测变量cno,ctan,cchl,a,cdrp,三个ODE 方程。(已经尝试用过MATLAB求解,结果不理想,所以想尝试下1stopt)。另外,4个待拟合未知参数 n,v,c,s的范围分别是n,v,c,s. 但不知该如何写成code。

我的版本运行不了下面的code,说有语法错误;请有高版本的朋友帮忙运行下,万分感谢。(参数不在范围也行,在已知范围最好)


Parameters n,v,c,s;
Variable t,cno,ctan,cchl,a,cdrp;
ODEFunction cno'=ctan*n - cno/10 - (88754884465587073*cchl*cdrp*cno*c*(1/exp(12/(5*exp((357*cchl)/20 + 6))) - 0.0907))/(22517998136852480*((119*cchl)/10 + 4)*(cdrp + 1/200)*(cno + ctan + 1/125));
ctan'=a - ctan*(n + v + 1/10) - (88754884465587073*cchl*cdrp*ctan*c*(1/exp(12/(5*exp((357*cchl)/20 + 6))) - 0.0907))/(22517998136852480*((119*cchl)/10 + 4)*(cdrp + 1/200)*(cno + ctan + 1/125));
cchl'=(88754884465587073*cchl*cdrp*(1/exp(12/(5*exp((357*cchl)/20 + 6))) - 0.0907)*(cno + ctan))/(22517998136852480*((119*cchl)/10 + 4)*(cdrp + 1/200)*(cno + ctan + 1/125)) - cchl*(s + 1/10);
Data;
1        0.023836178        1.378559603        0.048573        0.0150 0.0165114
15        0.031473861        2.249387417        0.127282        0.0618 0.03788735
29        0.084099887        3.435099338        0.209995        0.1301 0.09540635
43        0.114199706        3.299503311        0.239937        0.2092 0.1421958
57        0.114486588        4.05281457        0.254132        0.2908 0.29265445
71        0.234158178        4.090480132        0.329151        0.3694 0.3577899
85        1.116465505        5.032119205        0.554579        0.4410 0.36191775
99        1.647905083        6.546274834        0.786195        0.5038 0.4706321
页: [1]
查看完整版本: 帮忙运行下1stopt 程序?(MATLAB没有解出结果)