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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 非线性二阶微分方程组求解，matlab
131/2返回列表12下一页
查看: 4714  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 332求调剂 +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[考研] 297求调剂 +9 戏精丹丹丹 2026-03-17 9/450 2026-03-21 01:49 by JourneyLucky
[考研] 华东师范大学-071000生物学-293分-求调剂 +3 研究生何瑶明 2026-03-18 3/150 2026-03-21 01:30 by JourneyLucky
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 323求调剂 +3 洼小桶 2026-03-18 3/150 2026-03-20 22:54 by JourneyLucky
[考研] 一志愿苏州大学材料求调剂，总分315（英一） +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 一志愿 南京航空航天大学大学 ，080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考研] 0703化学调剂 ，六级已过，有科研经历 +13 曦熙兮 2026-03-15 13/650 2026-03-20 19:35 by Dream007008
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 环境工程调剂 +9 大可digkids 2026-03-16 9/450 2026-03-20 17:38 by 醉在风里
[考研] 0703化学调剂 +10 妮妮ninicgb 2026-03-15 14/700 2026-03-19 22:59 by 学员8dgXkO
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 0817调剂 +3 没有答案_ 2026-03-14 3/150 2026-03-19 09:51 by Xu de nuo
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 一志愿，福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 321求调剂 +5 大米饭！ 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
信息提示
关闭
请填处理意见
关闭