24小时热门版块排行榜    

查看: 2505  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:55 by OEbVnUOu01ol
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:43 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 03:40 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[基金申请] 2026国自然函评费到账 +11 羊腰板 2026-08-21 11/550 2026-08-22 22:03 by lsbin3733
[基金申请] 范进中举一文的中心思想 +3 zjuhero 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 8/400 2026-08-22 18:29 by 淀粉搬运工
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见