24小时热门版块排行榜    

查看: 3598  |  回复: 17
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

musejianglin

银虫 (小有名气)

[求助] 求助,用matlab解微分方程组,希望高手能给予指点

微分方程的形式如下所示:
d2y/dx2 =a*sinh(y)
边界条件为y(0)=b,dy/dx(x=0处)=c。
求y与x的曲线关系。
因为已经知道y 不存在解析解,只存在数值解,而且文献中提到是用matlab中的bvp4c路径做的,所以请matlab高手和数学高手驻足,帮忙解决一下,不胜感激。
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

somomo91

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答

xxz903

金虫 (小有名气)

用1stopt应该也是可以求解的~
2楼2013-05-14 23:20:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

arising2010

铜虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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
3楼2013-05-15 01:26:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

musejianglin

银虫 (小有名气)

引用回帖:
2楼: Originally posted by xxz903 at 2013-05-14 23:20:11
用1stopt应该也是可以求解的~

具体怎么做,可不可以指教?
4楼2013-05-15 09:03:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 06:46 by ZPa0EcMwuECS
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 06:43 by ZPa0EcMwuECS
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-30 05:48 by winsaint
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 13/650 2026-08-30 05:38 by ZPa0EcMwuECS
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 02:40 by ZPa0EcMwuECS
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:22 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:47 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 01:36 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 01:14 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 00:41 by ZPa0EcMwuECS
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 9/450 2026-08-29 09:59 by jklily
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
信息提示
请填处理意见