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

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (473)
>考研 (46)
>虫友互识 (37)
>硕博家园 (26)
>导师招生 (25)
>论文投稿 (8)
>休闲灌水 (6)
>考博 (2)
>航空航天 (1)
>机械 (1)
>微生物 (1)
>博后之家 (1)
>电化学 (1)
>基金申请 (1)
>找工作 (1)
>文献求助 (1)

返回列表

本帖产生 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. 谢谢。

方程组

回复此楼

» 猜你喜欢

366求调剂已经有9人回复
材料工程085601，270求调剂已经有37人回复
279学硕食品专业求调剂院校已经有18人回复
290调剂生物0860 已经有31人回复
一志愿085802 323分求调剂已经有13人回复
277求调剂已经有23人回复
322求调剂，08工科已经有4人回复
材料工程281还有调剂机会吗已经有30人回复
化学070300 求调剂已经有16人回复
化工学硕294分，求导师收留已经有12人回复

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

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 (院士)
金币: 146246
红花: 1374
帖子: 93091
在线: 7694.3小时
虫号: 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 (院士)
金币: 146246
红花: 1374
帖子: 93091
在线: 7694.3小时
虫号: 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 (小学生)
金币: 20572.5
散金: 4751
红花: 29
帖子: 7420
在线: 636.4小时
虫号: 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 (小学生)
金币: 20572.5
散金: 4751
红花: 29
帖子: 7420
在线: 636.4小时
虫号: 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 (院士)
金币: 146246
红花: 1374
帖子: 93091
在线: 7694.3小时
虫号: 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的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer

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

【答案】应助回帖

引用回帖:

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

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

【答案】应助回帖

对了，Γ（2）等于1. Γ（n）=(n-1)！

赞一下

回复此楼

10楼2013-11-07 17:05:15

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

相关版块跳转我要订阅楼主 zgchen9 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 生物学308求调剂 +5	相信必会光芒万� 2026-04-11	5/250	2026-04-12 18:14 by zhouxiaoyu
[考研] 求调剂，一志愿材料科学与工程985，365分， +8	材化李可 2026-04-11	10/500	2026-04-12 08:42 by 852137818
[考研] 347求调剂 +4	mhyqyy 2026-04-06	4/200	2026-04-12 02:27 by 秋豆菜芽
[考研] 药学305求调剂 +8	玛卡巴卡boom 2026-04-10	8/400	2026-04-12 00:07 by zhouwenxian
[考研] 086003调剂求助 +21	苏弋万 2026-04-09	22/1100	2026-04-11 20:25 by dongdian1
[考研] 085400 328分求调剂 +10	喂你一个大橙子 2026-04-09	14/700	2026-04-11 19:53 by lqspecial
[考研] 一志愿东北大学控制工程085406数二英二385，求调剂 +8	Ezra_Zhang 2026-04-09	8/400	2026-04-11 09:15 by 猪会飞
[考研] 一志愿211，化学310分，本科重点双非，求调剂 +23	努力奋斗112 2026-04-08	23/1150	2026-04-10 23:29 by 314126402
[考研] 22408 366分，本科211，一志愿西工大 +4	Rubt 2026-04-09	4/200	2026-04-10 19:51 by chemisry
[考研] 085404 298分求调剂 +10	呼啦呼啦呼呼呼 2026-04-10	11/550	2026-04-10 16:44 by wangy0907
[考研] 314求调剂 +23	wakeluofu 2026-04-09	24/1200	2026-04-10 15:31 by MOF_Catal
[考研] 江苏大学工科调剂捡漏 +3	Evan_Liu 2026-04-09	5/250	2026-04-10 10:22 by Evan_Liu
[考研] 08工学 309分求调剂 +6	Yin DY 2026-04-08	6/300	2026-04-10 09:18 by Delta2012
[考研] 085404，334分，求调剂 +5	sunjie8888 2026-04-08	8/400	2026-04-09 07:26 by sunjie8888
[考研] 327求调剂 +10	Xxjc1107. 2026-04-06	11/550	2026-04-09 01:21 by lature00
[考研] 材料调剂 +14	一样YWY 2026-04-06	14/700	2026-04-08 23:00 by 猪会飞
[考研] 293分求调剂，外语为俄语 +7	加一一九 2026-04-07	10/500	2026-04-08 20:14 by yutian743
[考研] 一志愿哈工大，初试329，求环境科学与工程调剂！ +11	余未辛 2026-04-06	11/550	2026-04-08 15:21 by screening
[考研] 化工学硕 285求调剂 +26	Wisjxn 2026-04-07	26/1300	2026-04-08 14:42 by screening
[考研] 287分求调剂有专利国奖一志愿哈工大085406 +6	白易辰 2026-04-06	7/350	2026-04-06 22:46 by 875465