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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 理论与计算方法 » 二阶微分方程组化为一阶微分方程组的方法
192/2返回列表上一页12
查看: 4590  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

伊景曦zyx

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

具体见附件中PPT，非常感谢您
赞一下(2人) 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：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的回帖

pdl9527

专家顾问 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
9楼: Originally posted by 伊景曦zyx at 2017-02-13 21:07:55
请问您在二阶变为一阶时，是怎么对方程进行化解的？...

令dx1/dt=y1   dx2/dt=y2   dx1^2/dt^2=dy1/dt=y3      dx2^2/dt^2=dy2/dt=y4

要使用ODE23求解，要保证方程里面右边都小于二阶，因此需要将第二个方程的二阶微分带入到第一个方程中，求出一个二阶的关系式，反过来将第一个方程的二阶微分带入到第二个方程中求出另一个二阶微分关系式，你在MATLAB里面doc ode23里面有比较简单例子，方法是一样的。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

伊景曦zyx
12楼2017-02-13 23:13:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WanderingHeart

铁杆木虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
11楼: Originally posted by 伊景曦zyx at 2017-02-13 21:26:44
具体见附件中PPT，非常感谢您...

意思是说matlab求解时候必须要求每个方程是dy/dt=f(x)这种形式了？
那你从两个方程里面把d2x1/dt2和d2x2/dt2当做未知数求解一下，相当于把原方程组变形为d2x1/dt2=……和d2x2/dt2=……，然后在变化成四个一阶方程组。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

伊景曦zyx
13楼2017-02-14 09:22:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

伊景曦zyx

木虫 (小有名气)
送红花一朵
引用回帖:
12楼: Originally posted by pdl9527 at 2017-02-13 23:13:58
令dx1/dt=y1   dx2/dt=y2   dx1^2/dt^2=dy1/dt=y3      dx2^2/dt^2=dy2/dt=y4

要使用ODE23求解，要保证方程里面右边都小于二阶，因此需要将第二个方程的二阶微分带入到第一个方程中，求出一个二阶的关系式，反过 ...

我明白您的意思了，非常感谢！我试试看，其实这样对方程的处理方式实质上仍然是交叉化解的，不知道这样的化为一阶微分方程组的方法对计算结果有没有影响？
赞一下 回复此楼
永远不要放弃真正想要的东西，殊不知：等待虽难，但后悔更甚！
14楼2017-02-14 09:41:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

伊景曦zyx

木虫 (小有名气)
送红花一朵
引用回帖:
13楼: Originally posted by WanderingHeart at 2017-02-14 09:22:27
意思是说matlab求解时候必须要求每个方程是dy/dt=f(x)这种形式了？
那你从两个方程里面把d2x1/dt2和d2x2/dt2当做未知数求解一下，相当于把原方程组变形为d2x1/dt2=……和d2x2/dt2=……，然后在变化成四个一阶方程 ...

恩，对的！问题就在这里，我之前的化解思路就是这样的，但是不同的化解办法得到的结果竟然是不一样的，所以我就在想：这样的交叉化解办法是否正确？  龙格库塔法在求解的时候，会不会因为方程的微小变化而导致结果的巨大变化呢？
赞一下 回复此楼
永远不要放弃真正想要的东西，殊不知：等待虽难，但后悔更甚！
15楼2017-02-14 09:47:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WanderingHeart

铁杆木虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
15楼: Originally posted by 伊景曦zyx at 2017-02-14 09:47:05
恩，对的！问题就在这里，我之前的化解思路就是这样的，但是不同的化解办法得到的结果竟然是不一样的，所以我就在想：这样的交叉化解办法是否正确？  龙格库塔法在求解的时候，会不会因为方程的微小变化而导致结果 ...

你这方程看起来也不是很复杂，应该不至于这样就引起数值上大的差异。
赞一下 回复此楼
16楼2017-02-14 09:58:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer
★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
伊景曦zyx: 金币+3 2017-02-14 11:01:06
刚刚看了楼主的附件中的二阶线性微分方程组。感觉完全可以用拉式变换求得解析解的，为何一定要用数值方法求数值解呢？有具体的边界条件吗？
赞一下 回复此楼
17楼2017-02-14 10:19:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

伊景曦zyx

木虫 (小有名气)
引用回帖:
16楼: Originally posted by WanderingHeart at 2017-02-14 09:58:18
你这方程看起来也不是很复杂，应该不至于这样就引起数值上大的差异。...

恩，是的。不同化解办法得到的解不同，之前的附件中【我的化解思路】有附过图
赞一下 回复此楼
永远不要放弃真正想要的东西，殊不知：等待虽难，但后悔更甚！
18楼2017-02-14 10:54:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

伊景曦zyx

木虫 (小有名气)
引用回帖:
17楼: Originally posted by peterflyer at 2017-02-14 10:19:45
刚刚看了楼主的附件中的二阶线性微分方程组。感觉完全可以用拉式变换求得解析解的，为何一定要用数值方法求数值解呢？有具体的边界条件吗？

之前计算的是一个二阶微分方程，并且是非线性的，只能用数值求解，所以一直在用龙格库塔法。后来扩充为两个二阶微分方程，先用线性微分方程组调一下程序，紧接着会加入非线性项，所以现在用数值方法求解。
给定的初始值
赞一下 回复此楼
永远不要放弃真正想要的东西，殊不知：等待虽难，但后悔更甚！
19楼2017-02-14 11:00:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 伊景曦zyx 的主题更新
192/2返回列表上一页12
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 310求调剂 +4 baibai1314 2026-03-16 4/200 2026-03-22 20:19 by edmund7
[考研] 材料与化工085600，总分304，本科有两篇sci参与，求调剂 +4 幸运的酱酱 2026-03-22 5/250 2026-03-22 20:15 by edmund7
[考研] 289材料与化工（085600）B区求调剂 +3 这么名字咋样 2026-03-22 4/200 2026-03-22 17:56 by 云民大李老师
[考研] 319求调剂 +4 小力气珂珂 2026-03-20 4/200 2026-03-22 15:53 by ColorlessPI
[考研] 269专硕求调剂 +6 金恩贝 2026-03-21 6/300 2026-03-22 14:31 by ColorlessPI
[考研] 303求调剂 +5 安忆灵 2026-03-22 6/300 2026-03-22 12:46 by 素颜倾城1988
[考研] 工科0856求调剂 +3 沐析汀汀 2026-03-21 3/150 2026-03-21 18:30 by 学员8dgXkO
[考研] 求助 +5 梦里的无言 2026-03-21 6/300 2026-03-21 17:51 by 学员8dgXkO
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 299求调剂 +5 shxchem 2026-03-20 7/350 2026-03-21 17:09 by ColorlessPI
[考研] 299求调剂 +4 某某某某位 2026-03-21 4/200 2026-03-21 16:30 by barlinike
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
[考研] 一志愿南昌大学，327分，材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 材料学求调剂 +4 Stella_Yao 2026-03-20 4/200 2026-03-20 20:28 by ms629
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
信息提示
关闭
请填处理意见
关闭