AM@微分方程相关概念@线性微分方程@一阶线性微分方程的通解

abstract

  • AM@微分方程相关概念@线性微分方程@一阶线性微分方程的通解

preface

  • 经验表明,获得微分方程的一般性数学理论是困难的
    • 有少数类型的微分方程,比如线性微分方程具有一般的求解理论
    • 还有一些一阶的简单的微分方程类型

一般的微分方程@ n n n阶微分方程

  • 含有未知函数,未知函数的导数与自变量之间的关系的方程,称为微分方程;

  • 未知函数导数的最高阶数称为该微分方程的

    • 方程 F ( x , y , y ′ , ⋯   , y n ) = 0 F(x,y,y',\cdots,y^{n})=0 F(x,y,y,,yn)=0(1)或方程 y ( n ) y^{(n)} y(n)= f ( x , y , y ′ , ⋯   , y ( n − 1 ) ) f(x,y,y',\cdots,y^{(n-1)}) f(x,y,y,,y(n1))(1-1)成为 n n n阶微分方程(形式(1-1)是一种常见的微分方程形式)

    • 其中 x , y , y ′ ⋯   , y ( n − 1 ) x,y,y'\cdots,y^{(n-1)} x,y,y,y(n1)可以没有,但必须有最高阶导数 y ( n ) y^{(n)} y(n),并且其系数为常数(通常化为常数1)

    • n = 1 n=1 n=1时,方程(1)称为一阶微分方程

常微分方程👺

  • 未知函数是一元函数的微分方程称为常微分方程

    • 例如 y = f ( x ) y=f(x) y=f(x)
  • 注意常微分方程和常系数微分方程的区别,两个是完全不同的概念

偏微分方程

  • 如果未知函数式多元函数,则成为偏微分方程
  • 例如 y = f ( x , y ) y=f(x,y) y=f(x,y)

微分方程的解

  • y = ϕ ( x ) y=\phi(x) y=ϕ(x)在区间 I = ( a , b ) I=(a,b) I=(a,b)上连续且 n n n阶可导,使得 F ( x , ϕ ( x ) , ϕ ′ ( x ) , ⋯   , ϕ ( n ) ( x ) ) = 0 F(x,\phi(x),\phi'(x),\cdots,\phi^{(n)}(x))=0 F(x,ϕ(x),ϕ(x),,ϕ(n)(x))=0,(2)即式(1)恒成立,则称 y = ϕ ( x ) y=\phi(x) y=ϕ(x)为该微分方程在区间 I I I上的一个

隐式解

  • 若关系 F ( x , y ) = 0 F(x,y)=0 F(x,y)=0(3)确定的隐函数 y = ϕ ( x ) y=\phi(x) y=ϕ(x)是(1)的解,则称(3)是(1)的隐式解(隐式通解)
  • 例如: x 2 + y 2 = 1 x^2+y^2=1 x2+y2=1(4)是一阶微分方程 y ′ = − x y y'=-\frac{x}{y} y=yx(4-1)的隐式解
    • 对(4)求导, 2 x + 2 y y ′ = 0 2x+2yy'=0 2x+2yy=0,变形得(4-1)式

通解和特解

  • 若含有 n n n独立任意常数的函数 y = ϕ ( x , C 1 , ⋯   , C n ) y=\phi(x,C_1,\cdots,C_n) y=ϕ(x,C1,,Cn), x ∈ I x\in{I} xI(5) n n n阶微分方程(1)的解,则称(5)是(1)的通解
  • 不含任意常数的解称为特解
  • 详见函数线性相关性

初始条件

  • 关系式: y ( i ) ( x 0 ) = y 0 ( i ) y^{(i)}(x_0)=y_0^{(i)} y(i)(x0)=y0(i), ( i = 0 , 1 , 2 , ⋯   , n ) (i=0,1,2,\cdots,n) (i=0,1,2,,n)称为 n n n阶微分方程的初始条件
    • 其中 y 0 i y_{0}^{i} y0i, ( i = 0 , 1 , 2 , ⋯   , n ) (i=0,1,2,\cdots,n) (i=0,1,2,,n) n n n个给定的数

初值问题

  • 一般地,由初始条件确定通解中的任意常数,就得到相应的一个特解

  • 上述方式确定特解的问题称为初值问题

  • 利用初始条件定义特解:微分方程中,满足初始条件的称为特解

微分方程的积分曲线

  • 微分方程的解(函数)的图形是一条曲线,称为微分方程的积分曲线
  • 一阶微分方程的初值问题
    • y ′ = f ( x , y ) y'=f(x,y) y=f(x,y); y ′ ∣ x = x 0 = y 0 y'|_{x=x_0}=y_0 yx=x0=y0的几何意义,就是求通过点 ( x 0 , y 0 ) (x_0,y_0) (x0,y0)的那条积分曲线

线性微分方程👺

  • 方程 ∑ i = 0 n a i ( x ) y ( i ) \sum_{i=0}^{n}a_{i}(x)y^{(i)} i=0nai(x)y(i)= f ( x ) f(x) f(x),(1);式(1)展开写为: y ( n ) + a 1 ( x ) y ( n − 1 ) + ⋯ + a n − 1 ( x ) y ′ + a n ( x ) y y^{(n)}+a_1(x)y^{(n-1)}+\cdots+a_{n-1}(x)y^{'}+a_n(x)y y(n)+a1(x)y(n1)++an1(x)y+an(x)y= f ( x ) f(x) f(x)(1-1);该方程称为** n n n阶线性微分方程**

  • 若式(1)中 f ( x ) = 0 f(x)=0 f(x)=0,此时方程(1)作 ∑ i = 1 n a i ( x ) y ( i ) = 0 \sum_{i=1}^{n}a_{i}(x)y^{(i)}=0 i=1nai(x)y(i)=0(2),称为** n n n阶线性齐次微分方程**,并且称(2)是(1)对应的齐次方程

  • 若式(1)中 f ( x ) ≠ 0 f(x)\neq{0} f(x)=0,其中系数 a i ( x ) a_i(x) ai(x)是已知函数,并假设 a i ( x ) a_i(x) ai(x), f ( x ) f(x) f(x)在某个区间 ( a , b ) (a,b) (a,b)连续,则式(1)为** n n n阶线性非齐次微分方程**, f ( x ) f(x) f(x)称为自由项

  • 等式左边没有常数项,常数项可归到右边的自由项

一阶线性微分方程👺

  • 形如 d y d x + P ( x ) y = Q ( x ) \frac{\mathrm{d}y}{\mathrm{d}x}+P(x)y=Q(x) dxdy+P(x)y=Q(x)(1)的方程称为一阶线性微分方程

通解

  • y = C e − ∫ P ( x ) d x + e − ∫ P ( x ) d x ⋅ ∫ Q ( x ) e ∫ P ( x ) d x d x y=Ce^{-\int{P(x)}\mathrm{d}x}+{e^{-\int{P(x)}\mathrm{d}x}}\cdot\int{Q(x)e^{\int{P(x)\mathrm{d}x}}\mathrm{d}x} y=CeP(x)dx+eP(x)dxQ(x)eP(x)dxdx

  • 另一种表示方式: y = exp ⁡ ( − ∫ P ( x ) d x ) ( ∫ [ Q ( x ) exp ⁡ ( ∫ P ( x )    d x )    d x ] + C ) y=\exp\left(-\int{P(x)\mathrm{d}x}\right) \left(\int{[Q(x)}\exp\left(\int{P(x)\;\mathrm{d}x}\right)\;\mathrm{d}x]+C \right) y=exp(P(x)dx)([Q(x)exp(P(x)dx)dx]+C)

推导

一阶齐次和非齐次线性微分方程

  • Q ( x ) ≡ 0 Q(x)\equiv{0} Q(x)0,则称 d y d x + P ( x ) y = 0 \frac{\mathrm{d}y}{\mathrm{d}x}+P(x)y=0 dxdy+P(x)y=0(2)为对应于(1)的齐次线性方程
  • Q ( x ) ≢ 0 Q(x)\not\equiv{0} Q(x)0(函数 Q ( x ) Q(x) Q(x)不总是取 0 0 0,这不同于函数 Q ( x ) ≠ 0 Q(x)\neq{0} Q(x)=0(不取0)),则称 d y d x + P ( x ) y = Q ( x ) ≢ 0 \frac{\mathrm{d}y}{\mathrm{d}x}+P(x)y=Q(x)\not\equiv{0} dxdy+P(x)y=Q(x)0(3)非齐次线性方程

一阶齐次线性微分方程的解

  • 方程(2)是可分离变量的: d y y \frac{\mathrm{d}y}{y} ydy= − P ( x ) d x -P(x)\mathrm{d}x P(x)dx(4),两边积分,得 ln ⁡ ∣ y ∣ \ln|y| lny= − ∫ P ( x ) d x + C 1 -\int{P(x)\mathrm{d}x}+C_1 P(x)dx+C1(4-1)
  • 两边取指数: ∣ y ∣ |y| y= e − ∫ P ( x ) d x + C 1 e^{-\int{P(x)\mathrm{d}x}+C_1} eP(x)dx+C1= e − ∫ P ( x ) d x ⋅ e C 1 e^{-\int{P(x)\mathrm{d}x}}\cdot{e^{C_1}} eP(x)dxeC1(5),为了便于书写复杂指数,使用 exp ⁡ x \exp{x} expx表示 e x e^{x} ex,则式(5)可以表示为 ∣ y ∣ |y| y= exp ⁡ ( − ∫ P ( x ) d x + C 1 ) \exp{(-\int{P(x)}\mathrm{d}x+C_1)} exp(P(x)dx+C1)= exp ⁡ ( C 1 ) ⋅ exp ⁡ ( − ∫ P ( x ) d x ) \exp(C_1)\cdot\exp{(-\int{P(x)}\mathrm{d}x)} exp(C1)exp(P(x)dx)
  • y = ± exp ⁡ ( C 1 ) ⋅ exp ⁡ ( − ∫ P ( x ) d x ) y=\pm{\exp(C_1)\cdot\exp{(-\int{P(x)}\mathrm{d}x)}} y=±exp(C1)exp(P(x)dx)= C exp ⁡ ( − ∫ P ( x ) d x ) C\exp{(-\int{P(x)}\mathrm{d}x)} Cexp(P(x)dx)(6),其中 C = ± exp ⁡ ( C 1 ) C=\pm{\exp(C_1)} C=±exp(C1)
  • 式(6)就是方程(3)的通解

一阶非齐次线性微分方程的解

  • 显然方程(2)是方程(1)的特殊情况,两者存在一定的联系
  • 对方程(3)进行变形: d y y \mathrm{d}y\over{y} ydy= ( − P ( x ) + 1 y Q ( x ) ) d x (-P(x)+\frac{1}{y}Q(x))\mathrm{d}x (P(x)+y1Q(x))dx(7),两边积分 ln ⁡ ∣ y ∣ \ln|y| lny= ∫ ( − P ( x ) + 1 y Q ( x ) ) d x + C 1 \int{(-P(x)+\frac{1}{y}Q(x))\mathrm{d}x}+C_1 (P(x)+y1Q(x))dx+C1= ∫ ( − P ( x ) + 1 y Q ( x ) ) d x + ln ⁡ ∣ C ∣ \int{(-P(x)+\frac{1}{y}Q(x))\mathrm{d}x}+\ln{|C|} (P(x)+y1Q(x))dx+lnC(7-1),其中 C 1 = ln ⁡ ∣ C ∣ C_1=\ln{|C|} C1=lnC, C 1 C_1 C1可以取任何常数,但为了得到 y y y,我们要对(7-1)两边取指数, exp ⁡ ln ⁡ ∣ C ∣ \exp{\ln|C|} explnC= ∣ C ∣ |C| C,是一个简单的值,继续展开(7-1),得 − ∫ P ( x ) d x + ∫ 1 y Q ( x ) d x + ln ⁡ ∣ C ∣ -\int{P(x)}\mathrm{d}x+\int{\frac{1}{y}Q(x)\mathrm{d}x}+\ln{|C|} P(x)dx+y1Q(x)dx+lnC(7-2)
  • 取指数, ∣ y ∣ |y| y= exp ⁡ ( − ∫ P ( x ) d x + ∫ 1 y Q ( x ) d x + ln ⁡ ∣ C ∣ ) \exp{(-\int{P(x)}\mathrm{d}x+\int{\frac{1}{y}Q(x)\mathrm{d}x}+\ln{|C|})} exp(P(x)dx+y1Q(x)dx+lnC)(8),即 y = ± ∣ C ∣ exp ⁡ ( − ∫ P ( x ) d x ) exp ⁡ ( ∫ 1 y Q ( x ) d x ) y=\pm{|C|\exp(-\int{P(x)\mathrm{d}x})}\exp{(\int{\frac{1}{y}Q(x)\mathrm{d}x})} y=±Cexp(P(x)dx)exp(y1Q(x)dx)= C exp ⁡ ( − ∫ P ( x ) d x ) exp ⁡ ( ∫ 1 y Q ( x ) d x ) C{\exp(-\int{P(x)\mathrm{d}x})}\exp{(\int{\frac{1}{y}Q(x)\mathrm{d}x})} Cexp(P(x)dx)exp(y1Q(x)dx)(9)
    • 其中 C = ± ∣ C ∣ C=\pm{|C|} C=±C
      • C exp ⁡ ( − ∫ P ( x ) d x ) C{\exp(-\int{P(x)\mathrm{d}x})} Cexp(P(x)dx),这就是式(6),即式(9)包含一阶齐次线性微分方程的通解
      • T = exp ⁡ ( ∫ 1 y Q ( x ) d x ) T=\exp{(\int{\frac{1}{y}Q(x)\mathrm{d}x})} T=exp(y1Q(x)dx)
  • T T T是关于 x x x的函数( y y y是关于 x x x的一元函数,所以 T T T x x x的函数,可以表示为 T ( x ) T(x) T(x)),
    • 式(9)也是关于 x x x的函数,但其表达式包含 y y y,下面的工作是化去式等号右边的 y y y,使之仅含有已知的关于 x x x的函数式
  • 因此,比较式(9),(6),利用常数变易法,将方程(6)中的 C C C变易为 x x x待定函数 C ( x ) C(x) C(x),使之满足方程(1),从而求出 C ( x ) C(x) C(x)(其表示的是 C T ( x ) CT(x) CT(x)
    • y = C ( x ) exp ⁡ ( − ∫ P ( x ) d x ) y=C(x)\exp{(-\int{P(x)}\mathrm{d}x)} y=C(x)exp(P(x)dx);(10),对其两边求导
    • d y d x \frac{\mathrm{d}y}{\mathrm{d}x} dxdy= C ′ ( x ) exp ⁡ ( − ∫ P ( x ) d x ) C'(x)\exp{(-\int{P(x)}\mathrm{d}x)} C(x)exp(P(x)dx)+ C ( x ) exp ⁡ ( − ∫ P ( x ) d x ) ( − P ( x ) ) C(x)\exp{(-\int{P(x)}\mathrm{d}x)}(-P(x)) C(x)exp(P(x)dx)(P(x))
      • = C ′ ( x ) exp ⁡ ( − ∫ P ( x ) d x ) C'(x)\exp{(-\int{P(x)}\mathrm{d}x)} C(x)exp(P(x)dx)- C ( x ) P ( x ) exp ⁡ ( − ∫ P ( x ) d x ) ) C(x)P(x)\exp{(-\int{P(x)}\mathrm{d}x)}) C(x)P(x)exp(P(x)dx))(11)
    • 将(10),(11)代入方程(1),得 C ′ ( x ) exp ⁡ ( − ∫ P ( x ) d x ) C'(x)\exp{(-\int{P(x)}\mathrm{d}x)} C(x)exp(P(x)dx)- C ( x ) P ( x ) exp ⁡ ( − ∫ P ( x ) d x ) ) C(x)P(x)\exp{(-\int{P(x)}\mathrm{d}x)}) C(x)P(x)exp(P(x)dx))+ P ( x ) C ( x ) exp ⁡ ( − ∫ P ( x ) d x ) P(x)C(x)\exp{(-\int{P(x)}\mathrm{d}x)} P(x)C(x)exp(P(x)dx)= Q ( x ) Q(x) Q(x);
    • 即得 C ′ ( x ) exp ⁡ ( − ∫ P ( x ) d x ) C'(x)\exp{(-\int{P(x)}\mathrm{d}x)} C(x)exp(P(x)dx)= Q ( x ) Q(x) Q(x),(11-1),整理可得 C ′ ( x ) C'(x) C(x)= Q ( x ) exp ⁡ ( ∫ ( P ( x ) d x ) ) Q(x)\exp{(\int(P(x)\mathrm{d}x))} Q(x)exp((P(x)dx))(11-2),两边积分,就求得函数 C ( x ) C(x) C(x)= ∫ [ Q ( x ) exp ⁡ ( ∫ ( P ( x ) d x ) ) ] d x + C \int{[Q(x)\exp{(\int(P(x)\mathrm{d}x))}]\mathrm{d}x}+C [Q(x)exp((P(x)dx))]dx+C(12)
    • 将式(12)代入到(10),得方程(3)的通解: y = [ ∫ [ Q ( x ) exp ⁡ ( ∫ ( P ( x ) d x ) ) ] d x + C ] y=[\int{[Q(x)\exp{(\int(P(x)\mathrm{d}x))}]\mathrm{d}x} +C] y=[[Q(x)exp((P(x)dx))]dx+C] ⋅ \cdot [ exp ⁡ ( − ∫ P ( x ) d x ) ] [\exp{(-\int{P(x)}\mathrm{d}x)]} [exp(P(x)dx)](13),习惯上把指数式放在前面,即 y y y= [ exp ⁡ ( − ∫ P ( x ) d x ) ] [\exp{(-\int{P(x)}\mathrm{d}x)]} [exp(P(x)dx)] ⋅ \cdot [ ∫ [ Q ( x ) exp ⁡ ( ∫ ( P ( x ) d x ) ) ] d x + C ] [\int{[Q(x)\exp{(\int(P(x)\mathrm{d}x))}]\mathrm{d}x} +C] [[Q(x)exp((P(x)dx))]dx+C](13-1),或者展开成两项和:
      • y y y= C [ exp ⁡ ( − ∫ P ( x ) d x ) ] C[\exp{(-\int{P(x)}\mathrm{d}x)]} C[exp(P(x)dx)]+ [ exp ⁡ ( − ∫ P ( x ) d x ) ] [\exp{(-\int{P(x)}\mathrm{d}x)]} [exp(P(x)dx)] [ ∫ [ Q ( x ) exp ⁡ ( ∫ ( P ( x ) d x ) ) ] d x ] [\int{[Q(x)\exp{(\int(P(x)\mathrm{d}x))}]\mathrm{d}x}] [[Q(x)exp((P(x)dx))]dx](14)
      • 式(14)中第一项(包含任意常数)就是齐次方程(2)的解,第二项(不包含任意常数)是非齐次方程(3)的一个特解
    • 总之,一阶非齐次线性微分方程的通解等于对应的齐次方程的通解加上非齐次方程的一个特解之和;这一点和线性代数中非齐次线性方程组解的结构结论相仿

小结

  • 一阶线性微分方程的解

    方程通解
    y ′ + p ( x ) y = 0 y'+p(x)y=0 y+p(x)y=0 y = C e − ∫ p ( x ) d x y=Ce^{-\int{p(x)}\mathrm{d}x} y=Cep(x)dx
    y ′ + p ( x ) y = q ( x ) y'+p(x)y=q(x) y+p(x)y=q(x) y = C e − ∫ p ( x ) d x + e − ∫ P ( x ) d x ⋅ ∫ q ( x ) e ∫ p ( x ) d x d x y=Ce^{-\int{p(x)}\mathrm{d}x}+ {e^{-\int{P(x)}\mathrm{d}x}}\cdot\int{q(x)e^{\int{p(x)\mathrm{d}x}}\mathrm{d}x} y=Cep(x)dx+eP(x)dxq(x)ep(x)dxdx= e − ∫ p ( x ) d x [ C + ∫ q ( x ) e ∫ p ( x ) d x d x ] e^{-\int{p(x)}\mathrm{d}x}[C+\int{q(x)e^{\int{p(x)\mathrm{d}x}}\mathrm{d}x}] ep(x)dx[C+q(x)ep(x)dxdx]
  • 步骤

    • 观察线性方程是否为一阶线性常微分方程(或变形),若是,则将其化为标准形,以便看出 p ( x ) , q ( x ) p(x),q(x) p(x),q(x)套用公式
    • 根据方程是齐次还是非齐次的,套用对应的公式
      • 无论是否为齐次的,都建议先计算: e − ∫ p ( x ) d x e^{-\int{p(x)}\mathrm{d}x} ep(x)dx= g ( x ) g(x) g(x)(1)
      • 从而齐次方程的解就是 y = C g ( x ) y=Cg(x) y=Cg(x)
      • 而非齐次方程的解就是 y = g ( x ) [ C + ∫ q ( x ) ⋅ 1 g ( x ) d x ] y=g(x)[C+\int{q(x)\cdot{\frac{1}{g(x)}}\mathrm{d}x}] y=g(x)[C+q(x)g(x)1dx]
    • 也可以先计算 ∫ p ( x ) d x \int{p(x)}\mathrm{d}x p(x)dx= t ( x ) t(x) t(x)`(2)
    • 齐次方程的解就是 y = C e − t ( x ) y=Ce^{-t(x)} y=Cet(x)
    • 非齐次线性方程的解就是 y = e − t ( x ) [ C + ∫ q ( x ) ⋅ e t ( x ) d x ] y=e^{-t(x)}[C+\int{q(x)\cdot{e^{t(x)}}}\mathrm{d}x] y=et(x)[C+q(x)et(x)dx]
  • y ′ = 2 x y + 2 x e x 2 y'=2xy+2xe^{x^2} y=2xy+2xex2
    • 分析可知这是一个一阶非齐次线性常微分方程,标准化为 y ′ − 2 x y y'-2xy y2xy= 2 x e x 2 2xe^{x^2} 2xex2
    • g ( x ) g(x) g(x)= e − ∫ ( − 2 x ) d x e^{-\int(-2x)\mathrm{d}x} e(2x)dx= e x 2 e^{x^2} ex2
    • 从而通解为: y = e x 2 ( C + ∫ 2 x e x 2 ⋅ e − x 2 d x ) y=e^{x^2}(C+\int{2xe^{x^2}}\cdot e^{-x^2}\mathrm{d}x) y=ex2(C+2xex2ex2dx)= ( x 2 + C ) e x 2 (x^2+C)e^{x^2} (x2+C)ex2
posted @   xuchaoxin1375  阅读(27)  评论(0编辑  收藏  举报  
相关博文:
阅读排行:
· 分享4款.NET开源、免费、实用的商城系统
· 全程不用写代码,我用AI程序员写了一个飞机大战
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· 记一次.NET内存居高不下排查解决与启示
· 白话解读 Dapr 1.15:你的「微服务管家」又秀新绝活了
历史上的今天:
2022-11-01 PT_随机变量函数的分布_随机变量线性函数的正态分布
2022-11-01 PT_二维随机变量:正态分布的可加性/一维随机变量函数与正态分布
2021-11-01 python_pandas入门(by offical document/reference)/loc和iloc操作/dataframe插入操作/pandas读取无表头的文件/查找某一列是否有某个值
点击右上角即可分享
微信分享提示