复数与复平面 Complex algebra and the complex plane
复数的定义及基本术语
−1=±i
数学中用 i,工程中常用 j,称作虚数 imaginary number
算术基本定理:n 阶的多项式有 n 个复数根(包括重根)
复数可表示为:
z=x+yi
其中 x,y 都是实数 real number
- x 称为实部 real part
- y 称为虚部 imaginary part,虚部不包括 i
- 复数的集合用 C 表示
复数的加法、减法、乘法、除法
共轭复数 complex conjugation:
x+iy=x−iy
有用的性质:若 z=x+iy,则
zz=(x+iy)(x−iy)=x2+y2
复数 z=x+iy 的大小 magnitude / 绝对值 absolute value / 模 modulus / 范式 norm:
∣z∣=x2+y2
复平面:x 轴为实轴,y 轴为虚轴
三角形不等式:
∣z1∣+∣z2∣≥∣z1+z2∣
当且仅当其中有一个为 0 或 arg(z1)=arg(z2) 时,取 =
极坐标与欧拉公式
极坐标 polar coordinates:
r=∣z∣,θ=arg(z)
欧拉公式 Euler's Formula:
eiθ=cosθ+isinθ
证明:
复数指数的极坐标:
z=x+iy=reiθ
优点:
- arg(z)=θ
- 乘法:z1z2=r1r2ei(θ1+θ2)
- 除法:z1/z2=r2r1ei(θ1−θ2)
![乘2i旋转90度,扩大2倍](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
自乘:(1+i)6=(2eiπ/4)6=−8i
复化或复变换:例如在计算以下积分时
∫excos(2x)dx=Re(∫exei2xdx)
第 N 个根:zN=c
r=R1/N,θ=Nϕ+N2πn
![n次根](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
逆欧拉函数:
cos(t)=2eit+e−it
sin(t)=2ieit−e−it
用欧拉公式很容易得到的 de Moivre's formula:
(cosθ+isinθ)n=cosnθ+isinnθ
复数z=a+bi的矩阵表示:
Z=[ab−ba]
极坐标形式(分解成伸展因子和旋转矩阵):
Z=[r00r][cosθsinθ−sinθcosθ]
指数函数与映射
复数的指数函数:
ez=ex(cosy+isiny)
eit 将实轴映射到单位圆上:
![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
复函数是二维对二维的,所以通常将其视为映射 mapping
![复数映射1](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![复数映射2](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
穿孔平面 punctured plane 是复平面减去原点,表示为 C−{0}
辐角
辐角主值 principal branch of arg(z) :−π<arg(z)≤π,记号为 Arg(z)
若想要幅角连续,则把定义域限制为平面减去割线 branch cut
对数函数
通过复指数定义 logz:elogz=z
也就是:
logz=log∣z∣+iarg(z)
注意:
- 因为 arg(z) 有多个可能的值,所以 logz 也有多个可能值
- log0 无定义
- 选择一组辐角,使得 logz 单值,则称我们选了 log 函数的一个分支 a branch of the log function
- log 的主分支 principal branch of log 来自辐角主值,即 −π<arg(z)≤π
![对数函数](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![对数函数2](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
幂函数
定义:
za=ealog(z)
是多值的,通常需要选择 log(z) 的一个分支
解析函数 Analytic functions
极限
在微积分中我们用极限定义导数,在复数中也一样
z0 周围半径为 r 的开圆盘 open disk是点 z 的集合,∣z−z0∣<r
z0 周围半径为 r 的**开去心圆盘 open deleted disk / 有孔圆盘 punctured disk **是点 z 的集合,0<∣z−z0∣<r
![开圆盘](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
开区域 open region 其中每个点都被开圆盘包围
极限的定义:如果 f(z) 定义在 z0 周围的一个有孔圆盘上,则
z→z0limf(z)=w0
如果有多条点序列可以靠近 z0,则所有的序列都必须趋向 w0
连续函数
复变量函数的一种写法:
f(z)=f(x+iy)=u(x,y)+iv(x,y)
当且仅当 u(x,y) 和 v(x,y) 都连续时,f(z) 连续
一些连续函数:
- 多项式函数在整个平面上连续
- 指数函数在整个平面上连续
- 辐角主值函数在平面减去非正实轴连续
- 对数函数在平面减去非正实轴连续
连续函数的性质:
- 加法连续
- 乘法连续
- 除法连续(除了分母为 0)
- 复合函数连续
无穷远
定义:延伸的复平面 extended complex plane = C∪{∞}
如果在 0 周围画一个大圆,则称在这个圆之外的区域为无穷远的领域 a neighborhood of infinity
![无穷远的领域](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
无穷远与极限:
z→∞limz1=0
黎曼球面的立体投影
![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
导数
f′(z0)=z→z0limz−z0f(z)−f(z0)
如果该极限存在,则称 f 在 z0 处可解析 analytic 或可导 differentiable
如果 f 在开区域 A 中所有点都可解析,则称 f 在 A 中可解析 analytic
求导法则:同实变量函数
定理:如果 f(z) 在开圆盘上有定义且可解析,f′(z)=0,则 f(z) 是一个常数
柯西-黎曼方程
柯西-黎曼方程 Cauchy-Riemann Equations:
如果 f(z)=u(x,y)+iv(x,y) 可解析,则
f′(z)=∂x∂u+i∂x∂v=∂y∂v−i∂y∂u
其中
∂x∂u=∂y∂v,∂y∂u=−∂x∂v
定理:如果 u 和 v 满足柯西-黎曼方程且有连续的偏导数,则 f(z) 在 A 可微
f′(z) 的矩阵表示:
[uxvxuyvy]
即 f(x,y) 的雅可比矩阵
定理:假设 u 和 v 的二阶偏导存在且连续,如果 f(z)=u+iv 可解析,则 f′(z) 也可解析
定义:如果一个函数在复平面上的每一点都可解析,则称为整函数 entire function
双曲正弦和余弦:
cosh(z)=2ez+e−z,sinh(z)=2ez−e−z
性质:
- dzdcosh(z)=sinh(z),dzdsinh(z)=cosh(z)
- cosh2(z)−sinh2(z)=1
- 对于实数 x,两者都是实数
- cosh(iz)=cos(z),sinh(z)=−isin(iz)
线积分和柯西定理 Line integrals and Cauchy's theorem
线积分
线积分也被称为路径 path 或轮廓 contour 积分,表示为:
∫γf(z)dz:=∫abf(γ(t))(γ′(t))dt
另一种写法:
∫γf(z)dz=∫γ(u+iv)(dx+idy)
复数线积分基本定理:如果 f(z) 是一个在开区域 A 上可解析的复函数,γ 是 A 中一条从 z0 到 z1 的曲线,则
∫γf′(z)dz=f(z1)−f(z0)
如果任何两条端点相同的路径的积分值相等,则称该线积分路径无关 path independence
定理:如果 f(z) 在开区域 A 上有不定积分,则路径积分 ∫γf(z)dz 在 A 的所有路径路径无关
两者等价:
- ∫γf(z)dz 路径无关
- ∫γf(z)dz 沿任何封闭路径为 0
柯西定理
类似于无旋场的格林定理
如果 A 是简单连通区域,f(z) 在 A 上可解析且 C 是 A 上的简单闭合曲线,则满足以下结论:
- ∫Cf(z)dz=0
- f 的路径积分路径无关
- f 有不定积分
扩展柯西定理:
![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
∫C1−C2f(z)dz=0
柯西积分公式 Cauchy's integral formula:假设 C 是一个简单闭合曲线且函数 f(z) 在包括 C 和它内部的区域可解析,C 是逆时针方向的,则对于 C 中任何 z0
![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
f(z0)=2πi1∫Cz−z0f(z)dz
导数的柯西积分公式:
f(n)(z)=2πin!∫Cw−zf(w)n+1dw,n=0,1,2,…
积分的三角不等式
三角不等式:
- ∣z1+z2∣≤∣z1∣+∣z2∣
- ∣z1∣−∣z2∣≤∣z1−z2∣
积分的三角不等式 triangle inequality for integrals:
∣∣∣∣∣∫abg(t)dt∣∣∣∣∣≤∫ab∣g(t)∣dt
三角不等式 2:
∣∣∣∣∣∫γf(z)dz∣∣∣∣∣≤∫γ∣f(z)∣∣dz∣
推论:如果 ∣f(z)<M,则
∣∣∣∣∣∫Cf(z)dz∣∣∣∣∣≤M⋅(C 的长度)
柯西定理的第二种扩展:如果 A 是简单连通区域,包括点 z0,假设函数 g 满足:
- 在 A−{z} 上可解析
- 在 A 上连续
则对于 A 中所有闭合曲线 C
∫Cg(z)dz=0
定理:如果 f(z) 在区域 A 上可解析,则 f 有任意阶导数
注意:观察导数柯西公式的证明,这个定理的因为 w−z1 在不包括 z 的区域有任意阶导数
柯西不等式 Cauchy's inequality:CR 是圆 ∣z−z0∣=R,假设 f(z) 在 CR 上和其内部可解析,MR=CR 上最大的 ∣f(z)∣,则
∣f(n)(z0)∣≤Rnn!MR
刘维尔定理 Liouville's theorem:假设 f(z) 是整函数且被约束在一个复平面内,即对于所有 z∈C,有 ∣f(z)∣<M,则 f(z) 是一个常数
推论:算术基本定理
中值性质 Mean value property:如果 f(z) 在开圆盘上可解析,则
f(z0)=2π1∫02πf(z0+reiθ)dθ
最大模原理 Maximum modulus principle:如果 f 在连通区域 A 上可解析,z0 是其中的一个点
- 如果 ∣f∣ 在 z0 上有极大值,则 f(z) 在 z0 的邻域是一个常数
- 如果 A 被围起来且 f 在 A 及其内部连续,则要么 f 是一个常数,要么 ∣f∣ 的最大值只出现在边界上
介绍调和函数 Introduction to harmonic functions
定义:函数 u(x,y) 满足以下条件:
∇2u=uxx+uyy=0
也称拉普拉斯方程 Laplace's equation
定理:如果 f(z)=u(x,y)+iv(x,y) 在区域 A 可解析,则 u 和 v 都是调和函数
定理:如果 u(x,y) 在简单闭区域调和,则 u 是解析函数 f(z)=u(x,y)+iv(x,y) 的实部
重要推论:u 无限可导
如果 u 和 v 是解析函数的实部和虚部,则它们调和共轭 harmonic conjugate
调和的第二种证明:柯西积分公式
与解析函数类似:
中值性质 Mean value property:如果 f(z) 在开圆盘上可解析,则
u(x0,y0)=2π1∫02πf(z0+reiθ)dθ
最大值原理 Maximum principle:如果 u(x,y) 在连通区域 A 上可解析,z0 是其中的一个点
- 如果 u 在 z0 上有极大值或极小值,则 u 在 z0 的邻域是一个常数
- 如果 A 被围起来且 u 在 A 及其内部连续,则 u 的最大值和最小值只出现在边界上
注意:解析函数使用模长,所以不谈最小值
引理:解析函数 f(z)=u(x,y)+iv(x,y) ,则其梯度的点积为零
∇u⋅∇v=0
定理:f(z)=u(x,y)+iv(x,y) 可解析,f′(z)=0,则 u 的等高线与 v 的正交
一些图形:
![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
二维流体力学和复势 Two dimensional hydrodynamics and complex potentials
流体
速度场 velocity field:F(x,y,z)=(u(x,y,t),v(x,y,t))
物理假设 |
数学结果 |
静止 stationary |
只和 x,y 有关,和 t 无关:F(x,y)=(u(x,y),v(x,y)) |
不可压缩 incompressible |
无散度:divF=ux+vy=0 |
无旋 irrotational |
无旋度:curlF=vx−uy=0 |
和柯西-黎曼方程相似
复势
ϕ(z) 是区域 A 上的解析函数,z=x+iy
ϕ(z)=ϕ(x,y)+iψ(x,y)
由此可以定义一个向量场
F=∇ϕ=(ϕx,ϕy)=:(u,v)
即此处 u 和 v 由 ϕx 和 ϕy 定义
根据我们对解析函数和调和函数的研究,得到这些函数有以下性质:
- ϕ 和 ψ 都是调和的
- ϕ 和 ψ 的等高线正交
- ϕ′=ϕx−iϕy
- F 是不可压缩、无旋的
一些术语:
- ϕ 被称为 F 的电势函数 potential function
- ϕ 称为 F 的复势函数 complex potential function$
- ψ:F 的流函数 stream function
- ϕ′:复速度 complex velocity
定理:假设 F=(u,v) 是一个在简单连通区域 A 上的不可压缩、无旋的场,则有一个解析函数 ϕ 是 F 的一个复势函数
定理:F 在每个地方都和 ψ 的等高线相切,ψ 的等高线叫做流线 streamlines
驻点 stagnation point 是速度场为 0 的点,即 ϕ′(z)=0
一些函数和图形:
线性涡旋:
ϕ(z)=ilog(z)
F=(r2y,−r2x)
![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
双源点:
ϕ(z)=log(z−1)+log(z+1)
F=(0,y2+12y)
![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
在均匀流中的源点:
ϕ(z)=z+2πQlog(z)
![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
源点 + 汇点
ϕ(z)=log(z−2)−log(z+2)
![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
泰勒和洛朗级数 Taylor and Laurent series
几何级数与幂级数
几何级数 geometric series:
S=aj=0∑∞rj
r 称为几何级数比率,若 ∣r∣<1,则几何级数收敛于
S=1−ra
和柯西积分公式的关系:柯西积分公式中有 w−z1,可以写成几何级数的形式:
如果 ∣z∣<∣w∣,则有
w−z1=w1⋅1−z/w1=w1(1+(wz)+(wz)2+⋯)
如果 ∣z∣>∣w∣,则有
w−z1=−z1⋅1−w/z1=−z1(1+(zw)+(zw)2+⋯)
幂级数 power series:
f(z)=n=0∑∞an(z−z0)n
存在 R≥0 使得:
-
若 R>0,对于 ∣z−z0∣<R,级数绝对收敛于一个解析函数
-
对于 ∣z−z0∣>R,级数发散
-
R 被称为收敛半径 radius of convergence,圆盘 ∣z−z0∣<R 称为收敛圆盘 disk of convergence
-
可以逐项求导:
f′(z)=n=0∑∞nan(z−z0)n−1
导数的收敛半径也是 R
-
如果 γ 是在收敛圆盘内的曲线,则可以逐项求积分:
∫γf(z)dz=n=0∑∞∫γan(z−z0)n
注意:
- 定理没有说 ∣z−z0∣=R 时的情况
- 如果 R=∞,则函数 f(z) 是整函数
- 如果 R=0,则级数仅在 z=z0 处收敛,不代表一个解析函数
比值法 ratio test
若 L=limn→∞∣cn+1/cn∣ 存在,则:
- L<1,则级数绝对收敛
- L>1,发散
- L=1,无法确定
根值法 root test
若 L=limn→∞∣cn∣1/n 存在,则:
- L<1,则级数绝对收敛
- L>1,发散
- L=1,无法确定
泰勒级数
泰勒定理:如果 f(z) 是一个在区域 A 上的解析函数,z0∈A,则
f(z)=n=0∑∞an(z−z0)n
级数收敛于任何在 A 中 ∣z−z0∣<r 的圆盘,系数为
an=n!f(n)(z0)=2πi1∫γ(z−z0)n+1f(z)dz
其中 γ 是 A 中 z0 周围的简单闭合曲线
这种级数称为在 z0 周围表示 f 的幂级数
定理:假设 f(z) 在圆盘 ∣z−z0∣<r 上可解析且 f 不为 0,则存在一个整数 k≥0,使得 ak=0 且 f 有在 z0 周围的泰勒级数:
f(z)=(z−z0)kn=k∑∞an(z−z0)n−k
其中的 k 称为 f 在 z0 处零的阶数,如果 k=1,则 z0 是一个简单零 simple zero
定理:零是孤立的
如果一个函数 f(z) 在 z0 处不可解析,则这个函数在该点是奇异的 singular
对于一个函数 f(z),如果 f 在可去圆盘 0<∣z−z0∣<r 上是可解析的,则奇点 z0 是孤立奇点 isolated singularity
洛朗级数
假设 f(z) 在圆环上可解析
r1<∣z−z0∣<r2
则 f(z) 可以表示为一个级数:
f(z)=n=1∑∞(z−z0)nbn+n=0∑∞an(z−z0)n
系数的公式:
an=2πi1∫γ(w−z0)n+1f(w)dw
bn=2πi1∫γf(w)(w−z0)n−1dw
其中 γ 是在圆环中的任何圆 ∣w−z0∣=r
- 级数 ∑n=0∞an(z−z0)n 对于 ∣z−z0∣<r2 收敛于一个解析函数
- 级数 ∑n=1∞(z−z0)nbn 对于 ∣z−z0∣>r1 收敛于一个解析函数
- 两者共同收敛于一个圆环
一些术语:
- 整个级数叫做 f 在 z0 周围的洛朗级数 Larent series
- 级数 ∑n=0∞an(z−z0)n 是洛朗级数的可解析或常规部分 analytic or regular part
- 级数 ∑n=1∞(z−z0)nbn 是洛朗级数的奇异或主要部分 singular or principal part
如果只有有限的系数 bn 非零,则称 z0 是 f 的有限极点 finite pole。如果 bi=0 且对于所有 n>k,bn=0,则 z0 是 k 阶极点
- 如果 z0 是一阶极点,则称为 f 的一个简单极点 simple pole
- 如果无限多个 bn 非零,则 z0 是本性奇点 essential singularity 或无限阶极点
- 如果所有的 bn 都为 0,则 z0 是一个可去奇点 removable singularity
留数定理 Residue Theorem
极点和零点
如果一个函数在区域 A 上可解析,则称其在 A 上全纯 holomorphic
如果一个函数在区域 A 上除了有限阶极点之外可解析,则其在 A 上亚纯 meromorphic
函数在零点处的表现:如果 f 在 z0 处有 n 阶零点,则在 z0 附近
f(z)≈an(z−z0)n
函数在极点处的表现:如果 f 在 z0 处有 n 阶极点,则在 z0 附近
f(z)≈(z−z0)nbn
皮卡大定理 Picard's theorem:如果 f(z) 在 z0 处有本性奇点,则在每个 z0 的领域内,f(z) 取所有可能的值,只有一个例外
函数的商:假设 f 在 z0 处有 m 阶零点,g 在 z0 处有 n$ 阶极点,$$h(z) = \frac{f(z)}{g(z)}$,则
- n>m,则 h(z) 在 z0 处有 n−m 阶极点
- n<m,则 h(z) 在 z0 处有 m−n 阶零点
- n=m,则 h(z) 在 z0 处可解析且没有零点
即 h(z) 在 z0 处有 n−m 阶极点,若 n−m 是负数,则极点是一个零点
f(z) 在 z0 处有一个孤立奇点,洛朗级数为
f(z)=n=1∑∞(z−z0)nbn+n=0∑∞an(z−z0)n
留数
f 在 z0 处的留数 residue 是 b1,表示为
Res(f,z0)=b1或Resz=z0f=b1
留数的意义:将洛朗级数逐项求积分时,唯一的非零积分来自 zb1
简单极点的留数性质:
-
f(z) 的洛朗级数有如下形式:
z−z0b1+a0+a1(z−z0)+⋯
则 f 在 z0 处有一个简单极点且 Res(f,z0)=b1
-
如果
g(z)=(z−z0)f(z)
在 z0 处可解析,则 Res(f,z0)=g(z0)
-
如果 f 在 z0 处有一个简单极点,则
z→z0lim(z−z0)f(z)=Res(f,z0)
-
如果 f 在 z0 处有一个简单极点,g(z) 在 z0 处可解析,则
Res(fg,z0)=g(z0)Res(f,z0)
如果 g(z0)=0,则
Res(f/g,z0)=g(z0)1Res(f,z0)
-
如果 g(z) 在 z0 处有一个简单零点,则
Res(g1,z0)=g′(z0)1
高阶极点:如果 f(z) 在 z0 处有 k 阶极点,则
g(z)=(z−z0)kf(z)
在 z0 处可解析,如果
g(z)=a0+a1(z−z0)+⋯
则
Res(f,z0)=ak−1=(k−1)!g(k−1)(z0)
柯西留数定理
柯西留数定理 Cauchy's residue theorem:假设 f(z) 在区域 A 上除了一些孤立奇点外可解析,C 是 A 中的简单闭合曲线,不经过任何 f 的奇点,方向为逆时针,则
∫Cf(z)dz=2πi∑C 中 f 的留数
无穷处的留数:
Res(f,∞)=−2πi1∫Cf(z)dz
定理:如果 f 在 C 中除了有限的奇点之外可解析,则
Res(f,∞)=−Res(w21f(w1),0)
使用留数定理的定积分 Definite integrals using the residue theorem
函数积分
定理:假设 f(z) 定义在上半平面,如果有 a>1 且 M>0 使得
∣f(z)∣<∣z∣aM
对于较大的 ∣z∣,则
R→∞lim∫CRf(z)dz=0
其中 CR 是如下所示的半圆:
![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
类似的,若 f(z) 定义在上半平面,则积分路径如下:
![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
定理 假设 f(z) 定义在上半平面,存在 M>0 使得
∣f(z)∣<∣z∣M
则
x1→∞,x2→∞lim∫C1+C2+C3f(z)eiazdz=0
其中 C1+C2+C3 如图所示:
![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
下半平面的积分路径为:
![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
三角积分:
- e−iθ=z1
- cosθ=2z+z1
- sinθ=2iz−z1
柯西主值 Cauchy principal value:f(x) 在实轴上的除了某个点 x1 外都连续
p.v.∫−∞∞f(x)dx=R→∞,r1→0lim∫−Rx1−r1f(x)dx+∫x1+r1Rf(x)dx
也可以推广到多个间断点
定理:如果 f(x) 有间断点 x1<x2<…<xn 且 ∫−∞∞f(x) 收敛,则 p.v.∫−∞∞f(x) 也收敛
定理:假设 f(z) 在 z0 处有简单极点,Cr 是半圆 γ(θ)=z0+reiθ,0≤θ≤π,则
r→0lim∫Crf(z)dz=πiRes(f,z0)
推广为扇形:假设 f(z) 在 z0 处有简单极点,Cr 是半圆 γ(θ)=z0+reiθ,θ0≤θ≤θ0+α,则
r→0lim∫Crf(z)dz=αiRes(f,z0)
傅里叶变换
函数 f(x) 的傅里叶变换 Fourier transform 定义为:
f^(ω)=∫−∞∞f(x)e−ixωdx
傅里叶逆变换:还原函数 f(x)
f(x)=2π1∫−∞∞f^(ω)eixωdω
用傅里叶变换解微分方程:
(P(D)f)(ω)=P(iω)f^
保角映射
如果一个角 ϕ 和一个伸缩 a>0,对于任何经过 z0 的曲线 γ(t),映射 f 旋转 z0 处的切向量 ϕ 并伸缩 a,则称 f(z) 在 z0 处保角 conformal
如果 f(z) 定义在区域 A 上,如果 A 中的每个点都是保角的,则称它是 A 上的保角映射 conformal map
注意:
用复数作为切线:
γ′(t)=x′(t)+iy′(t)
切向量乘以 c=aeiϕ 相当于伸缩 a,旋转 ϕ
解析函数是保角的,相当于将切向量乘以 f′(z0)
黎曼映射定理 Riemann mapping theorem:如果 A 简单连通且不是整个平面,则存在一个从 A 到单位圆盘的双射
推论:如果对于两个区域存在双射,则称它们保角等价
分式线性变换 fractional linear transformation 是以下形式的函数:
T(z)=cz+daz+b,ad−bc=0
也被称作莫比乌斯变换 Mobius transforms 或双线性变换 bilinear transforms,简称 FLT
如果 ad−bc=0,则 T(z) 是一个常函数
延伸到无穷:
T(∞)={ca,c=0∞,c=0
T(−cd)=∞,c=0
转置:T(z)=z1,将圆内部的点映射到圆的外部,反之亦然
![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
线和圆
定理:分式线性变换把线和圆映射到线和圆
找到把 z1,z2,z3 分别映射到 w1,w2,w3 的线性变换:
T1(z)=(z−z3)(z2−z1)(z−z1)(z2−z3)
T2(w)=(w−w3)(w2−w1)(w−w1)(w2−w3)
则 T(z)=T2−1∘T1(z) 就是所求的映射
变换的矩阵表示:
[acbd]
优点:
- 乘以 r 和 FLT 一样
- T(z)=A,S(z)=B,则 T∘S(z)=AB
- T(z)=A=[acbd],则 T−1(w)=A−1=[d−c−ba]
反射和对称
如果 z2 是 z1 关于线 S 的反射,则称两者关于线 S 对称 symmetric
如果 z1 和 z2 关于直线 S 对称,则任何穿过两者的圆都和 S 正交
把定义推广到线和圆:如果每条穿过两个点的线或圆都和线或圆 S 正交,则这两个点关于 S 对称
方法是通过 FLT 把圆变换为线,在使用线的对称定义
定理:如果 S 是线或圆,z1 不在 S 上,则存在唯一的 z2 使得 z1,z2 关于 S 对称
定理:反射 z=x+iy=reiθ 在单位圆上的反射是
zˉ1=∣z∣2z=x2+y2x+iy=reiθ
定理:圆 S 的中心 c 关于圆的对称是无穷
狄利克雷问题:在边界处求解微分方程
函数为:
u(x,0)=⎩⎪⎪⎪⎪⎨⎪⎪⎪⎪⎧c0,x<x1c1,x1<x<x2c2,x2<x<x3c3,x3<x
写成这种形式:
u(x,y)=c3+(c2−c3)πθ3+(c1−c2)πθ2+(c0−c1)πθ1
![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
单位圆上的调和函数:
u(eiθ)={1,−2π<θ<2π0,2π<θ<23π
![使用的变换如上](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
米尔恩-汤姆森圆定理 Milne-Thomson circle theorem:如果 f(z) 是复势,所有的奇点都在 ∣z∣=R 外,则
ϕ(z)=f(z)+f(zˉR)
是一个流线在 ∣z∣=R 的复势,和 f 相同的奇点都在区域 ∣z∣>R
辐角原理 Argument Principle
辐角原理
辐角原理 Argument Principle:γ 是简单闭合曲线,逆时针方向;f(z) 在 γ 上和内部可解析,除了一些 γ 内部的有限极点和零点,p1,…,pm 是 γ 内部的 f 的极点,z1,…,zn 是零点,mult(zk)=zk 处零点的阶数,mult(pk)=pk 处极点的阶数
∫γf(z)f′(z)dz=2πi(∑mult(zk)−mult(pk))
Zf,γ 为 f 在 γ 内部零点的乘数的和,Pf,γ 类似,则
∫γff′dz=2πi(Zf,γ−Pf,γ)
如果 γ 是一个闭合曲线,则可以通过柯西公式定义它在 z0 处的卷绕数 turning number 或指标 index 为:
Ind(γ,z0)=2πi1∫γz−z01dz
z=γ(t) 是一条曲线,w=f(z) 是一个函数,则 w=f∘γ(t) 是另一条曲线,称 f 把 γ 映射成 f∘γ
完整的辐角原理 Argument Principle:
∫γf(z)f′(z)dz=2πi(Zf,γ−Pf,γ)=2πiInd(f∘γ,0)
推论:如果 f∘γ 没有经过 −1,则
∫γf+1f′=2πiInd(f∘γ,−1)=2πi(Z1+f,γ−Pf,γ)
鲁歇定理 Rouche's theorem:γ 是简单闭合曲线,f,h 是 γ 上和内部的解析函数,除了一些有限极点,f,h 在 γ 上内有极点,∣h∣<∣f∣ 在 γ 上恒成立,则
Ind(f∘γ,0)=Ind((f+h)∘γ,0)
即
Zf,γ−Pf,γ=Zf+h,γ−Pf+h,γ
推论,如果 f,h 可解析(没有极点),则
Zf,γ=Zf+h,γ
奈奎斯特稳定判据
系统函数 system function 是一个复变量函数,通常用 G(s) 表示
通常用极点-零点图 pole-zero diagram 可视化 G(s),例如
![](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
一个不稳定的系统 G(s) 可以通过负反馈环路 negative feedback loop 变稳定,新系统叫做闭环系统 closed loop system,其系统函数通过布莱克公式给出:
GCL(s)=1+kG(s)G(s)
k 是反馈因子 feedback factor,G(s) 叫做开环系统函数 open loop system function
定理:GCL(s) 的极点是 1+kG(s) 的零点
曲线 γ 始终在虚轴上,即
s=γ(ω)=iω,∞<ω<∞
给定系统 G(s) 和 k 尼奎斯特曲线 Nyquist plot 如下:
w=kG∘γ(ω)=kG(iω)
即尼奎斯特曲线是虚轴在映射 w=kG(s) 的图像
奈奎斯特稳定判据 Nyquist stability criterion:GCL(s) 是稳定的 ⇔Ind(kG∘γ,−1)= 右半平面 G(s) 的极点数
信号与系统
信号 signal 是时间的函数
系统 system 是某个机器或程式,其接受一个信号作为输入 input,对其做一些处理,并产生其它信号作为输出 output
线性系统 linear system 接受输入 c1f1+c2f2,产生输出 c1y1+c2y2
如果对于输入 f(t−a) 产生输出 y(t−a),则该系统时间不变 time invariant,简称 LTI 系统
拉普拉斯变换
函数 f(t) 的拉普拉斯变换 Laplace transform:
L(f;s)=∫0∞e−stf(t)dt=F(s)
注意:拉普拉斯变换只关注 t≥0 的部分,所以一般要求 f(t)=0,t<0
和傅里叶变换的关系:(事实上,拉普拉斯变换又被称作傅里叶-拉普拉斯变换)
傅里叶变换:
f^(ω)=∫−∞∞f(t)e−iωtdt
假设 f(t)=0,t<0,变成
f^(ω)=∫0∞f(t)e−iωtdt
令拉普拉斯变换中的 s=iω,则
L(f;iω)=∫0∞f(t)e−iωtdt
两者相同
如果存在 M 使得对于所有 t≥0,∣f(t)∣<Meat,则称 f(t) 有指数类型 a
定理:如果 f 有指数类型 a,则对于 Re(s)>a,L(f) 绝对收敛
性质
L(af+bg)=aL(f)+bL(g)
- 拉普拉斯变换把 t 的导数变成 s 的倍数(加上一些细节),即
L(f′(t);s)=sF(s)−f(0),对Re(s)>a有效
进一步推广:
L(f′′(t);s)=s2F(s)−sf(0)−f′(0)
- 如果 f(t) 有指数类型 a,则 F(s) 对于 Re(s)>a 是解析函数且
F′(s)=−L(tf(t);s)
也称 s-导数法则 s-derivative rule,可延伸到更多阶导数:
L(tnf(t);s)=(−1)nF(n)(s)
L(f(t−a);s)=e−asF(s)
系统函数和拉普拉斯变换
对非齐次微分方程拉普拉斯变换得到:
P(s)Y(s)=F(s)
令 G(s)=P(s)1,也就是系统函数(也称为转换函数 transfer function,因为总是把输入转换成输出),则
Y(s)=G(s)⋅F(s)
该公式用文字表述:
输出 = 系统函数 × 输入
定理:f 和 g 连续,对于所有的 Re(s)>a 的 s,F(s)=G(s),则 f(t)=g(t),t≥0
拉普拉斯变换的逆:
- 假设 f 连续且有指数类型 a,则对于 c>a,有
f(t)=2πi1∫c−i∞c+i∞F(s)estds
- 假设 F(s) 有有限个数的极点且退化类似 s1,则
f(t)=∑Res(F(s)est,pk)
对于 h(t)=f(t−a),则称 h(t) 是 f(t) 的延时版本
H(s)=e−asF(s)
有延时的闭环系统为
GCL(s)=1+ke−asGG
拉普拉斯变换表
性质和规则:
函数 |
变换 |
名称 |
f(t) |
F(S)=∫0∞e−stf(t)dt |
定义 |
af(t)+bg(t) |
aF(s)+bG(s) |
线性 |
eatf(t) |
F(s−a) |
s-位移 |
f′(t) |
sF(s)−f(0) |
|
f′′(t) |
s2F(s)−sf(0)−f′(0) |
|
tf(t) |
−F′(s) |
|
tnf(t) |
(−1)nF(n)(s) |
|
f(t−a) |
e−asF(s) |
t-位移 |
∫0tf(τ)dτ |
sF(s) |
积分法则 |
tf(t) |
∫s∞f(σ)dσ |
|
函数表:
函数 |
变换 |
收敛域 |
1 |
s1 |
Re(s)>0 |
eat |
s−a1 |
Re(s)>Re(a) |
t |
s21 |
Re(s)>0 |
tn |
sn+1n! |
Re(s)>0 |
cos(ωt) |
s2+ω2s |
Re(s)>0 |
sin(ωt) |
s2+ω2ω |
Re(s)>0 |
δ(t) |
1 |
Re(s)>0 |