matlab里1stopt,请帮忙运行下1stopt 程序？（MATLAB没有解出结果）

最新推荐文章于 2021-03-24 14:34:16 发布

转载最新推荐文章于 2021-03-24 14:34:16 发布 · 407 阅读

文章标签：

博客内容涉及使用ODE方程进行数值模拟，尝试对四个参数n, v, c, s进行拟合，这些参数的取值范围分别为n[0,2], v[0,2], c[3,27], s[0,1]。已有的数据包括五个已测变量cno, ctan, cchl, a, cdrp。博主分享了尝试MATLAB求解但效果不佳的经历，并希望借助1stopt等其他工具寻求更好的解决方案。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

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