Typesetting math: 0%

微积分 A(1) 要点整理

期中考试前太鸽了就不补了,这里主要是期中考试之后的部分。

不定积分

不定积分的定义

不定积分的本质:找原函数。

对于两个函数 F,fF,f 而言,称函数 FfIR 上的原函数,当且 xI 都有 F(x)=f(x)

f(x)dxf 所有原函数的集合,称作 f 的不定积分。

可以证明,如果 F(x)f(x) 的原函数,那么 f(x)=F(x)+C,即两者只差一个常数。因为如果 G(x),F(x) 均为 f(x) 的原函数,那么 G(x)=F(x)(GF)(x)=0,因此 G(x)F(x) 为常值函数。也就是说,我们只用找到 f(x) 的一个原函数 F(x),那么 F(x)+C 就包含了 f(x) 的所有原函数。

容易证明积分满足线性性,即 f,g:IR,μ,λR,都有

(μf(x)+λg(x))dx=λf(x)dx+μg(x)dx

常用函数的不定积分

这里列举几个容易忘记的,太基础的就不在这里赘述了:

11+x2dx=arctanx+C

11x2dx=arcsinx+C

1sin2xdx=cotx+C

1cos2xdx=tanx+C

sinxcos2xdx=secx+C

cosxsin2xdx=cscx+C

不定积分的一些手段

1. 凑微分

如果被积函数可以写成 g(h(x))·h(x) 的形式,那么可以提一个 h(x) 变成 d(h(x)),这样换元令 y=h(x) 就变成了 g(y)dy。换句话说,g(h(x))h(x)dx=g(h(x))d(h(x)),其依据是复合函数求导的链式法则 (F(G(x)))=F(G(x))G(x)

2. 有理分式的积分

对于有理分式 P(x)Q(x),不妨假设 Q(x) 最高次项系数为 1,那么可以将 Q(x) 分解为 (x+a1)(x+a2)(x+an)(x2+b1x+c1)(x2+b2x+c2)(x2+bmx+cx) 的形式,这样 P(x)Q(x) 可以表示为

u1x+a1+u2x+a2++unx+an+v1x+w1x2+b1x+c1++vmx+wmx2+bmx+cm

使用待定系数法解出 u,v,w,此后根据

1x+adx=ln(|x+a|)+C

1(ax+b)2+1dx=1aarctan(ax+b)+C

将所有一次多项式变成带 ln 的项,二次多项式变成带 arctan 的项即可。

如果有重根(例如质因数分解形式中存在 (x+1)3 这种),那么假设 (x+ai) 这一项的次数为 k,那么待定系数法中应有 u1(x+a1),u2(x+a1)2,,uk(x+a1)k,对于二次式也同理。

3. 含三角函数的积分

  1. 凑微分:将 cosx 变成 dsinxsinx 变成 dcosx

  2. 二倍角公式 cos2x=1+cos2x2,sin2x=1cos2x2 来降次,但是不要把单个 cosx 变成 1sin2x,否则会引入更麻烦的根式。

  3. 一种能将含三角函数的积分全部转化为有理分式的方法:对于 cosxsinx 的多项式的积分,可以考虑万能公式:令 t=tan(x2),那么 sinx=2t1+t2,cosx=1t21+t2,这样

f(sinx,cosx)dx=f(2t1+t2,1t21+t2)d(2arctant)

还有一种比较特殊的处理方法:

tan4xdx=tan2x(sec2x1)dx=tan2xdx+tan2xdtanx=tan3x3(sec2x1)dx=tan3x3tanx+x

其原理是 tan2x=sec2x1,而 sec2x 刚好是 tanx 的导数,方便凑微分凑成 dtanx

4. 含根式的积分

  • 对于根号底下是关于 x 的一次多项式(如 x1)的根式,直接换元令 y=x1,这样 x=y2+1

  • 对于根号底下是关于 x 的二次多项式(如 1x2)的根式,采用三角换元:

    • 1x2 型,令 x=cosθ,这样 1x2=sinθ
    • 1+x2 型,令 x=tanθ,这样 1+x2=secθ
    • x21 型,令 x=secθ,这样 x21=tanθ

5. 分部积分

f(x)g(x)dx=f(x)g(x)f(x)g(x)dx

其依据是乘积的求导法则 (f(x)g(x))=f(x)g(x)+f(x)g(x)

6. 一类特殊的求解形如 eλxxksinx 的积分的方法

考虑不断将 eλxxksinx 求导会得到什么,显然你不论求多少次导,最终得到的式子都会形如 eλxsinxP(x)+eλxsinxQ(x),其中 P(x),Q(x) 都是最高项次数 xk 的多项式,也就是说我们可以将所有有用的式子写成 a0eλxsinx+a1eλxxsinx+a2x2eλxsinx++akeλxxksinx+ak+1eλxcosx+akeλxxksinx+ak+2eλxxcosx++a2k+1eλxxkcosx 的形式,那么一次求导相当于对这个系数序列 a 进行一次线性变换,即乘上一个矩阵,积分则是乘上该矩阵的逆矩阵,使用线性代数的知识求解即可。

定积分

定积分的定义

定积分的本质:求有向面积。

对于一闭区间 [a,b],定义其一组划分 P 为一组 a=x0<x1<x2<<xk=b,并定义这组划分的长度 P=max1ik{xixi1}

  • Darboux 积分的定义:对于一定义在 [a,b] 上的函数 f(x),考虑 [a,b] 的一组划分 P,定义这组划分的 Darboux 上和 ¯S(f,P)=ki=1supx[xi1,xi]f(x),下和 S_(f,P)=ki=1infx[xi1,xi]f(x),如果当划分长度趋近于 0 时,划分的 Darboux 上和与下和均存在极限且相等,那么定义 f(x)[a,b] 上的 Darboux 积分为这个相同的极限。
  • Riemann 积分的定义:对于一组 [a,b] 上的划分 P 和一组标志点集合 {ξi} 满足 xi1ξixi,定义其 Riemann 和 S(f,P,ξ)=ki=1f(ξi)(xixi1),如果存在 AR 满足 ϵ>0,存在 δ 使得对于所有 P<δ 的划分和任意标志点集合 (P,ξ),都有 |S(f,P,ξ)A|<ϵ,那么就称 Af(x)[a,b] 上的 Riemann 积分。

(Lebesgue 定理)以下结论等价:

  • f(x)[a,b] 上 Riemann 可积。
  • f(x)[a,b] 上有界,并且 Darboux 可积。
  • f(x)[a,b] 上有界,并且所有间断点组成了一个零测度集。(零测度集的定义,ϵ>0,存在一个长度之和不超过 ϵ 的区间集覆盖了这个区间上所有间断点)

根据 Lebesgue 定理,因为黎曼函数 R(x)={1q(x=pqQ,p,qZ)0(xQ) 在所有无理数处连续,有理数处间断,所以 R(x) 在任意有界区间上黎曼可积,而迪利克雷函数 D(x)={1(xQ)0(xQ) 则不是黎曼可积的。

定积分的性质

定积分满足线性性、保序性,证明是容易的。

Newton-Leibniz 公式

  1. 如果 f(x)[a,b] 上黎曼可积,那么其变上限积分函数 g(x)=xaf(x)dxf(x)[a,b] 上的原函数。
  2. 如果 f(x)[a,b] 上黎曼可积,且存在原函数 F(x),那么 baf(x)dx=F(b)F(a)

这意味这如果我们能求出被积函数的原函数,那么直接在原函数两端点处做差即可得到定积分的值

积分平均值定理

对于函数 f(x),g(x),如果 x[a,b],g(x)0,且 bag(x)dx>0f(x)[a,b] 上连续,那么 ξ(a,b),使得

f(ξ)bag(x)dx=baf(x)g(x)dx

其正确性可以由连续函数的介值性证得。

定积分的应用

注:在本节中,我们默认我们研究的曲线都是简单闭曲线,即曲线首尾相连、不自交且是正则的。(正则的定义:对于平面直角坐标系下由参数方程确定的曲线 x(x(t),y(t)):[α,β]R2 而言,该参数曲线是正则的当且仅当 t[α,β]x(t)0,其中 x(t)=(x(t),y(t)),在这个定义之下,参数曲线 {x=t3y=t2 就不是正则的,因为在 (0,0)x(t)=0)。

1. 计算封闭区域的面积

  1. 计算由 y1(x),y2(x)(x[α,β],y2(x)y1(x)) 为上、下边界围成的区域面积:

    βα(y2(x)y1(x))dx

  2. 平面直角坐标系下由参数方程确定的简单闭曲线:对于沿着自然正向(沿着 t 增加的方向走时左手边位于曲线封闭区域内)的曲线,有

S=ydx=βαy(t)x(t)dt

S=xdy=βαx(t)y(t)dt

S=12βαdet[x(t)y(t)x(t)y(t)]

​ 若是自然反向则需要加负号。

  1. 极坐标方程系下由参数方程确定的曲线:

    S=12θ2θ1r2dθ=12βαr(t)2θ(t)dt

    计算依据:S=12absinθ,且 θ 足够小时 θsinθ

2. 计算弧长

对于一条曲线而言,考虑在这条曲线上钉一些钉子,并顺次将钉的钉子连接起来形成一条折线,如果当钉子数量趋于无穷时,折线长的极限存在,那么就称这条曲线为可求长曲线,该极限称为曲线的弧长。计算弧长的方法:

  1. n 维平面直角坐标系下由参数方程确定的曲线:dl=dx21+dx22++dx2n=x1(t)2+x2(t)2++xn(t)2dt。据此可以积分得到

    L=βαx(t)dt

  2. 在极坐标系下由参数方程确定的曲线:dl=dr2+(rdθ)2=r(t)2+r(t)2θ(t)dt

3. 计算曲率

对于曲线上一点 M(x(t),y(t)) 而言,考虑曲线上另一点 M1,我们在 M,M1 处分别做曲线的切线,设两条切线形成的夹角为 θ(锐角),M,M1 在曲线上的距离为 l,那么我们定义曲线在 M 点处的曲率 κ=limM1M|θ|l。根据这个定义我们知道曲线在一个点处的曲率同时也是曲线在这个点处曲率圆半径的倒数。

对于曲线上一点 M(x(t),y(t)),计算曲线在此点处曲率的公式为 κ=x(t)×x(t)x(t)3。如果我们考虑记 l(t)=tαx(s)ds 的反函数为 t(l),并定义 ˜x(l)=x(t(l)),那么曲线在 M(x(t),y(t)) 处的曲率可以用更简洁的式子表示:κ=˜x(l(t))

如果函数以极坐标参数方程的方式求出,那么类似地可以变形得到 κ=|r(θ)2+2r(θ)2r(θ)r(θ)|(r(θ)2+r(θ)2)32

4. 计算旋转体体积

对于一曲线 y=f(x)(x[a,b]),如果我们要求其沿着 x 轴旋转一周得到的旋转体的体积,考虑微元法,将其沿着 x 轴切成若干薄片,每一片看作一圆柱体,那么圆柱体的体积就是 πy2dx,这样可以得到旋转体体积

V=πbaf2(x)dx

如果曲线以参数方程 x=x(t),y=y(t) 的形式给出,一样可以求得

V=πβαy2(t)x(t)dt

5. 计算旋转体表面积

对于一曲线 y=f(x)(x[a,b]),如果我们要求其沿着 x 轴旋转一周得到的旋转体的表面积,考虑微元法,将整段弧长切成若干段,如果将每一段视作与 x 轴平行的线段,那么其绕一周得到的旋转体表面积就是 2πydl,这样可以得到旋转体表面积

S=2πba|f(x)|dl

如果曲线以参数方程 x=x(t),y=y(t) 的形式给出,则有

S=2πβα|y(t)|x(t)2+y(t)2dt

广义积分

带无穷的积分

对于函数 f(x):[a,)R,如果 baf(x)[a,b] 上黎曼可积,并且

limb+baf(x)dx

存在,那么称

+af(x)dx=limb+baf(x)dx

同理定义左端点是负无穷的积分。

对于左右端点都是无穷的积分,如果

limb+baf(x)dx

limbabf(x)dx

都存在,那么称

+f(x)dx=limb+baf(x)dx+limbabf(x)dx

注意:

  1. lima+aaf(x)dx 来定义左右端点都是无穷的积分则是不正确的,比方说要判断 +x1+x2dx 的收敛性,因为 x1+x2是奇函数,所以用 lima+aaf(x)dx 的定义得到的极限是 0,但按照正确的定义这个积分是不收敛的。
  2. 如果 +af(x)dx 收敛,并不能推出 limx+f(x)=0,反例:f(x)={12n|xn|(x[n2n,n+2n])0(x[n2n,n+2n])

带瑕点的积分

对于函数 f(x):[a,b)R 无界,如果其在任意 [a,bϵ](ϵ>0) 上黎曼可积,并且

limϵ0+bϵaf(x)dx

存在,那么称

baf(x)dx=limϵ0+bϵaf(x)dx

并称 b 为该积分的瑕点。同理可以定义左端点是瑕点的积分。

以上两类积分统称为广义积分。

判断广义积分是否收敛的四种方法

1. Newton-Leibniz 公式

对于广义积分而言,需将常义积分的 Newton-Leibniz 公式进行一些推广。

f(x)[a,ω)(ω[a,){+}) 上任意子区间上可积,并且 F(x) 为其原函数,那么 ωaf(x)dx 存在当且仅当 limxωF(x) 存在,此时

ωaf(x)dx=limxωF(x)F(a)

Newton-Leibniz 公式适用于那些可以较容易用不定积分的手段求出原函数的函数的积分。

2. 柯西收敛准则

+af(x)dx 收敛,当且仅当 ϵ>0Nϵ 满足 A2>A1>Nϵ,都有 |A2A1f(x)dx|<ϵ

证明可以用类似于有理数的柯西列的方法,这里不再赘述。

3. 比较法

称广义积分 +af(x)dx 绝对收敛,当且仅当 +a|f(x)|dx 也收敛,反之其条件收敛。

容易证明绝对收敛的积分必然积分,其依据是三角不等式

|baf(x)dx|ba|f(x)|dx

对于广义积分 +af(x)dx 而言,如果存在函数 g 满足:

  • +ag(x)dx 绝对收敛。
  • x+ 时,f(x)=O(g(x))

那么 +af(x)dx 绝对收敛。

常用于做等价代换的 g(x) 有:

  • +0eλxdx。该广义积分收敛当且仅当 λ>0,收敛时极限为 1λ
  • +1xpdx。该广义积分收敛当且仅当 p<1,收敛时极限为 1p+1
  • 10xpdx。该广义积分收敛当且仅当 p>1,收敛时极限为 1p+1

此种方法适用于不存在原函数,但容易用泰勒展开求出其在无穷或者瑕点附近的阶的广义积分。

4. Dirichlet-Abel 判别法

  • Dirichlet 定理:对于在任意 [a,A] 上黎曼可积的函数 f(x) 和单调函数 g(x),如果 f 的变上限积分函数 baf(x)dx 有界,并且 limx+g(x)=0,那么 +af(x)g(x)dx 的极限存在。
  • Abel 定理:对于在任意 [a,A] 上黎曼可积的函数 f(x) 和单调函数 g(x),如果 +af(x)dx 存在,并且 limx+g(x)=A(AR),那么 +af(x)g(x)dx 的极限存在。

证明是使用分部积分。

比较重要的广义积分

Γ 函数

考虑积分

+0xα1exdx

首先,

+0xα1ex=10xα1ex+10xα1ex

对于前者,当 x0 时候,ex 趋近于 1,因此 xα1ex=O(xα1),根据比较法可知积分收敛当且仅当 α1>1,即 α>0

对于后者,因为指数函数的阶高于一切幂函数,所以存在 N 使得 x>N,xα1ex<ex/2,根据比较法可知不论 α 取何值加号后面的部分都是收敛的。

换句话说上述积分收敛当且仅当 α>0

这样我们定义

Γ(α)=+0xα1exdx

那么分部积分可得

Γ(α)=+0xα1exdx=(xα1ex|+0)+(α1)+0xα2exdx=(α1)Γ(α1)

换句话说,对于正整数 n,有 Γ(n+1)=n!

B 函数

考虑积分

10xα1(1x)β1dx

首先,

10xα1(1x)β1dx=120xα1(1x)β1dx+112xα1(1x)β1dx

对于前者,当 x0+ 时,xα1(1x)β1=O(xα1),其收敛当且仅当 α>0

对于后者,当 x1 时,xα1(1x)β1=O((1x)β1),其收敛当且仅当 β>0

换句话说上述积分收敛当且仅当 α,β>0

这样我们定义

B(α,β)=10xα1(1x)β1dx

那么分部积分可得

B(α,β)=10xα1(1x)β1dx=1βxα2(1x)β1|10+αβ10xα2(1x)βdx=αβB(α1,β+1)

又因为当 α=0 时,B(α,β)=10(1x)β1dx=(1x)ββ|10=1β

所以 B(α,β)=Γ(α)Γ(β)Γ(α+β)

B 函数和 Γ 函数有不少类似之处,因此前者被称为第一类欧拉积分,后者被称为第二类欧拉积分。

计算积分 π20ln(sinx)dx

比较经典也比较有意思的一类定积分。

首先收敛性的证明:

π20ln(sinx)dx=π20ln(x)dx+π20ln(sinxx)dx

前者等于 (xlnxx)|π20,因为 limx0xlnx=0,所以积分收敛。后者因为 limx0sinxx=1 也收敛。所以两者之和 π20ln(sinx)dx 收敛。

π20ln(sinx)dx=12π0ln(sinx)dx=π20ln(sin2x)dx=π20ln(2sinxcosx)dx=πln22+π20ln(sinx)dx+π20ln(cosx)dx=πln22+2π20ln(sinx)dx=πln22

有点类似于分部积分时将式子变成 f(x)dx=Af(x)dx+g(x)(A1) 的形式的套路,但是不完全一样。

+0sinxxαdx 的收敛性

x+ 时,sinxxα=O(xα),所以当 α>1 时,+1sinxxαdx 绝对收敛,但因为 sinx 可正可负,当 α1 的时候的收敛性并不能用比较法直接判断。

x0+ 时,sinxxα=x1α(1+o(1)),因此当 α<210sinxxαdx 收敛,反之 10sinxxαdx 发散。

0<α1 时,AB,都有

|BAsinxxαdx|=|BA1xαdcosx|=|cosxxα|BAαBAcosxxα+1dx|=|cosAAαcosBBααBAcosxxα+1dx|cosAAα+cosBBααBA1xα+1dx 1Aα+1Bα+1xα|BA3Aα

根据柯西准则可知,当 0<α1 时积分收敛但不绝对收敛。

α0 时,

2kπ+5π62kπ+π6sinxxαdx2kπ+5π62kπ+π612dxπ3

根据柯西准则可知积分不收敛。

因此:

  • α0α2 时积分发散。
  • 1<α<2 时积分绝对收敛。
  • 0<α1 时积分条件收敛。

微分方程

微分方程的定义

微分方程:含有未知函数及其导函数的方程或方程组。

常微分方程:微分方程中所有未知函数都是同一一元函数的若干阶导数的微分方程,即可以表示为 F(x,y,y,y,,y(n))=0 的微分方程。此时的 n 被称为微分方程的阶数。

满足上述微分方程的 y 称为微分方程的特解。含 n 个任意常数的解被称为微分方程的通解通解不一定包含特解。

一些经典一阶微分方程解法

1. y=f(x)

直接对 f(x) 做不定积分即可。

y(x)=f(x)dx+C

2. y=p(x)q(y)

将式子变形为

dyq(y)=p(x)dx

积分得到 p(x),1q(y) 的原函数 P(x),Q(y),这样微分方程解的隐式形式可以表示为

P(x)=Q(y)+C

当然,如果 Q(y) 存在反函数,那么可以进一步变形为 y=f(x) 的形式,这种情况一般具体题目具体分析。

3. y=q(yx)

u=yx,那么根据除法求导法则

dudx=yxxyx2

然后把 y=q(u) 带入可得

dudx=q(u)ux

右边有个 1x 的系数我们看着不爽,将其与 dx 结合起来变成 d(ln|x|) 可以化简为

dudln|x|=q(u)u

这样积分得到

ln|x|=1q(u)udu

不妨设 1q(u)u 存在原函数 P(u),那么上式等价于

x=CeP(yx)

因为 C 可正可负,这样我们左边的 x 就不用加绝对值了。这样我们得到了参数方程的隐式表达。

4. y=f(a1x+b1y+c1a2x+b2y+c2)

首先解方程组 {a1x+b1y+c1=0a2x+b2y+c2=0,会出现两种情况:

  • 方程组有解,假设解是 {x=x0y=y0,那么我们换元 p=xx0,q=yy0,这样式子变成 f(a1p+b1qa2p+b2q),这是关于 qp 的函数,使用 y=f(yx) 的解法求解即可。
  • 方程组无解,这时 a1a2=b1b2,微分方程可以写成 y=f(α+βa2x+b2y+c2) 的形式,此时再换元 u=a2x+b2y+c2u=b2y+a2,微分方程进一步变成 u=b2f(α+βu)+a2 的形式,令 p(x)=1,q(u)=b2f(α+βu)+a2,那么式子就是 u=p(x)q(u) 的形式,使用相应的方式求解即可。

5. y=f(ax+by+c)

换元令 u=ax+by+c,那么 dudx=a+by=a+bf(u),令 p(x)=1,q(u)=a+bf(u),那这个式子就是 u=p(x)q(u) 的形式,使用相应的方式求解即可。

6. 一阶线性方程

一阶线性方程是形如 y=a(x)y+b(x) 的微分方程,当 b(x)=0 时微分方程又被称为齐次方程

可以注意到:

  • 如果 f(x),g(x) 为齐次方程 y=a(x)y 的解,那么 f(x)+g(x) 也是齐次方程 y=a(x)y 的解。
  • 如果 f(x) 为齐次方程 y=a(x)y 的解,g(x) 为非齐次方程 y=a(x)y+b(x) 的解,那么 f(x)+g(x) 是非齐次方程 y=a(x)y+b(x) 的解。
  • 如果 f(x),g(x) 为非齐次方程 y=a(x)y+b(x) 的解,那么 f(x)g(x) 是齐次方程 y=a(x)y 的解。

换句话说,一阶线性方程的解形成了一组线性空间,这启发我们先求出对应的齐次方程的解 f(x),再求出非齐次方程的一组特解 g(x),这样非齐次方程的通解就可以表示为 Cf(x)+g(x) 的形式。

求解齐次方程是容易的,令 p(x)=a(x),q(y)=y,那么式子是 y=p(x)q(y) 的形式,使用分离变量的方法可以求得

f(x)=Cea(t)dt

求解非齐次方程的方法:常数变易法。设 y(x)=C(x)f(x),那么 C(x)f(x)+f(x)C(x)=a(x)C(x)f(x)+b(x),而因为 f(x) 是齐次方程 y=a(x)y 的解,因此 f(x)=a(x)f(x),因此化简可得 C(x)=b(x)f(x),积分可求得 C(x),进而求得非齐次方程的一组特解。

C(x) 积分得到 C(x) 的那一步中,C(x) 常常能够写成 eλxxk 或者 eλxxksinx 的形式,此时除了分部积分以外,还有一种使用线性代数知识求解的方法,可以参考「不定积分」那一部分的方法 6。

推广:伯努利方程:y=a(x)y+b(x)yα。令 z=y1α,那么 z=(1α)yαy=(1α)(b(x)+a(x)z),这是一个线性方程,可以解出 z 的通解,进而得到 y 的通解。

一些经典高阶微分方程的解法

1. 不显含 y 的微分方程

对于形如

F(x,y(k),y(k+1),,y(n))=0

的微分方程,考虑令 p=y(k),那么微分方程可以化作

F(x,p,p,,p(nk))=0

解出 p=u 以后对 u 积分 k 次即可得到 y

2. 不显含 x 的微分方程

对于形如

F(y,y,,y(n))=0

的微分方程,考虑 u(y)=y,那么微分方程可化为

F(y,u,u·dudy,(uddy)n1u)=0

特别地,当 n=2 时,微分方程可化为

F(y,u,uu)=0

此时 y 代替了 x 的地位,u 代替了 y 的地位,解出 u 以后再解一遍一阶微分方程 y=u(y) 即可。

3. 常系数线性齐次方程

形如 any(n)+an1y(n1)++a1y+a0y=0 的微分方程被称为常系数线性齐次方程。

对于此类方程,求解的通用方法是,先解出方程的 n特征根,即考虑方程 anxn+an1xn1++a1x+a0=0。有以下四种情形:

  • 特征根方程的解存在单个实根 λ,则微分方程的解里存在带自由元的项 Ceλx
  • 特征根方程的解里存在 k 重实根 λ,则微分方程的解里存在带自由元的项 eλx(Ck1xk1+Ck2xk2++C0)
  • 特征根方程的解里存在单个复根 a+bi,则微分方程的解里存在带自由元的项 eax(C1sinbx+C2cosbx)
  • 特征根方程的解里存在 k 重复根 a+bi,则微分方程的解里存在带自由元的项 eax((Ck1xk1+Ck2xk2+C0)sinbx+(Dk1xk1+Dk2xk2++D0)cosbx)

把这些带自由元的项相加即可得到常系数线性齐次方程的通解。

4. 常系数非齐次线性方程

形如 any(n)+an1y(n1)++a1y+a0y=f(x) 的微分方程被称为常系数线性非齐次方程。

因为常系数线性非齐次方程的解是一组线性空间,因此我们只需按照前文所述的方法求出对应的常系数线性齐次方程的通解以后,再求出一组常系数线性非齐次方程的特解即可。

而对于此种形式的微分方程而言,f(x) 通常是 eλx·P(x) 的形式,其中 P(x) 是关于 x 的多项式,因为形如 eλx·P(x) 的式子不管求多少遍导以后仍是这种形式,因此不妨假设非齐次方程的特解 Y(x)=eλxQ(x),然后通过待定系数法可以求出 Q(x) 的系数。

如果 f(x) 中带三角函数,比如 eλxsin(ωx),那么因为 eix=cosx+isinx,我们可以将 sin(ωx) 看成指数是复数的指数函数,使用类似的方法同样可以求解。

5. 欧拉方程

形如 pny(n)xn+pn1y(n1)xn1++p1yx+p0y=f(x) 的方程被称为欧拉方程。

对于欧拉方程而言,一般性的求解方法是,令 t=lnx,那么 y(n)xn=(ddtn+1)(ddtn+2)(ddt1)ddty,使用这种方法,我们设 y=f(t),那么式子可以转化为关于 y 的常系数高阶线性方程,求出 f(t) 以后根据 t=lnx 带回去即可。

posted @   tzc_wk  阅读(371)  评论(6编辑  收藏  举报
相关博文:
阅读排行:
· 在鹅厂做java开发是什么体验
· 百万级群聊的设计实践
· WPF到Web的无缝过渡:英雄联盟客户端的OpenSilver迁移实战
· 永远不要相信用户的输入:从 SQL 注入攻防看输入验证的重要性
· 全网最简单!3分钟用满血DeepSeek R1开发一款AI智能客服,零代码轻松接入微信、公众号、小程
历史上的今天:
2022-12-12 析合树建树
点击右上角即可分享
微信分享提示