24小时热门版块排行榜    

查看: 3182  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zhl717

铜虫 (小有名气)

[求助] matlab拟合微分方程组中的参数 已有2人参与

最近一直在看matlab求解微分方程组中参数问题,但是由于本人以前没有接触过matlab,所以还是不太懂,请懂的大神帮忙写一下,谢谢!!

有高版本1st opt的大神也可以帮忙跑一下,上回发了一个帖子,没人回应,谢谢各位!!

dx/dt=a*x*(1-x/b)*(1+s/c)^(-1)
dp/dt=d*dx/dt+e*[s/(s+f)]*x
ds/dt=g*dx/dt+h*dp/dt+i*x

t为自变量,x、p、s为因变量,a-i  均为参数

数据如下:
t,   X,   P,       S
0   0.35  0      99.8
12  2.75  0.55   88.96
24  4.6   3.72   73.49
36  6.27  8.19   61.57
48  8.07  14.12  47.38
60  9.73 19.37  33.63
72  10.41 22.54  24.65
84  10.7  27.61  11.19
96  10.53 32.49  3.2
108 10.56 35.62  0
120 10.59 37.58  0
132 10.4  38.44  0
144 10.77 39.88  0
回复此楼

» 本帖@通知

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

加油!!!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhl717

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by 月只蓝 at 2015-08-31 16:48:55
要拟合九个参数,这种问题还是推荐用1stOpt,参数个数太多,用MATLAB比较困难。
下面给出MATLAB参考代码,代码中初值k0可以试着调一下,我初步调了一下初值,拟合效果不好。

function k1k2k3
format long
cle ...

好的,非常感谢
加油!!!
3楼2015-08-31 17:03:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhl717 的主题更新
信息提示
请填处理意见