24小时热门版块排行榜    

查看: 2987  |  回复: 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的回帖
信息提示
请填处理意见