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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 请问如何求解二元一阶微分方程组
252/3返回列表上一页123下一页
查看: 7799  |  回复: 24
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 数学EPI ，点击这里进行查看

peterflyer

木虫之王 (文学泰斗)

peterflyer

【答案】应助回帖

引用回帖:
7楼: Originally posted by xxxfield at 2012-01-04 17:10:23
对第一式关于t再求一次导数，然后将y, y'用第一、二式代入，得到一个关于x的(一元)二阶常系数线性方程，这个方程的解有公式可用，解出x后马上就可求出y了。

简单的一阶方程就这样让你搞成了复杂的二阶方程。这样一来初始边界条件也不够了，那又将如何呢？
赞一下 回复此楼
11楼2013-11-07 19:14:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer

【答案】应助回帖

引用回帖:
10楼: Originally posted by peterflyer at 2013-11-07 17:05:15
对了，Γ（2）等于1. Γ（n）=(n-1)！

续前：
X(t)=[(B*F-C*E]/[A*E-B*D]+[(B*F-C*E]/[A*E-B*D]*
            exp{[(B*D- A*E)/(B-E)]*t}

Y(t)=(C*D-A*F)(A*E-B*D)*t+(F-C)/(B-E+A*E-B*D)+
         (F-C)/(B-E+A*E-B*D)*exp{[(A*E-B*D)/(B-E)]*t}
赞一下 回复此楼
12楼2013-11-07 19:23:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer

【答案】应助回帖

楼主能给个信箱，我把手算的计算过程用手机拍成图片给你发过去。
赞一下(1人) 回复此楼
13楼2013-11-07 19:29:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer

【答案】应助回帖

sweety: 应助指数+1 2013-11-10 12:49:18
sweety: 数学EPI+1, 耐心解答 2013-11-10 12:49:37
求解过程：对方程（1）、（2）两边求Laplace变换，记X(s) 、Y(s)分别为x=x(t)和y=y(t)关于t的Laplace变换。由拉氏变换的性质，有：
s*X(s)-k=A*X(s)+B*Y(s)+C/s               (3)
s*Y(s)-k=D*X(s)+E*Y(s)+F/s               (4)
联解方程（3）和（4），有：
X(s)=[k*s+(B*k-E*k+C)+(B*F-C*E)/s]/{[s-(A+E)/2]^2-[(A-E)^2/4
             +B*D]}
      =[k*s+(B*k-E*k+C)]/{[s-(A+E)/2]^2-[(A-E)^2/4+B*D]}
      +(B*F-C*E)/(A*E-B*D)*1/s
      +{[(C*E-B*F)/(A*E-B*D)]*s+(B*F-C*E)*(A+E)/(A*E-B*D)}/
         /{[s-(A+E)/2]^2-[(A-E)^2/4+B*D]}
Y(s)=[k*s+(F+D*k-A*k)+(C*D-A*F)/s]/{[s-(A+E)/2]^2-[(A-E)^2/4
             +B*D]}
      =[k*s+(F+D*k-A*k)]/{[s-(A+E)/2]^2-[(A-E)^2/4+B*D]}
       +[(C*D-A*F)/(A*E-B*D)]*1/s
      +{[A*F-C*D)/(A*E-B*D)]*s+(A+E)*(C*D-A*F)/(A*E-B*D)}/
          /{[s-(A+E)/2]^2-[(A-E)^2/4+B*D]}
对X(s)和Y(s)求拉氏反变换，得到x=x(t)和y=y(t)的表达式。查拉氏反变换表可知：1/s的拉氏反变换为1；s/[s^2+ω^2]的反变换为Cosωt;
1/[s^2+ω^2]的反变换为Sinωt; 1/[(s-δ)^2+ω^2]的反变换为
  exp(at)*Cosωt;  (s-δ)/[(s-δ)^2+ω^2]的反变换为exp(at)*Sinωt;
1/(s-a)的反变换为exp(at),同时并注意到拉氏变换与反变换均具有线性叠加的性质，故可得到如下结果：
(1) 若-B*D-(A-E)^2/4 ≥0, 令ω^2=-B*D-(A-E)^2/4 ,此处ω≥0
则上面的X(s)和Y(s)的表达式均可通过代数中的知识表达为由1/s、(s-γ)/[(s-γ)^2+ω^2]以及1/[(s-γ)^2+ω^2]等的线性表达式。因此， x(t)与y(t)均为1、exp(γ*t)*Cosωt、exp(γ*t)*Sinωt等的线性表达式，只不过各自的系数不同而已。由于表达式过于复杂，这里就不具体写出了。
（2）若-B*D-(A-E)^2/4 ≤0，令-ω^2=-B*D-(A-E)^2/4 ,此处ω≥0
此时，X(s)和Y(s)均可表达为1/s、1/(s-γ)、1/(s+γ)等的线性表达式，因此， x(t)与y(t)均为1、exp(γ*t）、exp(-γ*t)等的线性表达式，只不过各自的系数不同而已。由于表达式过于复杂，这里就不具体写出了。
赞一下(5人) 回复此楼
14楼2013-11-10 11:38:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

优雅的狮子

铜虫 (初入文坛)
牛，赞一个
回复此楼
www.tutu168.com
15楼2014-05-05 23:59:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

量子若水

银虫 (小有名气)
引用回帖:
14楼: Originally posted by peterflyer at 2013-11-10 11:38:50
求解过程：对方程（1）、（2）两边求Laplace变换，记X(s) 、Y(s)分别为x=x(t)和y=y(t)关于t的Laplace变换。由拉氏变换的性质，有：
s*X(s)-k=A*X(s)+B*Y(s)+C/s               (3)
s*Y(s)-k=D*X(s)+E*Y(s)+F/s    ...

这是用的变分法么？
赞一下 回复此楼
一切灵感，一切信仰，一切时代的结晶，一切人类天才的光华，都将随着太阳系的崩溃而毁灭，人类全部的成就神殿将不可避免的被埋藏在宇宙的废墟之中，所有这一切，
16楼2014-05-13 00:12:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer
引用回帖:
16楼: Originally posted by 量子若水 at 2014-05-13 00:12:39
这是用的变分法么？...

不是变分法，是拉普拉斯变换法。
赞一下 回复此楼
17楼2014-05-13 08:05:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

量子若水

银虫 (小有名气)
引用回帖:
17楼: Originally posted by peterflyer at 2014-05-13 08:05:28
不是变分法，是拉普拉斯变换法。...

虽然是变换到了s空间，但方程的解析式还是没有得出来。是不是变过去再转回来更麻烦了。

[ 发自小木虫客户端 ]
赞一下 回复此楼
一切灵感，一切信仰，一切时代的结晶，一切人类天才的光华，都将随着太阳系的崩溃而毁灭，人类全部的成就神殿将不可避免的被埋藏在宇宙的废墟之中，所有这一切，
18楼2014-05-13 15:14:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer
引用回帖:
18楼: Originally posted by 量子若水 at 2014-05-13 15:14:52
虽然是变换到了s空间，但方程的解析式还是没有得出来。是不是变过去再转回来更麻烦了。
...

对于像本题这样的情形，方程的解析式很容易求出来的，我已经写得比较明白了。当然，对于其他复杂的情形，拉氏变换和其逆变换会存在很大的难度的。这就是今天各种解析法发展不快而数值解法发展迅速的原因。
赞一下 回复此楼
19楼2014-05-13 16:35:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

量子若水

银虫 (小有名气)
引用回帖:
19楼: Originally posted by peterflyer at 2014-05-13 16:35:11
对于像本题这样的情形，方程的解析式很容易求出来的，我已经写得比较明白了。当然，对于其他复杂的情形，拉氏变换和其逆变换会存在很大的难度的。这就是今天各种解析法发展不快而数值解法发展迅速的原因。...

嗯！很佩服你的数学，只是很多的物理问题都用数值求解，失去了原来解析解的那种简洁，解析解中包含了很多的数学思想，数值求解的那种无限取近似的我一开始接触的时候真有点接受不了。

[ 发自小木虫客户端 ]
赞一下 回复此楼
一切灵感，一切信仰，一切时代的结晶，一切人类天才的光华，都将随着太阳系的崩溃而毁灭，人类全部的成就神殿将不可避免的被埋藏在宇宙的废墟之中，所有这一切，
20楼2014-05-13 19:01:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zgchen9 的主题更新
252/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +37 医学老男孩 2026-08-13 86/4300 2026-08-18 00:18 by longsantaizi
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 时间戳变了，能看出什么问题？ +4 基诺咪客 2026-08-17 4/200 2026-08-17 16:10 by Vivilian
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
信息提示
关闭
请填处理意见
关闭