[toc]
原章为微分方程,但事实上方程的两边不一定是微分
起始方程或微分方程解的过程中,方程两边几乎都会出现积分符号,而非只有微分(但主要用的是不定积分而不用定积分)
所以我这里改名为微分积分方程
,简称微积方程
,以体现这章是微分和积分的混合应用
在微分方程中,微分的写法往往比导数的写法更好用,不写y′写dxdy
比如:
- 经常会交换dx、dy(分离变量时)
- 高阶函数便于写清是谁导谁
- 便于两端同时积分,变为积分方程的形式
- 概念
- 【微分方程】:表示未知函数、未知函数导数、自变量之间的关系的方程,F(x,y,y′,⋯,y(n))=0
- 【微分方程的阶】:微分方程所出现的未知函数的最高阶导数的阶数
- 【微分方程的通解】:微分方程的解中含有任意常数,且其个数与微分方程的阶数相同(原因:每积分一次会有一个常数C出来)
- 【微分方程的特解】:确定通解中的任意常数(使用初值条件能求出常数)
- 【微分方程的积分曲线】:微分方程的解的图形,是一条曲线
定义
- 一般的一阶微分方程:y′=f(x,y),有时也写成对称形式:P(x,y)dx+Q(x,y)dy=0(dxdy=−Q(x,y)P(x,y))
- 可分离变量微分方程:如果一个一阶微分方程能写成g(y)dy=f(x)dx,那么原方程就称为可分离变量的微分方程
概念
- 隐式解:用隐式给出方程的解,即G(y)=F(x)+C,而非y=F(x)+C
- 隐式通解:含有任意常量的隐式解
- 齐次方程:如果一阶微分方程可化成dxdy=φ(xy)的形式,那么就称这方程为齐次方程
- dxdy=a2x+b2y+c2a1x+b1y+c1当c=c1=0时是齐次,否则不是齐次的当不是齐次时,可以使用变换{x=X+hy=Y+k,使方程组{a1h+b1k+c1=0a2h+b2k+c2=0那么方程组可化为齐次方程dXdY=a2X+b2Ya1X+b1Y
- 一阶非齐次线性微分方程
- dxdy+P(x)y=Q(x),(y、y′都是一次的方程,而不是说图像是线性的。而线性方程是未知数都是一次的方程,图像是线性的)
- 一阶 齐次 线性微分方程
- 当Q(x)≡0时,为齐次线性方程
- 一阶 齐次 线性微分方程求法
- 先分离变量,得ydy=−P(x)dx
- 再两端积分,得通解y=Ce−∫P(x)dx
- 一阶非齐次线性微分方程求法
- 用
常数变易法
求,即把齐次通解中的C换成x的未知函数u(x),得y=ue−∫P(x)dx - 用上式求dxdy,再代y和dxdy回原式,得通解
- 定义
- dxdy+P(x)y=Q(x)yn (n=0,1)
- 当n=0,1时,是齐次线性微分方程
- 可化为线性方程
- 伯努利方程是非线性的,但通过变换类的代换,可以把它化为线性
- 过程:较复杂、略
思路:通过代换把高阶微分方程化成低阶方程来求解
下面介绍三种容易降价的高阶微分方程
方法:两边同积分即可降一阶,∫dxndnydx=y(n−1)=∫f(x)dx+C
即F(x,y,y′,y′′)=0中没了变量y
方法:设y′=p,则y′′=p′=f(x,y)。算完后还原回来又得到一个一阶微分方程。即要算两次一阶微分方程
原理(核心思想):映射关系,y′′→p′y′→px →x
即F(x,y,y′,y′′)=0中没了变量x
方法:设y′=p,这里有两种转化方法(用哪种理解都行)
- y′′=dxdp=dydpdxdy=dydpp=f(y,p)
- y′′=dxdp=y′1dydp=dydpp=f(y,p)
- 有阻尼情况下的物体自由振动的微分方程
- 强迫振动的微分方程
- 串联电路的振荡方程
讨论二阶齐次线性方程y′′+P(x)y′+Q(x)y=0
概念
二阶线性微分方程:像这样的y′′+P(x)y′+Q(x)y=f(x)
齐次的:等式右端f(x)≡0
线性相关:设y1(x),y2(x),⋯,yn(x)为定义在区间I上的n个函数,如果存在n个不全为零的常数k1,k2,⋯,kn,使得当x∈I时有恒等式k1y1+k2y2+⋯+knyn≡0成立,那么就称这n个函数在区间I上线性相关
例如1、cos2、sin2线性相关(1−cos2−sin2=0)
线性无关:否则称线性无关
例如1,x,x2线性无关
判断是否线性相关:当函数只有两个时,只需判断它们的比是否为常数
定理
【定理1】:如果y1(x)与y2(x)是方程y′′+P(x)y′+Q(x)y=0的两个特解, 那么y=C1y1(x)+C2y2(x)也是其特解
【定理2】:如果y1(x)与y2(x)是方程y′′+P(x)y′+Q(x)y=0的两个线性无关特解,那么y=C1y1(x)+C1y2(x)就是方程的通解
【定理2推论】:如果y1(x),y2(x),⋯,yn(x)是n阶齐次线性方程y(n)+a1(x)y(n−1)+⋯+an−1(x)y′+an(x)y=0的n个线性无关的解,那么,此方程的通解为y=C1y1(x)+C2y2(x)+⋯+Cnyn(x)
【定理3】:设y∗(x)是二阶非齐次线性方程y′′+P(x)y′+Q(x)y=f(x)的一个特解Y(x)是与之对应的齐次方程的通解,则y=Y(x)+y∗(x)是二阶非齐次线性微分方程的通解
(可推广到n阶,略)
【定理4】解的叠加原理:设非齐次线性方程的右端f(x)是两个函数之和,即y′′+P(x)y′+Q(x)y=f1(x)+f2(x)而y1∗(x)与y2∗(x)分别是方程y′′+P(x)y′+Q(x)y=f1(x)与y′′+P(x)y′+Q(x)y=f2(x)的特解,则y1∗(x)+y2∗(x)就是原方程的特解
(可推广到n阶,略)
很麻烦,应该不考。方法比一阶非齐次线性方程的解法(常数变易法)更复杂:需要多变易一个常数,多处理一个未知函数,需要一些新手段
定义
- 二阶常系数齐次线性微分方程:P(x)、Q(x)位置均全为常数,即y′′+py′+qy=0
- 二阶变系数齐次线性微分方程:P(x)、Q(x)位置不全为常数
求解方法与思路,特征方程法的由来
- (1) 当r为常数时,y=erx与它的各阶导数都只相差一个常数因子。因此看能否令y=erx满足常系数齐次线性微分方程
- (2) y=erx、y′=rerx、y′′=r2erx代入到y′′+py′+qy=0得到(r2+pr+q)erx=0(erx=0)
- (3) 其中特征方程为r2+pr+q=0,特征根解:r1,2=2−p±p2−4q
- (4) 情况判断
- (1) 当Δ=p2−4q>0时,r1,r2是两个不相等的实根能得知两个特解根微分方程的通解为y=C1er1x+C2er2x
- (2) 当Δ=p2−4q=0时,r1,r2是两个 相等的实根能得知一个特解根,还缺一个y2,且要求y1y2不是常数。设y2=er1xu(x)代入得u′′=0,即可设u=x,y2为xer1x微分方程的通解为y=(C1+C2x)er1x
- (3) 当Δ=p2−4q<0时,r1,r2是一对 共轭复根需要通过欧拉公式能把结果替换为非复数的版本:y1=eaxcosβx,y2=eaxsinβx微分方程的通解为y=eax(C1cosβx+C2sinβx)
总结
特征方程r2+pr+q=0的两个根 | 微分方程y′′+py′+qy=0的通解 | 解算复杂度 |
---|
两个不相等的实根,r1,2=2−p±p2−4q=2−p±Δ | y=C1er1x+C2er2x | 简单 |
两个 相等的实根,r1=r2=−2p | y=C1er1x+C2xer1x y=(C1+C2x)er1x | 多一步,稍麻烦 但套公式简单 |
一对 共轭复根,r1,2=α±βi=2−p±Δ=2−p±i−Δ | y=eax(C1cosβx+C2sinβx) | 稍麻烦 |
补充:化复数特解为实数特解的方法
- (1) 欧拉公式,eiθ=cosθ+i⋅sinθ,y1,2=e(α±βi)x=eax(∗/)eβxi=eax(cosβx±isinβx)
- (2) 共轭关系 + 解叠加原理,y3=21(y1+y2)=eaxcosβx,y4=2i1(y1−y2)=eaxsinβx
- 其中α=−2p,β=2Δ
本节只介绍f(x)取两种常见形式时求y∗的方法
这种方法的特点是不用积分就可求出y∗来,它叫做待定系数法
变系数的线性微分方程比较难解,但有些特殊的变系数微分方程可以通过常量代换
化为常系数线性微分方程
例如欧拉方程就是其中的一种
- 注意:欧拉公式不等于欧拉方程,前者是eiθ=cosθ+isinθ
- 概念:形如xny(n)+p1xn−1y(n−1)+⋯+pn−1xy′+pny=f(x)
- 解算方法:略(备选章懒得看了)
(这章是应用层面上的一些举例)
微分方程组
常系数线性微分方程组
表达的写法有点多,这里捋一下
阶数 | 通用写法 |
---|
n 阶微分方程 | F(x,y,y′,⋯,y(n))=0 |
一阶微分方程 | F(x,y,y′)=0 |
高阶微分方程 | (n≥2) |
—————— | —————— |
可分离 | 表示二元、且可分离 |
线性 | 表示y(n)和y的次数都是1,且彼此间不相乘而分别均为单项式的因子 |
齐次 | 等式右端恒等于0,即f(x)不单独存在 |
常系数 | 等号左侧的P(x)、Q(x)等均为常数 |
排列顺序:从特殊到一般
微分方程 | 导数写法 | 微分写法 | 对称写法 / 分离变量 |
---|
一阶普通导数方程 | y′+P(x)=0 | dxdy+P(x)=0 | |
一阶微分方程 (可分离变量) | y′=g(y)f(x) | dxdy=g(x)f(y) | g(y)dy=f(x)dx |
一阶微分方程 (齐次方程) | y′=φ(xy)=φ2(xy)φ1(xy) | dxdy=φ(xy)=φ2(xy)φ1(xy) | ———— |
一阶微分方程 (可化为齐次) | y′=a2x+b2y+c2a1x+b1y+c1等 | dxdy=a2x+b2y+c2a1x+b1y+c1等 | ———— |
一阶线性微分方程 (齐次方程) | y′+P(x)y=0 | dxdy+P(x)y=0 | ydy=−P(x)dx |
一阶线性微分方程 (非齐次方程) | y′+P(x)y=Q(x) | dxdy+P(x)y=Q(x) | ———— |
一阶微分方程 (通用形式) | y′=Q(x,y)P(x,y)=ξ(x,y) | dxdy=Q(x,y)P(x,y)=ξ(x,y) | Q(x,y)dy=P(x,y)dx |
伯努利方程 (可化为线性) | y′+P(x)y=Q(x)yn(n=0,1) | dxdy+P(x)y=Q(x)yn(n=0,1) | |
———————— | ———————— | ———————— | ———————— |
二级普通导数方程 | y′′+P(x)=0 | dx2d2y+P(x)=0 | |
二阶微分方程 (可降阶的) | y′′=f(x,y′) | 略 | |
二阶微分方程 (可降阶的常系数) | y′′=f(y,y′) | 略 | |
二阶常系数线性微分 (齐次方程) | y′′+py′+qy=0 | 略 | |
二阶常系数线性微分 (非齐次方程) | f(x)=eλxPm(x)型 | 略 | |
二阶常系数线性微分 (非齐次方程) | f(x)=eλx[Pl(x)cosωx+Qn(x)sinωx]型 | 略 | |
二阶常系数线性微分 (非齐次方程) | y′′+py′+qy=f(x) 一般型 | 略 | |
二阶线性微分方程 (齐次方程) | y′′+P(x)y′+Q(x)y=0 | 略 | |
二阶线性微分方程 (非齐次方程) | y′′+P(x)y′+Q(x)y=f(x) | 略 | |
———————— | ———————— | ———————— | ———————— |
高阶微分方程 (可降阶的) | y(n)=f(x) | 略 | |
n阶常系数线性微分 (齐次方程) | y(n)+p1y(n−1)+⋯+pny=0 | 略 | |
n阶变系数线性微分 (欧拉方程) | xny(n)+p1xn−1y(n−1)+⋯+pn−1xy′+pny=f(x) | 略 | |
n阶线性微分方程 | y(n)+a1(x)y(n−1)+⋯+an(x)y=f(x) | 略 | |
排列顺序:从特殊到一般,从易解到困难
微分方程 | 微分写法 | 求解方法 | 通解形式 |
---|
一阶普通导数方程 | dxdy+P(x)=0 | 导数分离移项 + 同时积分(简单) | y=F(x)+C |
一阶微分方程 (可分离变量) | dxdy=g(x)f(y) | 分离变量移项 + 同时积分(简单) | G(y)=F(x)+C |
一阶微分方程 (齐次方程) | dxdy=φ(xy)=φ2(xy)φ1(xy) | 中间函数u(x)=xy (转化为可分离变量方程)(稍麻烦) | y=Cexy??? |
一阶微分方程 (可化为齐次) | dxdy=a2x+b2y+c2a1x+b1y+c1等 | 中间变量x=X+h,y=Y+k (转化为齐次方程)(稍稍麻烦) | (同上) |
一阶线性微分方程 (齐次必可分离) | dxdy+P(x)y=0 | 分离变量移项 + 同时积分(简单) 或,套公式(简单) | $\ln |
一阶线性微分方程 (非齐次方程) | dxdy+P(x)y=Q(x) | 齐次结果 + 常数变易法 + 求导 + 代回 或,套公式(都很麻烦) 或,求导公式法(推荐,较舒服) | y=e−∫P(x)dx[∫Q(x)e∫P(x)dxdx+C] 或,写成和的形式 |
一阶微分方程 (通用形式) | dxdy=Q(x,y)P(x,y)=ξ(x,y) | 无 | 无 |
伯努利方程 (可化为线性) | dxdy+P(x)y=Q(x)yn(n=0,1) | 两端除以yn再代换z=y1−n 可转换为线性方程(较麻烦) | 无 |
———————— | ———————— | ———————— | ———————— |
二级普通导数方程 | dx2d2y+P(x)=0 | 导数分离移项 + 同时积分(简单) | y=C1f(x)+C2g(x) |
二阶微分方程 (可降阶的) | y′′=f(x,y′),即少了y | y′′→p′,y′→p,x →x 变去变回后各算一次一阶(简单) | y=∫φ(x,C1)dx+C2 |
二阶微分方程 (可降阶的常系数) | y′′=f(y,y′),即少了x | y′′→y′dydy′ 算两次一阶(简单) | ∫φ(y,C1)dy=x+C2 |
二阶常系数线性微分 (齐次方程) | y′′+py′+qy=0 | 特征方程法(看情况,简单/简单/难) 或,套公式(麻烦) | y=C1er1x+C2er2x 或,y=C1er1x+C2xer1x 或,y=eax(C1cosβx+C2sinβx) |
二阶常系数线性微分 (非齐次方程) | f(x)=eλxPm(x)型 | | |
二阶常系数线性微分 (非齐次方程) | f(x)=eλx[Pl(x)cosωx+Qn(x)sinωx]型 | | |
二阶常系数线性微分 (非齐次方程) | y′′+py′+qy=f(x) 一般型 | 无 | 无 |
二阶线性微分方程 (齐次方程) | y′′+P(x)y′+Q(x)y=0 | 特解求法:两个线性无关特解套公式 (简单) | y=C1y1(x)+C2y2(x) |
二阶线性微分方程 (非齐次方程) | y′′+P(x)y′+Q(x)y=f(x) | 特解求法:齐次通解 加 非齐次特解 齐次结果 + 常数变易法 + 求导 + 代回 (特解法稍麻烦,而变易法很麻烦) | y=Y(x)+y∗(x) |
———————— | ———————— | ———————— | ———————— |
高阶微分方程 (可降阶的) | y(n)=f(x) | 正常两端积分多次即可(简单) | 无 |
n阶常系数线性微分 (齐次方程) | y(n)+p1y(n−1)+⋯+pny=0 | 特征方程法(看情况,简单/简单/难) 或,套公式(麻烦) | 略 |
n阶变系数线性微分 (欧拉方程) | xny(n)+p1xn−1y(n−1)+⋯+pn−1xy′+pny=f(x) | 常量代换化为常系数线性微分方程 | |
n阶线性微分方程 | y(n)+a1(x)y(n−1)+⋯+an(x)y=f(x) | 无 | 无 |
分离变量方法,可以直接两端加不定积分号的原因:
g(y)dyg(φ(x))φ′(x)dx∫g(y)dy=f(x)dx,设y=φ(x)=f(x)dx=∫f(x)dx
齐次微分方程,齐次方程中用中间函数u(x)替换xy而非用中间变量
dxdy=dxd(ux)=dxudx+xdu=u+xdxdu=u
非齐次线性微分方程,常数变易法详解
就是把常数变成f(x)的操作,从而从特殊到一般
该技巧并不是一定可行,只是在这里求齐次线性方程时恰好可行(是一种求解技巧)
书上并没有写真正意义上的求解方法,想必非常麻烦
补充说明: e∫P(x)dx亦称作积分因子
注意项:求伯努利方程时,dxdy1−n不是dydy1−n,要用分部积分法来算
解微分方程
- 验证函数是微分方程的解:从解开始求导会快点
- 非齐次线性微分方程:常数变易法