24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2301  |  回复: 10

1135725495

铁杆木虫 (著名写手)

[求助] 急!!!用matlab求解微分方程的问题 已有4人参与

各位虫友好,我对matlab不怎么了解,现在需要用这个软件进行计算,可是不知道是否可以使用这个软件,我遇到的问题是,一个微分方程,有一个自变量,两个因变量。我只会一个因变量,一个自变量的。不知道两个因变量该怎么编写程序。希望懂matlab帮忙看一下。不胜感激,
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

yooy

金虫 (初入文坛)

铁道出版社 钟尔杰写的一部数学实验方面的教材,里面有这样的例子。
可以参考一下。
555555
5楼2015-06-14 13:01:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1135725495

铁杆木虫 (著名写手)

引用回帖:
4楼: Originally posted by 月只蓝 at 2015-06-14 09:30:51
既然只有一个自变量,这种情况应该是常微分方程。
同时你说有两个因变量,就是应该是常微分方程组。
但你说只有一个方程,那说明这个方程组方程不完全,无法求解(解有去穷多个);包含两个因变量的方程,需要两个 ...

恩恩好滴!
c1'=-a*(c1^2-c2/b)/(1+d*c1+f*c2)^2
c2'=-0.5c1'
其中c1与c2都是关于t的变量。
我之前看过您的帖子,您之前好像说过可以把c2拟合成关于t的函数,在代入方程中,我这个方程也可以这样做吧?
6楼2015-06-14 21:39:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

amoroso_xmc

新虫 (初入文坛)

【答案】应助回帖


感谢参与,应助指数 +1
fegg7502: 金币+1, 鼓励交流 2015-06-15 11:17:45
只有一个自变量的是常微分方程,两个自变量的就是偏微分方程。上网查一下,资料很多的。
2楼2015-06-13 22:32:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1135725495

铁杆木虫 (著名写手)

引用回帖:
2楼: Originally posted by amoroso_xmc at 2015-06-13 22:32:14
只有一个自变量的是常微分方程,两个自变量的就是偏微分方程。上网查一下,资料很多的。

您好,您说的自变量的,我知道,我现在遇到的是两个因变量,一个自变量。不知道怎么解决

[ 发自小木虫客户端 ]
3楼2015-06-13 22:39:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖


感谢参与,应助指数 +1
fegg7502: 金币+1, 鼓励交流 2015-06-15 11:17:51
既然只有一个自变量,这种情况应该是常微分方程。
同时你说有两个因变量,就是应该是常微分方程组。
但你说只有一个方程,那说明这个方程组方程不完全,无法求解(解有去穷多个);包含两个因变量的方程,需要两个方程来定解。
说了半天,不如把方程给出来,让我们看看。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
4楼2015-06-14 09:30:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

changhw

铁虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
fegg7502: 金币+1, 鼓励交流 2015-06-15 11:17:57
c2=0.5c1+常数,带到①式里不就是一个常微分方程吗😊

[ 发自小木虫客户端 ]
7楼2015-06-15 08:05:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
6楼: Originally posted by 1135725495 at 2015-06-14 21:39:42
恩恩好滴!
c1'=-a*(c1^2-c2/b)/(1+d*c1+f*c2)^2
c2'=-0.5c1'
其中c1与c2都是关于t的变量。
我之前看过您的帖子,您之前好像说过可以把c2拟合成关于t的函数,在代入方程中,我这个方程也可以这样做吧?: ...

常微分方程组可直接求解了。

» 本帖已获得的红花(最新10朵)

MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
8楼2015-06-15 08:12:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1135725495

铁杆木虫 (著名写手)


送红花一朵
fegg7502: 金币+1, 鼓励交流 2015-06-15 11:18:05
引用回帖:
8楼: Originally posted by 月只蓝 at 2015-06-15 08:12:05
常微分方程组可直接求解了。...

您好!我不会做关于这个,我用1stOpt计算不出来结果。想用matlab进行计算,可以吗?您帮我看看数据。
c1'=-a*(c1^2-c2^2/b)/(1+d*c1+f*c2)
c2'=-0.5c1'
t         c1             c2
0        7.60599        0.94886
1        4.78957        4.01159
2        4.05789        4.92061
3        3.86145        5.13061
6        3.57636        5.52913
9        3.27073        5.9891
12        3.09418        6.23705
15        2.93546        6.38008
18        2.7715        6.54667
21        2.63196        6.65208
31        2.28437        6.98497
41        2.01135        7.20437
51        1.98294        7.21122
61        1.76631        7.40044
9楼2015-06-15 09:01:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖


感谢参与,应助指数 +1
fegg7502: 金币+1, 鼓励交流 2015-06-15 11:18:12
1stOpt计算:
CODE:
ConstStr s=-a*(c1^2-c2^2/b)/(1+d*c1+f*c2);
Variable t,c1,c2;
ODEFunction c1'=s;
            c2'=-0.5*s;
Data;
t         c1             c2
0        7.60599        0.94886
1        4.78957        4.01159
2        4.05789        4.92061
3        3.86145        5.13061
6        3.57636        5.52913
9        3.27073        5.9891
12        3.09418        6.23705
15        2.93546        6.38008
18        2.7715        6.54667
21        2.63196        6.65208
31        2.28437        6.98497
41        2.01135        7.20437
51        1.98294        7.21122
61        1.76631        7.40044

Root of Mean Square Error (RMSE):1.88748670832911
Sum of Square Error:92.6277579270955
Correlation Coef. (R): 0.983620664786593
R-Square: 0.967509612195219
Determination Coef. (DC): -2.93279440507291
F-Statistic: 1.46948652910799

Parameter                  Best Estimate
--------------------        -------------
a        6.51641293960292
b        333.976267430181
d        -7.40707089903982
f        58.3209857722069

公式似乎有问题,过拟合,参数组值a、d、f不唯一

» 本帖已获得的红花(最新10朵)

10楼2015-06-15 10:21:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 1135725495 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然面上和省基金B类撒花 +17 花田半亩~白 2026-04-21 17/850 2026-04-23 07:30 by john198197
[考研] 有没有学校收留 +3 蒋昌鹏qtj 2026-04-20 3/150 2026-04-22 20:25 by 学员JpLReM
[考研] 312求调剂 +3 山河似你温柔 2026-04-22 3/150 2026-04-22 20:17 by 学员JpLReM
[考博] 华师大读博 +3 xq83 2026-04-22 5/250 2026-04-22 10:42 by xq83
[论文投稿] 急需审稿人!!! +3 陆小果画大饼 2026-04-21 3/150 2026-04-21 23:54 by jzy_123456
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 8/400 2026-04-21 16:26 by 啃面包的小书虫
[考研] 085600材料与化工调剂 5+3 孜孜不倦2002 2026-04-19 6/300 2026-04-20 21:25 by babero
[论文投稿] 有没有接收比较快的sci期刊呀,最好在一个月之内的,研三孩子求毕业 20+4 之护着 2026-04-16 7/350 2026-04-20 15:45 by 豆豆7758
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
信息提示
请填处理意见