深圳建设网站哪家好,h5网站建设模板下载,做电脑网站起什么名字,去除wordpress版本号文章目录 abstract二阶常系数非齐次线性微分方程待定系数法可解类型类型1小结例 abstract
二阶常系数非齐次线性微分方程待定系数法可解决的经典类型1及其解法总结与应用 本文给出类型1为什么可以通过待定求出特解,并且待定函数要设成什么形式推理过程有一定工作量,而在应用中… 文章目录 abstract二阶常系数非齐次线性微分方程待定系数法可解类型类型1小结例 abstract
二阶常系数非齐次线性微分方程待定系数法可解决的经典类型1及其解法总结与应用 本文给出类型1为什么可以通过待定求出特解,并且待定函数要设成什么形式推理过程有一定工作量,而在应用中只需要记住可以用待定系数法求解,以及待定系数函数的形式公式以及公式中各部分的确定方法即可
二阶常系数非齐次线性微分方程 二阶常系数非齐次线性微分方程的一般形式为 y ′ ′ p y ′ q y f ( x ) ypyqyf(x) y′′py′qyf(x)(1),其中 p , q p,q p,q是常数 求方程(1)的通解,归结为求对应齐次方程: y ′ ′ p y ′ q y 0 ypyqy0 y′′py′qy0(2)的通解 Y ( x ) Y(x) Y(x)和一个(1)的特解 y ∗ ( x ) y^*(x) y∗(x);则 Y ( x ) y ∗ ( x ) Y(x)y^*(x) Y(x)y∗(x)为(1)的通解 对于非齐次的二阶常系数线性微分方程,仅有限的类型(以 f ( x ) f(x) f(x)的不同类型作区分)是容易解决的,这里介绍两种类型
待定系数法可解类型
当 f ( x ) f(x) f(x)取两种特殊类型的函数时,可以不用积分的方法求 y ∗ y^* y∗,而是通过待定系数法两种形式分别为: f ( x ) e λ x P m ( x ) f(x)e^{\lambda{x}}P_{m}(x) f(x)eλxPm(x), f ( x ) f(x) f(x) e λ x ( P l ( x ) cos ω x P n ( x ) sin ω x ) e^{\lambda{x}}(P_{l}(x)\cos{\omega}xP_{n}(x)\sin{\omega{x}}) eλx(Pl(x)cosωxPn(x)sinωx)
类型1 当 f ( x ) f(x) f(x) e λ x P m ( x ) e^{\lambda{x}}P_{m}(x) eλxPm(x)(2),其中 λ \lambda λ是常数, P m ( x ) P_{m}(x) Pm(x)为 x x x的一个** m m m次**多项式(设为 ∑ i 0 m a i x i \sum_{i0}^{m}a_ix^{i} ∑i0maixi) 此时方程(1)表示为 y ′ ′ p y ′ q y ypyqy y′′py′qy e λ x P m ( x ) e^{\lambda{x}}P_{m}(x) eλxPm(x)(2-1) 当 λ 0 \lambda0 λ0时, f ( x ) P m ( x ) f(x)P_m(x) f(x)Pm(x),方程(2-1)进一步改写为 y ′ ′ p y ′ q y ypyqy y′′py′qy P m ( x ) P_{m}(x) Pm(x)(2-2) 对式(2)求导 f ′ ( x ) f(x) f′(x) λ e λ x P m ( x ) \lambda e^{\lambda{x}}P_{m}(x) λeλxPm(x) e λ x P m ′ ( x ) e^{\lambda{x}}P_{m}(x) eλxPm′(x) e λ x ( λ P m ( x ) P m ′ ( x ) ) e^{\lambda{x}}(\lambda{P_{m}(x)}P_{m}(x)) eλx(λPm(x)Pm′(x))其中 Q m λ P m ( x ) P m ′ ( x ) Q_{m}\lambda{P_{m}(x)}P_{m}(x) QmλPm(x)Pm′(x)仍然是 m m m次多项式,继续求高阶导数,得到相仿的结论,即 f ( n ) ( x ) f^{(n)}(x) f(n)(x) e λ x Q m [ n ] ( x ) e^{\lambda{x}}Q_{m}^{[n]}(x) eλxQm[n](x),其中 Q m [ n ] ( x ) Q_{m}^{[n]}(x) Qm[n](x)表示 f ( x ) f(x) f(x)第 n n n阶导数包含的 m m m阶多项式 此类型的特点是,多项式函数和指数型函数的乘积的导数仍然是多项式函数和指数型函数的乘积,再官产方程(2-1),这恰好可以并为等号右端形式,即 e λ x P m ( x ) e^{\lambda{x}}P_{m}(x) eλxPm(x) 由此可以推测 y ∗ y^* y∗ Q ( x ) e λ x Q(x)e^{\lambda{x}} Q(x)eλx(3)可能是方程(1)的特解[待定系数法] 其中 Q ( x ) Q(x) Q(x)是某个多项式 对(3)求导, y ∗ ′ y^{*} y∗′ Q ′ ( x ) e λ x Q(x)e^{\lambda{x}} Q′(x)eλx Q ( x ) λ e λ x Q(x)\lambda{e^{\lambda{x}}} Q(x)λeλx e λ x ( Q ′ ( x ) Q ( x ) λ ) e^{\lambda{x}}(Q(x)Q(x)\lambda) eλx(Q′(x)Q(x)λ)(3-1) y ∗ ′ ′ y^{*} y∗′′ λ e λ x ( ( Q ′ ( x ) Q ( x ) λ ) \lambda e^{\lambda{x}}((Q(x)Q(x)\lambda) λeλx((Q′(x)Q(x)λ) e λ x ( Q ′ ′ ( x ) λ Q ′ ( x ) ) e^{\lambda{x}}(Q(x)\lambda Q(x)) eλx(Q′′(x)λQ′(x)) e λ x ( λ 2 Q ( x ) 2 λ Q ′ ( x ) Q ′ ′ ( x ) ) e^{\lambda{x}}(\lambda^2{Q(x)}2\lambda{Q(x)}Q(x)) eλx(λ2Q(x)2λQ′(x)Q′′(x))(3-2) 将(3,3-1,3-2)代入方程(2-1),得 e λ x ( λ 2 Q ( x ) 2 λ Q ′ ( x ) Q ′ ′ ( x ) ) e^{\lambda{x}}(\lambda^2{Q(x)}2\lambda{Q(x)}Q(x)) eλx(λ2Q(x)2λQ′(x)Q′′(x)) p e λ x ( Q ′ ( x ) Q ( x ) λ ) pe^{\lambda{x}}(Q(x)Q(x)\lambda) peλx(Q′(x)Q(x)λ) q Q ( x ) e λ x qQ(x)e^{\lambda{x}} qQ(x)eλx e λ x P m ( x ) e^{\lambda{x}}P_{m}(x) eλxPm(x),整理得 Q ′ ′ ( x ) ( 2 λ p ) Q ′ ( x ) ( λ 2 p λ q ) Q ( x ) Q(x)(2\lambdap)Q(x)(\lambda^2p\lambdaq)Q(x) Q′′(x)(2λp)Q′(x)(λ2pλq)Q(x) P m ( x ) P_{m}(x) Pm(x)(4) 根据 λ \lambda λ与方程(1)的特征方程 ( r 2 p r q 0 ) (r^2prq0) (r2prq0)(5)的根(特征根)关系,分为:不是特征根,单根,重根,这三种情形讨论 若 λ \lambda λ不是(5)的根,则 λ 2 p λ q ≠ 0 \lambda^2p\lambdaq\neq{0} λ2pλq0(5-1), 由于 P m ( x ) P_{m}(x) Pm(x)是 m m m次多项式,要使(4)式两边相等,则 Q ( x ) Q(x) Q(x)必须也是 m m m次多项式,记为 Q ( x ) Q m ( x ) Q(x)Q_{m}(x) Q(x)Qm(x) ∑ i 0 m b i x i \sum_{i0}^{m}b_ix^{i} ∑i0mbixi(6)将(6)代入方程(4),比较两边同次幂的系数,得到 0 ∼ m 0\sim{m} 0∼m次共 m 1 m1 m1个方程解这 m 1 m1 m1个方程,可得 b 0 , ⋯ , b m b_0,\cdots,b_m b0,⋯,bm代入(3),从而得到(2-1)的特解 若 λ \lambda λ是(5)的单根(两个互异根中的一个),则 λ 2 p λ q 0 \lambda^2p\lambdaq{0} λ2pλq0(6-1) 此处 λ \lambda λ是单根,设另一根是 λ ‾ \overline\lambda λ;由韦达定理: λ λ ‾ − p \lambda\overline{\lambda}-p λλ−p,而 λ λ ≠ λ λ ‾ \lambda\lambda\neq{\lambda\overline{\lambda}} λλλλ,所以 2 λ ≠ − p 2\lambda\neq{-p} 2λ−p,即 2 λ p ≠ 0 2\lambdap\neq{0} 2λp0(7) 此时方程(4)改写为 Q ′ ′ ( x ) ( 2 λ p ) Q ′ ( x ) Q(x)(2\lambdap)Q(x) Q′′(x)(2λp)Q′(x) P m ( x ) P_{m}(x) Pm(x);因此 Q ′ ( x ) Q(x) Q′(x)必须是 m m m次多项式;相应的, Q ( x ) Q(x) Q(x)就得是 m 1 m1 m1次多项式 令 Q ( x ) x Q m ( x ) Q(x)xQ_m(x) Q(x)xQm(x)(8);仍然可以用系数比较法确定出 Q m ( x ) Q_{m}(x) Qm(x)的系数 b 0 , b 1 , ⋯ , b m b_0,b_1,\cdots,b_m b0,b1,⋯,bm;就可以得出 Q m ( x ) Q_m(x) Qm(x),代入(8)得出 Q ( x ) Q(x) Q(x),再代入(3),得方程(2-1)的特解 若 λ \lambda λ是(5)的重根,此是也有(6-1)成立,并且 2 λ p 0 2\lambdap0 2λp0(8-1) 此时方程(4)改写为 Q ′ ′ ( x ) P m ( x ) Q(x)P_{m}(x) Q′′(x)Pm(x)要使方程(4)两端恒等,必有 Q ′ ′ ( x ) Q(x) Q′′(x)为 m m m次多项式,从而可以令 Q ( x ) x 2 Q m ( x ) Q(x)x^2Q_{m}(x) Q(x)x2Qm(x)同样使用系数比较法确定出 Q m ( x ) Q_{m}(x) Qm(x)的系数,从而得出 Q ( x ) Q(x) Q(x),最后代入(3)得出(2-1)的特解 上述类型1的结论可知,方程(1)确实可以使用待定系数法求特解,该特解可以设为 y ∗ y^* y∗ x k Q m ( x ) e λ x x^kQ_{m}(x)e^{\lambda{x}} xkQm(x)eλx,并且保证这个待定形式是可求解且正确的
小结
二阶常系数非齐次线性微分方程 y ′ ′ p y ′ q y ypyqy y′′py′qy P m ( x ) e λ x P_m(x)e^{\lambda{x}} Pm(x)eλx(1)具有形如 y ∗ x k Q m ( x ) e λ x y^*x^kQ_{m}(x)e^{\lambda{x}} y∗xkQm(x)eλx的特解, ( k 0 , 1 , 2 ) (k0,1,2) (k0,1,2) 其中 Q m ( x ) Q_{m}(x) Qm(x), P m ( x ) P_m(x) Pm(x)是同为 m m m次多项式, Q m ( x ) Q_{m}(x) Qm(x)的 m 1 m1 m1个系数由系数比较法构造 m 1 m1 m1个方程分别求出 k k k按照 λ \lambda λ是方程(1)的特征方程 r 2 p r q 0 r^2prq0 r2prq0(2)的根的重数决定的, 0 0 0重根, k 0 k0 k0;(表示 λ \lambda λ不是方程(2)的根的简称) 1 1 1重根, k 1 k1 k1 2 2 2重根, k 2 k2 k2
例
求 y ′ ′ − 2 y ′ − 3 y y-2y-3y y′′−2y′−3y 3 x 1 3x1 3x1(1)的一个特解 方程类型分析:方程(1)是常系数非齐次线性微分方程中的经典类型1 令 f ( x ) f(x) f(x) 3 x 1 3x1 3x1 e 0 x ( 3 x 1 ) e^{0x}(3x1) e0x(3x1),即对应类型1问题模型中有(1-1) λ 0 \lambda0 λ0; P m ( x ) P_m(x) Pm(x)是一次多项式 3 x 1 3x1 3x1 m 1 m1 m1 特征方程为 r 2 − 2 r − 3 0 r^2-2r-30 r2−2r−30(2)检查 λ \lambda λ和(2)的根的关系: λ 0 \lambda0 λ0不是(2)的根,应用类型1的待定系数法结论可知,方程(1)的特解可以设为 y ∗ y^* y∗ x k Q m ( x ) e λ x x^kQ_{m}(x)e^{\lambda{x}} xkQm(x)eλx,代入(1-1),可以具体为 y ∗ y^* y∗ Q 1 ( x ) Q_{1}(x) Q1(x) b 0 x b 1 b_0{x}b_1 b0xb1(3)( Q 1 ( x ) Q_1(x) Q1(x)是一次多项式)将(3)代入到方程(1): − 2 b 0 − 3 ( b 0 x b 1 ) -2b_0-3(b_0{x}b_1) −2b0−3(b0xb1) 3 x 1 3x1 3x1,整理得 − 3 b 0 x − 2 b 0 − 3 b 1 -3b_0x-2b_0-3b_1 −3b0x−2b0−3b1 3 x − 1 3x-1 3x−1(4)由系数比较法, − 3 b 0 3 -3b_03 −3b03; − 2 b 0 − 3 b 1 -2b_0-3b_1 −2b0−3b1 1 1 1,解得 b 0 − 1 b_0-1 b0−1; b 1 1 3 b_1\frac{1}{3} b131,从而得到特解为 y ∗ − x 1 3 y^*-x\frac{1}{3} y∗−x31
