应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 隐式的微分方程组求解问题。
81/1返回列表
查看: 1990  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

cyl19881108

铁虫 (初入文坛)

[求助] 隐式的微分方程组求解问题。

方程组如下：
dx/dt=a*x*(1-x/b)*(1/(1+s/c));
dp/dt=d*dx/dt+ex;
ds/dt=f*dx/dt+g*dp/dt+h*x;
其中a-h都是未知参数，解这个方程目的是得到x(t)、p（t）、s（t），进而使用实验数据拟合这3个方程并得到a-h参数的具体数值。
我就是想知道怎么通过上面的微分方程组求得x(t)、p（t）、s（t）这三个方程啊？
哪位大神可以教教我啊。编程这块我确实是个小白，要是哪里写的不清楚各位大神可以问我啊。
使用什么软件的什么程序能把这个微分方程组解出我想要的这三个方程呢？求各位大神指教啊。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-09-04 21:56:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
就是微分方程拟合，有数据吗？放上来看看。
赞一下 回复此楼
2楼2013-09-05 08:57:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cyl19881108

铁虫 (初入文坛)
引用回帖:
2楼: Originally posted by dingd at 2013-09-05 08:57:06
就是微分方程拟合，有数据吗？放上来看看。

t[0,6,12,18,24,30,36,42,48,54,60]
s[24,22,19.3,16.0,13.4,9.9,6.5,3.7,1.92,1.3,0.96]
x[0.12,0.27,1.321,1.264,1.22,1.33,1.321,1.245,1.264,1.151,1.151]
p[0,1.02,1.98,3.37,5.05,6.66,8.65,10.95,14.00,15.50,16.40.
赞一下 回复此楼
3楼2013-09-06 13:47:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
★ ★
csgt0: 金币+2, 谢谢 2013-09-06 15:41:28
1stOpt求解代码，x项的效果较差：
CODE:
ConstStr f1=a*x*(1-x/b)*(1/(1+s/c)),
         f2=d*f1+e*x,
         f3=f*f1+g*f2+h*x;
Variable t,x,p,s;
ODEFunction x'=f1;
            p'=f2;
            s'=f3;
Data;
t=[0,6,12,18,24,30,36,42,48,54,60];
x=[0.12,0.27,1.321,1.264,1.22,1.33,1.321,1.245,1.264,1.151,1.151];
p=[0,1.02,1.98,3.37,5.05,6.66,8.65,10.95,14.00,15.50,16.40];
s=[24,22,19.3,16.0,13.4,9.9,6.5,3.7,1.92,1.3,0.96];

均方差(RMSE):0.499058750578204
残差平方和(SSE):7.47178909586035
相关系数(R): 0.999326893708328
相关系数之平方(R^2): 0.998654240488736
确定系数(DC): 0.998768314212006
F统计(F-Statistic): 0.0522044372000502

参数                  最佳估算
--------------------        -------------
a        0.0446360761798174
b        -4203.33960099553
c        18491.7942131467
d        -14454.704777378
e        645.620163822298
f        738.371259987294
g        -2.60499188266666
h        -32.3637736613547
隐式的微分方程组求解问题。
c1.jpg


隐式的微分方程组求解问题。-1
c2.jpg


隐式的微分方程组求解问题。-2
c3.jpg
赞一下(2人) 回复此楼
4楼2013-09-06 14:45:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cyl19881108

铁虫 (初入文坛)
引用回帖:
4楼: Originally posted by dingd at 2013-09-06 14:45:43
1stOpt求解代码，x项的效果较差：

ConstStr f1=a*x*(1-x/b)*(1/(1+s/c)),
         f2=d*f1+e*x,
         f3=f*f1+g*f2+h*x;
Variable t,x,p,s;
ODEFunction x'=f1;
            p'=f2;
            s'=f ...

你用的是1stOpt 哪个版本的啊 ？ 第一个拟合的好像不好啊。。
赞一下 回复此楼
5楼2013-09-06 20:21:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cyl19881108

铁虫 (初入文坛)
引用回帖:
4楼: Originally posted by dingd at 2013-09-06 14:45:43
1stOpt求解代码，x项的效果较差：

ConstStr f1=a*x*(1-x/b)*(1/(1+s/c)),
         f2=d*f1+e*x,
         f3=f*f1+g*f2+h*x;
Variable t,x,p,s;
ODEFunction x'=f1;
            p'=f2;
            s'=f ...

我只有istOpt1.5哦，好像不能带着微分拟合啊。。所以想把微分方程组先解出来。
赞一下 回复此楼
6楼2013-09-06 20:26:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

怀疑你x的数据是否正确（看看图）。

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
7楼2013-09-06 22:11:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cyl19881108

铁虫 (初入文坛)
引用回帖:
7楼: Originally posted by dingd at 2013-09-06 22:11:08
怀疑你x的数据是否正确（看看图）。

X数据没错啊，我源数据的图也是那样的
赞一下 回复此楼
8楼2013-09-06 22:26:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cyl19881108 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿211，0703化学310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-21 22:21 by peike
[考研] 材料求调剂 +5 @taotao 2026-03-21 5/250 2026-03-21 20:55 by lbsjt
[考研] 297求调剂 +3 喜欢还是不甘心 2026-03-20 3/150 2026-03-21 18:33 by 学员8dgXkO
[考研] 265求调剂 +12 梁梁校校 2026-03-19 14/700 2026-03-21 13:38 by lature00
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 316求调剂 +6 梁茜雯 2026-03-19 6/300 2026-03-21 06:32 by Ecowxq666！
[考研] 一志愿天津大学化学工艺专业（081702）315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 085700资源与环境308求调剂 +12 墨墨漠 2026-03-18 13/650 2026-03-21 01:42 by JourneyLucky
[考研] 材料专硕英一数二306 +7 z1z2z3879 2026-03-18 7/350 2026-03-20 23:48 by JourneyLucky
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 5/250 2026-03-20 22:11 by 云游重阳
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +5 西南交通专材355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
[考研] 一志愿华中农业071010，总分320求调剂 +3 困困困困坤坤 2026-03-20 3/150 2026-03-20 20:38 by 学员8dgXkO
[考研] 289求调剂 +6 怀瑾握瑜l 2026-03-20 6/300 2026-03-20 20:30 by 学员8dgXkO
[考研] 一志愿 南京航空航天大学大学 ，080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 生物学071000 329分求调剂 +3 我爱生物生物爱� 2026-03-17 3/150 2026-03-18 10:12 by macy2011
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 一志愿，福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
信息提示
关闭
请填处理意见
关闭