版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

本帖产生 1 个数学EPI ，点击这里进行查看

zgchen9

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1175.8
散金: 1
帖子: 113
在线: 79小时
虫号: 288553
注册: 2006-10-21
性别: GG
专业: 无机材料化学

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

各位朋友新年好，我需要解一个二元一阶微分方程组，但是本人数学水平有限，特请交各位朋友。
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. 谢谢。

方程组

回复此楼

» 猜你喜欢

2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复
孩子确诊有中度注意力缺陷已经有6人回复
2026博士申请-功能高分子，水凝胶方向已经有6人回复
论文投稿，期刊推荐已经有4人回复
硕士和导师闹得不愉快已经有13人回复
请问2026国家基金面上项目会启动申2停1吗已经有5人回复
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？已经有3人回复

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

matlab数值求解边界条件微分方程组已经有7人回复
一个微分方程组，求解已经有18人回复
LABVIEW四阶龙格库塔法求解一阶微分方程组已经有12人回复
【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题已经有21人回复
【求助】用mathematica 5.0求解一个非线性方程组失败，特发帖求助！已经有5人回复
【求助】这个微分方程组simulink求解到底错在哪里？已经有3人回复
【求助】请教非齐次常微分方程组的解析解法已经有4人回复
【求助】积分微分方程matlab求解已经有6人回复
【求助】急请微分方程高手看下这个方程组！已经有14人回复
【求助】matlab怎么求解偏微分方程组啊，先谢谢了已经有13人回复
【求助】有限元法求解二阶偏微分方程组已经有6人回复
【求助】如何用Runge-Kutta迭代求解二阶常微分方程组【已解决】已经有9人回复

一份耕耘，一份收获

1楼2012-01-02 21:31:25

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

回帖支持 ( 显示支持度最高的前 50 名 )

peterflyer

木虫之王 (文学泰斗)

peterflyer

数学EPI: 10
应助: 20282 (院士)
金币: 145640
红花: 1374
帖子: 93062
在线: 7693.7小时
虫号: 1482829
注册: 2011-11-08
性别: GG
专业: 功能陶瓷

【答案】应助回帖

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

数学EPI: 10
应助: 20282 (院士)
金币: 145640
红花: 1374
帖子: 93062
在线: 7693.7小时
虫号: 1482829
注册: 2011-11-08
性别: GG
专业: 功能陶瓷

【答案】应助回帖

楼主能给个信箱，我把手算的计算过程用手机拍成图片给你发过去。

赞一下(1人)

回复此楼

13楼2013-11-07 19:29:05

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

普通回帖

一山

铁杆木虫 (正式写手)

应助: 17 (小学生)
金币: 5709
散金: 15
红花: 3
帖子: 988
在线: 287.6小时
虫号: 92730
注册: 2005-09-10
专业: 光学

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
soliton923(金币+2): 谢谢参与讨论~~~ 2012-01-02 22:31:57
zgchen9(金币+5): ★★★很有帮助 2012-01-02 22:34:49
zgchen9(金币+5): ★★★很有帮助 2012-01-02 22:35:21
zgchen9(金币+5): ★★★很有帮助 2012-01-05 20:58:34

用mathematica软件可以求解：
DSolve[{x'[t] == A x[t] + B y[t] + C, y'[t] == D x[t] + E y[t] + F,
x[0] == K, y[0] == K}, {x, y}, t] // Simplify

答案是：

x -> Function[{t}, (2 E^(-(1/
      2) (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) (-2 B C D E^(
      1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) +
      2 B C D E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) +
      2 B C D E^(
      Sqrt[A^2 + 4 B D - 2 A E + E^2] t +
         1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) -
      A C E^(1 + Sqrt[A^2 + 4 B D - 2 A E + E^2] t +
         1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) +
      C E^(2 + Sqrt[A^2 + 4 B D - 2 A E + E^2] t +
         1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) +
      A C E^(1 + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) -
      C E^(2 + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) -
      2 B C D E^(
      1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) +
      A C E^(1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) -
      C E^(2 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) -
      A C E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) +
      C E^(2 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) +
      C E^(1 + Sqrt[A^2 + 4 B D - 2 A E + E^2] t +
         1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[
      A^2 + 4 B D - 2 A E + E^2] +
      C E^(1 + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t)
         Sqrt[A^2 + 4 B D - 2 A E + E^2] -
      C E^(1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[
      A^2 + 4 B D - 2 A E + E^2] -
      C E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[
      A^2 + 4 B D - 2 A E + E^2] +
      A B E^(1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F -
      A B E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F -
      A B E^(Sqrt[A^2 + 4 B D - 2 A E + E^2] t +
         1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F -
      B E^(1 + Sqrt[A^2 + 4 B D - 2 A E + E^2] t +
         1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F +
      B E^(1 + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F +

      A B E^(1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F +
      B E^(1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F -
      B E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F -
      B E^(1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[
      A^2 + 4 B D - 2 A E + E^2] F +
      B E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[
      A^2 + 4 B D - 2 A E + E^2] F -
      B E^(Sqrt[A^2 + 4 B D - 2 A E + E^2] t +
         1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[
      A^2 + 4 B D - 2 A E + E^2] F +
      B E^(1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[
      A^2 + 4 B D - 2 A E + E^2] F +
      A B D E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K +
      2 B^2 D E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K -
      A B D E^(
      1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K -
      2 B^2 D E^(
      1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K +
      A^2 E^(1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K +
      2 A B E^(
      1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K +
      B D E^(1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K -
      A E^(2 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K -
      A^2 E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K -
      2 A B E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K -
      B D E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K +
      A E^(2 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K +
      B D E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[
      A^2 + 4 B D - 2 A E + E^2] K +
      B D E^(1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[
      A^2 + 4 B D - 2 A E + E^2] K -
      A E^(1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[
      A^2 + 4 B D - 2 A E + E^2] K -
      A E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[
      A^2 + 4 B D - 2 A E + E^2] K))/(Sqrt[
   A^2 + 4 B D - 2 A E +
      E^2] (-A - E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) (A + E + Sqrt[
      A^2 + 4 B D - 2 A E + E^2]))],
  y -> Function[{t}, -(2 E^(-(1/
      2) (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) (-A C D E^(
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) +
      A C D E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) +
      A C D E^(
         Sqrt[A^2 + 4 B D - 2 A E + E^2] t +
         1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) +
      C D E^(1 + Sqrt[A^2 + 4 B D - 2 A E + E^2] t +
         1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) -

      C D E^(1 +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) -
      A C D E^(
         1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) -
      C D E^(1 +
         1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) +
      C D E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) +
      C D E^(1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t)
         Sqrt[A^2 + 4 B D - 2 A E + E^2] -
      C D E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[
         A^2 + 4 B D - 2 A E + E^2] +
      C D E^(Sqrt[A^2 + 4 B D - 2 A E + E^2] t +
         1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[
         A^2 + 4 B D - 2 A E + E^2] -
      C D E^(1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[
         A^2 + 4 B D - 2 A E + E^2] +
      A^2 E^(1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F +
      2 B D E^(1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t)
         F -
      A^2 E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F -
      2 B D E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F -
      A^2 E^(Sqrt[A^2 + 4 B D - 2 A E + E^2] t +
         1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F -
      2 B D E^(
         Sqrt[A^2 + 4 B D - 2 A E + E^2] t +
         1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F +
      A E^(1 + Sqrt[A^2 + 4 B D - 2 A E + E^2] t +
         1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F -
      A E^(1 + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t)
         F + A^2 E^(
         1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F +
      2 B D E^(
         1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F -
      A E^(1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F +
      A E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F -
      A E^(1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[
         A^2 + 4 B D - 2 A E + E^2] F +
      A E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[
         A^2 + 4 B D - 2 A E + E^2] F -
      A E^(Sqrt[A^2 + 4 B D - 2 A E + E^2] t +
         1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[
         A^2 + 4 B D - 2 A E + E^2] F +
      A E^(1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[
         A^2 + 4 B D - 2 A E + E^2] F +
      A B D E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K -
      2 B D^2 E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K -
      A B D E^(
         1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K +
      2 B D^2 E^(
         1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K +
      A^2 E^(1 +
         1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K -
      2 A D E^(
         1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K +
      B D E^(1 +
         1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K -
      A E^(2 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K -
      A^2 E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K +
      2 A D E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t)
         K - B D E^(
         1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K +
      A E^(2 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K -
      B D E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[
         A^2 + 4 B D - 2 A E + E^2] K -
      B D E^(1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[
         A^2 + 4 B D - 2 A E + E^2] K +
      A E^(1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t +
         1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[
         A^2 + 4 B D - 2 A E + E^2] K +

      A E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[
         A^2 + 4 B D - 2 A E + E^2] K))/(Sqrt[
      A^2 + 4 B D - 2 A E +
      E^2] (-A - E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) (A + E +
      Sqrt[A^2 + 4 B D - 2 A E + E^2]))]}}

赞一下(1人)

回复此楼

忽悠王之俗家弟子

2楼2012-01-02 22:28:24

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

soliton923

铁杆木虫 (职业作家)

数学村村长

数学EPI: 1
应助: 36 (小学生)
贵宾: 0.565
金币: 4146.4
散金: 5275
红花: 72
沙发: 21
帖子: 3489
在线: 1501.2小时
虫号: 1005450
注册: 2010-04-25
性别: GG
专业: 数学物理

看图吧，希望对你有用

赞一下

回复此楼

soliton;sato-theory;algebre-geometry；Random-Matrices-Theory; Riemann-Hilbert method

3楼2012-01-02 22:29:45

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lilac_c

至尊木虫 (知名作家)

应助: 26 (小学生)
金币: 20544.9
散金: 4751
红花: 29
帖子: 7413
在线: 635.8小时
虫号: 1314366
注册: 2011-06-03
性别: GG
专业: 无机材料化学

【答案】应助回帖

感谢参与，应助指数 +1
zgchen9(金币+2): ★有帮助 2012-01-03 19:56:55

我们一般不直接求出解析解，太难了．

　都采用数值格式求解

离散，．．．．．．．龙格库塔

赞一下

回复此楼

我生活在一个经常爆发地震的年代

4楼2012-01-03 09:49:02

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

duotojh

5楼2012-01-03 15:34:35

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lilac_c

至尊木虫 (知名作家)

应助: 26 (小学生)
金币: 20544.9
散金: 4751
红花: 29
帖子: 7413
在线: 635.8小时
虫号: 1314366
注册: 2011-06-03
性别: GG
专业: 无机材料化学

【答案】应助回帖

★
soliton923(金币+1): 谢谢参与讨论~~~ 2012-01-04 22:20:52

dx/dt=Ax+By+C
dy/dt=Dx+Ey+F

基于　crank-Nicloson格式的求解程序编写思路

第一步
显示求解出　当前的x(n+1)与y(n+1)
以第一个方程为例：
x(n+1)=x(n)+dt*( A*x(n)+B*y(n)+C)

第二部校正开始
由于用显格式求解误差会越来越大，故此，要用
x(n+1)=x(n)+dt*( A/2*(x(n)+x(n+1))+B*(y(n)+y(n+1))/2. +c)
当两次求解误差在设置的误差范围内，结束迭代过程，否则，重新回到校正．

．．．．．．．．．．．．．．．．．

赞一下(1人)

回复此楼

我生活在一个经常爆发地震的年代

6楼2012-01-03 21:12:14

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

xxxfield

银虫 (小有名气)

数学EPI: 1
应助: 57 (初中生)
金币: 667.5
散金: 50
红花: 2
帖子: 230
在线: 53.5小时
虫号: 1496745
注册: 2011-11-17
专业: 函数论

【答案】应助回帖

★
感谢参与，应助指数 +1
soliton923(金币+1): 谢谢参与讨论~~~ 2012-01-04 22:21:03

对第一式关于t再求一次导数，然后将y, y'用第一、二式代入，得到一个关于x的(一元)二阶常系数线性方程，这个方程的解有公式可用，解出x后马上就可求出y了。

赞一下(1人)

回复此楼

7楼2012-01-04 17:10:23

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer

数学EPI: 10
应助: 20282 (院士)
金币: 145640
红花: 1374
帖子: 93062
在线: 7693.7小时
虫号: 1482829
注册: 2011-11-08
性别: GG
专业: 功能陶瓷

【答案】应助回帖

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