微分方程 —— 含导数或微分的方程称为微分方程。
微分方程的阶数 —— 微分方程中所含的导数或微分的最高阶数称为微分方程的阶数。
微分方程的解 —— 使得微分方程成立的函数称为微分方程的解。不含任意常数的解称为微分方程的特解;若微分方程的解所含有的相互独立的任意常数的个数与微分方程的阶数相等,称此解为微分方程的通解。
形如 d y / d x = f ( x , y ) dy/dx =f(x,y) dy/dx=f(x,y) 的一阶微分方程,称为可分离变量的微分方程。顾名思义,可以把 x x x 和 d x dx dx 放在一起,这样只需要两边同时积分就行。
举个例子: d y / d x = 1 + x + y 2 + x y 2 dy/dx =1+x+y^2+xy^2 dy/dx=1+x+y2+xy2 ,等式右端因式分解为 ( 1 + x ) ( 1 + y 2 ) (1+x)(1+y^2) (1+x)(1+y2) ,可分离变量为 d y / ( 1 + y 2 ) = ( 1 + x ) d x dy/(1+y^2)=(1+x)dx dy/(1+y2)=(1+x)dx ,两边同时积分,得到 a r c t a n y = x + 1 2 x 2 + C . arctany=x+\frac{1}{2}x^2+C. arctany=x+21x2+C.
形如 d y / d x = f ( x , y ) dy/dx =f(x,y) dy/dx=f(x,y) ,且 f ( x , y ) = φ ( y x ) f(x,y)=\varphi(\frac{y}{x}) f(x,y)=φ(xy) ,称为齐次微分方程。即 f ( t x , t y ) = f ( x , y ) . f(tx,ty)=f(x,y). f(tx,ty)=f(x,y).
令 u = y / x u=y/x u=y/x ,则 d y / d x = d u / d x + u dy/dx=du/dx+u dy/dx=du/dx+u,代入原方程得 u + x d u / d x = φ ( u ) u+xdu/dx=\varphi(u) u+xdu/dx=φ(u) ,分离变量即可。
举个例子,求微分方程 x d y / d x = y + x 2 + y 2 ( x > 0 ) xdy/dx=y+\sqrt{x^2+y^2}(x>0) xdy/dx=y+x2+y2(x>0) 满足 y ( 1 ) = 0 y(1)=0 y(1)=0 的特解。
解: 等式两边除以 x x x ,得 d y / d x = y / x + 1 + ( y / x ) 2 . dy/dx=y/x+\sqrt{1+(y/x)^2}. dy/dx=y/x+1+(y/x)2. 令 u = y / x u=y/x u=y/x ,整理可得 d u 1 + u 2 = d x x . \frac{du}{\sqrt{1+u^2}}=\frac{dx}{x}. 1+u2du=xdx. 积分得 u + 1 + u 2 = C x u+\sqrt{1+u^2}=Cx u+1+u2=Cx ,由 y ( 0 ) = 1 y(0)=1 y(0)=1 ,可得 C = 1 C=1 C=1 ,即 u + 1 + u 2 = x u+\sqrt{1+u^2}=x u+1+u2=x ,取倒数,则有 u − 1 + u 2 = 1 / x u-\sqrt{1+u^2}=1/x u−1+u2=1/x ,联立可得到特解 y = ( x 2 − 1 ) / 2. y=(x^2-1)/2. y=(x2−1)/2.
形如 d y / d x + P ( x ) y = 0 dy/dx+P(x)y=0 dy/dx+P(x)y=0 的方程,称为一阶齐次线性微分方程。其通解公式为 y = C e − ∫ P ( x ) d x . y=Ce^{-\int P(x)dx}. y=Ce−∫P(x)dx.
形如 d y / d x + P ( x ) y = Q ( x ) dy/dx+P(x)y=Q(x) dy/dx+P(x)y=Q(x) 的方程,称为一阶非齐次线性微分方程。其通解公式为 y = [ ∫ Q ( x ) e ∫ P ( x ) d x d x + C ] e − ∫ P ( x ) d x . y=[\int Q(x)e^{\int P(x)dx}dx+C ]e^{-\int P(x)dx}. y=[∫Q(x)e∫P(x)dxdx+C]e−∫P(x)dx.
通解公式推导第一步为等式两边同时乘 e ∫ P ( x ) d x e^{\int P(x)dx} e∫P(x)dx ,这公式用的相对蛮多的,还是记忆一下比较方便些。
形如 d y / d x + P ( x ) y = Q ( x ) y n ( n ≠ 0 , 1 ) dy/dx+P(x)y=Q(x)y^n(n \ne0,1) dy/dx+P(x)y=Q(x)yn(n=0,1) 的方程,称为伯努利方程。
令 z = y 1 − n z=y^{1-n} z=y1−n ,代入原方程得 d z / d x + ( 1 − n ) P ( x ) z = ( 1 − n ) Q ( x ) dz/dx+(1-n)P(x)z=(1-n)Q(x) dz/dx+(1−n)P(x)z=(1−n)Q(x) ,求解该关于 z , x z,x z,x 的一阶非齐次线性微分方程即可。
设 P ( x , y ) d x + Q ( x , y ) d y = 0 P(x,y)dx+Q(x,y)dy=0 P(x,y)dx+Q(x,y)dy=0 ,满足 ∂ Q / ∂ x = ∂ P / ∂ y \partial Q/\partial x=\partial P/\partial y ∂Q/∂x=∂P/∂y ,则称该方程为全微分方程。
第一种:形如 y ( n ) = f ( x ) y^{(n)}=f(x) y(n)=f(x) ,进行多次不定积分即可求解。
第二种:形如 f ( x , y ′ , y ′ ′ ) = 0 f(x,y',y'')=0 f(x,y′,y′′)=0(缺 y y y 型)。解法是令 y ′ = d y / d x = p ( x ) y'=dy/dx=p(x) y′=dy/dx=p(x) ,则 y ′ ′ = d p / d x . y''=dp/dx. y′′=dp/dx.
第三种:形如 f ( y , y ′ , y ′ ′ ) = 0 f(y,y',y'')=0 f(y,y′,y′′)=0(缺 x x x 型)。解法是令 y ′ = d y / d x = p ( y ) y'=dy/dx=p(y) y′=dy/dx=p(y) ,则 y ′ ′ = p d p / d y . y''=pdp/dy. y′′=pdp/dy.
y ( n ) + a 1 ( x ) y ( n − 1 ) + ⋯ + a n − 1 ( x ) y ′ + a n ( x ) y = 0 ( 4 − 1 ) y^{(n)}+a_1(x)y^{(n-1)}+\dots+a_{n-1}(x)y'+a_n(x)y=0(4-1) y(n)+a1(x)y(n−1)+⋯+an−1(x)y′+an(x)y=0(4−1) 称为 n n n 阶齐次线性微分方程,其中 a 1 ( x ) , … , a n − 1 ( x ) , a n ( x ) a_1(x),\dots,a_{n-1}(x),a_n(x) a1(x),…,an−1(x),an(x) 为 x x x 的函数。
y ( n ) + a 1 ( x ) y ( n − 1 ) + ⋯ + a n − 1 ( x ) y ′ + a n ( x ) y = f ( x ) ( 4 − 2 ) y^{(n)}+a_1(x)y^{(n-1)}+\dots+a_{n-1}(x)y'+a_n(x)y=f(x)(4-2) y(n)+a1(x)y(n−1)+⋯+an−1(x)y′+an(x)y=f(x)(4−2) 称为 n n n 阶非齐次线性微分方程,其中 a 1 ( x ) , … , a n − 1 ( x ) , a n ( x ) a_1(x),\dots,a_{n-1}(x),a_n(x) a1(x),…,an−1(x),an(x) 为 x x x 的函数。
若 f ( x ) = f 1 ( x ) + f 2 ( x ) f(x)=f_1(x)+f_2(x) f(x)=f1(x)+f2(x) ,则非齐次方程可分解为如下两个方程: y ( n ) + a 1 ( x ) y ( n − 1 ) + ⋯ + a n − 1 ( x ) y ′ + a n ( x ) y = f 1 ( x ) ( 4 − 3 ) y^{(n)}+a_1(x)y^{(n-1)}+\dots+a_{n-1}(x)y'+a_n(x)y=f_1(x)(4-3) y(n)+a1(x)y(n−1)+⋯+an−1(x)y′+an(x)y=f1(x)(4−3) y ( n ) + a 1 ( x ) y ( n − 1 ) + ⋯ + a n − 1 ( x ) y ′ + a n ( x ) y = f 2 ( x ) ( 4 − 4 ) y^{(n)}+a_1(x)y^{(n-1)}+\dots+a_{n-1}(x)y'+a_n(x)y=f_2(x) (4-4) y(n)+a1(x)y(n−1)+⋯+an−1(x)y′+an(x)y=f2(x)(4−4)
(1)设 φ 1 ( x ) , φ 2 ( x ) , … , φ s ( x ) \varphi_1(x),\varphi_2(x),\dots,\varphi_s(x) φ1(x),φ2(x),…,φs(x) 为方程(4-1)的一组解,则 k 1 φ 1 ( x ) + k 2 φ 2 ( x ) + ⋯ + k s φ s ( x ) k_1\varphi_1(x)+k_2\varphi_2(x)+\dots+k_s\varphi_s(x) k1φ1(x)+k2φ2(x)+⋯+ksφs(x) 也为方程(4-1)的解。
齐次解线性组合仍为齐次解。
(2)若 φ 1 ( x ) , φ 2 ( x ) \varphi_1(x),\varphi_2(x) φ1(x),φ2(x) 分别为方程(4-1)和(4-2)的解,则 φ 1 ( x ) + φ 2 ( x ) \varphi_1(x)+\varphi_2(x) φ1(x)+φ2(x) 为方程(4-2)的解。
齐次解+非齐次解等于非齐次解。
(3)若 φ 1 ( x ) , φ 2 ( x ) \varphi_1(x),\varphi_2(x) φ1(x),φ2(x) 为方程(4-2)的两个解,则 φ 1 ( x ) − φ 2 ( x ) \varphi_1(x)-\varphi_2(x) φ1(x)−φ2(x) 为方程(4-1)的解。
非齐次解相减为齐次解。
(4)若 φ 1 ( x ) , φ 2 ( x ) \varphi_1(x),\varphi_2(x) φ1(x),φ2(x) 分别为方程(4-3)和(4-4)的解,则 φ 1 ( x ) + φ 2 ( x ) \varphi_1(x)+\varphi_2(x) φ1(x)+φ2(x) 为方程(4-2)的解。
(5)设 φ 1 ( x ) , φ 2 ( x ) , … , φ s ( x ) \varphi_1(x),\varphi_2(x),\dots,\varphi_s(x) φ1(x),φ2(x),…,φs(x) 为方程(4-2)的一组解,则 k 1 φ 1 ( x ) + k 2 φ 2 ( x ) + ⋯ + k s φ s ( x ) k_1\varphi_1(x)+k_2\varphi_2(x)+\dots+k_s\varphi_s(x) k1φ1(x)+k2φ2(x)+⋯+ksφs(x) 为方程(4-2)的解的充分必要条件是 k 1 + k 2 + ⋯ + k s = 1. k_1+k_2+\dots+k_s=1. k1+k2+⋯+ks=1.
非齐次解线性组合,系数之和为 1 才还是非齐解。
(6)设 φ 1 ( x ) , φ 2 ( x ) , … , φ s ( x ) \varphi_1(x),\varphi_2(x),\dots,\varphi_s(x) φ1(x),φ2(x),…,φs(x) 为方程(4-2)的一组解,则 k 1 φ 1 ( x ) + k 2 φ 2 ( x ) + ⋯ + k s φ s ( x ) k_1\varphi_1(x)+k_2\varphi_2(x)+\dots+k_s\varphi_s(x) k1φ1(x)+k2φ2(x)+⋯+ksφs(x) 为方程(4-1)的解的充分必要条件是 k 1 + k 2 + ⋯ + k s = 0. k_1+k_2+\dots+k_s=0. k1+k2+⋯+ks=0.
非齐次解线性组合,系数之和为 0 才是齐次解。
(7)设 φ 1 ( x ) , φ 2 ( x ) , … , φ n ( x ) \varphi_1(x),\varphi_2(x),\dots,\varphi_n(x) φ1(x),φ2(x),…,φn(x) 为方程(4-1)的 n n n 个线性无关解,则 k 1 φ 1 ( x ) + k 2 φ 2 ( x ) + ⋯ + k n φ n ( x ) k_1\varphi_1(x)+k_2\varphi_2(x)+\dots+k_n\varphi_n(x) k1φ1(x)+k2φ2(x)+⋯+knφn(x) 为方程(4-1)的通解。
(8)设 φ 1 ( x ) , φ 2 ( x ) , … , φ n ( x ) \varphi_1(x),\varphi_2(x),\dots,\varphi_n(x) φ1(x),φ2(x),…,φn(x) 为方程(4-1)的 n n n 个线性无关解, φ 0 ( x ) \varphi_0(x) φ0(x) 为方程(4-2)的一个特解,则 k 1 φ 1 ( x ) + k 2 φ 2 ( x ) + ⋯ + k n φ n ( x ) + φ 0 ( x ) k_1\varphi_1(x)+k_2\varphi_2(x)+\dots+k_n\varphi_n(x)+\varphi_0(x) k1φ1(x)+k2φ2(x)+⋯+knφn(x)+φ0(x) 也为方程(4-2)的通解。
非齐通解为齐次通解+非齐特解。
(一)二阶常系数齐次线性微分方程的解法
形如 y ′ ′ + p y ′ + q y = 0 y''+py'+qy=0 y′′+py′+qy=0(其中 p , q p,q p,q 为常数)的方程称为二阶常系数齐次线性微分方程。其求解步骤如下:
(1)求解方程 y ′ ′ + p y ′ + q y = 0 y''+py'+qy=0 y′′+py′+qy=0 的特征方程 λ 2 + p λ + q = 0. \lambda^2+p\lambda+q=0. λ2+pλ+q=0.
(2)根据特征方程根的不同,分为如下三种情形:
(2.1)当判别式 Δ = p 2 − 4 q > 0 \Delta =p^2-4q>0 Δ=p2−4q>0 时,两特征根 λ 1 ≠ λ 2 \lambda_1 \ne \lambda_2 λ1=λ2 ,则原方程的通解为 y = C 1 e λ 1 x + C 2 e λ 2 x ; y=C_1e^{\lambda_1x}+C_2e^{\lambda_2x}; y=C1eλ1x+C2eλ2x;(2.2)当判别式 Δ = p 2 − 4 q = 0 \Delta =p^2-4q=0 Δ=p2−4q=0 时,方程有两个相同的特征根 λ 1 = λ 2 \lambda_1 = \lambda_2 λ1=λ2 ,则原方程的通解为 y = ( C 1 + C 2 x ) e λ 1 x ; y=(C_1+C_2x)e^{\lambda_1x}; y=(C1+C2x)eλ1x; (2.3)当判别式 Δ = p 2 − 4 q < 0 \Delta =p^2-4q<0 Δ=p2−4q<0 时,特征方程有两个共轭虚根 α ± i β \alpha \pm i\beta α±iβ ,则原方程的通解为 y = e α x ( C 1 c o s β x + C 2 s i n β x ) . y=e^{\alpha x}(C_1cos\beta x+C_2sin\beta x). y=eαx(C1cosβx+C2sinβx).
对于三阶常系数齐次线性微分方程 y ′ ′ ′ + p y ′ ′ + q y ′ + r y = 0 y'''+py''+qy'+ry=0 y′′′+py′′+qy′+ry=0 ,解法类似二阶,不过根的种类组合较多。其特征方程 λ 3 + p λ 2 + q λ + r = 0 \lambda^3+p\lambda^2+q\lambda+r=0 λ3+pλ2+qλ+r=0 有三个根 λ 1 , λ 2 , λ 3 \lambda_1,\lambda_2,\lambda_3 λ1,λ2,λ3 ,可分为如下几种情形:
(1) λ 1 , λ 2 , λ 3 ∈ R \lambda_1,\lambda_2,\lambda_3 \in R λ1,λ2,λ3∈R 且两两不相等,类比判别式大于 0 的情况,其通解为 y = C 1 e λ 1 x + C 2 e λ 2 x + C 3 e λ 3 x ; y=C_1e^{\lambda_1x}+C_2e^{\lambda_2x}+C_3e^{\lambda_3x}; y=C1eλ1x+C2eλ2x+C3eλ3x; (2) λ 1 , λ 2 , λ 3 ∈ R \lambda_1,\lambda_2,\lambda_3 \in R λ1,λ2,λ3∈R 且有一对根相等 λ 1 = λ 2 ≠ λ 3 \lambda_1=\lambda_2\ne \lambda_3 λ1=λ2=λ3 ,相等的那一对根类比判别式等于 0 的情况,剩下的一根类比大于 0 的情况,则其通解为 y = ( C 1 + C 2 x ) e λ 1 x + C 3 e λ 3 x ; y=(C_1+C_2x)e^{\lambda_1x}+C_3e^{\lambda_3x}; y=(C1+C2x)eλ1x+C3eλ3x; (3) λ 1 , λ 2 , λ 3 ∈ R \lambda_1,\lambda_2,\lambda_3 \in R λ1,λ2,λ3∈R 且 λ 1 = λ 2 = λ 3 \lambda_1=\lambda_2= \lambda_3 λ1=λ2=λ3 ,类比判别式为 0 的情况,则其通解为 y = ( C 1 + C 2 x + C 3 x 2 ) e λ 1 x ; y=(C_1+C_2x+C_3x^2)e^{\lambda_1x}; y=(C1+C2x+C3x2)eλ1x; (4) λ 1 ∈ R , λ 2 , 3 = α ± i β ( β ≠ 0 ) \lambda_1 \in R,\lambda_{2,3}=\alpha \pm i\beta(\beta \ne 0) λ1∈R,λ2,3=α±iβ(β=0) ,共轭虚跟类比判别式小于 0 的情况,剩余的实根类别判别式大于 0 的情况,则其通解为 y = C 1 e λ 1 x + e α x ( C 2 c o s β x + C 3 s i n β x ) . y=C_1e^{\lambda_1x}+e^{\alpha x}(C_2cos\beta x+C_3sin\beta x). y=C1eλ1x+eαx(C2cosβx+C3sinβx).
(二)二阶常系数非齐次线性微分方程的特解
形如 y ′ ′ + p y ′ + q y = f ( x ) y''+py'+qy=f(x) y′′+py′+qy=f(x)(其中 p , q p,q p,q 为常数)的方程称为二阶常系数非齐次线性微分方程。根据 f ( x ) f(x) f(x) 的不同形式,求特解方程分为如下两种情况:
(1) f ( x ) = P n ( x ) e k x . f(x)=P_n(x)e^{kx}. f(x)=Pn(x)ekx.
C a s e I Case I CaseI —— 若 k k k 非特征根,可令特解 y 0 = ( a 0 + a 1 x + ⋯ + a n x n ) e k x y_0=(a_0+a_1x+\dots +a_nx^n)e^{kx} y0=(a0+a1x+⋯+anxn)ekx ,即 P n ( x ) P_n(x) Pn(x) 是几次,特解中的多项式就几次。如 y ′ ′ − y ′ − 2 y = ( x + 1 ) e x y''-y'-2y=(x+1)e^x y′′−y′−2y=(x+1)ex ,其特解方程应设为 y 0 = ( a x + b ) e x ; y_0=(ax+b)e^x; y0=(ax+b)ex;
C a s e I I Case II CaseII —— 若 k k k 与其中一个特征根相同,可令 y 0 = x ( a 0 + a 1 x + ⋯ + a n x n ) e k x y_0=x(a_0+a_1x+\dots +a_nx^n)e^{kx} y0=x(a0+a1x+⋯+anxn)ekx ,即比(1.1)多乘了一个 x x x 。
C a s e I I I Case III CaseIII —— 若 k k k 与两个特征根均相同,可令 y 0 = x 2 ( a 0 + a 1 x + ⋯ + a n x n ) e k x y_0=x^2(a_0+a_1x+\dots +a_nx^n)e^{kx} y0=x2(a0+a1x+⋯+anxn)ekx ,即比(1.1)多乘了一个 x 2 x^2 x2 。
得出特解形式后,将特解代入原方程进行整理,可得出相关系数。
(2) f ( x ) = e α x ( P l ( x ) c o s β x + P s ( x ) s i n β x ) . f(x)=e^{\alpha x}(P_l(x)cos\beta x+P_s(x)sin\beta x). f(x)=eαx(Pl(x)cosβx+Ps(x)sinβx).
令 n = m a x n=max n=max { l , s l,s l,s} 。
C a s e I Case I CaseI —— 若 α ± i β \alpha \pm i\beta α±iβ 不是特征根,则令 y 0 ( x ) = e α x ( Q n ( 1 ) ( x ) c o s β x + Q n ( 2 ) ( x ) s i n β x ) ; y_0(x)=e^{\alpha x}(Q_n^{(1)}(x)cos\beta x+Q^{(2)}_n(x)sin\beta x); y0(x)=eαx(Qn(1)(x)cosβx+Qn(2)(x)sinβx);
C a s e I I Case II CaseII —— 若 α ± i β \alpha \pm i\beta α±iβ 是特征根,则令 y 0 ( x ) = x e α x ( Q n ( 1 ) ( x ) c o s β x + Q n ( 2 ) ( x ) s i n β x ) ; y_0(x)=xe^{\alpha x}(Q_n^{(1)}(x)cos\beta x+Q^{(2)}_n(x)sin\beta x); y0(x)=xeαx(Qn(1)(x)cosβx+Qn(2)(x)sinβx);
举个例子,设 y ′ ′ − 2 y ′ + 2 y = x e x c o s x , y''-2y'+2y=xe^xcosx, y′′−2y′+2y=xexcosx, 求该方程特解形式。
解: 由特征方程 λ 2 − 2 λ + 2 = 0 \lambda^2-2\lambda+2=0 λ2−2λ+2=0 解得 λ 1 , 2 = 1 ± i \lambda_{1,2}=1 \pm i λ1,2=1±i 。 x e x c o s x = e x ( x c o s x + 0 s i n x ) xe^xcosx=e^x(xcosx+0sinx) xexcosx=ex(xcosx+0sinx) ,故 α = 1 , β = 1 , n = 1 \alpha=1,\beta=1,n=1 α=1,β=1,n=1 ,所以 α ± i β \alpha \pm i\beta α±iβ 是特征根,故特解形式为 y 0 ( x ) = x e x ( ( a x + b ) c o s x + ( c x + d ) s i n x ) . y_0(x)=xe^{ x}((ax+b)cos x+(cx+d)sin x). y0(x)=xex((ax+b)cosx+(cx+d)sinx).
假设特解有两个要点:
(1)根据右边形式进行假设,且正弦余弦都需要,多项式取最高次。
(2)如果 α ± i β \alpha \pm i\beta α±iβ 是特征根,需要多乘一个 x . x. x.
形如 x n y ( n ) + a n − 1 x n − 1 y ( n − 1 ) + a 1 x y ′ + a 0 y = f ( x ) x^ny^{(n)}+a_{n-1}x^{n-1}y^{(n-1)}+a_1xy'+a_0y=f(x) xny(n)+an−1xn−1y(n−1)+a1xy′+a0y=f(x) 的方程称为欧拉方程,其中 a 0 , a 1 , … , a n − 1 a_0,a_1,\dots,a_{n-1} a0,a1,…,an−1 为常数。
令 x = e t , d x / d t = e t = x x=e^t,dx/dt=e^t=x x=et,dx/dt=et=x ,则有 x y ′ = x d y / d t d x / d t = d y / d t xy'=x\frac{dy/dt}{dx/dt}=dy/dt xy′=xdx/dtdy/dt=dy/dt x 2 y ′ ′ = x 2 ( d ( d y d x ) d x ) = x 2 d ( d y / d t d x / d t ) d x = x 2 d ( d y / d t x ) / d t d x / d t = x 2 d 2 y / d t 2 − d y / d t x x = d 2 y d t 2 − d y d t x^2y''=x^2(\frac{d(\frac{dy}{dx})}{dx})=x^2\frac{d(\frac{dy/dt}{dx/dt})}{dx}=x^2\frac{d(\frac{dy/dt}{x})/dt}{dx/dt}=x^2\frac{\frac{d^2y/dt^2-dy/dt}{x}}{x}=\frac{d^2y}{dt^2}-\frac{dy}{dt} x2y′′=x2(dxd(dxdy))=x2dxd(dx/dtdy/dt)=x2dx/dtd(xdy/dt)/dt=x2xxd2y/dt2−dy/dt=dt2d2y−dtdy 记 d d t = D \frac{d}{dt}=D dtd=D ,有 D y = d y / d t , D 2 y = d 2 y / d t 2 Dy=dy/dt,D^2y=d^2y/dt^2 Dy=dy/dt,D2y=d2y/dt2 ,则 x y ′ = D , x 2 y ′ ′ = D 2 y − D = [ D ( D − 1 ) ] y xy'=D,x^2y''=D^2y-D=[D(D-1)]y xy′=D,x2y′′=D2y−D=[D(D−1)]y ,故 x n y ( n ) = D ( D − 1 ) ( D − n + 1 ) y . x^ny^{(n)} =D(D-1)(D-n+1)y. xny(n)=D(D−1)(D−n+1)y. 代入原方程,得到一个关于 y , t y,t y,t 的高阶常系数非齐次线性微分方程。
举个例子,求微分方程 x 2 y ′ ′ + 2 x y ′ − 2 y = 3 x − 2 ( x > 0 ) x^2y''+2xy'-2y=3x-2(x>0) x2y′′+2xy′−2y=3x−2(x>0) 的通解。
解: 令 x = e t x=e^t x=et ,则有 x y ′ = d y / d t , x 2 y ′ ′ = d 2 y / d t 2 − d y / d t xy'=dy/dt,x^2y''=d^2y/dt^2-dy/dt xy′=dy/dt,x2y′′=d2y/dt2−dy/dt ,代入原方程,得 d 2 y d t 2 + d y d t − 2 y = 3 e t − 2 ( 1 ) \frac{d^2y}{dt^2}+\frac{dy}{dt}-2y=3e^t-2(1) dt2d2y+dtdy−2y=3et−2(1) 根据解的结构性质(4),将其拆分为两个方程: d 2 y d t 2 + d y d t − 2 y = 3 e t ( 1.1 ) \frac{d^2y}{dt^2}+\frac{dy}{dt}-2y=3e^t(1.1) dt2d2y+dtdy−2y=3et(1.1) 和 d 2 y d t 2 + d y d t − 2 y = − 2 ( 1.2 ) \frac{d^2y}{dt^2}+\frac{dy}{dt}-2y=-2 (1.2) dt2d2y+dtdy−2y=−2(1.2) 方程(1)的特征方程 λ 2 + λ − 2 = 0 \lambda^2+\lambda-2=0 λ2+λ−2=0 ,解得 λ 1 = 1 , λ 2 = − 2 \lambda_1=1,\lambda_2=-2 λ1=1,λ2=−2 ,故方程(1)对应的齐次方程通解为 y = C 1 e t + C 2 e − 2 t . y=C_1e^t+C_2e^{-2t}. y=C1et+C2e−2t. 设 y 1 ( t ) = a t e t y_1(t)=ate^t y1(t)=atet 为方程(1.1)的特解,代入可解得 a = 1 a=1 a=1 ,故 y 1 ( t ) = t e t y_1(t)=te^t y1(t)=tet 。 y 2 ( t ) = 1 y_2(t)=1 y2(t)=1 显然为方程(1.2)的特解,可得方程(1)的特解为 t e t + 1 te^t+1 tet+1 ,故方程(1)的通解为 y = C 1 e t + C 2 e − 2 t + t e t + 1 y=C_1e^t+C_2e^{-2t}+te^t+1 y=C1et+C2e−2t+tet+1,把 t t t 换回 x x x ,可得原方程的通解为 y = C 1 x + C 2 x 2 + x l n x + 1. y=C_1x+\frac{C_2}{x^2}+xlnx+1. y=C1x+x2C2+xlnx+1.
呼~,内容真多,不过记忆的较多,难度不太大。中间还漏了全微分方程的解法,因为曲线积分还没学(惭愧)。接下来就还剩下级数、曲线曲面积分和空间解析几何以及前面落下的应用部分了,稳扎稳打吧。