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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 请问如何求解二元一阶微分方程组
121/2返回列表12下一页
查看: 6812  |  回复: 24
查看全部回帖@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 数学EPI ，点击这里进行查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zgchen9

金虫 (小有名气)

[求助] 请问如何求解二元一阶微分方程组

各位朋友新年好，我需要解一个二元一阶微分方程组，但是本人数学水平有限，特请交各位朋友。
A,B,C,D,E,F,K为常数，x 和y为t 的函数。dx/dt和dy/dt为导数，二元一阶微分方程组如下：
dx/dt=Ax+By+C
dy/dt=Dx+Ey+F
边界条件为t=0时，x=y=K.


请问如何得到x 和y. 谢谢。

方程组
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
一份耕耘，一份收获
1楼2012-01-02 21:31:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer

【答案】应助回帖

设X(t)和Y(t)的Laplace变换分别记为F1（s）和F2（s）。分别对两个方程的两边取Laplace变换，得：
  s*F1(s)=A*F1(s)+B*F2(s)+C/s
   s*F2(s)=D*F1(s)+E*F2(s)+F/s
联解以上二元一次方程求出F1(s)和F2(s)
F1(s)=(B*F-C*E)/{s*[(B-E)*s+a*e-b*d]}
       =[(B*F-C*E]/[A*E-B*D]/s+[(B*F-C*E]/[A*E-B*D]/{s-[B*D-A*E]/(B-E)}
F2(s)=[(F-C)*s+C*D-A*F]/{s^2*[(B-E)*s+a*e-b*d]}
       =(C*D-A*F)(A*E-B*D)/s^2+(F-C)/(B-E+A*E-B*D)/s+(F-C)/(B-E+A*E-
        B*D)/[s-(A*E-B*D)/(B-E)]
对F1(s)和F2(s)分别求反变换，并注意到变换与反变换具有线性叠加性质，且反变换公式：1/s的反变换为1；1/ (s-a)的反变换为exp(at),1/s^2的反变换为t/Γ(2),由此便可求的X(t)和Y(t)。

完毕。
赞一下 回复此楼
8楼2013-11-07 16:53:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer

【答案】应助回帖

引用回帖:
2楼: Originally posted by 一山 at 2012-01-02 22:28:24
用mathematica软件可以求解：
DSolve // Simplify

答案是：

x -> Function,
  y -> Function}}...

那要这么麻烦？很简单的问题让楼上太忙了。呵呵。
赞一下 回复此楼
9楼2013-11-07 16:54:41
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer

【答案】应助回帖

对了，Γ（2）等于1. Γ（n）=(n-1)！
赞一下 回复此楼
10楼2013-11-07 17:05:15
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

peterflyer

木虫之王 (文学泰斗)

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

不是变分法，是拉普拉斯变换法。
赞一下 回复此楼
17楼2014-05-13 08:05:28
 已阅   回复此楼   关注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的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer
引用回帖:
20楼: Originally posted by 量子若水 at 2014-05-13 19:01:33
嗯！很佩服你的数学，只是很多的物理问题都用数值求解，失去了原来解析解的那种简洁，解析解中包含了很多的数学思想，数值求解的那种无限取近似的我一开始接触的时候真有点接受不了。
...

我们不是搞纯数学的，而是学数学为应用而已。只要能满足工程应用的计算精度，采用数值法就行了。比如现在在各行业大行其道的有限差分法、有限单元法、边界单元法等就是这样的。
赞一下 回复此楼
21楼2014-05-13 20:14:39
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zgchen9 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭