[学习笔记] - 复变函数与积分变换

复数

z=x+iyz = x + iy

称为复数,其中ii为虚数单位,并规定i2=i×i=1i^2 = i \times i = -1,或i=1i = \sqrt{-1}x,yRx, y \in \mathbb{R}

  • 实部:Rez=x(Real Part)\text{Re}z = x(\text{Real Part})
  • 虚部:Lmz=x(Imaginary Part)\text{Lm}z = x(\text{Imaginary Part})

此笔记仅作为参考手册,不进行证明等操作

性质

  • z1=z2x1=x2,y1=y2z_1 = z_2 \Leftrightarrow x_1 = x_2, y_1 = y_2
  • z=0x=y=0z = 0 \Leftrightarrow x = y = 0

复数不能比较大小

复数四则运算

加法

z1+z2=(x1±x2)+i(y1±y2)z_1 + z_2 = (x_1 \pm x_2) + i(y_1 \pm y_2)

乘法

z1z2=(x1x2y1y2)+i(x1y2+x2y1)z_1z_2 = (x_1x_2 - y_1y_2) + i(x_1y_2 + x_2y_1)

除法

z1z2=x1x2+y1y2x22+y22+ix2y1x1y2x22+y22\frac{z_1}{z_2} = \frac{x_1x_2 + y_1y_2}{x_2^2 + y_2^2} + i\frac{x_2y_1 - x_1y_2}{x_2^2 + y_2^2}

除法的本质是将分母实体化

复数几何表示

复数类似向量可以在二维平面上表示

共轭复数z=x+iyz = x + iy的共轭复数为xiyx - iy,它们以x=0x = 0对称
复数的长度: z=x2+y2|z| = x^2 + y^2

  • z=z|z| = |\overline{z}|
  • zz=z2z\overline{z} = |z|^2
  • z=z\overline{\overline{z}} = z
  • 如果zz是实数:z=zz = \overline{z}

Rez=z+z2Imz=zz2i\text{Re}z = \frac{z + \overline{z}}{2} \quad \text{Im}z = \frac{z - \overline{z}}{2i}

z1±z2=z1±z2z1z2=z1z2z2z1=z2z1\overline{z_1 \pm z_2} = \overline{z_1} \pm \overline{z_2} \quad \overline{z_1\cdot z_2} = \overline{z_1} \cdot \overline{z_2} \quad \overline{\frac{z_2}{z_1}} = \frac{\overline{z_2}}{\overline{z_1}}

复数三角表示

z=x+iy{x=rcosθy=rsinθ=r(cosθ+isinθ)\begin{aligned} z &= x + iy \rightarrow \begin{cases} x = r\cos\theta \\ y = r\sin\theta \end{cases} \\ &= r(\cos\theta + i\sin\theta) \end{aligned}

  • z1z2=r1r2[cos(θ1+θ2)+isin(θ1+θ2)]z_1 \cdot z_2 = r_1r_2[\cos(\theta_1 + \theta_2) + i\sin(\theta_1 + \theta_2)]
  • z1z2=r1r2[cos(θ1θ2)+isin(θ1θ2)]\frac{z_1}{z_2} = \frac{r_1}{r_2}[\cos(\theta_1 - \theta_2) + i\sin(\theta_1 - \theta_2)]

复数指数形式

欧拉公式

eix=cosx+isinxe^{ix} = \cos x + i \sin x

从而可以得到

z=r(cosθ+isinθ)=reiθ\begin{aligned} z &= r(\cos\theta + i\sin\theta) \\ &= r \cdot e^{i\theta} \end{aligned}

θ\theta的值为zz的辐角,记为argz\arg zargz=θ+2nπ\arg z = \theta + 2n\pi

  • z1z2=r1r2ei(θ1+θ2)z_1 \cdot z_2 = r_1r_2e^{i(\theta_1 + \theta_2)}
  • z1z2=r1r2ei(θ1θ2)\frac{z_1}{z_2} = \frac{r_1}{r_2}e^{i(\theta_1 - \theta_2)}

乘幂公式

zn=rn(cosnθ+isinnθ)=rneinθzn=1zn=1rn(cosnθisinnθ)nz=r1n(cosθn+isinθn)=nr(cosθ+2kπn+isinθ+2kπn)\begin{aligned} z^n &= r^n(\cos n\theta + i\sin n\theta) \\ &= r^ne^{in\theta} \\ z^{-n} &= \frac{1}{z^n} = \frac{1}{r^n} (\cos n\theta - i\sin n\theta) \\ ^n\sqrt{z} &= r^{\frac{1}{n}}(\cos \frac{\theta}{n} + i\sin\frac{\theta}{n}) \\ &= ^n\sqrt{r} (\cos \frac{\theta + 2k\pi}{n} + i\sin\frac{\theta + 2k\pi}{n}) \end{aligned}

n个根的几何意义为以原点为中心,nn^n\sqrt{n}为半径的内接正nn边形的顶点

无穷远点

全部复数中,没有一个数和00相乘其积异于0,所以不存在一个属,是被0所除的商,引入\infty,称为无穷大

=10\infty = \frac{1}{0}

  • a+=+a=a + \infty = \infty + a = \infty
  • a=\infty - a = \infty
  • a=a - \infty = -\infty
  • a=a=a \cdot \infty = \infty \cdot a = \infty
  • a=0\frac{a}{\infty} = 0
  • a=\frac{\infty}{a} = \infty

在复平面上。 没有一个点和\infty对应,但是可以设想有一个理想点和它对应,称为无穷远点,平面上所有直线都通过无穷远点(黎曼球面

平面点集

领域

平面上以z0z_0为中心,δ\delta为半径的圆,表示为

zz0<δ|z - z_0| < \delta

称它为z0z_0的领域,记作U(z0,δ)U(z_0, \delta),称0<zz0<δ0 < |z - z_0| < \delta所确定的点集为z0z_0的去心领域
z0=z_0 = \inftykk为任意正数,满足z>k|z|>k的点称为无穷远点的领域

  • 内点:如果δ>0\exist \delta > 0,使得U(z0,δ)GU(z_0, \delta) \subset G,那么称z0z_0为集GG的内点
  • 开集:如果集GG中的每个点都是它的内点,那么称集GG为开集
  • 余集:平面上不属于GG的点的全体称为GG的余集,记为GcG^c
  • 内集:开集的余集称为闭集
  • 边界点z0z_0是一点,若在z0z_0的任意领域既有GG的点也有GcG^c的点,则为边界点
  • 边界GG的边界点全体称为GG的边界
  • 孤立点z0Gz_0 \in G,若在z0z_0的领域内除z0z_0外不含G的点,则称z0z_0是G的一个孤立点
    • 孤立点一定是一个边界点
  • 有界集:如果存在一个以z=0z = 0为中心的圆盘包含GG
  • 无界集:如果不存在一个以z=0z = 0为中心的圆盘包含GG
  • 连同集:设DD为一点集,对DD中任意两点可用一条全部属于D的折现连接起来,则D称为联通集
  • 区域
    • DD是开集
    • DD是连同的
    • 区域是开集
  • 闭区域:区域DD与它的边界一起构成的

复平面上的曲线

z=a(cost+isint)(0t2π)z = a(\cos t + i\sin t) \quad (0 \leq t \leq 2\pi)

圆周

z=a|z| = a

连通z1z_1z2z_2的直线段

z=z1+(z2z1)t(0t1)z = z_1 + (z_2 - z_1)t \quad (0 \leq t \leq 1)

  • 光滑曲线:如果在区间atba \leq t \leq bx(t)x'(t)y(t)y'(t)都是连续的,且对于tt的每一个值,有x(t)2+y(t)20x'(t)^2 + y'(t)^2 \neq 0
  • 重点:满足a<t1<b,a<t1<ba < t_1 < b,a < t_1 < b,当t1t2t_1 \neq t_2z(t1)=z(t2)z(t_1) = z(t_2)时,点z(t1)z(t_1)称为曲线CC的重点
  • 简单曲线:没有重点的连续曲线
  • 简单闭合曲线:闭合的简单曲线

复变函数

设有一复数z=x+iyz = x + iy的集合GG,如果有一个确定的法则存在,对于集合GG中的每一个复数zz,按照这一法则,复数w=u+ivw = u + iv就随之而定,那么称ww是复变数zz的函数,简称复变函数,记为w=f(z)w = f(z)

w=u+iv=f(z)(z=x+iy){u=u(x,y)v=v(x,y)\begin{gathered} w = u + iv = f(z) \quad (z = x + iy) \\ \Downarrow \\ \begin{cases} u = u(x, y) \\ v = v(x, y) \end{cases} \end{gathered}

把函数w=f(z)w = f(z)的定义集合GG看成是zz平面上的一个点集合,而把函数值集合DD看成是ww平面上的一个点集合,那么函数w=f(z)w = f(z)在集合上就可用看成是把集合GG变到DD的一个映射

f(z)f(z)将G中的点zz映射到DD中的ww,即集GG映射为集DD,则称点ww为点zz的像,点zz为点ww的原像,同样称DDGG的像,GGDD的原像

复变函数就是把点集进行映射为好研究的形状进行研究

复变函数的极限

设函数w=f(z)w = f(z)定义在z0z_0的领域0<zz0<ρ0 < |z - z_0| < \rho内,总有一确定数AA存在,ε>0\forall \varepsilon > 0,相应地必有整数δ(ε)\delta(\varepsilon),使得0<zz0<ρ0 < |z - z_0| < \rho时,f(z)A<ε|f(z) - A| < \varepsilon,则称AAf(z)f(z)zz趋向z0z_0时的极限,记为

limzz0f(z)=A\lim_{z \rightarrow z_0} f(z) = A


f(z)=u(x,y)+iv(x,y)f(z) = u(x, y) + iv(x, y)A=a+bi,z0=x0+iy0A = a + bi, z_0 = x_0 + iy_0limzz0f(z)=A\lim_{z \rightarrow z_0}f(z) = A的充要条件为

limxx0,yy0u(x,y)=alimxx0,yy0v(x,y)=b\lim_{x \rightarrow x_0, y \rightarrow y_0} u(x, y)= a \quad \lim_{x \rightarrow x_0, y \rightarrow y_0} v(x, y) = b

可以把求复变函数的极限转化为求两个二元实变量函数的极限

  • limzz0f(z)=A\lim_{z \rightarrow z_0} f(z) = A
  • limzz0g(z)=B\lim_{z \rightarrow z_0} g(z) = B
    • limzz0[f(z)+g(z)]=A+B\lim_{z \rightarrow z_0}[f(z) + g(z)] = A + B
    • limzz0f(z)g(z)=AB\lim_{z \rightarrow z_0}f(z)g(z) = AB
    • limzz0f(z)g(z)=AB\lim_{z \rightarrow z_0}\frac{f(z)}{g(z)} = \frac{A}{B}

复变函数的连续性

limzz0f(z)=f(z0)\lim_{z \rightarrow z_0} f(z) = f(z_0),则称f(z)f(z)z0z_0处连续
f(z)f(z)在区域DD内处处连续,则称f(z)f(z)在D内连续

  • z0z_0这点要有定义

  • z0z_0这点的极限值要存在

  • z0z_0这点的极限值要等于在这一点的函数值

  • 定理:函数f(z)=u(x,y)+iv(x,y)f(z) = u(x, y) + iv(x, y)z0=x0+iy0z_0 = x_0 +iy_0处连续的充要条件为u(x,y)u(x, y)v(x,y)v(x, y)(x0,y0)(x_0, y_0)处连续

  • 定理:连续函数的和,差,积,商,仍然时连续函数,连续函数的复合函数还是连续的

复变函数的导数

w=f(z)w = f(z)定义于区域DDz0,z0+ΔZDz_0, z_0 + \Delta Z \in D,如果

limΔz0ΔWΔZ=limΔz0f(z0+Δz)f(z0)Δz\lim_{\Delta z \rightarrow 0}\frac{\Delta W}{\Delta Z} = \lim_{\Delta z \rightarrow 0}\frac{f(z_0 + \Delta z) - f(z_0)}{\Delta z}

存在,则称f(z)f(z)z0z_0处可导,极限值称为f(z)f(z)z0z_0的导数

在一点可导和在一点可微是等价的

  • 连续不一定可导
  • 可导一定连续
  • C=0C' = 0
  • (f(z)±g(z))=f(z)±g(z)(f(z) \pm g(z))' = f'(z) \pm g'(z)
  • (f(z)g(z))=f(z)g(z)+f(z)g(z)(f(z)g(z))' = f'(z)g(z) + f(z)g'(z)
  • f(z)g(z)=f(z)g(z)f(z)g(z)g2(z)\displaystyle \frac{f(z)}{g(z)}' = \displaystyle\frac{f'(z)g(z) - f(z)g'(z)}{g^2(z)}
  • f(g(z))=f(w)g(z)f(g(z))' = f'(w)g'(z)
  • f(z)=1φ(w)f'(z) = \frac{1}{\varphi'(w)},其中w=f(z)w = f(z)z=φ(w)z = \varphi(w)互为单只反函数

解析函数

  • f(z)f(z)z0z_0z0z_0的领域内处处可导,则称f(z)f(z)z0z_0解析
  • f(z)f(z)在区域DD内每一点解析,则称f(z)f(z)在区域DD内解析,或称f(z)f(z)是区域DD内的一个解析函数
  • f(z)f(z)z0z_0处不解析,则称z0z_0f(z)f(z)的奇点

解析 > 可导 > 连续 > 极限,解析要求这一点以及某个领域可导

  • 函数f(z)f(z)在点zz解析就是函数f(z)f(z)在以点zz为圆心的某个领域内可微
  • 两个解析函数的和,差,积,商都是解析函数,解析函数的复合函数仍是解析
  • 所有多项式在复平面上是处处解析的
  • 任意的一个有理函数p(z)q(z)\frac{p(z)}{q(z)}(均为多项式)在不含分母为零的点的区域是解析函数

解析函数充要条件

w=f(z)=u(x,y)=iv(x,y)w = f(z) = u(x, y) = iv(x, y),柯西-黎曼方程为

ux=vyuy=vx\frac{\partial u}{\partial x} = \frac{\partial v}{\partial y} \quad \frac{\partial u}{\partial y} = -\frac{\partial v}{\partial x}

函数w=f(z)=u(x,y)=iv(x,y)w = f(z) = u(x, y) = iv(x, y)z=x+iyz = x + iy处可导的充分必要条件是

  • u(x,y)u(x, y)v(x,y)v(x, y)在点(x,y)(x, y)处可微
  • 满足柯西黎曼方程

函数w=f(z)=u(x,y)=iv(x,y)w = f(z) = u(x, y) = iv(x, y)在区域DD内解析的充要条件是

  • u(x,y)u(x, y)v(x,y)v(x, y)DD内任一点z=x+iyz = x + iy处可微
  • 满足柯西黎曼方程

u(x,y)u(x, y)是一个二元函数,说它可微,指的是它的增量函数可以表示为两部分

  • 自变量的线性部分
  • 自变量的高阶无穷小部分

Δu=uxΔx+uyΔy+ρ(Δx2+Δy2)limΔx0,Δy0ρ=0\Delta u = \frac{\partial u}{\partial x}\Delta x + \frac{\partial u}{\partial y}\Delta y + \rho(\sqrt{\Delta x^2 + \Delta y^2}) \quad \lim_{\Delta x \rightarrow 0, \Delta y \rightarrow 0} \rho = 0

则称u(x,y)u(x, y)可微

调和函数

如果二元实函数φ(x,y)\varphi(x, y)在区域DD内有二阶连续偏导数,且满足拉普拉斯方程

2φx2+2φy2=0\frac{\partial^2\varphi}{\partial x^2} + \frac{\partial^2 \varphi}{\partial y^2} = 0

则称φ=φ(x,y)\varphi = \varphi(x, y)为区域DD内的调和函数,或称函数φ(x,y)\varphi(x, y)在区域DD内调和

  • 定理:如果f(z)=u+ivf(z) = u + iv在区域DD内解析,则uuvvDD内都是调和函数

φ(x,y)\varphi(x, y)ψ(x,y)\psi(x, y)均为区域DD内的调和函数,且满足柯西黎曼方程

φx=ψyψx=ψy\frac{\partial \varphi}{\partial x} = \frac{\partial\psi}{\partial y} \quad \frac{\partial\psi}{\partial x} = -\frac{\partial \psi}{\partial y}

则称ψ(x,y)\psi(x, y)φ(x,y)\varphi(x, y)的共轭调和函数

f(z)=x+iyf(z) = x + iy如果在整个复平面上是解析的,也知道它的虚部yy也是实部xx的共轭调和

  • 定理:复变函数f(z)=u(x,y)+iv(x,y)f(z) = u(x, y) + iv(x, y)在区域DD内解析的充分必要条件是:在区域DD内,f(z)f(z)的虚部v(x,y)v(x, y)是实部u(x,y)u(x, y)的共轭调和

指数函数

ez=expz=ex(cosy+isiny)e^z = \exp z = e^x(\cos y + i\sin y)

为复指数函数

此处的ee是复指数函数的专用记号

  • eze^z是单值函数
  • z,ez0\forall z,e^z \neq 0
  • z1=x1+iy1,z2=x2+iy2\forall z_1 = x_1 + iy_1, z_2 = x_2 + iy_2,有ez1+z2=ez1ez2e^{z_1 + z_2} = e^{z_1}e^{z_2}
  • ez+2kπi=eze^{z + 2k\pi i} = e^z
  • eze^z是处处解析的

对数函数

expw=ew=z\exp w = e^w = z,则称wwzz的对数函数,记为

w=f(z)=Lnzw = f(z) = \text{Ln}z

  • 多值性:设z=reiθz = re^{i\theta}w=u+ivw = u + iv,有ew=reiθe^w = re^{i\theta},即eu+iv=eueiv=rei(argz+2kπ)=rei(θ+2kπ)e^{u + iv} = e^ue^{iv} = re^{i(\arg z + 2k\pi)} = re^{i(\theta + 2k\pi)}
    • w=Lnz=lnr+i(θ+2kπ)w = \text{Ln}z = \ln r + i(\theta + 2k\pi)
  • Lnz1z2=Lnz1+Lnz2\text{Ln}z_1z_2 = \text{Ln}z_1 + \text{Ln}z_2
  • Lnz1z2=Lnz1Lnz2\text{Ln}\displaystyle\frac{z_1}{z_2} = \text{Ln}z_1 - \text{Ln}z_2

三角函数

cosz=eiz+eiz2sinz=eizeiz2i\cos z = \frac{e^{iz} + e^{-iz}}{2} \quad \sin z = \frac{e^{iz} - e^{-iz}}{2i}

  • 在复平面解析
    • sinz=cosz\sin z' = \cos z
    • cosz=sinz\cos z' = -\sin z
  • 奇偶性实函数
  • 周期性以2π2\pi为周期
  • 无界性
  • 恒等关系
    • sin2z+cos2z=1\sin^2z + \cos^2z = 1
    • sin(z1±z2)=sinz1cosz2±cosz1sinz2\sin(z_1 \pm z_2) = \sin z_1\cos z_2 \pm \cos z_1 \sin z_2
    • cos(z1±z2)=cosz1cosz2sinz1sinz2\cos(z_1 \pm z_2) = \cos z_1\cos z_2 \mp \sin z_1 \sin z_2

双曲函数

shz=ezez2chz=ez+ez2thz=shzchz=ezezez+ez\sh z = \frac{e^z - e^{-z}}{2} \quad \ch z = \frac{e^z + e^{-z}}{2} \quad \th z = \frac{\sh z}{\ch z} = \frac{e^z - e^{-z}}{e^z + e^{-z}}

  • shz=isiniz\sh z = -i \sin iz
    • shz=chzsh z' = \ch z
    • 周期:2πi2\pi i
  • chz=cosiz\ch z = \cos iz
    • chz=shzch z' = \sh z
    • 周期:2πi2\pi i
  • thz=itaniz\th z = -i \tan iz
    • 周期:πi\pi i

复积分

Cf(z)dz=limλ0k=1nf(ςk)Δzk\int_C f(z)\mathrm{d}z = \lim_{\lambda \rightarrow 0}\sum^n_{k=1}f(\varsigma_k)\Delta z_k

  • Cf(z)dz\int_{C^-}f(z)\mathrm{d}z:沿曲线CC的负方向积分
  • Γf(z)dz\oint_\Gamma f(z)\mathrm{d}z:沿闭曲线Γ\Gamma的积分

复积分的性质

  • C[af(z)+βg(z)]dz=aCf(z)dz+βCg(z)dz\int_C [af(z) + \beta g(z)]\mathrm{d}z = a \int_C f(z)\mathrm{d}z + \beta \int_C g(z)\mathrm{d}z
  • Cf(z)dz=Cf(z)dz\int_C f(z)\mathrm{d}z = -\int_{C^-}f(z)\mathrm{d}z
  • Cf(z)dz=c1f(z)dz+C2f(z)dz\int_C f(z)\mathrm{d}z = \int_{c_1}f(z)\mathrm{d}z + \int_{C_2}f(z)\mathrm{d}z
    • C=C1+C2C = C_1 + C_2
  • Cf(z)dzCf(z)dz=Cf(z)dsML|\int_C f(z)\mathrm{d}z| \leq \int_C |f(z)||\mathrm{d}z| = \int_C |f(z)|\mathrm{d}s \leq ML
    • M=maxzCf(z)M = \max_{z \in C}|f(z)|

例子:

I=Cdz(zz0)nCzz0=rI = \oint_C \frac{\mathrm{d}z}{(z - z_0)^n} \quad C \Rightarrow |z - z_0| = r

曲线CC的方程为z=z0+reiθz = z_0 + re^{i\theta}

I=02πreiθi(reiθ)ndθ=irn102πei(1n)θdθ\begin{aligned} I &= \int^{2\pi}_0 \frac{re^{i\theta}i}{(re^{i\theta})^n}\mathrm{d}\theta \\ &= \frac{i}{r^{n-1}}\int^{2\pi}_0 e^{i(1-n)\theta}\mathrm{d}\theta \end{aligned}

  • n=1n = 1I=2πiI = 2\pi i
  • n1n \neq 1I=0I = 0

柯西积分定理

设函数f(z)f(z)在单连通域DD内解析,Γ\GammaDD内的任意一条简单闭曲线,则有

Γf(z)dz=0\oint_\Gamma f(z)\mathrm{d}z = 0

设单连域DD的边界为CC,函数f(z)f(z)DD内解析,在D=D+C\overline{D} = D + C上连续,则有

Γf(z)dz=0\oint_\Gamma f(z)\mathrm{d}z = 0

设二连域DD的边界为C=C1+C2C = C_1 + C_2^-,函数f(z)f(z)DD内解析,在D\overline{D}上连续,则Cf(z)dz=0\oint_C f(z)\mathrm{d}z = 0

C1f(z)dz=C2f(z)dz\oint_{C_1} f(z)\mathrm{d}z = \oint_{C_2} f(z)\mathrm{d}z

设多连域的边界为C=C0+C1+C2++CnC = C_0 + C_1^- + C_2^- + \cdots + C_n^-,函数f(n)f(n)DD内解析,在D\overline{D}上连续,则

Cf(z)dz=0\oint_C f(z)\mathrm{d}z = 0

或者

C0f(z)dz=C1f(z)dz+C2f(z)dz++Cnf(z)dz\oint_{C_0} f(z)\mathrm{d}z = \oint_{C_1} f(z)\mathrm{d}z + \oint_{C_2} f(z)\mathrm{d}z + \cdots + \oint_{C_n} f(z)\mathrm{d}z

一个闭区域的复积分为里面的奇点函数值之和

路径无关性

设函数f(z)f(z)在单连通域DD内解析,C1,C2C_1, C_2DD内的任意两条从z0z_0z1z_1的简单曲线,则有

C1f(z)dz=C2f(z)dz\int_{C_1} f(z)\mathrm{d}z = \int_{C_2}f(z)\mathrm{d}z

z0z_0出发,沿C1C_1正向走到z0z_0

原函数

设在单连域DD内,函数F(z)F(z)恒满足条件F(z)=f(z)F'(z) = f(z),则F(z)F(z)称为f(z)f(z)DD内的一个原函数

f(z)dz=F(z)+c\int f(z)\mathrm{d}z = F(z) + c

Newton-Leibniz公式

f(z)f(z)在单连域DD内处处解析,G(z)G(z)f(z)f(z)的原函数,则

z0z1f(z)dz=G(z1)G(z0)\int_{z_0}^{z_1}f(z)\mathrm{d}z = G(z_1) - G(z_0)

柯西积分公式

如果函数f(z)f(z)在连续区域DD内解析,在DD的闭包上连续,z0Dz_0 \in D,则

f(z0)=12πiCf(z)zz0dzf(z_0) = \frac{1}{2\pi i} \oint_C \frac{f(z)}{z - z_0}\mathrm{d}z

解析函数在其解析区域的值完全由边界上的值确定,换句话说,解析函数可用其解析区域边界上的值以一种特定的积分形式表现出来

应用:反过来计算积分:Cf(z)zz0dz=2πif(z0)\oint_C \frac{f(z)}{z - z_0}\mathrm{d}z = 2\pi i f(z_0)

例子:

I=C2z1z2zdzI = \oint_C \frac{2z - 1}{z^2 - z}\mathrm{d}z

f(z)=2z1z2zdzf(z) = \frac{2z - 1}{z^2 - z}\mathrm{d}z,解得

I=C1f(z)dz+C2f(z)dz=C12z1z1zdz+C12z1zz1dz=2πi(2z1z1z=0+2z1zz=1)=4πi\begin{aligned} I &= \oint_{C_1}f(z)\mathrm{d}z + \oint_{C_2}f(z)\mathrm{d}z \\ &= \oint_{C_1}\frac{\frac{2z - 1}{z - 1}}{z}\mathrm{d}z + \oint_{C_1}\frac{\frac{2z - 1}{z}}{z - 1}\mathrm{d}z \\ &= 2\pi i (\frac{2z- 1}{z -1}|_{z = 0} + \frac{2z - 1}{z}|_{z = 1}) \\ &= 4\pi i \end{aligned}

平均值公式

如果函数f(z)f(z)zz0<R|z-z_0| < R内解析,在zz0R|z - z_0| \leq R上连续,则有

f(z0)=12π02πf(z0+Reiθ)dzf(z_0) = \frac{1}{2\pi} \int^{2\pi}_0 f(z_0 + Re^{i\theta})\mathrm{d}z

函数f(z)f(z)在圆心处的函数值 = 函数在圆周上的算数平均值

最大模定理

如果函数f(z)f(z)DD内解析,且不为常数,则在DDf(z)|f(z)|没有最大值

  • 在区域DD内解析的函数,如果其模在DD内达到最大值,则此函数比恒为常数
  • f(z)f(z)在有界区域DD内解析,在D\overline{D}上连续,则f(z)|f(z)|DD的边界上必能达到最大值

高阶导数

如果函数f(z)f(z)在区域DD内解析,在D=D+C\overline{D} = D + C上连续,则f(z)f(z)的各阶导数均在DD上解析,且

f(n)(z)=n!2πiCf(ς)(ςz)n+1dz(zD)f^{(n)}(z) = \frac{n!}{2\pi i}\oint_C \frac{f(\varsigma)}{(\varsigma - z)^{n+1}}\mathrm{d}z \quad (z \in D)

解析函数的导数仍解析

应用:反过来计算积分

Cf(z)(zz0)n+1dz=2πin!f(n)(z0)\oint_C \frac{f(z)}{(z - z_0)^{n+1}}\mathrm{d}z = \frac{2\pi i}{n!}f^{(n)}(z_0)

柯西不等式

设函数f(z)f(z)zz0<R|z - z_0| < R内解析,且f(z)<M|f(z)| < M,则

f(n)(z0)n!MRn|f^{(n)}(z_0)| \leq \frac{n!M}{R^n}

刘维尔定理

设函数f(z)f(z)在全平面上解析且有界,则f(z)f(z)为以常数

复数序列

znz_n为复数,称{Zn}n=1,2,3\lbrace Z_n \rbrace_n = 1,2,3为复数序列

  • limn+zn=a\lim_{n\rightarrow+\infty}z_n = a: 收敛
  • limn+zn=\lim_{n\rightarrow+\infty}z_n = \infty: 发散

zn=xn+iyn,a=a+iβz_n = x_n + iy_n, a = a + i\beta,则limn+zn=a\lim_{n\rightarrow+\infty}z_n = a的充要条件是

{limn+xn=alimn+yn=β\begin{cases} \lim_{n\rightarrow+\infty}x_n = a \\ \lim_{n\rightarrow+\infty}y_n = \beta \end{cases}

复数项级数

n=1+zn=z1+z2,+\sum_{n = 1}^{+\infty} z_n = z_1 + z_2,+ \cdots

复数项级数,简记为zn\sum z_n,称

k=1nzk=z1+z2++zk\sum_{k = 1}^{n} z_k = z_1 + z_2 + \cdots + z_k

为级数的部分和

  • limn+sn=s\lim_{n\rightarrow +\infty} s_n = s: 级数收敛

zn=xn+iyn,a=a+iβz_n = x_n + iy_n, a = a + i\beta,则limn+sn=s\lim_{n\rightarrow +\infty} s_n = s的充要条件是

{xnyn\begin{cases} \sum x_n \\ \sum y_n \end{cases}

都收敛

  • zn\sum |z_n|收敛,则称zn\sum z_n绝对收敛
  • zn\sum |z_n|发散,zn\sum z_n收敛,则称zn\sum z_n条件收敛
  • zn\sum |z_n|收敛,则zn\sum z_n必收敛

复变函数项级数

  • 复变函数序列:{fn(z)}n=1,2,\lbrace f_n(z) \rbrace_{n=1,2,\cdots}
  • 复变函数项级数:n=1+fn(z)=f1(z)+f2(z)+\sum_{n=1}^{+\infty} f_n(z) = f_1(z) + f_2(z) + \cdots
  • 级数fn(z)\sum f_n(z)的部分和:sn(z)=k=1nfk(z)s_n(z) = \sum_{k=1}^n f_k(z)
  • 收敛:limn+sn(z0)=s(z0)\lim_{n \rightarrow +\infty} s_n(z_0) = s(z_0)
  • 在区域DD内收敛,称s(z)s(z)为和函数,DD为收敛域

幂级数

n=0+an(za)n=a0+a1(za)+a2(za)2+\sum_{n=0}^{+\infty}a_n(z - a)^n = a_0 + a_1(z - a) + a_2(z - a)^2 + \cdots

其中an,aa_n, a为复常数,当a=0a = 0时有

n=0+anzn=a0+a1z+a2z2+\sum_{n=0}^{+\infty}a_nz^n = a_0 + a_1z+a_2z^2+\cdots

对于幂级数anzn\sum a_n z^n

  • 在0点必收敛
  • 如果级数在z0z_0点收敛,则它在z<z0|z| < |z_0|上绝对收敛
  • 如果级数在z1z_1点发散,则它在z>z1|z| > |z_1|上发散

  • n=0+anzn±n=0+bnzn=n=0+(an±bn)zn\sum_{n = 0}^{+\infty}a_nz^n \pm \sum_{n = 0}^{+\infty}b_nz^n = \sum_{n = 0}^{+\infty}(a_n \pm b_n)z^n
  • n=0+anznn=0+bnzn=n=0+(a0bn+a1bn1++anb0)zn\sum_{n = 0}^{+\infty}a_nz^n \cdot \sum_{n = 0}^{+\infty}b_nz^n = \sum_{n = 0}^{+\infty}(a_0b_n + a_1b_{n-1} + \cdot + a_nb_0)z^n

收敛圆

CRC_R的半径为RR

  • 称圆域z<R|z| < R为收敛圆
  • RR为收敛半径

级数在收敛圆的边界上各点的收敛情况时不一定的

收敛半径

  • 比值法:如果limn+an+1an=λ\lim_{n \rightarrow +\infty}\frac{|a_{n+1}|}{a_n} = \lambda
    • 收敛半径:R=1λR = \frac{1}{\lambda}
  • 根植法:如果limn+ncn=ρ\lim_{n\rightarrow +\infty}^n\sqrt{|c_n|} = \rho
    • 收敛半径:R=1ρR = \frac{1}{\rho}

幂级数性质

f(z)=n=0+an(zz0)n,zz0<Rf(z) = \sum^{+\infty}_{n=0}a_n(z - z_0)^n, |z - z_0 < R|,则

  • 函数f(z)f(z)在收敛圆zz0<R|z - z_0| < R内解析
  • 函数f(z)f(z)的导数可由其幂级数逐项求导得到
    • f(z)=n=1+nan(zz0)n1f'(z) = \sum^{+\infty}_{n=1} na_n(z - z_0)^{n-1}
  • 在收敛圆内可以逐项积分,即
    • F(z)=z0zf(z)dz=n=0+ann+1(zz0)n+1F(z) = \int^z_{z_0} f(z)\mathrm{d}z = \sum^{+\infty}_{n=0}\frac{a_n}{n + 1}(z - z_0)^{n+1}

设级数n=0+anzn\sum^{+\infty}_{n=0}a_nz^nz<R|z| < R内收敛,和函数为f(z)=n=0+anznf(z) = \sum^{+\infty}_{n=0}a_nz^n,又设函数g(n)g(n)z<r|z| < r内解析,且满足g(z)<R|g(z)| < R,则当z<r|z| < r时,有f(g(z))=n=0+an(g(z))nf(g(z)) = \sum^{+\infty}_{n=0}a_n(g(z))^n

泰勒级数

泰勒定理:设函数f(z)f(z)在区域DD内解析,CCDD的边界,z0Dz_0 \in DR=minzCzz0R = \min_{z \in C}|z - z_0|,则当zz0<R|z - z_0| < R,有

f(z)=n=0+an(zz0)nan=1n!f(n)(z0)=12πiLf(z)(zz0)n+1dzf(z) = \sum^{+\infty}_{n=0} a^n(z - z_0)^n \quad a_n = \frac{1}{n!}f^{(n)}(z_0) = \frac{1}{2\pi i}\oint_L \frac{f(z)}{(z - z_0)^{n+1}}\mathrm{d}z

  • 幂级数的收敛域必须是圆域
  • 幂级数一旦收敛,其和函数一定解析
  • 对于一个给定的函数,用任何方法展开为幂级数,其结果唯一
  • 函数f(z)f(z)z0z_0点展开为泰勒级数,其收敛半径等于从z0z_0点到f(z)f(z)的最近一个奇点zz的距离

在比较小的范围内,泰勒级数无限接近原函数,可以用它近似原函数


  • 直接计算

an=1n!f(n)(0)a_n = \frac{1}{n!}f^{(n)}(0)

常见展开式

11z=n=0+zn=1+z+z2+z3+z<1\frac{1}{1 - z} = \sum^{+\infty}_{n=0}z^n = 1 + z + z^2 + z^3 + \cdots \quad |z| < 1

ez=n=0+znn!=1+z+z22!+z33!+z<+e^z = \sum^{+\infty}_{n=0}\frac{z^n}{n!} = 1 + z + \frac{z^2}{2!} + \frac{z^3}{3!} + \cdots \quad |z| < +\infty

洛朗级数

引入负幂次项,就可能将一个函数在整个复平面上展开(除了奇点所在的圆周)

设函数f(z)f(z)在圆环域D:R1<zz0<R2D: R_1 < |z - z_0| < R_2内解析,则f(z)f(z)一定能再次圆环域中展开为

f(z)=n=+an(zz0)nan=12πiCf(ζ)(ζz0)n+1dζf(z) = \sum^{+\infty}_{n = -\infty}a_n(z - z_0)^n \quad a_n = \frac{1}{2\pi i}\oint_C \frac{f(\zeta)}{(\zeta - z_0)^{n+1}}\mathrm{d}\zeta

CC是圆环内绕z0z_0的任何一条简单闭曲线

  • 正幂次项:解析部分
  • 负幂次项:主要部分
  • 展开项唯一
  • 若在0zz0<R0 \leq |z - z_0| < R内解析,则f(z)f(z)在此圆环内的洛朗展开式就是泰勒展开式

洛朗级数的展开

在求展开式之前,都需要根据函数的奇点位置,将复平面分为若干个解析环
比如奇点为z1,z2,z3z_1, z_2, z_3,展开点为z0z_0则复平面被分为

  • 0zz0<r10 \leq |z - z_0| < r_1
  • r1<zz0<r2r_1 < |z - z_0| < r_2
  • r2<zz0<r3r_2 < |z - z_0| < r_3
  • r3<zz0<+r_3 < |z - z_0| < +\infty

复闭路积分

  • f(z)f(z)DD内解析,在Γ\Gamma上连续,由柯西积分定理Γf(z)dz=0\oint_\Gamma f(z)\mathrm{d}z = 0
  • f(z)f(z)DD内由唯一的奇点z0z_0,由闭路变形原理:Γf(z)dz=Cf(z)dz\oint_\Gamma f(z)\mathrm{d}z = \oint_C f(z)\mathrm{d}z
    • Cdz(zz0)n={2πin=10n1\oint_C \frac{\mathrm{d}z}{(z - z_0)^n} = \begin{cases} 2\pi i & n = 1\\0 & n\neq 1\end{cases}
    • 则积分Γf(z)dz=2πic1\oint_\Gamma f(z)\mathrm{d}z = 2\pi i c_{-1}c1c_{-1}为洛朗展开负一次项系数

孤立奇点

z0z_0为函数f(z)f(z)的奇点,且存在δ>0\delta > 0使得f(z)f(z)z0z_0的去心领域0<zz0<δ0 < |z - z_0| < \delta内解析,则称z0z_0f(z)f(z)孤立奇点

  • 可去奇点:如果在洛朗级数中不含zz0z - z_0的负幂项,则孤立奇点z0z_0称为f(z)f(z)的可去奇点
  • 极点:如果在洛朗级数中只有有限多个zz0z - z_0的负幂项,且其中关于(zz0)1(z - z_0)^{-1}的最高负幂次为mm,则为mm阶极点
  • 本性奇点:洛朗级数中含有无穷多zz0z - z_0的负幂项,则孤立奇点z0z_0称为f(z)f(z)的本性奇点
    • limzz0f(z)\lim_{z \rightarrow z_0}f(z)不存在且不为\infty

零点

  • f(z0)=0f(z_0) = 0,则称z=z0z = z_0f(z)f(z)的零点
  • f(z)=(zz0)mφ(z)f(z) = (z - z_0)^m\varphi(z)φ(z)\varphi(z)z0z_0处解析且φ(z0)0\varphi(z_0) \neq 0,则称为mm阶零点

对于不恒为0的解析函数,其零点是孤立的
即在零点的一个小领域内,函数无其他零点

  • 定理:如果z0z_0f(z)f(z)mm阶零点,则z0z_0就是1/f(z)1/f(z)mm阶极点,反过来也成立

设函数f(z)f(z)z0z_0处解析,则下列条件是等价的

  • z0z_0f(z)f(z)mm阶零点
  • f(k)(z0)=0,k=0,1,2,,m1;f(m)(z0)0f^{(k)}(z_0) = 0, k =0,1,2,\cdots,m-1; f^{(m)}(z_0) \neq 0

无穷远点是复平面外的理想点,故无穷远点总是函数f(z)f(z)的奇点
如果函数f(z)f(z)在无穷远点z=z = \infty的去心领域R<z<R < |z| < \infty内解析,称点\inftyf(z)f(z)的孤立奇点
做变换w=1/zw = 1/z把扩充zz平面上\infty的去心领域R<z<+R < |z| < +\infty映射称扩展ww平面上原点的去心领域0<w<1/R0 < |w| < 1/R
若记

f(z)=f(1w)=φ(w)limzf(z)=limw0φ(w)f(z) = f(\frac{1}{w}) = \varphi(w) \quad \lim_{z \rightarrow \infty}f(z) = \lim_{w \rightarrow 0}\varphi(w)

函数f(z)f(z)在无穷远点z=z = \infty的性态可由函数φ(w)\varphi(w)在原点w=0w = 0的性态来刻画

  • limzf(z)\lim_{z \rightarrow \infty}f(z)存在且有限,称\inftyf(z)f(z)的可去奇点
    • limzf(z)<\lim_{z \rightarrow \infty}f(z) < \infty
  • limzf(z)\lim_{z \rightarrow \infty}f(z),称\inftyf(z)f(z)的极点
    • f=g(zz0)m,limzf(z)=f = \frac{g}{(z - z_0)^m}, \lim_{z \rightarrow \infty}f(z) = \infty
  • limzf(z)\lim_{z \rightarrow \infty}f(z)不存在且不为无穷,称\inftyf(z)f(z)的本性奇点

留数

z0z_0为函数f(z)f(z)的孤立奇点,将f(z)f(z)z0z_0的去心领域内展开为洛朗级数,称c1c_{-1}是积分过程中唯一残留下来的洛朗系数,称c1c_{-1}f(z)f(z)z0z_0处的留数,记作

Res[f(z),z0]=c1=12πiCf(z)dzRes[f(z), z_0] = c_{-1} = \frac{1}{2\pi i} \oint_C f(z)\mathrm{d}z

留数,就是被留下来的数,c1c_{-1}被留下了,这就是留数了

留数定理

设函数f(z)f(z)在区域DD内除有限个孤立奇点z1,z2,,znz_1, z_2, \cdots, z_n外处处解析,CCDD内包围诸奇点的一条正向简单闭曲线,则

Cf(z)dz=2π1k=1nRes[f(z),zk]\oint_C f(z)\mathrm{d}z = 2\pi 1 \sum^n_{k = 1}Res[f(z), z_k]

  • 如果z0z_0f(z)f(z)的可去奇点,则Res[f(z),z0]=0Res[f(z), z_0] = 0
  • 如果z0z_0f(z)f(z)的本性奇点,则需要洛朗级数展开
  • 如果z0z_0f(z)f(z)的极点,则有对c1c_{-1}的法则

法则:如果z0z_0f(z)f(z)的一阶极点,则

Res[f(z),z0]=limzz0(zz0)f(z)Res[f(z), z_0] = \lim_{z \rightarrow z_0}(z - z_0)f(z)

法则:如果z0z_0f(z)f(z)mm阶极点,则

Res[f(z),z0]=1(m1)!limzz0dm1dzm1{(zz0)mf(z)}Res[f(z), z_0] = \frac{1}{(m - 1)!}\lim_{z \rightarrow z_0} \frac{\mathrm{d}^{m-1}}{\mathrm{d}z^{m-1}} \lbrace (z - z_0)^mf(z) \rbrace

法则:若f(z)=P(z)Q(z)f(z) = \frac{P(z)}{Q(z)}Q(z0)=0Q(z_0) = 0Q(z0)0Q'(z_0) \neq 0P(z0)0P(z_0) \neq 0,且P(z),Q(z)P(z), Q(z)z0z_0点解析,则z0z_0f(z)f(z)的简单极点

Res[f(z),z0]=P(z)Q(z)Res[f(z), z_0] = \frac{P(z)}{Q'(z)}

无穷远点的留数

设函数f(z)f(z)在圆环域R<z<R < |z| < \infty内解析,CC为该圆环域内绕原点的任何一条简单闭曲线,则积分

12πiCf(z)dz\frac{1}{2\pi i} \oint_{C^-} f(z)\mathrm{d}z

f(z)f(z)\infty点的留数,记作

Res[f(z),]=12πiCf(z)dzRes[f(z), \infty] = \frac{1}{2\pi i}\oint_{C^-} f(z)\mathrm{d}z

定理:如果f(z)f(z)在扩充复平面内只有有限个孤立奇点,那么f(z)f(z)在所有各奇点(包括\infty)的留数总和必等于零
法则:

Res[f(z),]=Res[f(1z)1z2,0]Res[f(z), \infty] = -Res[f(\frac{1}{z})\cdot\frac{1}{z^2}, 0]

留数应用

z=eiθz = e^{i\theta}

θ[0,2π]\theta \in [0, 2\pi]为复平面上的单位圆周,dz=ieiθdθ=izdθdθ=dziz\mathrm{d}z = ie^{i\theta}\mathrm{d}\theta = iz\mathrm{d}\theta \Rightarrow \mathrm{d}\theta = \frac{\mathrm{d}z}{iz}

cosθ=eiθ+eiθ2=z+z12=z2+12zsinθ=eiθeiθ2i=zz12i=z212iz\begin{aligned} \cos\theta = \frac{e^{i\theta} + e^{-i\theta}}{2} = \frac{z + z^{-1}}{2} = \frac{z^2 + 1}{2z} \\ \sin\theta = \frac{e^{i\theta} - e^{-i\theta}}{2i} = \frac{z - z^{-1}}{2i} = \frac{z^2 - 1}{2iz} \end{aligned}

即可转换积分

02πR(cosθ,sinθ)dθ=z=1R(z2+12z,z212iz)1izdz=z=1f(z)dz\begin{aligned} \int^{2\pi}_0 R(\cos\theta, \sin\theta)\mathrm{d}\theta &= \oint_{|z| = 1}R(\frac{z^2 + 1}{2z}, \frac{z^2 - 1}{2iz})\frac{1}{iz}\mathrm{d}z \\ &= \oint_{|z| = 1}f(z)\mathrm{d}z \end{aligned}

其中f(z)f(z)zz的有理函数,且在单位圆周z=1|z| = 1上分母不为零,根据留数定数有

02πR(cosθ,sinθ)dθ=z=1f(z)dz=2πik=1nRes[f(z),zk]\int^{2\pi}_0 R(\cos\theta, \sin\theta)\mathrm{d}\theta = \oint_{|z| = 1}f(z)\mathrm{d}z = 2\pi i \sum^{n}_{k=1}Res[f(z), z_k]

其中zk(k=1,2,,n)z_k(k=1,2,\cdots,n)为单位圆z=1|z| = 1内的f(z)f(z)的孤立奇点


形如

+R(x)dz\int^{+\infty}_{-\infty}R(x)\mathrm{d}z

  • R(x)=P(x)Q(x)R(x) = \frac{P(x)}{Q(x)},其中P,QP, Q为多项式
  • 分母Q(x)Q(x)的次数比分子P(x)P(x)的次数至少高2次
  • 分母Q(x)Q(x)无实零点

+R(x)d=2πikRes[R(z),zk]\int^{+\infty}_{-\infty}R(x)\mathrm{d} = 2\pi i \sum_{k}Res[R(z), z_k]

这里zkz_kR(z)R(z)在上半平面的所有孤立奇点


形如

+R(x)eiaxdx(a>0)\int^{+\infty}_{-\infty}R(x)e^{iax}\mathrm{d}x(a > 0)

  • R(x)=P(x)Q(x)R(x) = \frac{P(x)}{Q(x)},其中P,QP, Q为多项式
  • 分母Q(x)Q(x)的次数比分子P(x)P(x)的次数至少高1次
  • 分母Q(x)Q(x)无实

+R(x)eiaxdx=2πikRes[R(z)eiaz,zk]\int^{+\infty}_{-\infty}R(x)e^{iax}\mathrm{d}x = 2\pi i \sum_{k}Res[R(z)e^{iaz}, z_k]

这里zkz_kR(z)R(z)在上半平面的所有孤立奇点

共形映射

伸缩率:经过z0z_0点的曲线C0C_0w=f(z)w = f(z)映射后,变成了过w0w_0点的曲线Γ0\Gamma_0,曲线被伸缩和选中

\lim_{\substack{z \rightarrow z_0 \\ C_0}} \frac{|w - w_0|}{|z - z_0|} = \lim_{\substack{\Delta z \rightarrow 0 \\ C_0}} \frac{|\Delta w|}{|\Delta z|}

为曲线C0C_0w=f(z)w = f(z)映射后,在z0z_0点的伸缩率

\lim_{\substack{z \rightarrow z_0 \\ C_0}} (\varphi - \theta) = \varphi_0 - \theta_0

为曲线C0C_0w=f(z)w = f(z)映射后,在z0z_0点的旋转角

导数的几何意义

  • f(z0)|f'(z_0)|为曲线C0C_0z0z_0点的伸缩率
  • argf(z0)\arg f'(z_0)为曲线C0C_0z0z_0点的旋转角

  • 伸缩率不变性:任何一条经过z0z_0点的曲线的伸缩率均为f(z0)|f'(z_0)|
  • 旋转角不变性:任何一条经过z0z_0点的曲线的旋转角均为argf(z0)\arg f'(z_0)
  • 保角性:即w=f(z)w = f(z)保持了两条曲线的夹角的大小与方向不变

  • 第一类保角映射
    • 函数w=f(z)w = f(z)在区域DD内满足
      • 保角性
      • 伸缩率不变性
      • f(z)0f'(z) \neq 0
  • 第二类保角映射
    • 函数w=f(z)w = f(z)在区域DD内满足
      • 能保持两条曲线的交角的大小不变,但方向相反
      • 伸缩率不变性

共性映射:若函数w=f(z)w = f(z)在区域DD的第一类保角映射,且当z1z2z_1 \neq z_2时,f(z1)f(z2)f(z_1) \neq f(z_2),则称w=f(z)w = f(z)为区域DD内的共形映射,要求函数必须时一一映射(双方单值)

解析 + f(z)0+f'(z) \neq 0 + 一一映射 = 共性映射

  • 保域性定理:设函数w=f(z)w = f(z)在区域DD内解析,且不恒为常数,则其像集合G=f(D)G = f(D)仍位区域
  • 边界对应原理:设区域DD的边界位简单闭曲线CC,函数w=f(z)w = f(z)的闭区域D=D+C\overline{D} = D + C上解析,且将曲线C双方单值地映射位简单闭曲线Γ\Gamma
    • zz沿CC的正向绕行时,相应的ww的绕行方向定为Γ\Gamma的正向,并令GG是以Γ\Gamma为边界的区域,则w=f(z)w = f(z)DD共性映射为GG
  • 黎曼存在唯一性定理
    • DDGG是任给的两个单连域,在其各自的边界上至少含有两个点,则一定存在解析函数w=f(z)w = f(z)将区域DD共形映射为GG
    • 如果在区域DDGG内再分别任意指定一点z0z_0w0w_0,并任给一实数θ0(π,π]\theta_0 \in (-\pi, \pi],使得argf(z0)=θ0\arg f'(z_0) = \theta_0f(z0)=w0f(z_0) = w_0,则映射w=f(z)w = f(z)是唯一的

对给定的单连域DD,求共性映射,使得DD映射为单位圆域,事实上,由此即可求得任意两个单连域之间的共性映射(单位圆域为中介)

D(z)ξ=g(z)单位圆域ξ=h(w)G(w)w=h1(g(z))=f(z)\begin{gathered} D(z) \xrightarrow{\xi = g(z)} \text{单位圆域} \xleftarrow{\xi = h(w)} G(w) \\ \Downarrow \\ w = h^{-1}(g(z)) = f(z) \end{gathered}

分式线性映射

由分式线性函数w=az+bcz+d,(acbd)w = \frac{az + b}{cz + d}, (\frac{a}{c} \neq \frac{b}{d})构成的映射,若c=0c = 0,则为整式线性映射

  • 两个分式线性映射的复合,仍是一个分式线性映射
  • 分式线性映射的逆隐射也是一个分式线性映射

w=2zz+i=2z+2i2iz+i=2+2iz+i=2+2eπ2i1z+izz+iz11z1z2eπ2iz2z32z3z4z4+2w\begin{gathered} w = \frac{2z}{z + i} = \frac{2z + 2i - 2i}{z + i} = 2 + \frac{-2i}{z + i} = 2 + 2e^{-\frac{\pi}{2}i}\cdot \frac{1}{z + i} \\ z \xrightarrow{z + i} z_1 \xrightarrow{\frac{1}{z_1}} z_2\xrightarrow{e^{-\frac{\pi}{2}i}z_2} z_3 \xrightarrow{2z_3} z_4 \xrightarrow{z_4 + 2} w \end{gathered}

一个一般形式的分式线性映射可以由几种最简单发分式线性映射复合而成

  • 整式线性映射
    • w=z+b(bZ)w = z + b(b \in \mathbb{Z})
    • w=eiθ0z(θ0R)w = e^{i\theta_0}z(\theta_0 \in \mathbb{R})
    • w=rz(r>0)w = rz(r > 0)
  • 分式线性映射
    • w=eiθ0z(θ0R)w = e^{i\theta_0}z(\theta_0 \in \mathbb{R})
    • w=rz(r>0)w = rz(r > 0)
    • w=1zw = \frac{1}{z}

  • 平移映射:w=z+bw = z + b
    • 将点集沿着向量b\overline{b}的方向平移一段距离b|\overline{b}|
  • 选择映射:w=eiθ0zw = e^{i\theta_0}z
    • 将点集绕原点旋转一个角度θ0\theta_0
      • θ0>0\theta_0 > 0:逆时针旋转
      • θ0<0\theta_0 < 0:顺时针旋转
  • 相似映射:w=rzw = rz
    • 将点集沿极径方向相似地扩大(缩小)rr
  • 反演映射:w=1zw = \frac{1}{z}
    • 将单位圆内(外)的点映射到单位圆外(内)的点
    • 它们的模互为导数,辐角反号

定义:若OAOB=R2\overline{OA}\cdot\overline{OB} = R^2,则称AA点和BB点关于圆周CC对称
保对称点性:设点z1,z2z_1,z_2关于圆周CC对称,则在分式线性映射下,它们的像点w1,w2w_1,w_2也关于像曲线CC'对称

分式线性映射特性

  • 保形性
  • 解析性
  • 保圆性

分式线性映射再扩充复平面上是共形映射
导数映射将圆变为圆或直线,将直线变为圆或直线

  • 线性映射将圆变成圆,将直线变成直线
  • 将直线看作是半径为无穷大的圆

在扩充复平面上,分式线性映射能把圆变成圆(可能是直线)

  • 如果给定的圆上没有点映射成无穷远点,则映射为半径有限的圆
  • 如果给定的圆上有一点映射成无穷远点,则映射为直线
  • 对于圆弧段,如果其中一个端点映射成无穷远点,则映射为射线

给定三个条件,就能决定一个分式线性映射

定理:在zz平面上任给三个不同的点z1,z2,z3z_1,z_2,z_3,在ww平面上也任给三个不同的点w1,w2,w3w_1,w_2,w_3,则存在唯一的分式线性映射,将z1,z2,z3z_1,z_2,z_3分别依次映射为w1,w2,w3w_1,w_2,w_3

ww1ww2:w3w1w3w2=zz1zz2:z3z1z3z2\frac{w - w_1}{w - w_2}:\frac{w_3 - w_1}{w_3 - w_2} = \frac{z - z_1}{z - z_2}:\frac{z_3 - z_1}{z_3 - z_2}

称为对应点公式,由于分式线性映射具有保圆性,因此该公式通常用于:把过z1,z2,z3z_1,z_2,z_3三点的圆(或圆弧)映射为过w1,w2,w3w_1,w_2,w_3的圆(或圆弧)
如果z1,z2,z3z_1,z_2,z_3w1,w2,w3w_1,w_2,w_3中有一个为无穷远点\infty,则只需将对应点公式中含有\infty的项换成11


w=f(z)w = f(z)为分式线性映射,且f(z1)=w1,f(z2)=w2f(z_1) = w_1, f(z_2) = w_2,则它可表示为:ww1www=kzz1zz2\frac{w - w_1}{w - w_w} = k\frac{z - z_1}{z - z_2},特别的若f(z1)=0,f(z2)=f(z_1) = 0, f(z_2) = \infty,则w=kzz1zz2w = k\frac{z - z_1}{z - z_2}

把过z1,z2z_1, z_2的点的弧映射成过原点的直线,若左式作为中间步骤,则kk可设为11

其他函数映射

  • 幂函数:w=znw = z^n
    • 扩大顶点在原点的角形域
  • 根式函数:w=nzw = ^n\sqrt{z}
    • 缩小顶点在原点的角形域

幂函数w=znw = z^nzz平面上除原点外是第一类保角映射
在角形域或扇形域(0<θ<θ0)(0 < \theta < \theta_0)上,如果θ02πn\theta_0 \leq \frac{2\pi}{n},则幂函数w=znw = z^n是共形映射

  • 指数函数:w=ezw = e^z
    • 将水平带形域(半带形域)变为角形域(或扇形域)

幂函数w=znw = z^nzz平面上是第一类保角映射
在水平带形域(0<y<h)(0 < y < h)上,如果h2πh \leq 2\pi,则指数函数w=ezw = e^z是共形映射

傅里叶变换

  • 对于给定的函数,其傅里叶级数展开式是唯一的
    • 傅里叶系数是固定的
  • 在计算展开系数cnc_n时,可在任意一个长度为TT的区间上计算其中的积分
  • 采用周期延拓级数,可以将结论应用到仅仅定义在某个有限区间上的函数
    • 对于定于在有限区间的函数,同样可以展开为傅里叶级数

任何一个周期性函数,都可以分解为正余弦函数的和

傅里叶级数的三角形式

fT(t)f_T(t)是以TT为周期的实值函数,且在[T2,T2][-\frac{T}{2}, \frac{T}{2}]上满足条件

  • 连续或只有有限个第一类间断点
  • 只有有限个极值点

ω0=2π/T\omega_0 = 2\pi /T称为基频,则在fT(t)f_T(t)的连续点处,有

fT(t)=a02+n=1+(ancosnω0t+bnsinnω0t)f_T(t) = \frac{a_0}{2} + \sum^{+\infty}_{n=1}(a_n\cos n\omega_0 t + b_n \sin n\omega_0 t)

在间断点处,上式左端为12[fT(t+0)+fT(t0)]\frac{1}{2}[f_T(t + 0) + f_T(t - 0)],其中

an=2TT/2T/2fT(t)cosnω0tdt(n=1,2,)bn=2TT/2T/2fT(t)sinnω0tdt(n=1,2,)\begin{aligned} a_n &= \frac{2}{T}\int^{T/2}_{-T/2}f_T(t)\cos n\omega_0 t\mathrm{d}t (n = 1,2,\cdots) \\ b_n &= \frac{2}{T}\int^{T/2}_{-T/2}f_T(t)\sin n\omega_0 t\mathrm{d}t (n = 1,2,\cdots) \end{aligned}

由欧拉公式

cosnω0t=einω0t+einω0t2sinnω0t=einω0t+einω0t2\begin{aligned} \cos n\omega_0 t &= \frac{e^{in\omega_0t} + e^{-in\omega_0t}}{2} \\ \sin n\omega_0 t &= \frac{-e^{in\omega_0t} + e^{-in\omega_0t}}{2} \end{aligned}

带入上式整理可得

fT(t)=a02+n=1+(anibn2einω0t+an+jbn2ejnω0t)f_T(t) = \frac{a_0}{2} + \sum^{+\infty}_{n=1}(\frac{a_n - ib_n}{2}e^{in\omega_0t} + \frac{a_n + jb_n}{2}e^{-jn\omega_0t})

A0=a02,An=an2+bn2,cosθn=anAn,sinθn=bnAnA_0 = \frac{a_0}{2}, A_n = \sqrt{a^2_n + b^2_n}, \cos\theta_n = \frac{a_n}{A_n}, \sin\theta_n = \frac{-b_n}{A_n},则

fT(t)=A0+n=1+Ancos(nω0t+θn)\color{red}f_T(t) = A_0 + \sum^{+\infty}_{n=1}A_n\cos(n\omega_0t + \theta_n)

  • 周期信号可以分解为一系列固定的简谐波之和,这些简谐波的角频率分别为一个基频ω0\omega_0的倍数
  • 任何一个周期为TT的周期信号fT(t)f_T(t)并不包含所有的频率成分,其频率是以基频ω0\omega_0为间隔离散取值的

这是周期信号的一个非常重要的特点

  • 振幅AnA_n:反映了在信号fT(t)f_T(t)中频率为nω0n\omega_0的简谐波所占有的份额
  • 相位θn\theta_n:反映了在信号fT(t)f_T(t)中频率为nω0n\omega_0的简谐波沿时间轴移动的大小

这2各指标完全的定理的刻画了信号的频率特性

傅里叶级数的指数形式

c0=a02,cn=anibn2,cn=an+jbn2c_0 = \frac{a_0}{2}, c_n = \frac{a_n - ib_n}{2}, c_{-n} = \frac{a_n + jb_n}{2},则有

fT(t)=n=+cneinω0tcn=1TT/2T/2fT(t)einω0tdtf_T(t) = \sum^{+\infty}_{n = -\infty}c_ne^{in\omega_0t} \quad c_n = \frac{1}{T}\int^{T/2}_{-T/2}f_T(t)e^{-in\omega_0t}\mathrm{d}t

称系数cnc_n为(离谱)频谱,记为F(nω0)=cnF(n\omega_0) = c_n,它的模cn|c_n|恰好反映了振幅,而它的辐角argcn\arg c_n就是相位

  • 振幅谱:cn|c_n|
  • 相位谱:argcn\arg c_n
  • 频谱:cnc_n

einω0te^{in\omega_0t}可以看作一组正交基

傅里叶积分公式

借助傅里叶级数的展开,可以了解信号的频率特性,认清信号的本质,这样的分析手段叫做频谱分析,但是傅里叶级数要求展开函数必须时周期函数,而实际中大部分时非周期函数,所以需要进行转换

没有周期的函数可以看作为一个周期无穷大的一个周期函数

设函数f(t)f(t)满足

  • (,+)(-\infty,+\infty)上的任一有限区间内满足Dirichlet条件
    • 在一周期内,连续或只有有限个第一类间断点;
    • 在一周期内,极大值和极小值的数目应是有限个;
    • 在一周期内,信号是绝对可积的。
  • 绝对可积,即:+f(t)dt<+\int^{+\infty}_{-\infty}|f(t)|\mathrm{d}t < +\infty

则在f(t)f(t)的连续点处,有

f(t)=12π+[+f(t)eiωtdt]eiωtdωf(t) = \frac{1}{2\pi}\int^{+\infty}_{-\infty}\left [\int^{+\infty}_{-\infty}f(t)e^{-i\omega t}\mathrm{d}t\right ]e^{i\omega t}\mathrm{d}\omega

在间断点处,公式的左端应为12[f(t+0)+f(t0)]\frac{1}{2}[f(t + 0) + f(t - 0)]

傅里叶正变换

F(ω)=+f(t)eiωtdt=记为F(f(t))\color{red}F(\omega) = \int^{+\infty}_{-\infty}f(t)e^{-i\omega t}\mathrm{d}t \xlongequal{\text{记为}} \mathfrak{F}(f(t))

傅里叶逆变换

f(t)=12π+F(ω)eiωtdω=记为F1(f(ω))\color{red}f(t) = \frac{1}{2\pi}\int^{+\infty}_{-\infty}F(\omega)e^{i\omega t}\mathrm{d}\omega \xlongequal{\text{记为}} \mathfrak{F}^{-1}(f(\omega))

  • 象函数:F(ω)F(\omega)
  • 象原函数:f(t)f(t)
  • 傅里叶变换对:f(t)F(ω)f(t) \leftrightarrow F(\omega)

广义积分取柯西主值

傅里叶变换的物理意义

傅里叶变换刻画了非周期函数的频谱特性,非周期函数的频谱是连续取值的,像函数F(w)F(w)反映的是函数f(t)f(t)中各频率分量的分布密度,它是复值函数,故可表示为

F(ω)=F(ω)eiargF(ω)F(\omega) = |F(\omega)|e^{i\arg F(\omega)}

F(ω)F(\omega)为频谱密度函数

  • 振幅谱:F(ω)|F(\omega)|
  • 相位谱:argF(ω)\arg F(\omega)

傅里叶常用公式

+sinaxxdx=π(a>0)0+sinaxxdx=π2(a>0)\begin{aligned} &\int^{+\infty}_{-\infty} \frac{\sin ax}{x}\mathrm{d}x = \pi \quad (a > 0) \\ &\int^{+\infty}_{0} \frac{\sin ax}{x}\mathrm{d}x = \frac{\pi}{2} \quad (a > 0) \end{aligned}


1π+sinωtωdω={1t>00t=01t<0==sgntsgnt2jω\begin{gathered} \frac{1}{\pi}\int^{+\infty}_{-\infty} \frac{\sin \omega t}{\omega}\mathrm{d}\omega = \begin{cases} 1 & t > 0 \\ 0 & t = 0 \\ -1 & t < 0 \end{cases} == \text{sgn} t \\ \text{sgn} t \leftrightarrow \frac{2}{j\omega} \end{gathered}

单位冲激函数

一些常用函数不能进行傅里叶变换,还有许多瞬时物理量不能用通常的函数形式来描述
比如,对线密度函数进行积分可以得到质量

Pa(x)={maoxa0其他+Pa(x)dx=m\begin{gathered} P_a(x) = \begin{cases} \frac{m}{a} & o \leq x \leq a \\ 0 & \text{其他} \end{cases} \\ \int^{+\infty}_{-\infty} P_a(x)\mathrm{d}x = m \end{gathered}

单位冲激函数δ(t)\delta(t)满足

  • t0t \neq 0时,δ(t)=0\delta(t) = 0
  • +δ(t)dt=1\int^{+\infty}_{-\infty} \delta(t)\mathrm{d}t = 1

单位冲激函数δ(t)\delta(t)不是经典意义下的函数,通常称其为广义函数
不能通过常规意义下的值的对应关系来理解和使用,而总是通过它的性质来使用它

单位冲激函数性质

筛选性质:设函数f(t)f(t)是定义在(,+)(-\infty, +\infty)上的有界函数,且在t=0t = 0处连续,则有

+δ(t)f(t)dt=f(0)\int^{+\infty}_{-\infty} \delta(t)f(t)\mathrm{d}t = f(0)

若在t=t0t = t_0处连续,则有

+δ(tt0)f(t)dt=f(t0)\int^{+\infty}_{-\infty} \delta(t - t_0)f(t)\mathrm{d}t = f(t_0)

对称函数:单位冲激函数为偶函数,即

δ(t)=δ(t)\delta(t) = \delta(-t)

积分性质:设

u(t)={1t>00t<0u(t) = \begin{cases} 1 & t > 0 \\ 0 & t < 0 \end{cases}

则有

tδ(t)dt=u(t)u(t)=δ(t)\int^t_{-\infty} \delta(t) \mathrm{d}t = u(t) \quad u'(t) = \delta(t)

称函数u(t)u(t)为单位阶跃函数(Heaviside)函数

单位冲激函数的图形采用一条从原点出发且长度为1的有向线段来表示,线段长度称为冲激强度

单位冲激函数的傅里叶变换

F(δ(t))=+eiωtdt=eiωtt=0=1δ(t)1\begin{gathered} \mathfrak{F}(\delta(t)) = \int^{+\infty}_{-\infty} e^{-i\omega t}\mathrm{d}t = e^{-i\omega t}|_{t=0} = 1 \\ \delta(t) \leftrightarrow 1 \end{gathered}

单位冲激函数包含所有频率成分,且它们具有相等的幅度,称此为均匀频谱

F1(1)=12π+1eiωtdω=δ(t)\mathfrak{F}^{-1}(1) = \frac{1}{2\pi} \int^{+\infty}_{-\infty} 1 \cdot e^{i\omega t}\mathrm{d}\omega = \delta(t)

由此可以得到重要公式

+eiωtdω=2πδ(t)+teiωtdω=2πiδ(t)\begin{aligned} \int^{+\infty}_{-\infty} e^{i\omega t}\mathrm{d}\omega &= 2\pi \delta(t) \\ \int^{+\infty}_{-\infty} t \cdot e^{i\omega t}\mathrm{d}\omega &= 2\pi i \delta'(t) \end{aligned}

δ\delta函数的傅里叶变换中,其广义积分是根据δ\delta函数的性质给出,所以无法直接计算出积分,称这种方式的傅里叶变换为广义傅里叶变换

傅里叶变换性质

线性性质:设a,ba, b为常数,则

F(af(t)+bg(t))=aF(ω)+bG(ω)\mathfrak{F}(af(t) + bg(t)) = aF(\omega) + bG(\omega)

位移性质:设t0,ω0t_0, \omega_0为常数,则

F(f(tt0))=eiωt0F(ω)时移性质F1(F(ωω0))=eiω0tf(t)频移性质\begin{aligned} \mathfrak{F}(f(t - t_0)) = e^{-i\omega t_0}F(\omega) \quad \text{时移性质} \\ \mathfrak{F}^{-1}(F(\omega - \omega_0)) = e^{i\omega_0 t}f(t) \quad \text{频移性质} \end{aligned}

相似性质:设aa为非零常数,则

F(f(at))=1aF(ωa)\mathfrak{F}(f(at)) = \frac{1}{|a|}F(\frac{\omega}{a})

微分性质:若limt+f(t)=0\lim_{|t|\rightarrow +\infty}f(t) = 0,则

F(f(t))=iωF(ω)\mathfrak{F}(f'(t)) = i\omega F(\omega)

一般的,若limt+f(k)(t)=0\lim_{|t|\rightarrow +\infty}f^{(k)}(t) = 0,则

F(f(k)(t))=(iω)kF(ω)\mathfrak{F}(f^{(k)}(t)) = (i\omega)^k F(\omega)

同理,可得到像函数的导数公式

F1(F(ω))=itf(t)F1(F(n)(ω))=(it)nf(t)\begin{aligned} \mathfrak{F}^{-1}(F'(\omega)) &= -itf(t) \\ \mathfrak{F}^{-1}(F^{(n)}(\omega)) &= -(it)^nf(t) \end{aligned}

积分性质

F(tf(t)dt)=1iωF(ω)\mathfrak{F}(\int^t_{-\infty}f(t)\mathrm{d}t) = \frac{1}{i\omega}F(\omega)

帕塞瓦尔等式

+f2(t)dt=12π+F(ω)2dω\int^{+\infty}_{-\infty} f^2(t)\mathrm{d}t = \frac{1}{2\pi}\int^{+\infty}_{-\infty}|F(\omega)|^2\mathrm{d}\omega

卷积与卷积定理

设函数f1(t)f_1(t)f2(t)f_2(t)在区间(,+)(-\infty, +\infty)上有定义,如果广义积分+f1(τ)f2(tτ)dτ\int^{+\infty}_{-\infty}f_1(\tau)f_2(t - \tau)\mathrm{d}\tau对任何实数tt都收敛,则他在(,+)(-\infty, +\infty)上定义了一个自变量为tt的函数,称此函数为f1(t)f_1(t)f2(t)f_2(t)的卷积,记为f1(t)f2(t)f_1(t)*f_2(t),即

卷积是一种新的运算,两函数相乘并积分

f1(t)f2(t)=+f1(τ)f2(tτ)dτf_1(t)*f_2(t) = \int^{+\infty}_{-\infty}f_1(\tau)f_2(t - \tau)\mathrm{d}\tau

  • 交换律:f1(t)f2(t)=f2(t)f1(t)f_1(t)*f_2(t) = f_2(t)*f_1(t)
  • 结合律:f1(t)[f2(t)f3(t)]=[f1(t)f2(t)]f3(t)f_1(t)*[f_2(t)*f_3(t)] = [f_1(t)*f_2(t)]*f_3(t)
  • 分配律:f1(t)[f2(t)+f3(t)]=f1(t)f2(t)+f1(t)f3(t)f_1(t)*[f_2(t)+f_3(t)] = f_1(t)*f_2(t) + f_1(t)*f_3(t)

采用图形方式则比较容易确定积分限
卷积有反褶,平移,相乘,积分四个步骤组成

  • 首先将函数g(τ)g(\tau)反褶并平移到tt,得到g(tτ)g(t - \tau)
  • 再与函数f(τ)f(\tau)相乘后求积分,得到卷积f1(t)f2(t)f_1(t)*f_2(t)

F(f1(t))=F1(ω)\mathfrak{F}(f_1(t)) = F_1(\omega)F(f2(t))=F2(ω)\mathfrak{F}(f_2(t)) = F_2(\omega),则有

F(f1(t)f2(t))=F1(ω)F2(ω)F1(F1(ω)F2(ω))=2πf1(t)f2(t)\begin{aligned} \mathfrak{F}(f_1(t) * f_2(t)) = F_1(\omega) \cdot F_2(\omega) \\ \mathfrak{F}^{-1}(F_1(\omega) * F_2(\omega)) = 2\pi f_1(t) \cdot f_2(t) \end{aligned}

拉普拉斯变换

  • 傅里叶变换具有局限性,只有函数f(t)f(t)满足Dirichlet条件,且在(,+)(-\infty,+\infty)上绝对可积时,才能进行古典意义下的傅里叶变换
  • 绝对可积限制过大
  • 许多函数在自变量维负时为0或无意义
  • 没有必要或者不可能在这个实轴上进行

改善:

  • 将函数f(t)f(t)乘以一个单位阶跃函数u(t)u(t),使得函数在t<0t < 0的部分补零
  • 将函数再乘上一个衰减指数函数eβt(β>0)e^{-\beta t}(\beta > 0),使得函数在t>0t > 0的部分尽快地衰减下来

就有希望使得函数f(t)u(t)eβtf(t)u(t)\cdot e^{\beta t}满足傅里叶变换条件

F(f(t)u(t)eβt)=+f(t)u(t)eβteiωtdt=0+f(t)e(β+iω)tdt\begin{aligned} \mathfrak{F}(f(t)u(t)e^{-\beta t}) &= \int^{+\infty}_{-\infty} f(t)u(t)e^{-\beta t}e^{-i\omega t}\mathrm{d}t \\ &= \int^{+\infty}_{0} f(t) e^{-(\beta + i\omega)t}\mathrm{d}t \end{aligned}

将上式中的β+iω\beta + i\omega记为ss,就得到一种新的变换

F(s)=0+f(t)estdtF(s) = \int^{+\infty}_{0} f(t) e^{-st}\mathrm{d}t

上述广义积分存在的关键:变量ss的实部Res=β\text{Res} = \beta足够大


设函数f(t)f(t)是定义在(0,+)(0, +\infty)上的实值函数,如果对于复参数s=β+iωs = \beta + i\omega,积分F(s)=0+f(t)estdtF(s) = \int^{+\infty}_{0} f(t) e^{-st}\mathrm{d}t在复平面ss的某一区域内收敛,则称F(s)F(s)f(t)f(t)拉普拉斯变换像函数,记为F(s)=LF(s) = \mathfrak{L},即

F(s)=L(f(t))=0+f(t)estdtF(s) = \mathfrak{L}(f(t)) = \int^{+\infty}_{0} f(t) e^{-st}\mathrm{d}t

相应的,称f(t)f(t)F(s)F(s)的拉普拉斯逆变换或像原函数,记为

f(t)=L1(F(s))f(t) = \mathfrak{L}^{-1}(F(s))

存在性定理:
设函数f(t)f(t)t0t \geq 0时,满足

  • 在任何有限区间上分段连续
  • 具有有限的增长性
    • 即存在常数ccM>0M > 0,使得f(t)Mect|f(t)| \geq Me^{ct}
    • cc称为增长指数

则象函数F(s)F(s)在半平面Res>cRes > c上一定存在且解析

  • 象函数F(s)F(s)的存在域一般是一个右半平面Res>cRes > c,即只要复数ss的实部足够大就可以了
  • 在拉普拉斯变换中的函数一般均约定在t<0t < 0时为零,即函数f(t)f(t)等价于f(t)u(t)f(t)u(t)

常见拉普拉斯变换

L(1)=L(u(t))=1sL(tm)=m!sm+1=Γ(m+1)sm+1L(eat)=1saL(cosat)=ss2+a2L(sinat)=as2+a2\begin{aligned} \mathfrak{L}(1) &= \mathfrak{L}(u(t)) = \frac{1}{s} \\ \mathfrak{L}(t^m) &= \frac{m!}{s^{m+1}} = \frac{\Gamma(m+1)}{s^{m+1}} \\ \mathfrak{L}(e^{at}) &= \frac{1}{s - a} \\ \mathfrak{L}(\cos at) &= \frac{s}{s^2 + a^2} \\ \mathfrak{L}(\sin at) &= \frac{a}{s^2 + a^2} \end{aligned}

拉普拉斯变换性质

线性性质:设a,ba,b为常数,则

L(af(t)+bg(t))=aF(s)+bG(s)L1(aF(t)+bG(t))=af(s)+bg(t)\begin{aligned} \mathfrak{L}(af(t) + bg(t)) &= aF(s) + bG(s) \\ \mathfrak{L}^{-1}(aF(t) + bG(t)) &= af(s) + bg(t) \end{aligned}

相似性质:设aa为正实数,则

L(f(at))=1aF(sa)\mathfrak{L}(f(at)) = \frac{1}{a}F(\frac{s}{a})

延迟性质:设当t<0t < 0时,f(t)=0f(t) = 0,则对任一非负实数τ\tau

L(f(tτ))=esτF(s)\mathfrak{L}(f(t - \tau)) = e^{s\tau}F(s)

强调了t<0t < 0时,f(t)=0f(t) = 0

L(f(tτ)u(tτ))=esτF(s)L1(esτF(s))=f(tτ)u(tτ)\begin{aligned} \mathfrak{L}(f(t - \tau)u(t - \tau)) &= e^{s\tau}F(s) \\ \mathfrak{L}^{-1}(e^{-s\tau}F(s)) &= f(t \tau)u(t - \tau) \end{aligned}

位移性质:设aa为任一复常数,则

L(eatf(t))=F(sa)\mathfrak{L}(e^{at}f(t)) = F(s - a)

微分性质:导数的像函数

L(f(t))=sF(s)f(0)\mathfrak{L}(f'(t)) = sF(s) - f(0)

可用于求解微分方程组的初值问题

像函数的导数

F(s)=L(tf(t))一般的F(n)(s)=(1)nL(tnf(t))\begin{gathered} F'(s) = -\mathfrak{L}(tf(t)) \\ \text{一般的} \\ F^{(n)}(s) = (-1)^n\mathfrak{L}(t^nf(t)) \end{gathered}

积分性质
一个函数的象函数与它积分之后的象函数的关系

L(0tf(t)dt)=1sF(s)\mathfrak{L}(\int^t_0 f(t)\mathrm{d}t) = \frac{1}{s}F(s)

一般的

L(0tdt0tdt0tf(t)dt)=1snF(s)\mathfrak{L}(\int^t_0 \mathrm{d}t\int^t_0 \mathrm{d}t\cdots \int^t_0 f(t)\mathrm{d}t) = \frac{1}{s^n}F(s)

象函数的原函数与它积分以后的原函数的关系

sF(s)ds=L(f(t)t)\int^{\infty}_s F(s)\mathrm{d}s = \mathfrak{L}(\frac{f(t)}{t})

一般的

sdssdssF(s)ds=L(f(t)tn)\int^{\infty}_s\mathrm{d}s\int^{\infty}_s\mathrm{d}s\cdots\int^{\infty}_s F(s)\mathrm{d}s = \mathfrak{L}(\frac{f(t)}{t^n})

拉普拉斯卷积与卷积定理

t<0t < 0时,f1(t)=f2(t)=0f_1(t) = f_2(t) = 0,则有

f1(t)f2(t)=0tf1(τ)f2(tτ)dτf_1(t)*f_2(t) = \int^{t}_{0}f_1(\tau)f_2(t - \tau)\mathrm{d}\tau

上式仍满足交换律,结合律以及分配律

L(f1(t)f2(t))=F1(s)F2(s)\begin{aligned} \mathfrak{L}(f_1(t) * f_2(t)) = F_1(s) \cdot F_2(s) \end{aligned}

拉普拉斯逆变换

也称反演积分公式,由拉普拉斯变换与傅里叶变换的关系可知:函数f(t)f(t)的拉普拉斯变换F(s)=F(β+iω)F(s) = F(\beta + i\omega)就是函数f(t)u(t)eβtf(t)u(t)e^{-\beta t}的傅里叶变换

F(s)=F(β+iω)=+(f(t)u(t)eβt)eiωtdtf(t)u(t)eβt=12π+F(β+iω)eiωtdωf(t)u(t)=12πjβjωβ+iωF(s)estds\begin{aligned} F(s) &= F(\beta + i\omega) = \int^{+\infty}_{-\infty}(f(t)u(t)e^{-\beta t})e^{-i\omega t}\mathrm{d}t \\ f(t)u(t)e^{-\beta t} &= \frac{1}{2\pi}\int^{+\infty}_{-\infty}F(\beta + i\omega)e^{i\omega t}\mathrm{d}\omega \\ &\Downarrow \\ f(t)u(t) &= \frac{1}{2\pi j} \int^{\beta + i\omega}_{\beta - j\omega} F(s)e^{st}\mathrm{d}s \end{aligned}

即可得到

f(t)=12πjβjωβ+iωF(s)estds(t>0)f(t) = \frac{1}{2\pi j} \int^{\beta + i\omega}_{\beta - j\omega} F(s)e^{st}\mathrm{d}s \quad (t > 0)

我们就有了如下的拉普拉斯变换对

{F(s)=0+f(t)eβtestdtf(t)=12πjβjωβ+iωF(s)estds(t>0)\begin{cases} F(s) = \int^{+\infty}_{0}f(t)e^{-\beta t}e^{-st}\mathrm{d}t \\ f(t) = \frac{1}{2\pi j} \int^{\beta + i\omega}_{\beta - j\omega} F(s)e^{st}\mathrm{d}s \quad (t > 0) \end{cases}

利用留数计算反演积分:设函数F(s)F(s)除在半平面Res<cRes < c内有有限个孤立奇点s1,s2,,sns_1, s_2, \cdots, s_n外是解析的,且当ss \rightarrow \infty时,F(s)0F(s) \rightarrow 0,则

f(t)=12πiβjωβ+iωF(s)estds=k=1nRes[F(s)est,sk](t>0)\begin{aligned} f(t) &= \frac{1}{2\pi i}\int^{\beta + i\omega}_{\beta - j\omega} F(s)e^{st}\mathrm{d}s \\ &= \sum^n_{k=1} Res[F(s)e^{st}, s_k] \quad (t > 0) \end{aligned}

或者利用查表来计算

L1(1s)=1L1(m!sm+1)=tmL1(ss2+b2)=cosbtL1(bs2+b2)=sinbtL1(1)=δ(t)L1(1sa)=eatL1(m!(sa)m+1)=eattmL1(sa(sa)2+b2)=eatcosbtL1(b(sa)2+b2)=eatsinbt\begin{aligned} \mathfrak{L}^{-1}(\frac{1}{s}) &= 1 \\ \mathfrak{L}^{-1}(\frac{m!}{s^{m+1}}) &= t^m \\ \mathfrak{L}^{-1}(\frac{s}{s^2 + b^2}) &= \cos bt \\ \mathfrak{L}^{-1}(\frac{b}{s^2 + b^2}) &= \sin bt \\ \mathfrak{L}^{-1}(1) &= \delta(t) \\ \mathfrak{L}^{-1}(\frac{1}{s - a}) &= e^{at} \\ \mathfrak{L}^{-1}(\frac{m!}{(s - a)^{m+1}}) &= e^{at}t^m \\ \mathfrak{L}^{-1}(\frac{s-a}{(s-a)^2 + b^2}) &= e^{at}\cos bt \\ \mathfrak{L}^{-1}(\frac{b}{(s-a)^2 + b^2}) &= e^{at}\sin bt \end{aligned}

求解常微分方程组

            graph LR
            微分方程组--拉普拉斯正变换-->象函数的代数方程组
象函数的代数方程组--求解-->得到象函数
得到象函数--拉普拉斯逆变换-->微分方程组的解
          

参考