| 查看: 3491 | 回复: 17 | ||||
musejianglin银虫 (小有名气)
|
[求助]
求助，用matlab解微分方程组，希望高手能给予指点
|
|||
|
微分方程的形式如下所示： d2y/dx2 =a*sinh(y) 边界条件为y(0)=b,dy/dx(x=0处)=c。 求y与x的曲线关系。 因为已经知道y 不存在解析解，只存在数值解，而且文献中提到是用matlab中的bvp4c路径做的，所以请matlab高手和数学高手驻足，帮忙解决一下，不胜感激。 |
回复此楼» 猜你喜欢
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有1人回复
-
如何从铁电相到顺电相。
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有53人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有41人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有24人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有1人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请问这个4阶的非线性偏微分方程组在PDE中怎么解？等待高手解答。
已经有17人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- 如何用matlab求解矩阵系数的二阶微分方程
已经有13人回复
- matlab求解偏微分方程组
已经有3人回复
- 求助 MATLAB解方程组-fslove
已经有7人回复
- 求助matlab---fsolve解非线性方程组
已经有6人回复
- 高手进来用matlab求解常微分方程，十分感谢！
已经有17人回复
- matlab有限差分法求解复杂的偏微分方程组
已经有19人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- matlab解微分方程
已经有10人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】有限元法求解二阶偏微分方程组
已经有6人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
xxz903
金虫 (小有名气)
- 应助: 33 (小学生)
- 金币: 1045.6
- 散金: 190
- 红花: 4
- 帖子: 253
- 在线: 66.7小时
- 虫号: 1360336
- 注册: 2011-08-03
- 性别: GG
- 专业: 化工系统工程
2楼2013-05-14 23:20:11
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68163.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.3小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
csgt0: 金币+2, 谢谢 2013-05-15 17:44:25
musejianglin: 金币+15 2013-05-16 21:19:25
musejianglin: 金币+3 2013-05-20 09:05:12
感谢参与，应助指数 +1
csgt0: 金币+2, 谢谢 2013-05-15 17:44:25
musejianglin: 金币+15 2013-05-16 21:19:25
musejianglin: 金币+3 2013-05-20 09:05:12
|
你的方程少个边界条件吧，是不是少个"dy/dx(x=某数值)= ”的边界条件？ 不然按你目前的边界条件，都是x=0处的初始值，这个微分方程已经不是边值问题了，而是初值问题了。在现有条件下，解该问题用ode可解，程序如下，复制进一个m文件，F5运行即可，根据需要输入a，b，c的数值。a=b=c=1的结果见附图1，该方程组有刚性，不知道你希望x什么多大的区间内求解。补充完边界条件后，可再发上来。 function solvebvp333 clear all;clc global a b c a=input('参数a='); b=input('参数b='); c=input('参数c='); y0=[1 1]; xspan=0:0.02:1; [x,y]=ode15s(@fun332,xspan,y0); figure plot(x,y(:,1),'bo-'); function dydx=fun332(x,y) global a dydx=[y(2);a*sinh(y(1))];  附图1.jpg |
6楼2013-05-15 09:24:36
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68163.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.3小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
7楼2013-05-15 09:32:00
arising2010
铜虫 (正式写手)
- 应助: 32 (小学生)
- 金币: 1433.4
- 红花: 11
- 帖子: 996
- 在线: 412.6小时
- 虫号: 2464689
- 注册: 2013-05-15
- 专业: 可再生与替代能源利用中的
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
musejianglin: 金币+2, ★有帮助 2013-05-15 09:06:33
fegg7502: 金币+2, 应助指数+1, 3ks 2013-05-18 14:31:46
musejianglin: 金币+3 2013-05-20 09:04:56
感谢参与，应助指数 +1
musejianglin: 金币+2, ★有帮助 2013-05-15 09:06:33
fegg7502: 金币+2, 应助指数+1, 3ks 2013-05-18 14:31:46
musejianglin: 金币+3 2013-05-20 09:04:56
|
http://wenku.baidu.com/view/953dce370912a21614792967.html http://wenku.baidu.com/view/690deff2770bf78a65295436.html 可以参考一下上面两个连接 |
3楼2013-05-15 01:26:56
musejianglin
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1301
- 帖子: 91
- 在线: 16.6小时
- 虫号: 1225342
- 注册: 2011-03-08
- 性别: MM
- 专业: 胶体与界面化学
4楼2013-05-15 09:03:36
musejianglin
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1301
- 帖子: 91
- 在线: 16.6小时
- 虫号: 1225342
- 注册: 2011-03-08
- 性别: MM
- 专业: 胶体与界面化学
5楼2013-05-15 09:06:08
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
csgt0: 金币+1, 谢谢 2013-05-16 17:09:32
musejianglin: 金币+15, ★★★很有帮助 2013-05-16 21:19:06
感谢参与，应助指数 +1
csgt0: 金币+1, 谢谢 2013-05-16 17:09:32
musejianglin: 金币+15, ★★★很有帮助 2013-05-16 21:19:06
|
不必用1stopt， MATLAB 完全可以做到 这里是程序 ： =================================== % main program close all; clear all; clc global a a = 0.1; t0 = 0; tend=5000; tspan = t0:0.5:tend; v_y_ini = [1 0]; [t,Y] = ode45(@ode_sys, tspan, v_y_ini); v = Y(:, 1); y = Y(:, 2); =================================== % ode system function [ dy ] = ode_sys( ~, Y ) global a v = Y(1); y = Y(2); dy = zeros(2,1); dy(1) = a * sinh( y ); dy(2) = v ; end =================================== |
8楼2013-05-15 18:28:25
9楼2013-05-15 18:29:29
musejianglin
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1301
- 帖子: 91
- 在线: 16.6小时
- 虫号: 1225342
- 注册: 2011-03-08
- 性别: MM
- 专业: 胶体与界面化学
10楼2013-05-16 21:46:54
50