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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 隐式的微分方程组求解问题。
81/1返回列表
查看: 2175  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] filecode=后面第一个是大写字母 +7 wangze12014 2026-08-14 8/400 2026-08-17 12:45 by a441914426
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +12 Ldrop2023 2026-08-13 13/650 2026-08-17 12:33 by gltch
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
信息提示
关闭
请填处理意见
关闭