应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 非线性二阶微分方程组求解，matlab
131/2返回列表12下一页
查看: 4703  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

lylpiao

木虫 (正式写手)

[求助] 非线性二阶微分方程组求解，matlab 已有1人参与

大家好，最近在计算过程中遇见一个矩阵求解问题，方程组为非线性二阶微分方程组，我不会求解，请大家指点下，问题在图片里面~~谢谢~

非线性二阶微分方程组求解，matlab
)8F_VRMW$Z)KN)A)~EZNYPH.jpg
回复此楼

» 收录本帖的淘帖专辑推荐

博士发论文的思考体会

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
把一切交给时间！
1楼 2014-03-25 15:58:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

george1990

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lylpiao: 金币+10, ★有帮助, 谢谢 2014-03-30 11:47:01
这种问题你可以参考机械动力学里面的4阶龙格库塔解法  它对于解这种方程很ok的
赞一下(1人) 回复此楼
天天多学一点
5楼2014-03-30 01:11:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

george1990

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
lylpiao: 金币+5, 哈，那我明天也去试一下，有结果再告诉你 2014-04-01 23:11:10
这个我当时试了的 消不消去都没什么问题  两种结果是一样的 至于多方程的情况我没试过
赞一下(1人) 回复此楼
天天多学一点
11楼2014-04-01 17:37:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

WWJ1991

木虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
参考：http://zhidao.baidu.com/link?url ... le8kHjgTQAnHcUx-EMa。
百度知道给的是一阶的求法，你换成二阶就行了。可以将每一阶微分的分量写在一起，组成矩阵形式。
赞一下 回复此楼
2楼2014-03-25 16:48:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

daiyouming

新虫 (初入文坛)

【答案】应助回帖

感谢参与，应助指数 +1
引用回帖:
2楼: Originally posted by WWJ1991 at 2014-03-25 16:48:07
参考：http://zhidao.baidu.com/link?url=n-6qzdYoFHcIxmVsFUFbKrknOkbLdZstCBx1IgbFeDhTH3rdU8gqbP1xQ5hPe-FAZ6Wle8kHjgTQAnHcUx-EMa。
百度知道给的是一阶的求法，你换成二阶就行了。可以将每一阶微分的分量写在 ...

您好，请问你的意思是变成24X24矩阵吗,然后一个个一起求解，有这种参考案例吗::
赞一下 回复此楼
3楼2014-03-25 18:39:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WWJ1991

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
lylpiao: 金币+5, ★有帮助, 谢谢 2014-03-30 11:47:20
引用回帖:
3楼: Originally posted by daiyouming at 2014-03-25 18:39:05
您好，请问你的意思是变成24X24矩阵吗,然后一个个一起求解，有这种参考案例吗::...

不是，还是12X12的。二阶微分用d2y表示啊。
赞一下 回复此楼
4楼2014-03-25 20:44:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lylpiao

木虫 (正式写手)
ben_ladeng: 求助不能申请程序强帖，应助和精彩讲座可以申请程序强帖。欢迎发表高质量作品 2014-03-31 06:57:26
引用回帖:
5楼: Originally posted by george1990 at 2014-03-30 01:11:47
这种问题你可以参考机械动力学里面的4阶龙格库塔解法  它对于解这种方程很ok的

请问下，该如何求解啊，送你一朵鲜花哈
非线性二阶微分方程组求解，matlab-1
H{W}]M5NY@[C48~~W@42B}6.jpg
赞一下 回复此楼
把一切交给时间！
6楼2014-03-30 11:43:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

george1990

金虫 (小有名气)

【答案】应助回帖

function dy=lyl(t,y)
dy=zeros(4,1);
dy(1)=y(2);
dy(3)=y(4);
dy(2)=(30*(300*sin(t)-(0.5*y(2)+10*y(4)+1000*y(1)+600*y(4).*y(3)+300*y(2).*y(3))))/2999;
dy(4)=(100*cos(t)-(dy(2)+y(4)+10*y(2)+2*y(3)+3*y(2).*y(3)+20*y(1).*y(3)))/30;


close all;clc;
y0=[0 0 0 0];  %x初值和x导数的初值
tspan=[0 3];
options = odeset('RelTol',1e-8,'AbsTol',1e-2);
[t y]=ode23s(@lyl,tspan,y0,options);
%plot(t,y(:,2),'r');
我这里设置的初值是随意的 导致计算到2.7s就不能计算了 你看看你那边的初值是多少 因为没给这个 我是随便设置的
赞一下 回复此楼
天天多学一点
7楼2014-03-31 18:46:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

george1990

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lylpiao: 金币+5, ★有帮助, 谢谢~ 2014-04-01 10:37:09
lylpiao: 金币+5, ★有帮助, 谢谢~ 2014-04-01 10:40:34
额 第一段程序发错了 应该是这样的：
function dy=lyl(t,y)
dy=zeros(4,1);
dy(1)=y(2);
dy(3)=y(4);
dy(2)=(300*sin(t)-(dy(4)+0.5*y(2)+10*y(4)+1000*y(1)+600*y(4).*y(3)+300*y(2).*y(3)))/100;
dy(4)=(100*cos(t)-(dy(2)+y(4)+10*y(2)+2*y(3)+3*y(2).*y(3)+20*y(1).*y(3)))/30;
赞一下 回复此楼
天天多学一点
8楼2014-03-31 18:53:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lylpiao

木虫 (正式写手)
引用回帖:
8楼: Originally posted by george1990 at 2014-03-31 18:53:11
额 第一段程序发错了 应该是这样的：
function dy=lyl(t,y)
dy=zeros(4,1);
dy(1)=y(2);
dy(3)=y(4);
dy(2)=(300*sin(t)-(dy(4)+0.5*y(2)+10*y(4)+1000*y(1)+600*y(4).*y(3)+300*y(2).*y(3)))/100;
dy(4)=(10 ...

你好，我有个疑惑，按照你的思路，一个方程里面同时含有二阶量，dy(2)和dy(4)，这可以吗
dy(2)=(300*sin(t)-(dy(4)+0.5*y(2)+10*y(4)+1000*y(1)+600*y(4).*y(3)+300*y(2).*y(3)))/100;
dy(4)=(100*cos(t)-(dy(2)+y(4)+10*y(2)+2*y(3)+3*y(2).*y(3)+20*y(1).*y(3)))/30;
赞一下 回复此楼
把一切交给时间！
9楼2014-04-01 10:32:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lylpiao

木虫 (正式写手)
引用回帖:
9楼: Originally posted by lylpiao at 2014-04-01 10:32:53
你好，我有个疑惑，按照你的思路，一个方程里面同时含有二阶量，dy(2)和dy(4)，这可以吗
dy(2)=(300*sin(t)-(dy(4)+0.5*y(2)+10*y(4)+1000*y(1)+600*y(4).*y(3)+300*y(2).*y(3)))/100;
dy(4)=(100*cos(t)-(dy(2) ...

查阅相关资料：若两个高阶微分方程同时含有最高价导数项，需要先进行相应处理，先消去某个高阶导数。但是方程组少的话计算还行，多的话消去就太麻烦了。所以还是不知道求解。不知道我理解的对不对呢
赞一下 回复此楼
把一切交给时间！
10楼2014-04-01 11:13:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lylpiao 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 307求调剂 +4 73372112 2026-02-28 6/300 2026-03-01 00:04 by ll247
[考研] 304求调剂 +3 52hz~~ 2026-02-28 5/250 2026-03-01 00:00 by 52hz~~
[考研] 317一志愿华南理工电气工程求调剂 +5 Soliloquy_Q 2026-02-28 8/400 2026-02-28 23:36 by xyx2012xyx
[考研] 292求调剂 +3 yhk_819 2026-02-28 3/150 2026-02-28 21:57 by gaoxiaoniuma
[考研] 290求调剂 +5 材料专硕调剂； 2026-02-28 6/300 2026-02-28 21:40 by gaoxiaoniuma
[考博] 26申博 +4 想申博！ 2026-02-26 4/200 2026-02-28 21:37 by limorning
[考研] 材料类求调剂 +6 wana_kiko 2026-02-28 6/300 2026-02-28 21:20 by gaoxiaoniuma
[考研] 284求调剂 +4 天下熯 2026-02-28 4/200 2026-02-28 21:13 by gaoxiaoniuma
[考研] 高分子化学与物理调剂 +4 好好好1233 2026-02-28 7/350 2026-02-28 20:42 by 好好好1233
[教师之家] 版面费该交吗 +15 苹果在哪里 2026-02-22 18/900 2026-02-28 18:20 by mibaomingg
[考研] 285求调剂 +5 满头大汗的学生 2026-02-28 5/250 2026-02-28 18:10 by 材料专硕调剂；
[考博] 博士自荐 +3 kkluvs 2026-02-28 3/150 2026-02-28 16:59 by StarAura
[考研] 265分求调剂不调专业和学校有行学上就 +4 礼堂丁真258 2026-02-28 6/300 2026-02-28 16:18 by 求调剂zz
[考研] 0856调剂 +3 刘梦微 2026-02-28 3/150 2026-02-28 13:22 by houyaoxu
[考研] 寻找调剂 +3 LYidhsjabdj 2026-02-28 3/150 2026-02-28 12:59 by miniwendy
[硕博家园] 博士自荐 +6 科研狗111 2026-02-26 9/450 2026-02-28 12:32 by seaskyy
[考研] 272求调剂 +3 田智友 2026-02-28 3/150 2026-02-28 12:31 by 王加浩to
[考研] 298求调剂 +4 axyz3 2026-02-28 4/200 2026-02-28 11:21 by wang_dand
[基金申请] 面上可以超过30页吧？ +12 阿拉贡aragon 2026-02-22 13/650 2026-02-26 22:09 by Hahaxia
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
关闭
请填处理意见
关闭