| 查看: 1092 | 回复: 2 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
m76acxb铜虫 (小有名气)
|
[求助]
耦合方程组的四阶龙格库塔编程 已有1人参与
|
|
|
对于耦合方程组：1.目的是要画plot（X,H)，但是在边界条件中出现了含T的函数； 2.边界条件：T为任意时刻时，X=0，H=H0; X范围为【0 1】 3.H0=0.03*sin(30*T)+1; 4.方程组中参数：ε=0.2，U1=0.8，第三个式子求导已经写在第二个图片中。 希望能够用四阶龙格库塔编程，不胜感激！   耦合方程组  相关式子.gif |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求用龙格库塔法求解下列方程
已经有4人回复
- 用四阶龙格库塔法求解常微分方程组
已经有15人回复
- 求助！！对于调用龙格库塔法如何对结果判断反馈
已经有6人回复
- 求助，帮忙求解下面方程组，谢谢！
已经有13人回复
- 急！求大神用Matlab四阶龙格库塔解个方程！
已经有5人回复
- 请问什么是分段龙格库塔的方法？需要用matlab软件采用这种方法求二阶微分方程的解
已经有3人回复
- mathcad怎样用四阶龙格库塔法求解二阶微分方程？
已经有8人回复
- 求解含高次项的一阶微分方程
已经有5人回复
- 四阶非线性偏微分方程组求解
已经有6人回复
- 非线性微分耦合方程的解法
已经有9人回复
- Matlab用龙格库塔法求二阶非线性方程
已经有13人回复
- 【求助】非线性方程组的求解问题
已经有6人回复
- 二阶非线性微分方程求解
已经有14人回复
- 请教mathematic问题？关于解微分方程组的和画图的一些问题。
已经有19人回复
- 龙格库塔法 matlab
已经有6人回复
- Matlab求解偏微分方程组
已经有22人回复
- LABVIEW四阶龙格库塔法求解一阶微分方程组
已经有12人回复
- 【求助】这个微分方程组simulink求解到底错在哪里？
已经有3人回复
- 【求助】Matlab中利用四阶龙格-库塔法求解微分方程！！！！
已经有9人回复
- 【求助】四阶龙格-库塔方法！！！！
已经有17人回复
m76acxb
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 209.5
- 散金: 26
- 帖子: 62
- 在线: 27.9小时
- 虫号: 3628589
- 注册: 2015-01-05
- 性别: GG
- 专业: 机械动力学
3楼2015-05-23 11:37:28
eidie
木虫之王 (文坛精英)
- 应助: 15 (小学生)
- 贵宾: 1.016
- 金币: 71294.4
- 散金: 8785
- 红花: 240
- 沙发: 1
- 帖子: 13888
- 在线: 1583.2小时
- 虫号: 2259480
- 注册: 2013-01-24
- 专业: 应用光学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
xiegangmai: 金币+2, 鼓励讨论交流！ 2015-05-22 17:41:21
m76acxb: 金币+5, ★有帮助 2015-05-24 09:12:28
感谢参与，应助指数 +1
xiegangmai: 金币+2, 鼓励讨论交流！ 2015-05-22 17:41:21
m76acxb: 金币+5, ★有帮助 2015-05-24 09:12:28
|
实在没看懂，方程组第一个式子对H和X影响没多大，就影响了H0初始值，我想H0初始值由你的三角函数T决定。 其它都像是一阶微分方程。 clc; clear; epsilon=0.2; U1=0.8; T=0; H0=0.03*sin(30*T)+1; dHX=@(X,H)[H/U1*(2*(1-epsilon)+2*epsilon*(1-X))*0.64*(1-X)/(1.6+0.4*(1-X)^2)^2]; [X,H]=ode45(dHX,[0 1],H0); plot(X,H,'r'); |
2楼2015-05-22 16:10:14
