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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 非线性二阶微分方程组求解，matlab
131/2返回列表12下一页
查看: 4713  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 301求调剂 +10 yy要上岸呀 2026-03-17 10/500 2026-03-21 03:14 by JourneyLucky
[考研] 303求调剂 +5 睿08 2026-03-17 7/350 2026-03-21 03:11 by JourneyLucky
[考研] 299求调剂 +6 △小透明* 2026-03-17 6/300 2026-03-21 02:42 by JourneyLucky
[考研] 083200学硕321分一志愿暨南大学求调剂 +3 innocenceF 2026-03-17 3/150 2026-03-21 02:35 by JourneyLucky
[考研] 265求调剂 +9 梁梁校校 2026-03-17 9/450 2026-03-21 02:17 by JourneyLucky
[考研] 一志愿西南交大，求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 311求调剂 +5 冬十三 2026-03-18 5/250 2026-03-21 00:16 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 求调剂，一志愿:南京航空航天大学大学 ，080500材料科学与工程学硕，总分289分 +4 @taotao 2026-03-19 4/200 2026-03-20 22:14 by JourneyLucky
[考研] 材料与化工 322求调剂 +4 然11 2026-03-19 4/200 2026-03-20 22:12 by luoyongfeng
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 5/250 2026-03-20 22:11 by 云游重阳
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 材料学求调剂 +4 Stella_Yao 2026-03-20 4/200 2026-03-20 20:28 by ms629
[考研] 0703化学调剂 ，六级已过，有科研经历 +13 曦熙兮 2026-03-15 13/650 2026-03-20 19:35 by Dream007008
[考研] 286分人工智能专业请求调剂愿意跨考！ +3 lemonzzn 2026-03-17 4/200 2026-03-20 11:04 by lemonzzn
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +5 枫桥ZL 2026-03-18 7/350 2026-03-19 14:52 by 功夫疯狂
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 321求调剂 +5 大米饭！ 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
信息提示
关闭
请填处理意见
关闭