24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 2278  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +9 小聂爱学习 2026-04-11 13/650 2026-04-14 19:18 by Art1977
[考研] 327求调剂 +23 Xxjc1107. 2026-04-13 25/1250 2026-04-14 18:48 by 逍遥三郎
[考研] 26药学专硕105500求调剂 +6 喽哈加油 2026-04-13 6/300 2026-04-14 16:40 by zhouxiaoyu
[考研] 245求调剂 +6 冰糖橘?汽水 2026-04-13 10/500 2026-04-14 10:49 by jyl0317
[考研] 085600材料与化工349分求调剂 +16 李木子啊哈哈 2026-04-12 17/850 2026-04-14 09:11 by fenglj492
[考研] 一志愿华南理工大学331分材料求调剂 +10 天下ww 2026-04-09 11/550 2026-04-13 23:25 by pies112
[基金申请] 2026 WR青拔 +3 冬日阳光CAS 2026-04-09 6/300 2026-04-13 18:40 by liuchb715
[考研] 生物学调剂,一志愿西南大学348,Top期刊一区二作、二区三作,三等奖学金三次 +5 candyyyi 2026-04-09 5/250 2026-04-13 09:02 by 可淡不可忘
[考研] 求调剂288 +7 ioodiiij 2026-04-10 9/450 2026-04-13 08:33 by Hayaay
[考研] 296求调剂 +14 汪!?! 2026-04-10 16/800 2026-04-12 10:48 by zhouyuwinner
[找工作] 山东高校教师考核超级无底线,员工过不下去啦 +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[考研] 283求调剂 +22 那个噜子 2026-04-09 22/1100 2026-04-11 10:41 by 逆水乘风
[考研] 机械专硕270求调剂,接受跨专业 +12 老师看看我吧aba 2026-04-09 14/700 2026-04-11 10:21 by laoshidan
[考研] 22408 366分,本科211,一志愿西工大 +4 Rubt 2026-04-09 4/200 2026-04-10 19:51 by chemisry
[考研] 298求调剂 +13 钉叮咚冬瓜 2026-04-09 13/650 2026-04-10 15:49 by jiajinhpu
[考研] 332,085601求调剂 +12 ydfyh 2026-04-09 14/700 2026-04-09 17:28 by wp06
[考研] 083200 初试305分 求调剂 暂不考虑跨专业 +15 Claireyyyy 2026-04-09 15/750 2026-04-09 16:11 by zhuimr
[考研] 334求调剂 +16 Riot2025 2026-04-08 17/850 2026-04-09 09:28 by wdyheheeh
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[考研] 求调剂,现在还能填的 +3 上岸小莹加油 2026-04-08 3/150 2026-04-08 14:30 by zhq0425
信息提示
请填处理意见