待定系数法在解非齐次线性微分方程组上的应用.docVIP

待定系数法在解非齐次线性微分方程组上的应用 长期以来，人们在探讨这样一个问题，即当f(t)具有特殊类型时，怎样用待定系数法求出非齐次线性微分方程组x’=ax+f(t)的特解，但是没有取得成功。本文将给出关于这种方法一个比较圆满的回答，它不仅为特殊类型非齐线性微分方程组求解拓宽了渠道，又给解决实际问题带来了方便。 非齐线性微分方程组 x’=ax+f(t)(1) 满足初始条件φ（t0）=η的解，由公式 φ（t）=exp［(t-t0)a］η+∫tt0exp［(t-s)a］·f(s)ds(2) 给出，这里 如果我们知道方程组（1）的一个特解φ（t），则方程组（1）满足初始条件的解就可以写成 φ（t）=exp［（t-t0）a］η+φ（t）.(3) 下面介绍当f(t)具有某些特殊型式时怎样求出特解。 引理［1］若方程（1）中矩阵a的互异特征根为λ1，λ2，…，λl,重数分别为n1,n2,…,nl,(n1+n2+…+n1=n),则有非奇异矩阵t（t为n×n阶矩阵），使得 exp［（t-s）a］=t［exp(t-s)j］t-1，其中j具有约当标准型，即有 矩阵空白处元素均为零。 定理1：方程组（1）当f(t)=(bmtm+…+b1t+b0)eat时有型如 φ（t）=∑m+k-1i=0citieat(4) 的特解。其中bj,ci,(j=0,1,2,…,m;i=0,1,…,m+k-1)分别为维列向量， k=max(n1,n2,…,nl) 证明：显然φ（t）=∫tt0exp［（t-s）a］·f(s)ds是方程组（1）的特解。由引理及f(t)的形式有exp［（t-s）a］·f(s)=t［exp（t-s）j］t-1·f(s),即 其中dnik为ni维列向量(i=1,2…,l)所以 定理2：方程组（1）当f(t)=［a(t)cosβt+b(t)sinβt］eat时有型如 φ（t）=eat［p(t)cosβt+q(t)sinβt］. (5) 的特解［其中p(t),q(t)是次数为k+m-1的实系数多项式，αβ为常数，a（t）,b(t)是多项式，次数一个m次，另一个不超过m次］. 证明：由定理1的证明知道，当α不是实数，而是复数时有关结论仍然正确。现将f(t)表示为指数形式 f(t)=a(t)-ib(t)2e(a+iβ)t+a(t)+ib(t)2e(a-iβ)t, 则根据非齐线性方程的叠加原理，直接利用定理1的结果有 φ（t）=∑k+m-1i=0citie(a+iβ)t eat［∑k+m-1i=0（ci+bi）ticosβt+i∑k+m-1i=0(ci-bi)sinβt］. 令∑k+m-1i=0（ci+bi）ti=p(t) i∑k+m-1i=0(ci-bi)ti=q(t)，则φ（t）=eat［p（t）cosβt+q(t)sinβt］ 由于a（t）-ib(t)与a(t)+ib(t)共轭，根据定理1的证明过程可知∑k+m-1i=0citi与∑k+m-1i=0biti共轭，所以p(t),q（t）为实系数多项式，可见定理2成立。 例1：求方程组x’=0 1 -11 1 0 1 1x+１-tt-t2-1+2t-t2的一个特解 解:∵ |λi-a|=λ(λ-1)2,m=2 ∴k=2,k+m-1=3 设φ（t）=α1α2α3+β1β2β3t+γ1γ2γ3t2+ω1ω2ω3t3 代入方程组得 整理方程组得α1α2α3=000=β1β2β3=10-1， γ1γ2γ3=011=ω1ω2ω3=000 故方程组的一个特解为φ（t）=tt2t2-t 例2：求方程组 x’=3 5-5 3x+e-t0的特解 解：∵|λi-a|=(λ-3-5i)(λ-3+5i) ∴k=1,m=0,k+m-1=0 ∴设φ（t）=αβe-t代入方程组得-α-β=3α+5β+1-5α+3β解得β=-541，α=-441 故方程组的解是φ（t）=-441 541e-t 参考文献 ［1］ 叶彦谦.常微分方程讲义,人民教育出版社. ［2］ 中山大学数力系.常微分方程，人民教育出版社.