声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 3510|回复: 0

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

[复制链接]
发表于 2011-10-19 09:29 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
问题描述: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
回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-26 03:30 , Processed in 0.139712 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表