24小时热门版块排行榜    

查看: 4590  |  回复: 18

伊景曦zyx

木虫 (小有名气)

[交流] 二阶微分方程组化为一阶微分方程组的方法

有哪位大神可以帮我看一下,附件中的二阶微分方程组(相互之间有关系)在利用龙格库塔法求解时,如何化为一阶微分方程组???我尝试过多种化解办法,但最终的计算结果都是错的,深深的在怀疑自己的变化方法是不是错误的???
         比较着急,希望大神们可以帮帮我,本人将不胜感激
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 二阶微分方程组.docx
  • 2017-02-13 10:25:18, 14.64 K

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

永远不要放弃真正想要的东西,殊不知:等待虽难,但后悔更甚!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

伊景曦zyx

木虫 (小有名气)

引用回帖:
10楼: Originally posted by WanderingHeart at 2017-02-13 21:11:55
本来就是未知待求的函数,不太懂你说没有定义过是什么意思……...

具体见附件中PPT,非常感谢您

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 2017.02.13_-_副本.pptx
  • 2017-02-13 21:26:30, 55.6 K
永远不要放弃真正想要的东西,殊不知:等待虽难,但后悔更甚!
11楼2017-02-13 21:26:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

伊景曦zyx

木虫 (小有名气)

引用回帖:
2楼: Originally posted by WanderingHeart at 2017-02-13 11:52:45
应该是化成四个方程的方程组咯?

恩,是的,是应该化解成为四个一阶微分方程。但是因为原来的两个微分方程之间是有关系的,导致化解的时候,第一个二阶变为的一阶要代入第二个原始的二阶微分方程,这样会对结果产生影响吧?
永远不要放弃真正想要的东西,殊不知:等待虽难,但后悔更甚!
3楼2017-02-13 15:06:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WanderingHeart

铁杆木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
不懂你为什么那么搞……

其实这样似乎就可以啊……

引入两个新的未知变量y1=dx1/dt, y2=dx2/dt,这两个就是两个一阶方程,然后原来的两个方程利用新引入的变量可以化简成一阶的。就完事儿了呀。
6楼2017-02-13 16:52:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

WanderingHeart

铁杆木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
应该是化成四个方程的方程组咯?
2楼2017-02-13 11:52:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WanderingHeart

铁杆木虫 (著名写手)

★ ★ ★ ★ ★
伊景曦zyx: 金币+5 2017-02-13 16:41:48
你是怎么化的?

发自小木虫Android客户端
4楼2017-02-13 16:05:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

伊景曦zyx

木虫 (小有名气)

引用回帖:
4楼: Originally posted by WanderingHeart at 2017-02-13 16:05:40
你是怎么化的?

我的具体化解思路请看附件

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 我的化解思路.pptx
  • 2017-02-13 16:29:53, 90.16 K
永远不要放弃真正想要的东西,殊不知:等待虽难,但后悔更甚!
5楼2017-02-13 16:30:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pdl9527

专家顾问 (小有名气)

★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
伊景曦zyx: 金币+5 2017-02-13 21:08:05
根据你的问题改变U, f, L y0, tspan的值,看看得到的结果是不是你想要的,方程的处理你可以doc ode23里面说的很详细了。
CODE:
function question_24
clear;clc
%2017-2-13
U=100;
f=[1 4];
L=[1 2];
y0=[1 2 3 4];
tspan=[0 5];
[t,y]=ode23(@fun,tspan,y0);
plot(t,y(:,1),'-o',t,y(:,2),'-.')
function dydt=fun(t,y)
dydt=[y(1);
    y(2);
    1/2*(U-(2*y(3)-y(4)+4*pi^2*(2*f(1)^2*L(1)*y(1)-f(2)^2*L(2)^2*y(2))));
    1/2*(U-2*y(4)-4*pi^2*f(2)^2*L(2)*y(2))];
end
end

7楼2017-02-13 18:55:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

伊景曦zyx

木虫 (小有名气)

引用回帖:
6楼: Originally posted by WanderingHeart at 2017-02-13 16:52:34
不懂你为什么那么搞……

其实这样似乎就可以啊……

引入两个新的未知变量y1=dx1/dt, y2=dx2/dt,这两个就是两个一阶方程,然后原来的两个方程利用新引入的变量可以化简成一阶的。就完事儿了呀。

如果是两个互相没有关系的微分方程的话,你的思路是完全没有问题的。但是关键是现在两个二阶微分方程之间有关系,就这样单纯的计算的话,化解为的第三个一阶微分方程中就会包含dy2,这样就会导致计算机认为这是一个没有定义过的未知量。
永远不要放弃真正想要的东西,殊不知:等待虽难,但后悔更甚!
8楼2017-02-13 20:59:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

伊景曦zyx

木虫 (小有名气)

引用回帖:
7楼: Originally posted by pdl9527 at 2017-02-13 18:55:16
根据你的问题改变U, f, L y0, tspan的值,看看得到的结果是不是你想要的,方程的处理你可以doc ode23里面说的很详细了。


function question_24
clear;clc
%2017-2-13
U=100;
f=;
L=;
y0=;
tspan=;
=od ...

请问您在二阶变为一阶时,是怎么对方程进行化解的?
永远不要放弃真正想要的东西,殊不知:等待虽难,但后悔更甚!
9楼2017-02-13 21:07:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WanderingHeart

铁杆木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
8楼: Originally posted by 伊景曦zyx at 2017-02-13 20:59:52
如果是两个互相没有关系的微分方程的话,你的思路是完全没有问题的。但是关键是现在两个二阶微分方程之间有关系,就这样单纯的计算的话,化解为的第三个一阶微分方程中就会包含dy2,这样就会导致计算机认为这是一个 ...

本来就是未知待求的函数,不太懂你说没有定义过是什么意思……
10楼2017-02-13 21:11:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 伊景曦zyx 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 328求调剂,英语六级551,有科研经历 +6 生物工程调剂 2026-03-17 10/500 2026-03-22 20:22 by edmund7
[考研] 考研调剂 +3 呼呼?~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
[考研] 一志愿深大,0703化学,总分302,求调剂 +4 七月-七七 2026-03-21 4/200 2026-03-21 18:20 by 学员8dgXkO
[考研] 材料 271求调剂 +5 展信悦_ 2026-03-21 5/250 2026-03-21 17:29 by 学员8dgXkO
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 22 350 本科985求调剂,求老登收留 +3 李轶男003 2026-03-20 3/150 2026-03-21 13:28 by 搏击518
[考研] 307求调剂 +3 wyyyqx 2026-03-17 3/150 2026-03-21 03:20 by JourneyLucky
[考研] 二本跨考郑大材料306英一数二 +3 z1z2z3879 2026-03-17 3/150 2026-03-21 02:29 by JourneyLucky
[考研] 华东师范大学-071000生物学-293分-求调剂 +3 研究生何瑶明 2026-03-18 3/150 2026-03-21 01:30 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-18 3/150 2026-03-21 00:46 by JourneyLucky
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 08工学调剂 +5 用户573181 2026-03-20 5/250 2026-03-20 15:47 by xia_2003
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 286求调剂 +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
信息提示
请填处理意见