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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 非线性二阶微分方程组求解，matlab
131/2返回列表12下一页
查看: 4701  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] (出分啦!!!)26 考研调剂 20w+调剂信息火速更新中！【捡漏】qq 群：1087823172 +3 db58 2026-02-28 133/6650 2026-02-28 18:25 by nbf68
[基金申请] 面上模板改不了页边距吧？ +3 ieewxg 2026-02-25 3/150 2026-02-28 18:25 by addressing
[考研] 0856材料求调剂 +9 hyf hyf hyf 2026-02-28 10/500 2026-02-28 18:22 by liqiongjy
[考研] 高分子化学与物理调剂 +4 好好好1233 2026-02-28 5/250 2026-02-28 18:20 by liqiongjy
[教师之家] 版面费该交吗 +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 材料专硕调剂；
[考研] 材料类求调剂 +4 wana_kiko 2026-02-28 4/200 2026-02-28 18:08 by djennjx
[考研] 298求调剂 +4 人间唯你是清欢 2026-02-28 4/200 2026-02-28 18:08 by njzyff
[考研] 材料调剂 +3 爱擦汗的可乐冰 2026-02-28 3/150 2026-02-28 18:06 by houyaoxu
[考研] 276求调剂 +3 路lyh123 2026-02-28 3/150 2026-02-28 18:02 by houyaoxu
[考研] 求调剂 +3 repeatt?t 2026-02-28 3/150 2026-02-28 18:00 by houyaoxu
[考研] 311求调剂 +6 南迦720 2026-02-28 6/300 2026-02-28 17:17 by 刘茉莉
[考研] 化工专硕348，一志愿985求调剂 +3 弗格个 2026-02-28 5/250 2026-02-28 17:04 by sandychj
[考博] 博士自荐 +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
[考博] 26申博 +3 想申博！ 2026-02-26 3/150 2026-02-28 16:07 by nxgogo
[考研] 290求调剂 +4 材料专硕调剂； 2026-02-28 5/250 2026-02-28 13:32 by houyaoxu
[硕博家园] 博士自荐 +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
[基金申请] 面上可以超过30页吧？ +12 阿拉贡aragon 2026-02-22 13/650 2026-02-26 22:09 by Hahaxia
信息提示
关闭
请填处理意见
关闭