差分方程与微分方程的求解

合集下载

微分方程与差分方程

微分方程与差分方程

λ = −1± i, 则齐次方程的通解为 y = e−x (C1 cos x + C2 sin x). 因 −1+ i 是单特征根,故设原非齐次方程的特解为
y* = xe−x[( A0 x + A1) cos x + (B0 x + B1) sin x].
402
把它代入原非齐次方程得
4B0 x cos x + 2(A0+B1) cos x − 4A0 x sin x + 2(B0−A1) sin x = x cos x + 3sin x,
解 将特解 y = e2x + (1+ x)ex 代入原非齐次微分方程得 (4 + 2 p + q)e2x + (3 + 2 p + q)ex + (1+ p + q)xex = rex.
比较系数,得方程组
⎧2 p + q = −4, ⎧ p = −3;
⎪⎨2 p + q − r = −3,⇒ ⎪⎨q = 2;
tan y
tan x

1 tan
y
d
tan
y
=
−∫
1 tan
x
d
tan
x,
ln(tan y) = − ln(tan x) + ln C, 故通解为 tan x tan y = C. 例3 求微方程 cos ydx + (1+ e−x ) sin ydy 在 y(0) = π 下的特解.
4
解 原方程变形为 (1+ e−x ) sin ydy = − cos ydx, 分离变量,得
过程,只要对所给通解求若干次导数,以消去所有任意常数即可.

微分方程与差分方程

微分方程与差分方程

微分方程与差分方程一、 微分方程一种机理分析方法研究两个变量之间的变化规律。

1.1.微分方程的建立变化率微元法(彭放等,《数学建模方法》,第章)1.2. 方程的求解和结果分析 1.2.1. 解析法一些常系数的或特殊函数形式的微分方程1.2.2. 数值解法大多数变系数的、非线性函数形式的微分方程一般只能求得微分方程的近似解。

给定等间距自变量点列{x n }。

1) 欧拉方法用差商代替导数,结合初始条件,推出计算{y n }的迭代公式 ⎪⎩⎪⎨⎧==00)(),(y x y y x g dx dy 的第一个方程变为))(,()()(1n n n n x y x g h x y x y ≈-+,于是 ))(,()()(1n n n n x y x hg x y x y +≈+——显式欧拉式))(,()()(111++++=n n n n x y x hg x y x y ——隐式欧拉式欧拉方法计算精度低,收敛速度慢。

))(,(2))(,(2)()(111+++++=n n n n n n x y x g h x y x g h x y x y ——梯形公式梯形公式比欧拉公式精度高,收敛速度快。

改进的欧拉方法第一步,由显式欧拉式计算1+n y 的预测值1+n y 第二步,将1+n y 代入梯形公式进行校正,即⎪⎩⎪⎨⎧++=+=++++),(2))(,(2)()())(,()(1111n n n n n n n n n n y x g h x y x g h x y x y x y x hg x y y ——改进的欧拉公式 (彭放等,《数学建模方法》,第4章)2) 龙格-库塔法(简称R-K 法)泰勒公式MATLAB 中数值求解的系统函数的实现原理就是龙格-库塔法(彭放等,《数学建模方法》,第4章)1.2.3. 图解法可以将微分方程解的全局信息直观地、形象地展现出来。

斜率场(彭放等,《数学建模方法》,第4章)1.2.4. 定性分析(徐全智等,《数学建模》,第6章)1.2.5. 稳定性分析平衡点及其稳定性的概念只针对自治方程有意义。

第5章微分方程与差分方程

第5章微分方程与差分方程

两边积分,得 故
dy = − p( x) d x , ( y ≠ 0) , y y = 0 对应于 ln | y | = − ∫ p ( x) d x + C1 , C= 。 0
y = ±e ⋅ e ∫
C1 − p( x)d x

记 C = ± eC1,得一阶齐线性方程 的通解为 y = Ce ∫
− p( x)d x
2d y = d x, 2 y −1
对上式两边积分, 对上式两边积分,得原方程的通解 y −1 ln = x + C1 。 y +1 经初等运算可得到原方程的通解为 隐函数形式
1 + Ce x y= 。 (C = ± eC1 ) 1 − Ce x 你认为做完了没有? 你认为做完了没有?
代入原方程可知: 令 y 2 − 1 = 0 ,得出 y = ±1,代入原方程可知:
5、初值条件: 给定微分方程的解所满足的条件. 初值条件: 给定微分方程的解所满足的条件. 初值问题: 求微分方程满足初始条件的解的问题. 初值问题: 求微分方程满足初始条件的解的问题.
y′ = f ( x , y ) 一阶: 一阶 y x = x0 = y 0
过定点的积分曲线; 过定点的积分曲线
dx = t2 dt
d2 y dy +b + cy = sin x 2 dx dx d x − x2 = t3 dt
2
一阶 线性 二阶 线性 一阶 非线性
微分方程的一般表示形式
n 阶微分方程的一般形式 为
F ( x, y′, y′′, L , y ( n ) ) = 0 。
dN = rN (1 例1、 ) dt N ( 0) = N 0

数的微分方程与差分方程

数的微分方程与差分方程

数的微分方程与差分方程微分方程和差分方程是数学中重要的研究对象,用于描述数学模型中的变化规律。

微分方程关注连续变化的问题,而差分方程则研究离散变化的情况。

本文将对数的微分方程和差分方程进行介绍,并比较它们之间的异同点。

一、数的微分方程微分方程是描述自变量与因变量之间的关系的方程,其基本形式为:dy/dx = f(x, y)其中dy/dx表示y对x的导数,f(x, y)是给定函数。

微分方程可分为常微分方程和偏微分方程。

常微分方程中只涉及一个自变量,而偏微分方程中涉及多个自变量。

解微分方程的方法有解析解和数值解两种。

解析解是指通过变量分离、恰当的变量换元等方法得到的精确解,而数值解则是利用数值方法进行近似计算得到的解。

二、数的差分方程差分方程是用差商表示的离散形式的方程,其基本形式为:Δy/Δx = f(x, y)其中Δy/Δx表示y对x的差商,f(x, y)是给定的函数。

差分方程可以用于描述离散的时间序列或空间序列中的变化规律。

与微分方程类似,差分方程也可分为常差分方程和偏差分方程。

解差分方程的方法主要有迭代法、插值法和递推法等。

通过这些方法,我们可以逐步逼近差分方程的解。

三、微分方程与差分方程的联系与区别微分方程和差分方程有很多共同之处,同时也存在一些区别。

首先,微分方程和差分方程都是用来描述变化规律的数学工具,它们都需要给定的函数和初始条件。

而微分方程描述的是连续变化,差分方程描述的是离散变化。

其次,微分方程和差分方程的解法也有相似之处。

两者都可以通过符号计算、数值方法等途径求解。

然而,由于微分方程是连续的,其解法更为灵活和复杂,常常需要应用高级的数学工具,而差分方程在求解过程中则更注重离散的计算方法。

最后,微分方程和差分方程在应用中具有不同的优势。

微分方程主要用于描述连续变化的物理、化学和工程等领域的问题,而差分方程则更适用于计算机科学、经济学和生物学等领域的离散模型。

总之,微分方程和差分方程是数学中研究变化规律的重要工具。

微分方程与差分方程简介

微分方程与差分方程简介

差分方程的分类
一阶差分方程
只包含一个差分的方程,如 (y(n+1) - y(n) = f(n))。
高阶差分方程
包含多个差分的方程,如 (y(n+2) - 2y(n+1) + y(n) = 0)。
线性差分方程
差分项之间线性关系的方程,如 (y(n+1) - y(n) = a + by(n))。
非线性差分方程
05
微分方程与差分方程的 稳定性分析
李雅普诺夫稳定性分析
李雅普诺夫稳定性分析是一种判断动 态系统稳定性的方法,通过分析系统 状态的变化趋势,判断系统是否具有 稳定性。
李雅普诺夫第二方法通过构造一个正 定的李雅普诺夫函数,来研究非线性 系统的稳定性,这种方法适用于非线 性系统的稳定性分析。
线性稳定性分析
经济问题
描述市场供需关系、价格变动、经 济增长等。
03
02
工程问题
控制工程、航空航天、机械工程等 领域。
生物医学问题
描述生理过程、药物动力学、流行 病传播等。
04
02
差分方程简介
差分方程的定义
差分方程是描述离散变量变化规律的数学模型,通常表示为离散变量的函数及其差分之间的关系式。
它与微分方程类似,但时间或空间变量是离散的,而不是连续的。
微分方程与差分方程 简介
目 录
• 微分方程简介 • 差分方程简介 • 微分方程与差分方程的联系与区别 • 微分方程与差分方程的数值解法 • 微分方程与差分方程的稳定性分析
01
微分方程简介
微分方程的定义
1
微分方程是包含一个或多个未知函数的导数的方 程。
2
它描述了某一函数随时间或其他变量的变化规律。

微分方程差分方程

微分方程差分方程

微分方程差分方程微分方程和差分方程是数学中非常重要的两个分支,它们分别研究连续和离散的变化规律。

在物理学、工程学和经济学等领域中,微分方程和差分方程被广泛应用于解决实际问题,为研究和预测系统的行为提供了有效的工具。

首先我们来了解一下什么是微分方程。

微分方程是描述变量及其变化率之间关系的数学方程。

它可以用来描述自然界中的许多现象,如物体的运动、弹簧的振动、电路中的电流变化等。

简单来说,微分方程可以帮助我们理解和预测系统的演变过程。

微分方程的求解通常需要一些基本的数学知识,如微积分、代数等。

通过运用这些知识,我们可以用解析方法或数值方法求解微分方程。

解析方法是通过数学推导来得到方程的解析解,这种方法常用于求解简单的微分方程。

而数值方法通常用于求解复杂的微分方程,它将微分方程转化为差分方程,通过迭代计算逼近方程的解。

差分方程与微分方程有着密切的关系,它是离散时间下的微分方程。

差分方程用来描述离散的变化过程,所以在计算机科学、电子工程和信号处理等领域中得到了广泛的应用。

差分方程的求解方法与微分方程类似,同样可以用解析方法或数值方法进行求解。

在实际应用中,我们常常需要根据已知条件建立微分方程或差分方程,然后利用方程求解问题。

例如,在物理学中,我们可以根据牛顿第二定律建立物体的微分方程,从而描述物体在力的作用下的运动规律。

在经济学中,我们可以根据供需关系建立市场的差分方程,从而研究市场的波动和价格变动。

微分方程和差分方程在实际问题中的应用非常广泛,它们不仅可以帮助我们理解自然和社会现象,还可以用来解决工程和科学中的实际问题。

通过对微分方程和差分方程的研究,我们可以预测和控制系统的行为,优化系统的性能,并取得更好的研究和应用效果。

总结来说,微分方程和差分方程是数学中的两个重要分支,它们分别研究连续和离散的变化规律。

微分方程用于描述变量及其变化率之间的关系,差分方程用于描述离散的变化过程。

通过对微分方程和差分方程的求解,我们可以预测和控制系统的行为,解决实际问题,为科学和工程提供有力的支持。

微分方程差分方程

微分方程差分方程

微分方程差分方程摘要:1.微分方程与差分方程的定义及区别2.微分方程的应用领域3.差分方程的应用领域4.求解微分方程和差分方程的方法5.两者在实际问题中的结合与转化正文:微分方程与差分方程是数学中的两种重要方程类型,它们在许多实际问题中有广泛的应用。

尽管它们具有一定的相似性,但它们之间仍然存在着明显的区别。

本文将对微分方程和差分方程进行简要介绍,并探讨它们在实际问题中的求解方法及应用领域。

一、微分方程与差分方程的定义及区别1.微分方程微分方程是一种描述变量随时间变化的数学方程。

它包含一个或多个未知函数及其导数,要求求解该未知函数在某一区间内的解。

微分方程可以分为线性和非线性两类。

2.差分方程差分方程是一种离散时间模型,它描述了变量在离散时间点上的关系。

差分方程包含一个或多个未知数,并要求求解这些未知数在离散时间点上的取值。

与微分方程类似,差分方程也可以分为线性和非线性两类。

二、微分方程的应用领域1.物理:微分方程在物理学中被广泛应用于描述力学、电磁学、热力学等领域中的现象。

2.生物学:微分方程在生物学中可以用于描述生物种群的数量变化、生长速率等。

3.经济学:微分方程在经济学中可以用于描述物价、产量等经济指标的变化。

4.工程:微分方程在工程领域中可以用于分析结构的动态特性、控制系统的稳定性等。

三、差分方程的应用领域1.计算机科学:差分方程在计算机科学中可以用于数值计算、图像处理等领域。

2.生物学:差分方程在生物学中可以用于模拟生物种群的动态行为。

3.社会科学:差分方程在社会科学中可以用于研究人口统计、经济学模型等。

4.工程:差分方程在工程领域中可以用于分析系统的稳定性、预测发展趋势等。

四、求解微分方程和差分方程的方法1.数值方法:对于微分方程和差分方程,可以采用数值方法求解,如欧拉法、龙格-库塔法等。

2.解析方法:对于一些简单的微分方程和差分方程,可以尝试通过解析方法求解,如分离变量法、常数变易法等。

有限差分解微分方程

有限差分解微分方程

有限差分解微分方程
解:
首先,我们需要将微分方程转化为差分方程。

对于一阶常系数线性微分方程$y' + p(x)y = q(x)$,其对应的差分方程为$y_{n+1} - y_n + p(x_n)y_n = q(x_n)$。

对于题目中的微分方程$y' - 2xy = 1$,我们可以将其转化为差分方程。

首先,我们计算微分方程的导数,得到$y' = 2xy - 1$,然后将其代入差分方程的右边,得到$y_{n+1} - y_n = 2x_ny_n - 1$。

然后,我们需要求解这个差分方程。

我们可以将差分方程转化为递推关系式,即$y_{n+1} = y_n + 2x_ny_n - 1$。

然后我们可以使用迭代法求解这个递推关系式,得到$y_{n+1} = (2x_n - 1)y_n + 1$。

最后,我们需要将解从离散点扩展到整个实数域。

我们可以将解表示为$y(x) = \lim_{n \to \infty} y_n$,然后利用递推关系式求解这个极限,得到$y(x) = \frac{1}{2x - 1}$。

所以,原微分方程的解为$y(x) = \frac{1}{2x - 1}$。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

求解 1. 求差分方程满足初值问题之解:11232133123123(1)3()()()(1)2()()(1)()()2()(1)(1)1,(1)0x n x n x n x n x n x n x n x n x n x n x n x x x +=-+⎧⎪+=+⎪⎨+=-+⎪⎪===⎩ 解:原差分方程组可化为:112233(1)311()(1)201()(1)112()x n x n x n x n x n x n +-⎛⎫⎛⎫⎛⎫⎪ ⎪⎪+= ⎪ ⎪⎪ ⎪ ⎪⎪+-⎝⎭⎝⎭⎝⎭则令311201112-⎛⎫⎪= ⎪ ⎪-⎝⎭A ,求矩阵A 的特征值及特征向量 设特征值分别为123,,λλλ,对应的特征向量分别为123β,β,β.则231121(2)(1)0112λλλλλλ---=-=--=--A E可解得1232,2,1λλλ===设1λ对应的特征向量1111a b c ⎛⎫ ⎪= ⎪ ⎪⎝⎭β,则满足111111022101100a b c -⎛⎫⎛⎫⎛⎫⎪⎪ ⎪-= ⎪⎪ ⎪ ⎪⎪ ⎪-⎝⎭⎝⎭⎝⎭可化简为11100a b c -=⎧⎨=⎩,令111a b ==可以得到特征向量1110⎛⎫⎪= ⎪ ⎪⎝⎭β同理可得到特征向量2110-⎛⎫ ⎪=- ⎪ ⎪⎝⎭β,3011⎛⎫ ⎪= ⎪ ⎪⎝⎭β设方程组的通解为:111222333()nnnx n c c c λλλ=++βββ代入特征值、特征向量,可得到方程组的通解为:123110()21211001n n x n c c c -⎛⎫⎛⎫⎛⎫⎪ ⎪ ⎪=+-+ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭代入初值条件:123(1)(1)1,(1)0x x x ===得到12123322122110n n n n n c c c c c c ⎛⎫-⎛⎫ ⎪ ⎪--= ⎪ ⎪⎪ ⎪⎝⎭⎝⎭可得123120c c c ⎧-=⎪⎨⎪=⎩,可以令11c =,所以212c =;综上所述,满足方程初值方程组的解为:11()210n x n -⎛⎫⎪= ⎪ ⎪⎝⎭2. 求差分方程之通解:2(4)2(2)()32nx n x n x n n n+-++=-+ 解:原方程的特征方程为:42210λλ-+= 即22(1)0λ-=从而求得特征根为11λ=-(二重),21λ=(二重) 因此原方程所对应的齐次方程的通解为:()(1)()1()n n xn A Bn C Dn =-+++ 即 ()(1)()nxn A Bn C Dn =-+++ 而原方程的特解为2(4)2(2)()3x n x n x n n +-++=-的特解1()x n与(4)2(2)()2n x n x n x n n +-++=的特解2()x n 之和.从而原方程具有如下的特解形式:221201201()()()()2()n x n x n x n n A n A n A B n B =+=++++将特解形式代入原方程,可得0010120014811922402244883914890A A A A A AB B B =⎧⎪+=⎪⎪++=-⎨⎪=⎪⎪+=⎩,从而0120114816124194881A A A B B ⎧=⎪⎪⎪=-⎪⎪⎪=⎨⎪⎪=⎪⎪⎪=-⎪⎩综上,原方程的通解为22111148()()()(1)()()2()48624981n n x n xn x n A Bn C Dn n n n n =+=-++++-++- 3. 求微分方程满足初值问题之解:211212212121120d d d 320d d d d d 20d d d (0)1,-1,(0)0d t x x x x x tt t xx x x t t x x x t =⎧++++=⎪⎪⎪++-=⎨⎪⎪===⎪⎩解:方法一:降阶法令13d d x x t =,则原方程组可表示为:13323122312d d d d 320d d d 20d x x t xx x x x tt x x x x t ⎧=⎪⎪⎪++++=⎨⎪⎪++-=⎪⎩化简得:132123323d d d 2d d 22d x x t xx x x t x x x t ⎧=⎪⎪⎪=-+-⎨⎪⎪=--⎪⎩它的系数矩阵为001211022⎛⎫ ⎪=-- ⎪ ⎪--⎝⎭A ,特征方程是01211(2)(2)(1)0022λλλλλλλ--=---=+-++=---A E特征根为1232,2,1λλλ=-==-求得特征根所对应的特征向量分别为1102⎛⎫ ⎪= ⎪ ⎪-⎝⎭T ,21221⎛⎫⎪ ⎪=- ⎪ ⎪⎪⎝⎭T ,31121⎛⎫ ⎪ ⎪= ⎪ ⎪-⎝⎭T .故方程组的通解为1222123311()121()e 0e 2e 221()1t t t x t x t C C C x t --⎛⎫⎛⎫ ⎪⎛⎫⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪=+-+ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪- ⎪⎝⎭⎝⎭- ⎪⎝⎭⎝⎭根据初值1120d (0)1,-1,(0)0d t x x x t====得12312323112211202C C C C C C C C ⎧++=⎪⎪-+-=-⎨⎪⎪-+=⎩解得123112,,463C C C === 则原方程组的解为:22122112()e e e 412311()e e 33t t t t tx t x t ---⎧=++⎪⎪⎨⎪=-+⎪⎩方法二:消元法设dd t λ=,则原方程组可化为21212(32)(1)0(1)(2)(1)0(2)x x x x λλλλλ⎧++++=⎨++-=⎩(1)(2)λ-得21(2)(21)0(3)x λλλ++--= (2)(3)-得22(2)0x λλ--=解得两个特征根为122,1λλ==- 则2x 可表示为:2212e e ttx C C -=+ 根据初值2(0)0x =得22e e ttx C C -=- 将2x 代入(2)得212e 2e ttx C C λ-+=+ 即211d 2e 2e (4)d t t x x C C t-+=+ 下面用常数变易法求解(4) 先解对应齐次方程11d 20d x x t+=得齐次通解211e t x C -= 由常数变易法,令211(t)etx C -=为非齐次方程(4)的解,代入后得221()e e 2e t t t C t C C --'=+积分得41()e 2e 4tt C C t C =+ 则(4)的通解为2211e e 2e 4t tt C x C C --=++ 根据初值110d (0)0,-1d t x x t===得112142212C C C C C C ⎧++=⎪⎪⎨⎪-+-=-⎪⎩解得11314C C ⎧=⎪⎪⎨⎪=⎪⎩ 则221112()e e e 4123t t tx t --=++ 将13C =代入22e e t tx C C -=-得方程组的解为 22122112()e e e 412311()e e 33t t t t tx t x t ---⎧=++⎪⎪⎨⎪=-+⎪⎩4. 利用待定系数法求解下列初值问题之解:Td (),(0)(0,1)d xA x f t x t=+= 其中TT 1235(,),,()(e ,0)53t x x x A f t -⎛⎫===⎪-⎝⎭解:方法一:待定系数法原方程组所对应的齐次方程组为112212d 35d d 53d x x x tx x xt⎧=+⎪⎪⎨⎪=-+⎪⎩特征方程235(3)25053λλλλ--==-+=--A E求得特征根为1,235i λ=±下求135i λ=+所对应的特征向量,设112αα⎛⎫=⎪⎝⎭ξ 则111225i 50()55i 0ααλαα-⎛⎫⎛⎫⎛⎫⎛⎫-==⎪ ⎪ ⎪ ⎪--⎝⎭⎝⎭⎝⎭⎝⎭A E 从而可取11α=,则2i α= 于是由132()1e (cos5isin 5)()i t x t t t x t ⎛⎫⎛⎫=+⎪ ⎪⎝⎭⎝⎭得到齐次方程的通解为:11322()cos5sin 5e ()sin 5cos5t xt C t t x t C t t ⎛⎫⎛⎫⎛⎫= ⎪ ⎪ ⎪-⎝⎭⎝⎭⎝⎭下求非齐次方程的特解利用待定系数法,可设特解为12()e ()e t t x t A x t B --⎛⎫⎛⎫= ⎪ ⎪⎝⎭⎝⎭将其代入原方程组,可得e 3e 5e ee 5e 3et t t tt t tA AB B A B -------⎧-=++⎪⎨-=-+⎪⎩ 即451540A B A B +=-⎧⎨-=⎩,从而求得441541A B ⎧=-⎪⎪⎨⎪=-⎪⎩ 因此原方程的通解为113224()cos5sin 541e e ()sin 5cos5541t t x t C t t x t C t t -⎛⎫-⎪⎛⎫⎛⎫⎛⎫=+ ⎪ ⎪ ⎪ ⎪-⎝⎭ ⎪⎝⎭⎝⎭- ⎪⎝⎭ 代入初值条件T(0)(0,1)x =得到1240415141C C ⎧-=⎪⎪⎨⎪-=⎪⎩,从而124414641C C ⎧=⎪⎪⎨⎪=⎪⎩.综上,原方程组满足初值条件的解为:13244()cos5sin 54141e e ()sin 5cos54654141t t x t t t x t t t -⎛⎫⎛⎫-⎪ ⎪⎛⎫⎛⎫=+ ⎪ ⎪ ⎪ ⎪-⎝⎭⎪ ⎪⎝⎭- ⎪ ⎪⎝⎭⎝⎭方法二:常数变易法利用常数变易法,可设特解为11322()()cos5sin 5e ()()sin 5cos5t x t C t t t x t C t t t ⎛⎫⎛⎫⎛⎫=⎪ ⎪ ⎪-⎝⎭⎝⎭⎝⎭ 带回到原方程,可得到132()cos5sin 5e e ()sin 5cos50t tC t t t C t t t -'⎛⎫⎛⎫⎛⎫= ⎪ ⎪ ⎪'-⎝⎭⎝⎭⎝⎭从而1132()cos5sin 5e e cos5e ()sin 5cos50e sin 5t t t t C t t t t C t t t t ----'⎛⎫⎛⎫⎛⎫⎛⎫== ⎪⎪ ⎪ ⎪'-⎝⎭⎝⎭⎝⎭⎝⎭进而4142()e cos5()e sin 5t tC t t C t t --'⎛⎫⎛⎫= ⎪ ⎪'⎝⎭⎝⎭两边积分可得414254()e (sin 5cos5)414145()e (sin 5cos5)4141t t C t t t C t t t --⎧=-⎪⎪⎨⎪=--⎪⎩因此原方程组的通解为111222()()()()()()x t xt x t x t x t x t ⎛⎫⎛⎫⎛⎫=+ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭13254sin 5cos5cos5sin 5cos5sin 54141e e sin 5cos5sin 5cos545sin 5cos54141t t t t C t t t t C t t t t t t -⎛⎫- ⎪⎛⎫⎛⎫⎛⎫=+⎪ ⎪ ⎪ ⎪--⎝⎭⎝⎭⎪⎝⎭-- ⎪⎝⎭344cos5sin 54141e e sin 5cos54654141t t t t t t -⎛⎫⎛⎫-⎪ ⎪⎛⎫=+ ⎪⎪ ⎪-⎝⎭⎪ ⎪- ⎪ ⎪⎝⎭⎝⎭代入初值条件T(0)(0,1)x =得到1240415141C C ⎧-=⎪⎪⎨⎪-=⎪⎩,从而124414641C C ⎧=⎪⎪⎨⎪=⎪⎩.综上,原方程组满足初值条件的解为13244()cos5sin 54141e e ()sin 5cos54654141t t x t t t x t t t -⎛⎫⎛⎫-⎪ ⎪⎛⎫⎛⎫=+ ⎪ ⎪ ⎪ ⎪-⎝⎭⎪ ⎪⎝⎭- ⎪ ⎪⎝⎭⎝⎭.。

相关文档
最新文档