F(x,ϕ(x),ϕ′(x),ϕ′′(x),...,ϕ(n)(x))≡0
最高几阶导,称为微分方程的阶
考核范围仅限线性微分方程
微分算子复习
考研范围大概只用掌握以下几类微分方程
P(x)dx=P(y)dy
解法是显然的
dxdy=ϕ(xy)
dxdy=a1x+b1y+c1ax+by+c
通过换元 x′=x+A,y′=y+B 来将其转化为齐次
a1x+b1y+c1ax+by+c=a1x′+b1y′+(c1−a1A−b1B)ax′+by′+(c−aA−bB)
因此只需通过
{c−aA−bB=0c1−a1A−b1B=0
来唯一确定 A,B
这个二元方程有四种情况,唯一解,一维任意解,二维任意解,无解
- 唯一解对应 A,B 系数矩阵可逆,解空间为整个平面
- 一维任意解对应系数矩阵不可逆,但非零矩阵,同时 c1c=a1a,(c,c1) 在一维的解空间中
- 二维任意解对应系数矩阵为零矩阵,同时 c1=c=0,解空间为零点
- 无解是除以上以外的任意情况,(c,c1) 不在 (a,a1) 和 (b,b1) 张成的向量空间中
只有 1. 需要换元,2. 的微分方程右侧为常数,3. 不可能出现(分母为 0),4. 可以分类讨论以下
- 系数矩阵秩为 0,此时微分方程右侧为常数,或者 c1=0 时不可能
- 系数矩阵秩为 1,此时微分方程右侧可以化为 a1a+a1x+b1y+c1d
具体例子来说
dxdy=1+x+y+11
此时设 u=x+y+1,则 dy/dx=d(u−x−1)/dx=du/dx−1
于是转化为可分离变量的一阶线性微分方程
dxdu=2+u1
a0(x)y+a1(x)y′+a2(x)y′′+...+an(x)y(n)=b(x)
【同济】形如以下形式的方程称为一阶线性微分方程
dxdy+P(x)y=Q(x)
ey+y′=0
柯西-尤拉方程是形式如
x2y′′+bxy′+cy=0
(其中 b,c 是常數)的二階變係數常微分方程。
dxdy=rxr−1dx2d2y=r(r−1)xr−2
于是方程转化为
(r(r−1)+br+c)y=0
wiki 对于这个二次方程的三个解做了分别讨论
- 两个不同解时通解为 y=c1xm1+c2xm2
- 两个相同解时通解为 y=c1xmln(x)+c2xm
- 比较复杂,懒得写了,总之是复数
为什么这里好像假设了 x>0
哦,t 可以是复数,那没事了
微分算子复习
D→dtd
Dy=xy′,D(D−1)y=x2y′′
原方程转化为 (D2−D+bD+c)y=0
于是变成 y 对 t 的二阶齐次线性微分方程
dxdy+P(x)y=Q(x)yn
积 n 次即可
设 y′=p,即有 dxdp=f(x,p)
TODO: y′=exy 似乎也有这种形式,常数变易法能做吗
设 y′=p,则有
dxdp=f(y,p)=dydpdxdy=dydpp
转换为 dydp=p1f(y,p)
a1f′′(x)+a2f′(x)+a3f(x)=ecx(b1x3+b2x2+b3x+b4)
写为矩阵乘法:
(a1a2a3)f′′(x)f′(x)f(x)=ecx(b1b2b3b4)x3x2x1
设特解为
g(x)=ex(p1x3+p2x2+p3x+p4)=ecx(p1p2p3p4)x3x2x1
则 dg(x)/dx 可求为
g′(x)=cecx(p1p2p3p4)x3x2x1+ecx(p1p2p3p4)0000300002000010x3x2xx=ecx(p1p2p3p4)[c+0000300002000010]x3x2x1
递归可得
g′′(x)=ecx(p1p2p3p4){c[c+0000300002000010]+0000300002000010}x3x2x1
我草,写起来好麻烦,记
A=0000300002000010
于是上面都简化为:
g(x)=ecx(p1p2p3p4)x3x2x1g(x)=ecx(p1p2p3p4)(c+A)x3x2x1g′′(x)=ecx(p1p2p3p4)(c(c+A)+A)x3x2x1
于是
(a1a2a3)g′′(x)g′(x)g(x)=ecx(p1p2p3p4)(a1a2a3)c(c+A)+Ac+A1x3x2x1
故
(p1p2p3p4)(a1a2a3)c(c+A)+Ac+A1=(b1b2b3b4)