复变函数
一、复变函数与解析函数
高等数学中常用的函数是”反对幂指三”,其中反函数和对数函数是后面的反函数,主要研究的是”幂指三”,用幂函数逼近函数的是泰勒展开,用三角函数逼近的是傅里叶变换,那么这里主要研究利用欧拉公式e i θ = cos θ + i s in θ ,用指数函数逼近一个函数。
欧拉公式的证明可以用泰勒展开:
e i x = 1 + ( i x ) + 2 ! ( i x ) 2 + 3 ! ( i x ) 3 + 4 ! ( i x ) 4 + 5 ! ( i x ) 5 + ...
= ( 1 − 2 ! x 2 + 4 ! x 4 + ... ) + i ( x − 3 ! x 3 + 5 ! x 5 )
= cos x + i sin x
1. 复数的概念、复数的四则运算
概念
设i 2 = 1 ,称i = − 1 为虚数单位,并称形如x + i y 或x + y i 的表达式为复数,其中x 和y 是任意两个实数,把这里的x 和y 分别称为复数z = x + i y (或z = x + y i )的实部和虚部,并记做x = Re z ,y = Im z
当复数的虚部为零、实部不为零(即y = 0 , x = 0 )时,复数z = x + i y 等于x + 0 i 为实数x ,而虚部不为零(即y = 0 )的复数称为虚数。在虚数中,实部为零(即x = 0 , y = 0 )的称为纯虚数.
共轭:若z = a + ib ,则z = a − ib
复数的模=Re z 2 + Im z 2 = z z
四则运算
1.z 1 + z 2 = ( x 1 + x 2 ) + i ( y 1 + y 2 )
2.z 1 − z 2 = ( x 1 − x 2 ) + i ( y 1 − y 2 )
3.z 1 z 2 = ( x 1 + i y 1 ) ( x 2 + i y 2 ) = ( x 1 x 2 − y 1 y 2 ) + i ( x 1 y 2 + x 2 y 1 )
4.z 2 z 1 = x 1 + i y 1 x 2 + i y 2 = x 2 2 + i y 2 2 x 1 x 2 + y 1 y 2 + i x 2 2 + i y 2 2 x 2 y 1 − x 1 y 2
运算律
交换律z 1 + z 2 = z 2 + z 1 ;
结合律z 1 + ( z 2 + z 3 ) = ( z 1 + z 2 ) + z 3 ,z 1 ( z 2 z 3 ) = ( z 1 z 2 ) z 3 ;
分配律z 1 ( z 2 + z 3 ) = z 1 z 2 + z 1 z 3 ;
z 1 ± z 2 = z 1 ± z 2 ;
z 1 z 2 = z 1 z 2 ;
( z 1 z 2 ) = z 1 z 2 ;
z = z ;
z + z = 2 x = 2 Re z ,z − z = 2 i Im z ;
z 1 z 2 = x 2 + y 2 = ( Re z ) 2 + ( Im z ) 2 ;
2. 复平面
复平面
用∣ OP ∣ 表示复数z 时,这个向量在x 轴和y 轴上的投影分别是x和y。把向量OP 的长度r 称为复数z 的模或者称为z 的绝对值,并记为∣ z ∣ ,显然∣ z ∣= r = x 2 + y 2 ,∣ z ∣≤∣ x ∣ + ∣ y ∣ ,∣ x ∣≤∣ z ∣ ,∣ y ∣≤∣ z ∣ .
辐角
如果点P 不是原点,那么把x 轴的正向与向量OP 的夹角θ 称为复数z 的辐角,记作Arg z ,对每个z = 0 ,都有无穷多个辐角,因为θ o 表示复数z 的一个辐角时,θ = θ 0 + 2 kπ ( k = 0 , ± 1 , ± 2 , ... ) 就是z 的辐角的一般表达式。
主辐角 arg z ∈ ( − π , π ]
主辐角
Arg z = arg z + 2 kπ ( k = 0 , ± 1 , ± 2 , ... )
通过tan ( Arg z ) = x y 算出角度后,还需判断所在象限,通过加减周期π 变换到( − π , π ] 之间,规律是Ⅰ、Ⅳ不变,Ⅱ、Ⅲ加减π
指数式表示和三角表示式
指数形式:z = r ( cos θ + i sin θ )
三角形式:z = r e i θ (由Euler公式:e i θ = cos θ + i sin θ 得,其中r =∣ z ∣ ,θ = Arg z )
Arg z ˉ = − Arg z (z = 0 );z ˉ = r e − i θ .
3. 乘幂与方根
乘幂
三角形式
z 1 z 2 ⋯ z n = r 1 r 2 ⋯ r n [ cos ( θ 1 + θ 2 + ... + θ n ) + i sin ( θ 1 + θ 2 + ... + θ n )]
特别地,当z 1 = z 2 = ... = z n = r [ cos ( n θ ) + i sin ( n θ )] 时,z n = r n [ cos ( n θ ) + i sin ( n θ )]
于是 ∣ z 1 z 2 ⋯ z n ∣=∣ z 1 ∣∣ z 2 ∣ ⋯ ∣ z n ∣
Arg ( z 1 z 2 ⋯ z n ) = Arg z 1 + Arg z 2 + ⋯ + Arg z n
当z 1 = z 2 = ... = z n 时
∣ z ∣ n =∣ z n ∣ , Arg ( z n ) = n Arg z
指数形式
z 1 z 2 ⋯ z n = r 1 r 2 ⋯ r n e i ( θ 1 + θ 2 + ⋯ + θ n ) , z n = r n e in θ
特别地,当∣ z ∣= r = 1 时,则( cos θ + i sin θ ) n = ( cos n θ + i sin n θ ) De Moivre公式
当用z − n = z n 1 定义负整数幂时,De
Moivre公式仍成立。
除法
指数形式
即当z 1 = r 1 e i θ 1 , z 2 = r 2 e i θ 2 时,z 2 z 1 = r 2 r 1 e i ( θ 1 − θ 2 )
于是
z 2 z 1 = ∣ z 2 ∣ ∣ z 1 ∣ , Arg z 2 z 1 = Arg z 1 − Arg z 2
方根
θ 取主辐角,若用指数表示式,则当z = r ⋅ e i θ 时,
z n 1 = r n 1 e n i ( θ + 2 kπ ) ( k = 0 , ± 1 , ± 2 , ⋯ , ± ( n − 1 ))
n z = z n 1 的n个根就是以原点为中心的、半径为r n 1 的圆的内接正多边形的n个顶点所表示的复数。
例:z 3 = − 27 ,求z
解:z 有三个值,如果是实数,z = − 3 ,此时可以写成z = 3 e πi
那么另外两个根是z = 3 e 3 i ( π + 2 kπ ) ,k = 0 , 1 , − 2
∴ z 1 = 3 , z 2 = 3 ( 2 1 + 2 3 i ) , z 3 = 3 ( − 2 1 − 2 3 i )
4. 图形的复数表示
x = 2 z + z ,y = 2 i z − z
比如直线y = 2 x + 3 ,代入得2 i z − z = z + z + 3
再比如∣ z − i + i ∣ = 2 可以看作是∣ z − ( 1 − i ) ∣ = 2 ,而∣ z 1 − z 2 ∣ 表示两点间的距离,1 − i 是复平面上的点( 1 , − 1 ) ,离一定点距离不变的曲线是个圆。这里也可以展开算,变成( x − 1 ) 2 + ( y + 1 ) 2 = 2 ,也是一个圆的表达式。
5. Jordan曲线、连通性
单连通区域 和多连通区域
正方向:沿着曲线走,区域在左手边。外圈一般是逆时针,内圈一般是顺时针。
6. 初等复变函数
指数函数e x + i y = e x ( cos y + i sin y )
对数函数ln z = ln ∣ z ∣ + i arg z ,Ln z = ln z + i 2 kπ ( k = 0 , ± 1 , ± 2 , ⋯ )
幂函数z a = e a Ln z ( z = 0 )
三角函数cos z = 2 e i z + e − i z ,sin z = 2 i e i z − e − i z
反三角函数
Arcsin z = − i Ln ( i z + 1 − z 2 )
Arccos z = − i Ln ( z + z 2 − 1 )
Arctan z = − 2 i Ln 1 − i z 1 + i z
双曲与反双曲函数
sh z = 2 e z − e − z ,ch z = 2 e z + e − z ,th z = ch z sh z
sh ( − z ) = − sh z ,ch ( − z ) = ch z
ch 2 z − sh 2 z = 1
sh ( z 1 + z 2 ) = sh z 1 ch z 2 + ch z 1 sh z 2
ch ( z 1 + z 2 ) = ch z 1 ch z 2 + sh z 1 sh z 2
sh ( i z ) = i sin z ,ch ( i z ) = cos z
sin z = sin x ch y + i cos x sh y ,cos z = cos x ch y − i sin x sh y
反双曲函数
Arcsh z = Ln ( z + z 2 + 1 )
Arcch z = Ln ( z + z 2 − 1 )
Arcth z = 2 1 Ln 1 − z 1 + z
7. 连续与导数
连续
若f ( z ) = u ( x , y ) + i v ( x , y ) ,z 0 = x 0 + y 0 ,A = a + ib ,则
z → z 0 lim f ( z ) = A ⟺ ( x , y ) → ( x 0 , y 0 ) lim u ( x , y ) = a , ( x , y ) → ( x 0 , y 0 ) lim v ( x , y ) = b
连续性的判别:lim z → z 0 f ( z ) = f ( z 0 ) ,Δ z → 0 , Δ f ( z ) → 0
连续函数之间加减乘除复合后,函数仍然连续
与实函数求导相似,可导一定连续,连续不一定可导
导数
定义:f ′ ( z ) = lim Δ z → 0 Δ z f ( z + Δ z ) − f ( z )
导数存在的充分必要条件(C-R方程):∂ x ∂ u = ∂ y ∂ v ,∂ y ∂ u = − ∂ x ∂ v 复合函数求导:g ′ [ f ( z 0 )] = g ′ ( h 0 ) f ′ ( z 0 )
反函数的导数:f ′ ( z ) = φ ′ ( w ) 1
微分:d w = f ′ ( z ) d z
例:f ( z ) = x + 2 y i ,求f ′ ( z )
解:f ′ ( z ) = Δ z → 0 lim Δ z f ( z + Δ z ) − f ( z )
注意其中z = x + i y ,Δ z = Δ x + i Δ y
∴ f ′ ( z ) = Δ z → 0 lim Δ x + i Δ y ( x + Δ x ) + 2 i ( y + Δ y ) − ( x + 2 i y ) = Δ z → 0 lim Δ x + i Δ y Δ x + 2 i Δ y
Δ z → 0 意味着Δ x → 0 且Δ y → 0 ,这和多元函数求极限一样,要求沿不同路径趋近于0极限都一样,才代表极限存在
如果令Δ y = k Δ x ,则f ′ ( x ) = Δ z → 0 lim Δ x + k Δ y Δ x + 2 ki Δ = 1 + ki 1 + 2 ki
显然,如果k取值不同,所得的极限不同,∴ 极限不存在,即不可导.
8. 解析函数与奇点
解析
设f ( z ) 在区域D有定义
设z 0 ∈ D ,若存在z 0 的一个邻域,使得f ( z ) 在此邻域内处处可导,则称f ( z ) 在z 0 处解析,也称z 0 是f ( z ) 的解析点.
若f ( z ) 在区域D上每一点都解析,则称f ( z ) 在区域D内解析,或称f ( z ) 是区域D内的解析函数;
若G是一个区域,若闭区域D ⊂ G ,且f ( z ) 在G上解析,则称f ( z ) 在闭区域D上解析.
根据定义,f ( z ) 在区域D内解析和在区域D内可导是等价的;但在z 0 处解析和在z 0 处可导意义不同,前者指的是z 0 的某一邻域内可导,但后者只要求z 0 处可导;z 0 处解析和z 0 的某一个邻域内解析是一个意思;而闭区域上D解析,说明属于D的每一点都是f ( z ) 的解析点.
解析的判别:
实部与虚部可微(导数存在且连续)
满足C-R方程:∂ x ∂ u = ∂ y ∂ v , ∂ y ∂ u = − ∂ x ∂ v
奇点和孤立奇点
若f ( z ) 在z 0 处不解析,则称z 0 是f ( z ) 的奇点. 若z 0 是f ( z ) 的奇点,但z 0 的某邻域内除z 0 之外,再没有其他的奇点,则称z 0 是f ( z ) 的孤立奇点. 根据导数的性质,如果f ( z ) ,g ( z ) 在区域D上解析,则f ( z ) ± g ( z ) 和f ( z ) g ( z ) 也在D内解析,当g ( z 0 ) = 0 时,z 0 是g ( z ) f ( z ) 的解析点.特别地, 多项式P ( z ) 是全平面内的解析函数, 而有理分式R ( z ) = Q ( z ) P ( z ) 在复平面内除分母为零的点之外,
处处解析,分母为零的点是R ( z ) 的孤立奇点.
找奇点一般找分母为零的点
9. 可导、解析的充要条件
可导充要条件
复变函数f ( z ) = u ( x , y ) + i v ( x , y ) 在点z 0 = x 0 + i y 0 处可微(即可导)的充分必要条件是二元函数u ( x , y ) , v ( x , y ) 在( x 0 , y 0 ) 处都可微,并且满足
CauchyRiemann(柯西-黎曼)方程
∂ x ∂ u = ∂ y ∂ v , ∂ y ∂ u = − ∂ x ∂ v
此时f ′ ( z ) = ∂ x ∂ u + i ∂ x ∂ v = ∂ x ∂ u − i ∂ y ∂ u = ∂ y ∂ v + i ∂ x ∂ v = ∂ y ∂ v − i ∂ y ∂ u
解析充要条件
复变函数f ( z ) = u ( x , y ) + i v ( x , y ) 在区域D内解析的充分必要条件是u ( x , y ) , v ( x , y ) 都在区域D内可微,且在D内满足
Cauchy-Riemann方程
∂ x ∂ u = ∂ y ∂ v , ∂ y ∂ u = − ∂ x ∂ v
在区域D内
f ′ ( z ) = ∂ x ∂ u + i ∂ x ∂ v = ∂ x ∂ u − i ∂ y ∂ u = ∂ y ∂ v + i ∂ x ∂ v = ∂ y ∂ v − i ∂ y ∂ u
注:判断一个复变函数是否可微?受否解析?首先,看它的实部与虚部是否可导(可微);之后分别求实部与虚部对x和y的偏导,判断是否符合Cauchy-Riemann方程;若符合Cauchy-Riemann方程,则说明复变函数解析;若部分符合Cauchy-Riemann方程,那么将在符合方程的条件下函数可微,而在任何一点的领域内都有不可微的点,故函数处处不解析。
柯西-黎曼方程(C-R方程)的证明
由f ′ ( z ) = Δ z → 0 lim Δ z f ( z + Δ z ) − f ( z )
得出f ( z + Δ z ) − f ( z ) = f ′ ( z ) Δ z + o ( ∣Δ z ∣ )
⇒ Δ u + i Δ v = ( a + ib ) ( Δ x + i Δ y ) + o ( ∣Δ z ∣ )
⇒ Δ u + i Δ v = ( a Δ x − b Δ y ) + i ( b Δ x + a Δ y ) + o ( ∣Δ z ∣ )
令实部等于实部,虚部等于虚部,得
{ Δ u = ( a Δ x − b Δ y ) + o ( ∣Δ z ∣ ) ( 1 ) Δ v = ( b Δ x + a Δ y ) + o ( ∣Δ z ∣ ) ( 2 )
u和v可以写成增量加高阶无穷小的形式,意味着要求他们都是可微的
由(1)式得a = ∂ x ∂ u ,b = − ∂ y ∂ u
由(2)式得a = ∂ y ∂ v ,b = ∂ x ∂ v
∴ ∂ x ∂ u = ∂ y ∂ v ,∂ y ∂ u = − ∂ x ∂ v ,而f ′ ( z ) = a + ib = ∂ x ∂ u + i ∂ x ∂ v
该过程反过来证也成立,所以是充要条件
C-R方程简单地理解就是,虚部的偏导可以由实部给出,即实部决定了虚部,这样就把一个虚数的问题变成了实数的问题
f ( z ) = u + i v 解析的条件
对于f ( z ) = u ( x , y ) + v ( x , y ) i ,(其中z = x + i y )
可以写成f ( z , z ) = u ( 2 z + z , 2 i z − z ) + v ( 2 z + z , 2 i z − z ) i ,其实是关于z 和z 的函数
如果要使f ( z ) 解析,则要求化简后不含z ,即∂ z ∂ f ( z , z ) = 0
证明:
∂ z ∂ f ( x , y ) = ∂ z ∂ u ( x , y ) + ∂ z ∂ v ( x , y ) i
而x ( z , z ) = 2 z + z ,y ( z , z ) = 2 i z − z ,根据全微分,有
∴ ∂ z ∂ f ( u , v ) = ∂ x ∂ u ( x , y ) ∂ z ∂ x + ∂ y ∂ u ( x , y ) ∂ z ∂ y + ( ∂ x ∂ v ( x , y ) ∂ z ∂ x + ∂ y ∂ v ( x , y ) ∂ z ∂ y ) i = 2 1 ∂ x ∂ u ( x , y ) − 2 i 1 ∂ y ∂ u ( x , y ) + ( 2 1 ∂ x ∂ v ( x , y ) − 2 i 1 ∂ y ∂ v ( x , y ) ) i = 2 1 ( ∂ x ∂ u ( x , y ) − ∂ y ∂ v ( x , y ) ) + 2 1 ( ∂ y ∂ u ( x , y ) + ∂ y ∂ v ( x , y ) ) i = 0
令实部和虚部分别等于0,得∂ x ∂ u = ∂ y ∂ v ,∂ y ∂ u = − ∂ x ∂ v
即,对于f ( z ) = u ( x , y ) + v ( x , y ) i ,∂ z ∂ f ( z , z ) = 0 等价于柯西-黎曼方程
10. 初等解析函数
复指数函数
假设z = x + i y ,则由e x ( cos y + i sin y ) 定义复指数函数,记为或简记为e z = e x ( cos y + i sin y ) 。 显然 Re ( e z ) = e x cos y ,Im ( e z ) = e x sin y
∣ e z ∣= e x , Arg ( e z ) = y + 2 kπ ( k = 0 , ± 1 , ± 2 , ⋯ ) .
性质
e z 为指数函数,则e z 在全平面解析,( e z ) ′ = e z , 且
e z + w = e x e w 对所有z , w ∈ C 成立,所以( e z ) n = e n z ;
e z = 0 , 如果z = x 为实数,当x > 0 , e x > 1 , 当x < 0 , e x < 1 ;
e z 是周期函数,其周期T = 2 nπi (其中n 为整数);
e 2 π i = i , e πi = − 1 , e 2 3 πi = − i , e 2 πi = 1 ;
e z = 1 的充要条件是z = 2 nπi (其中n 为整数);
注:复指数函数在全平面解析
复对数函数
指数函数的反函数称为对数函数,即把满足方程e w = z ( z = 0 ) 的函数w = f ( z ) 称为z 的对数函数,记作w = ln z .
令w = u + i v , z = r e i θ ,则由e w = z ( z = 0 ) 可得e u + i v = r e i θ ,从而由复数的相等的定义知e u = r , v = θ + 2 kπ , 即u = ln r , v = θ + 2 kπ (k 为整数) ,
或u = ln ∣ z ∣ , v = Arg z
所以
w = Ln z = ln ∣ z ∣ + i Arg z = ln ∣ z ∣ + i ( arg z + 2 kπ ) ( k = 0 , ± 1 , ⋯ )
记ln z = ln ∣ z ∣ + i arg z ; 则对数函数可以写为Ln z = ln z + 2 ikπ (其中k 为整数),
对应某个确定的k , 称为对数函数的第k 个分支,对应k = 0 的那个分支,则称为对数函数主支,ln z 称为对数函数的主值。
Ln ( z 1 z 2 ) = Ln z 1 + Ln z 2 , Ln ( z 2 z 1 ) = Ln z 1 − Ln z 2 ( z 1 z 2 = 0 )
注意:Ln z 2 = ln z 2 + 2 kπi ,2 Ln z = 2 ln z + 4 kπi ,两者并不相等,差了周期
注:对数主支ln z =∣ z ∣ + i arg z 在去掉原点和负实轴的平面上解析,且d z d l n z = z 1 .
因为ln x 在0处没有定义,而对于arg z ,取到π 后又变成了− π ,在负实轴处不连续
复幂函数
设z 为不等于零的复变数,μ 为任意一个复数,定义复幂函数z μ =e μ Ln z .
当z 为正实变数,μ 为实数时,与乘幂一致;当z 为复变函数,μ 为复数时
z μ = e μ Ln z = e μ ( l n z + i 2 kπ ) = e μ l n z ⋅ e 2 kπ μ i ( k = 0 , ± 1 , ± 2 , ⋯ )
注:复幂函数在去掉原点和负实轴的平面上解析
因为Ln z 在去掉原点和负实轴的平面上解析,所以复幂函数在该区域上也为解析函数,且由复合函数求导公式可得( z μ ) ′ = μ z μ − 1 .
在实数里ln 1 = 0 ,但是在复数里Ln 1 = ln 1 + 2 kπi = 2 kπi ,仍然有虚部,这就导致了1 n = 1
例1:1 2 = e 2 Ln 1 = e 2 2 kπi = cos 2 2 kπ + i sin 2 2 kπ
例2:( 1 + 3 i ) i = e i Ln ( 1 + 3 i ) = e i ( l n 2 + 3 π i + 2 kπi ) = e − 3 π − 2 kπ + i l n 2 = e − 3 π − 2 kπ ( cos ln 2 + i sin ln 2 )
三角函数和双曲函数
复三角函数与指数函数相关联(即与双曲函数相关联)
定义三角函数与双曲函数如下:
正弦函数sin z = 2 i e i z − e − i z ;
余弦函数cos z = 2 e i z + e − i z ;
双曲正弦函数sh z = 2 e z − e − z ;
双曲余弦函数ch z = 2 e z + e − z ;
注:由于e ± z , e ± i z 在C 复平面上是解析的,所以上述四个函数在整个复平面上解析
性质
( sin z ) ′ = cos z ,( cos z ) ′ = − sin z ,( sh z ) ′ = ch z ,( ch z ) ′ = sh z ,( tan z ) ′ = sec 2 z = ( c o s z 1 ) 2
sin z 、cos z 是以2 π 为周期的周期函数;sh z 、ch z 是以2 πi 为周期的周期函数,这是因为e ± i z 是以2 π 为周期的函数;e ± z 是以2 πi 为周期的函数.
sin z 、sh z 为奇函数;cos z 、ch z 为偶函数.
一些恒等式关系仍成立.
角度关系仍成立。比如sin z = 0 ⇒ 2 i e i z − e − i z = 0 ⇒ e i z = e − i z ⇒ e 2 i z = 1 = e 2 kπi ⇒ z = kπ ,cos z = 0 ⇒ 2 e i z + e − i z = 0 ⇒ e i z = − e − i z ⇒ e 2 i z = − 1 = e ( − π + 2 kπ ) i ⇒ z = − 2 π + kπ
三角函数与双曲函数函数之间满足关系式cos ( i z ) = ch z ,sin ( i z ) = i sh z ,ch ( i z ) = cos z ,sh ( i z ) = i sin z .
sin z , cos z 不是有界函数。
因为sin ( z ) = sin ( x + i y ) = sin x cos ( i y ) + cos x sin ( i y ) = sin x ch y + i cos x sh y
所以∣ sin z ∣ 2 = sin 2 x ch 2 y + cos 2 x sh 2 y = sin 2 x ( ch 2 y − sh 2 y ) + sh 2 y = sin 2 x + sh 2 y
虽然0 ≤ sin 2 x ≤ 1 , 但是当y → ∞ 时,sh 2 y → ∞ , 所以当y → ∞ 时,∣ sin z ∣→ ∞ , 即sin z 是无界函数。
(简单证明:cos ( i z ) = ch z ,而ch z 无界)
反三角反双曲函数
Arcsin z = − i Ln ( i z + 1 − z 2 ) ,Arccos z = − i Ln ( z + z 2 − 1 ) ,Arctan z = − 2 i Ln 1 − i z 1 + i z Arcsh z = Ln ( z + z 2 + 1 ) ,Arcch z = Ln ( z + z 2 − 1 ) ,Arcth z = 2 1 Ln 1 − z 1 + z
证明:
令z = sin w ,则w = Arcsin z
z = 2 i e i w − e − i w ⇒ 2 i z = e i w − e − i w
两边同乘e i w 得,( e i w ) 2 − 1 − 2 i z e i w = 0
令t = e i w ,t 2 − 2 i z t − 1 = 0
∴ t = 2 2 i z ± − 4 z + 4 = i z ± 1 − z 2 = e i w
两边取对数得w = − i Ln ( i z ± 1 − z 2 )
其他证明同理
11.调和函数
调和函数:(1)具有二阶偏导数(2)满足拉普拉斯方程h xx ( x , y ) + h yy ( x , y ) = 0
共轭调和的概念:若u 调和,∃ v 使得u + i v 是解析函数,则称v 是u 的共轭调和函数
共轭调和的性质:(1)一阶导数满足C-R方程 (2)实部和虚部是调和函数
f = u + i v 解析,则u 、v 调和
证明:∂ 2 x ∂ 2 u + ∂ 2 y ∂ 2 u = C-R 方程 ∂ y ∂ x ∂ 2 v + ( − ∂ x ∂ y ∂ 2 v ) = 二阶导连续 0 ,v 的部分也可以这样证明出来
解析与调和的关系:解析的充分必要条件为虚部是实部的共轭调和函数,或者说,负实部是虚部的共轭调和函数
找共轭调和函数的三种方法
线积分法
d v = 全微分 ∂ x ∂ v d x + ∂ y ∂ v d y = C-R 方程 − ∂ y ∂ u d x + ∂ x ∂ u d y ,然后v = ∫ ( x 0 , y 0 ) ( x , y ) d v + C
不定积分法
f ′ ( z ) = ∂ x ∂ u + i ∂ x ∂ v = ∂ x ∂ u − i ∂ y ∂ u
f ( z ) = ∫ f ′ ( z ) d z = ∫ ( ∂ x ∂ u − i ∂ y ∂ u ) d z ,v 就是f ( z ) 的虚部
偏积分法
v = ∫ ∂ y ∂ v d y + C 1 ( x ) = ∫ ∂ x ∂ u d y + C 1 ( x ) ,这里对y 积分x 是常数,然后∂ x ∂ v = [ ∫ ∂ x ∂ u d y ] ′ + C 1 ′ ( x ) ,观察C-R方程中的∂ x ∂ v ,可得出C 1 ′ ( x ) ,再让C 1 ′ ( x ) 对x 积分即可
或者再算出v = ∫ ∂ x ∂ v d x + C 2 ( y ) = − ∫ ∂ y ∂ u d x + C 2 ( y ) ,后对比得出C 1 ( x ) 和C 2 ( y )
注意:找v 之前,要先验证u 是否调和
例:已知u = y 3 − 3 x 3 y ,求共轭调和函数v
先验证u 是调和的:∂ x ∂ u = − 6 x y = ∂ y ∂ v ,∂ y ∂ u = 3 y 3 − 3 x 2 = − ∂ x ∂ v
∂ 2 x ∂ 2 u = ∂ 2 y ∂ 2 u = − 6 y ,∴ u 是调和的
下面分别用三种不同方法求v :
线积分法:
为方便( x 0 , y 0 ) 取( 0 , 0 )
v = ∫ ( x 0 , y 0 ) ( x , y ) d v + C = 全微分 ∫ ( x 0 , y 0 ) ( x , y ) ∂ x ∂ v d x + ∂ y ∂ v d y = C-R 方程 ∫ ( x 0 , y 0 ) ( x , y ) ( 3 x 2 − 3 y 2 ) d x + ( − 6 x y ) d y = ∫ x 0 x ( 3 x 2 − 3 y 2 ) d x − ∫ y 0 y 6 x y d y + C = x 3 − 3 x y 2 + C
不定积分法
f ( z ) = ∫ f ′ ( z ) d z = ∫ ∂ x ∂ u + i ∂ x ∂ v d z = ∫ − 6 x y + i ( 3 x 2 − 3 y 2 ) d z = ∫ 3 i ( x 2 + 2 i x y − y 2 ) d z = ∫ 3 i z 2 d z = i z 3 + C
v = i f ( z ) − u = x 3 − 3 x y 2 + C
之前说过,任何复变函数都可以写成f ( z , z ) 的形式,而且如果f = u + i v 解析,化简后一定没有z ,所以总能凑出z 来
可以用这些式子:z = x + i y ,z = x − i y ,z ⋅ z = x 2 + y 2 .
比如( x 2 + y 2 ) 2 x 2 − y 2 − 2 x y i = ( z z ) 2 z 2 = z 2 1 .
偏积分法
v = ∫ ∂ y ∂ v d y + C 1 ( x ) = ∫ ( − 6 x y ) d y + C 1 ( x ) = − 3 x y 2 + C 1 ( x )
∂ x ∂ v = − 3 y 2 + C 1 ′ ( x ) ,而C-R方程中∂ x ∂ v = 3 x 2 − 3 y 2 ,∴ C 1 ′ ( x ) = 3 x 2
∴ v = x 3 − 3 x y 2 + C
二、复变函数的积分
12. 积分的概念及性质
积分的概念
复变函数积分
设C是以z 0 为起点,Z 为终点的有向连续曲线,f ( z ) 是定义在C 上的复变函数. 对于曲线C 上任意的分割点0 z 0 , z 1 , z 2 , ⋯ , z n = Z , 在每个小弧段z k − 1 z k 上任取一点ς k ( k = 1 , 2 , ⋯ , n ) ,令Δ z k = z k − z k − 1 ( k = 1 , 2 , ⋯ , n ) , 做和数S n = ∑ k = 1 n f ( ς k ) Δ z k
设λ = 1 ≤ k ≤ n ma x { ∣ Δ z k ∣ } , 如果极限λ → 0 l im ∑ k = 1 n f ( ς k ) Δ z k 存在,那么称这个极限值为f ( z ) 在曲线C 上的积分 , 并记作C ∫ f ( z ) d z
即C ∫ f ( z ) d z = λ → 0 lim ∑ k = 1 n f ( ς k ) Δ z k
如果C 是一条封闭曲线,往往记作C ∮ f ( z ) d z 而当C 是实轴上的区间[ a , b ] 方向a → b ,f ( z ) 为实值函数,则这个积分就是定积分.
此积分满足:
线性:C ∫ f ( z ) + g ( z ) d z = C ∫ f ( z ) d z + C ∫ g ( z ) d z
可加性:C ∫ = C 1 ∫ + C 2 ∫
估值:∣ C ∫ f ( z ) d z ∣ ≤ C ∫ ∣ f ( z ) ∣ d z ,若M是∣ f ( z ) ∣ 在C 处的最大值,则∣ C ∫ f ( z ) d z ∣≤ M C ∫ 1 d z
方向: 若C 表示正方向,C − 表示反方向,则C ∫ f ( z ) d z = − C − ∫ f ( z ) d z ,
估值定理的应用
例:C 为从0 到3 + 4 i 的直线段,求C ∫ z − i 1 d z 的模长的上界
∣ C ∫ z − i 1 d z ∣≤ M L ,其中M是∣ z − i 1 ∣ 在C上的最大值,L是曲线长度,由于是直线,L = 5
∣ z − i 1 ∣=∣ x + ( y − 1 ) i 1 ∣= x 2 + ( y − 1 ) 2 1 ,而y = 3 4 x ,其中x 从0到3,代入后得知∣ z − i 1 ∣≤ 3 5
∴ ∣ C ∫ z − i 1 d z ∣≤ 3 15
积分的性质
定理 1
设C 是按段光滑的有向曲线,f ( z ) 在C 上连续,则∫ C f ( z ) d z 存在 , 并且
∫ C f ( z ) d z = ∫ C u d x − v d y + i ∫ C v d x + u d y
证明:
f ( z ) = u + i v ,z = x + i y
根据全微分d z = ∂ x ∂ z d x + ∂ y ∂ z d y = 1 d x + i d y
∴ C ∫ f ( z ) d z = C ∫ ( u + i v ) ( d x + i d y ) = C ∫ u d x − v d y + i C ∫ v d x + u d y
定理 2(z为参数方程时)
设C : z = z ( t ) = x ( t ) + i y ( t ) ( α ≤ t ≤ β ) 是一条光滑曲线,z ( α ) 是C 的起点,z ( β ) 是终点,则
∫ C f ( z ) d z = ∫ α β f [ z ( t )] z ′ ( t ) d t
注:∫ C z d z 与积分路径无关 ,∫ C Re z d z 与积分路径相关
例1:重要积分C ∮ ( z − z 0 ) n + 1 1 d z
其中C 是圆周∣ z − z 0 ∣= r ( r > 0 ) 的正向
∣ z − z 0 ∣= r ⇒ ( x − x 0 ) 2 + ( y − y 0 ) 2 = r 2 ⇒ { x = x 0 + r cos θ y = y 0 + r sin θ ⇒ z = x + i y = ( x 0 + i y 0 ) + r cos θ + i r sin θ ⇒ z = z 0 + r e i θ
∴ z − z 0 = r e i θ
∴ I = ∣ z − z 0 ∣ = r ∮ ( z − z 0 ) n + 1 1 d z = ∫ 0 2 π ( r e i θ ) n + 1 1 d r e i θ = ∫ 0 2 π ( r e i θ ) n + 1 i r e i θ d θ = ∫ 0 2 π r n i e − in θ d θ = r n i ∫ 0 2 π e − in θ d θ
当n = 0 时,I = i ∫ 0 2 π 1 d θ = i θ ∣ 0 2 π = 2 πi
当n = 0 时,I = r n i − in 1 e − in θ ∣ 0 2 π = − n r n 1 [ cos ( − n 2 π ) + i sin ( − n 2 π ) − cos 0 − i sin 0 ] = 0
于是得到了重要积分 :
C ∮ ( z − z 0 ) n + 1 1 d z = { 2 πi 0 , n = 0 , n = 0
注: 积分值与圆周的半径无关,与n 的取值有关,且C 是∣ z − z 0 ∣ = r 要和分母处保持一致(其实由后面的复合闭路定理知,只要C 能包住z 0 点就行,不一定要以z 0 为圆心)
例2:I = ∣ z ∣ = 2 ∮ ∣ z ∣ z d z
由于z z = ∣ z ∣ 2 ,∴ z = z ∣ z ∣ 2 ,代入后利用例1结论,得I = ∣ z ∣ = 2 ∮ z ∣ z ∣ = ∣ z ∣ = 2 ∮ z 2 = 4 πi
例3:C ∫ z 2 d z ,其中
(1)C是从(0,0)到(2,1)的直线段 (2)C是从(0,0)到(2,0)再到(2,1)的折线段
解:(1)积分路径为y = 2 1 x ,而z = x + i y = x + 2 1 i x ,∴ d z = ( 1 + 2 i ) d x
C ∫ z 2 d z = ∫ 0 2 [( 1 + 2 i ) x ] 2 ( 1 + 2 i ) d x = ( 1 + 2 i ) 3 ∫ 0 2 x 2 d x = 3 1 ( 2 + 11 i )
(2)利用可加性,C 1 表示从(0,0)到(2,0)的直线,C 2 表示从(2,0)到(2,1)的直线
对于C 1 ,z = x + i y = x ,其中x 从0到2,∴ d z = d x ,C 1 ∫ x 2 d x = 3 1 x 3 ∣ 0 2 = 3 8
对于C 2 ,z = x + i y = 2 + i y ,其中y 从0到1,∴ d z = i d y ,C 2 ∫ ( 2 + i y ) 2 i d y = − 2 + 3 11 i
C ∫ = C 1 ∫ + C 2 ∫ = 3 8 − 2 + 3 11 i = 3 1 ( 2 + 11 i )
可见有时积分值是与积分路径无关的
13. Cauchy积分定理
又叫柯西积分定理
设f ( z ) 是单连通区域D 上的解析函数,则对D 内的任何按段光滑的J or d an 曲线C ,
都有
C ∮ f ( z ) d z = 0
利用Green公式进行证明:
Riemann通过加强条件”f(z)在区域D内连续”,使证明更简单(Goursat的证明不需要该条件):
若f ( z ) 在单连通区域D内解析且连续,则∀ D 内封闭曲线C ,C ∮ f ( z ) d z = 0
I = C ∫ f ( z ) d z = C ∫ ( u + i v ) ( d x + i d y ) = C ∫ u d x − v d y + i C ∫ v d x + u d y
根据格林公式,I = G ∬ ( − ∂ x ∂ v − ∂ y ∂ u ) d x d y + i G ∬ ( ∂ x ∂ u − ∂ y ∂ v ) d x d y
由于f ( z ) 在D 内解析,满足C-R方程,∴ I = 0 .
推论1. 若f ( z ) 在D 内解析,积分与路径无关.
推论2. 曲线C 围成的是一个单连通区域D,若f ( z ) 在在D 内解析,则C ∮ f ( z ) d z = 0 。也就是说,含有限个孤立奇点闭区域内的积分只和奇点有关,只要包含的奇点不变,积分路径可以随意变换;若内部没有奇点,则积分直接等于0,
例:∣ z ∣ = 2 1 ∮ ( z − 2 ) ( z + 1 ) z 2 d z
解:∣ z ∣ = 2 1 内没有奇点,积分为0
Cauchy积分定理反过来也成立:若f ( z ) 在单连通区域D 内连续,且∀ 简单闭曲线积分为0,则f ( z ) 在D 内解析.(Morera定理)
14. 复合闭路定理
设 f(z) 是多连通区域 D 上的解析函数,C , C 1 , C 2 , … , C n 是D 内按段光滑的Jordan曲线,且满足:
C 1 , C 2 , … , C n 中的任意曲线都在其余曲线的外部(即互不包含也不相交),
而它们都在C的内部;
位于C 的内部和C 1 , C 2 , … , C n 外部的区域G 包含于D ,其边界C , C 1 , C 2 , … , C n 也包含于D ,即G⊂ D ,则
C ∮ f ( z ) d z = k = 1 ∑ n C k ∮ f ( z ) d z
其中C 和C k 均取逆时针方向
证明
复合闭路定理
任何多连通区域都可以“剪开”成单连通区域,然后根据Cauchy积分定理,C + L 1 + C 1 − + L 1 − ∮ f ( z ) d z = 0 ,此时L 1 和L 1 − 是方向相反的同一条线,于是L 1 ∫ + L 1 − ∫ = 0 ,根据可加性C + L 1 + C 1 − + L 1 − ∮ f ( z ) d z = C ∫ + L 1 ∫ + C 1 − ∫ + L 1 − ∫ = C ∫ + C 1 − ∫ = 0 ,即C ∫ = C 1 ∫ ,其中C 1 取逆时针方向,如果有C 1 ... C n ,也可以这样逐个“剪开”,就变成了C ∮ f ( z ) d z = ∑ k = 1 n C k ∮ f ( z ) d z
复合闭路定理说明,单连通闭区域内的积分可以转化为用逆时针闭曲线圈住奇点后,对这些曲线进行积分再求和
进一步的推论,若外圈取逆时针,包含奇点的内圈全部取顺时针,外圈可以转化为在这些内圈上取逆时针再求和,于是内圈的顺时针和逆时针方向积分相互抵消,积分为0。简单地说就是,多连通区域若全体都是正方向且解析,那么积分为0。
例1:C ∮ z cosz d z ,C 由∣ z ∣ = 2 的逆时针与∣ z ∣ = 1 的顺时针围成
解:奇点只有( 0 , 0 ) ,∣ z ∣ = 2 逆时针的曲线可以转化成∣ z ∣ = 1 逆时针曲线,此时与∣ z ∣ = 1 顺时针抵消,积分为0。
例2:I = ∣ z ∣ = 2 1 ∮ z 2 − z 1 d z
解:裂项得I = ∣ z ∣ = 2 1 ∮ ( z − 1 1 − z 1 ) d z
第一项的奇点是1,不在∣ z ∣ = 2 1 的圆内,由柯西积分定理知,积分为0
第二项根据重要积分∣ z − z 0 ∣ = r ∮ z − z 0 1 d z = 2 πi ,可知积分为2 πi
∴ I = − 2 πi
例3:I = C ∮ z 2 − z 1 d z ,其中C 是把∣ z ∣ = 1 围住的简单正向闭曲线
解:裂项得I = C ∮ ( z − 1 1 − z 1 ) d z ,把C 分成C 1 、C 2 ,C 1 是只把奇点( 0 , 0 ) 围住的一个小圆,C 2 是只把奇点( 0 , 1 ) 围住的一个小圆
于是I = C 1 ∮ ( z − 1 1 − z 1 ) d z + C 2 ∮ ( z − 1 1 − z 1 ) d z = ( 0 − 2 πi ) + ( 2 πi − 0 ) = 0
15. Cauchy积分公式
设f ( z ) 是在单连通区域 D
上的解析函数,z 0 是 D 内的一个定点,C
是任意一条把z 0 含在其内部区域的按段光滑的 Jordan
曲线,则
f ( z 0 ) = 2 πi 1 C ∮ z − z 0 f ( z ) d z
证明:C ∮ z − z 0 f ( z ) d z = C ∮ z − z 0 f ( z ) − f ( z 0 ) + f ( z 0 ) d z = ∣ z − z 0 ∣ = r ∮ z − z 0 f ( z ) − f ( z 0 ) d z + f ( z 0 ) ⋅ ∣ z − z 0 ∣ = r ∮ z − z 0 1 d z
根据估值定理,∣ ∣ z − z 0 ∣ = r ∮ z − z 0 f ( z ) − f ( z 0 ) d z ∣≤ ∣ z − z 0 ∣ = r ∮ ∣ z − z 0 f ( z ) − f ( z 0 ) ∣ d z
令∣ f ( z ) − f ( z 0 ) ∣ = ε ,当r 足够小,ε → 0 ,而∣ z − z 0 ∣ = r ,代入得
∣ z − z 0 ∣ = r ∮ ∣ z − z 0 f ( z ) − f ( z 0 ) ∣ d z = r ε ⋅ ∣ z − z 0 ∣ = r ∮ 1 d z = r ε ⋅ 2 π r = 2 π ε → 0
∴ C ∮ z − z 0 f ( z ) d z = 0 + f ( z 0 ) ⋅ ∣ z − z 0 ∣ = r ∮ z − z 0 1 d z = 2 πi f ( z 0 )
不严谨的推导就是,换路径后,当r 足够小,f ( z ) → f ( z 0 ) ,可以当作常数提出来,剩下的用重要积分算得2 πi
如果换元把∣ z − z 0 ∣ = r 变成z = z 0 + r e i θ ,有f ( z 0 ) = 2 πi 1 ∫ 0 2 π r e i θ f ( z 0 + r e i θ ) r e i θ i d θ = 2 π 1 ∫ 0 2 π f ( z 0 + r e i θ ) d θ ,说明解析函数在圆心处的取值等于圆周的平均值
如果让z 0 动起来,有f ( z ) = 2 πi 1 C ∮ ς − z f ( ς ) d ς ,说明如果边界确定了,内部每一点的取值都可以由边界给出,其中f ( ς ) 是边界上的任意一点的值,ς − z 1 叫做核函数
例1:∣ z ∣ = 2 ∮ z − 2 π i e i z d z ,奇点是( 0 , 2 π ) ,在∣ z ∣ = 2 里面,∴ 2 πi e i ( 0 + 2 π i ) = 2 πi e − 2 π .
例2:前面的题目也可以用柯西积分公式来做,C ∮ z 2 − z 1 d z = C 1 ∮ z z − 1 1 d z + C 2 ∮ z − 1 z 1 d z = 2 πi 0 − 1 1 + 2 πi 1 1 = 0 ,这样就不用再裂项了
16. Cauchy导数公式
设f ( z ) 是在单连通区域D 上的解析函数,C 是D 内按段光滑的Jordan曲线,z 0 是C 内部的任意点,则f ( z ) 在z 0 处存在各阶导数,且
f ( n ) ( z 0 ) = 2 πi n ! C ∮ ( z − z 0 ) n + 1 f ( z ) d z
此处C 取正向; 如果用ς 表示积分变量,即表示C 上的点,z 表示C 内部的任意点,则Cauchy导数公式可写为 f ( n ) ( z ) = 2 πi n ! C ∮ ( ς − z ) n + 1 f ( ς ) d ς
或者
C ∮ ( z − z 0 ) n + 1 f ( z ) d z = n ! 2 πi f ( n ) ( z 0 )
例1:I = ∣ z ∣ = 2 ∮ z 3 ( z + 1 ) 1 d z ,用复合闭路定理分成C 1 圈住( 0 , 0 ) ,C 2 圈住( − 1 , 0 )
I = C 1 ∮ z 3 z + 1 1 d z + C 2 ∮ z + 1 z 3 1 d z = 2 ! 2 πi ( z + 1 1 ) ′′ ∣ z = 0 + 2 πi z 3 1 ∣ z = − 1 = 2 ! 2 πi ⋅ ( z + 1 ) 3 2 ∣ z = 0 − 2 πi = 0 .
例2:I = ∣ z ∣ = 2 ∮ ( z 2 + 1 ) 2 e z d z ,奇点为± i ,∴ I = ∣ z ∣ = 2 ∮ ( z + i ) 2 ( z − i ) 2 e z d z ,用C 1 圈住( 0 , 1 ) ,C 2 圈住( 0 , − 1 ) ,得
I = C 1 ∮ ( z − i ) 2 ( z + i ) 2 e z d z + C 2 ∮ ( z + i ) 2 ( z − i ) 2 e z d z = 2 πi [ ( z + i ) 2 e z ] ′ ∣ z = i + 2 πi [ ( z − i ) 2 e z ] ′ ∣ z = − i = πi ( sin 1 − cos 1 )
17. 解析函数的原函数
设f ( z ) 是定义在区域D 上的复变函数,若存在D 上的解析函数F ( z ) ,使F ′ ( z ) = f ( z ) 在D 上成立,则称F ( z ) 是f ( z ) 在区域D 上的原函数.
f ( z ) 在区域D 上存在原函数F ( z ) ,则f ( z ) 必须是解析函数因为解析函数的导函数仍是解析函数 (所以解析函数有任意阶导数).
定理一
设F ( z ) 和G ( z ) 都是f ( z ) 在D 上的原函数,则F ( z ) − G ( z ) ≡ C (常数).
定理二
设f ( z ) 是在单连通区域D 上的解析函数,z 0 是D 内的一个定点,C 是D 内以z 0 为起点,z 为终点的按段光滑曲线,则积分C ∫ f ( ς ) d ς 只依赖于z 与z 0 ,而与路径C 无关.(解析区域上的积分只与起点和终点有关,与路径无关 ,此时才可以用原函数求积分)
定理三
f ( z ) 是在单连通区域D 上的解析函数,z 0 ∈ D 是定点,z ∈ D 是动点,则F ( z ) = ∫ z 0 z f ( ς ) d ς 是f ( z ) 在D 上的原函数.
定理四
设f ( z ) 是在单连通区域D 上的解析函数,F ( z ) 是f ( z ) 在区域D 上的原函数,z 1 , z 2 是D 内的任意两点,则
∫ z 1 z 2 f ( z ) d z = F ( z 2 ) − F ( z 1 )
三、复变函数的级数
18. 复数项级数
复数列
设α n = a n + i b n (n = 1 , 2 , ⋯ )是一列复数,其中a n b n ,
均是实数,则称它为复数列,简称为数列,并记为{α n } .
复数项级数
设α n = a n + i b n 是一复数列(其中a n , b n 均是实数列), 则称
∑ n = 1 ∞ α n = α 1 + α 2 + ⋯ + α n + ⋯ 为复数项级数,把S n = ∑ k = 1 n α k = α 1 + α 2 + ⋯ + α n 称为级数的前n项部分和.
级数∑ n = 1 ∞ α n 收敛的充分必要条件是实数项级数a n 和b n 都收敛。并且当∑ n = 1 ∞ α n 收敛时,∑ n = 1 ∞ α n = ∑ n = 1 ∞ a n + i ∑ n = 1 ∞ b n 。若n → ∞ lim a n = a ,n → ∞ lim b n = b ,则n → ∞ lim ( a n + i b n ) = a + ib .
收敛的性质
n → ∞ lim α n = α 的充要条件是n → ∞ lim a n = a ,n → ∞ lim b n = b ,即实部虚部都收敛.
级数∑ n = 1 ∞ α n 收敛,则n → ∞ lim α n = 0 .(通项趋于零是收敛的必要条件)
若级数∑ n = 1 ∞ α n 绝对收敛,则它一定是收敛级数,此时满足不等式∣ ∑ n = 1 ∞ α n ∣≤ ∑ n = 1 ∞ ∣ α n ∣
设α n 是收敛数列,则α n 有界, 即存在M > 0 ,对任意正整数n ,有α n ≤ M .
设∑ n = 1 ∞ α n 与∑ n = 1 ∞ β n 都是绝对收敛级数,令γ n = α 1 β n + α 2 β n − 1 + ⋯ + α n β 1 ( n = 1 , 2 , ⋯ ) ,则∑ n = 1 ∞ γ n 收敛且∑ n = 1 ∞ γ n = ( ∑ n = 1 ∞ α n ) ( ∑ n = 1 ∞ β n )
性质2的证明:
∣ α n ∣ = a n 2 + b n 2 ,∑ n = 1 ∞ ∣ a n ∣ ≤ ∑ n = 1 ∞ a n 2 + b n 2 ,∑ n = 1 ∞ ∣ b n ∣ ≤ ∑ n = 1 ∞ a n 2 + b n 2
由于∑ n = 1 ∞ a n 2 + b n 2 收敛,∴ ∑ n = 1 ∞ a n 和∑ n = 1 ∞ b n 收敛且绝对收敛,∑ n = 1 ∞ α n 也收敛
审敛法
和高等数学一样
正数项级数
比较判别法,参照物∑ n = 1 ∞ n p 1 = { 收敛 发散 , p > 1 , p ≤ 1 ,∑ n = 1 ∞ q n = { 收敛 发散 , ∣ q ∣ < 1 , ∣ q ∣ ≥ 1
比较判别法极限形式,n → ∞ lim v n u n = l ,若0 < l < + ∞ ,则u n 和v n 同敛散性
比值判别法,n → ∞ lim u n u n + 1 = l ⎩ ⎨ ⎧ 0 ≤ l < 1 l = 1 l > 1 , 收敛 , 不清楚 , 发散
根值判别法,n → ∞ lim n u n = ⎩ ⎨ ⎧ q < 1 q = 1 q > 1 , 收敛 , 不清楚 , 发散
交错级数,莱布尼茨判别法:交错级数∑ n = 1 ∞ ( − 1 ) n u n ,若(1)u n + 1 < u n ,(2)n → ∞ lim u n = 0 ,则收敛,且0 < ∑ n = 1 ∞ ( − 1 ) n u n < u 1
其他级数:u n 取绝对值后,变成正数项级数,判断是否绝对收敛
例1:判断∑ n = 1 ∞ n i n 是否收敛
展开得∑ n = 1 ∞ n i n = 1 i − 2 1 − 3 i + 4 1 ...
把第1、2项看作一对,3、4项看作一对…,得到∑ k = 1 ∞ ( 2 k − 1 i 2 k − 1 + 2 k i 2 k ) = ∑ k = 1 ∞ ( 2 k ( − 1 ) k + i 2 k − 1 ( − 1 ) k − 1 )
根据莱布尼茨判别法,虚部实部都收敛,∴ 级数收敛
例2:判断∑ n = 1 ∞ 5 n ( 2 + 3 i ) n 是否收敛
∑ n = 1 ∞ ∣ 5 n ( 2 + 3 i ) n ∣= ∑ n = 1 ∞ ( ∣5∣ ∣2 + 3 i ∣ ) n ,而5 13 < 1 ,属于∑ n = 1 ∞ q n ,∣ q ∣ < 1 ,∴ 收敛
19. 幂级数的概念
幂级数:∑ n = 0 ∞ C n ( z − z 0 ) n
Abel定理:对于∑ n = 0 ∞ C n z n ,若当z = z 0 时收敛(z 0 = 0 ),则对∀∣ z ∣ < ∣ z 0 ∣ 的z 都收敛
收敛半径:若n → ∞ lim ∣ C n C n + 1 ∣= ρ ,则收敛半径r = ρ 1 。若n → ∞ lim n ∣ C n ∣ = ρ ,则收敛半径r = ρ 1
20. 幂级数的性质
两个幂级数可以在公共收敛区域内作加、减、乘运算
若∑ n = 0 ∞ C n ( z − z 0 ) n ,在∣ z − z 0 ∣ < R 上收敛,则
和函数f ( z ) = ∑ n = 0 ∞ C n ( z − z 0 ) n 在该区域内解析
和函数可以在该区域内逐项求导、逐项积分
与高等数学中的幂级数类似
21. Taylor展开定理
Taylor 展开定理 设f ( z ) 是区域 D
上的解析函数,z 0 是 D 内的一个定点,
而z − z 0 < R 是包含在 D 内的z 0 的最大邻域(R就是z 0 点和 D 的边界点间的距离的最小值(若 D
是全平面时,R = + ∞ ),则f ( z ) 在z − z 0 < R 上可展开成为
f ( z ) = n = 0 ∑ ∞ c n ( z − z 0 ) n
其中c n = n ! 1 f ( n ) ( z 0 ) ( n = 0 , 1 , 2 , ⋯ ) .
用柯西导数公式c n 可以写成2 πi 1 C ∮ ( z − z 0 ) n + 1 f ( z ) d z ,其中C 是D内任意一个以z 0 为圆心的圆∣ z − z 0 ∣ = ρ < R
如何求Taylor级数的收敛半径
设α 是f ( z ) 的所有奇点中离z 0 最近的奇点,则R =∣ z 0 − α ∣ (展开的区域不能跨过奇点或边界的最小值)
Taylor级数收敛半径
例:求( z − 1 ) ( z − 2 ) 1 在z = 1 + i 处展开成幂级数之后的收敛半径
解:奇点有( 1 , 0 ) 、( 2 , 0 ) ,离( 1 , 1 ) 最近的点是( 1 , 0 ) ,距离为1
∴ 收敛半径是1
22. Taylor级数展开的唯一性
设f ( z ) 是 D 上的解析函数,z 0 是 D 内的定点,且在z − z 0 < R 可展成幂级数f ( z ) = ∑ n = 0 ∞ c n ′ ( z − z 0 ) n 则它一定是f ( z ) 在z 0 点的Taylor级数,即c n ′ = n ! f ( n ) ( z ) ( n = 1 , 2 , ⋯ ) .
Taylor级数展开的唯一性定理为间接展开方法提供了理论基础。
将函数展开为Taylor级数 的方法
直接方法:由Taylor展开定理计算级数的系数c n = n ! f ( n ) ( z 0 ) ( n = 0 , 1 , 2 , ⋯ ) ,然后将函数f ( z ) 在z 0 展开成幂级数。
间接方法:借助一些已知函数的展开式,结合解析函数的性质,幂级数运算性质,求函数的Taylor展开式。不需要求各阶导数,通常比直接展开更为简洁,使用范围更广泛。
23. 函数的零点
零点
零点 设函数f ( z ) 在解析区域D 内的一点z 0 的值为零,则称z 0 为解析函数f ( z ) 的零点.
孤立零点 如果函数f ( z ) 在点z 0 的某个邻域B ( z 0 , δ ) = { z ∣ ∣ z − z 0 ∣< δ } 内解析,f ( z 0 ) = 0 ,
且除了点z 0 外,在B ( z 0 , δ ) 内f ( z ) 处处不为零,则称z 0 为f ( z ) 孤立零点.
m级零点 如果解析函数f ( z ) 在点z 0 的邻域内可以表示为f ( z ) = ( z − z 0 ) m ψ ( z ) , 其中ψ ( z ) 在点z 0 解析,且ψ ( z 0 ) = 0 , m ≥ 1 则称z 0 为f ( z ) 的m级零点.
性质
定理_1 不恒为零的解析函数的零点必是孤立零点
推论一:设f ( z ) 在区域D 内解析,z n ( n = 1 , 2 , ⋯ ) 是f ( z ) 在D 内的一列零点 ,z m = z n ( m = n ) ,且z 0 → z n ( n → ∞ ) ,z 0 ∈ D ,则f ( z ) 在D 中恒为零.
推论二:(解析函数的惟一性定理) 设函数f ( z ) 与g ( z ) 在区域D 内解析,z n ( n = 1 , 2 , ⋯ ) 是D 内的点列,z m = z n ( m = n ) ,且z 0 → z n ( n → ∞ ) ,z 0 ∈ D . 若对一切n , 都有f ( z n ) = g ( z n ) , 则在D 内,恒有f ( z ) = g ( z ) 。
注:推论二指出定义在区域D 内的两个解析函数,只要在D 内的某一部分(子区域或孤段)上的值相等,则它们在整个区域D 上的值相等.
定理_2 不恒为零的解析函数f ( z ) 以z 0 为其m级零点的充分必要条件是f ( z 0 ) = f ′ ( z 0 ) = ⋯ = f ( m − 1 ) ( z 0 ) = 0 ,但f ( m ) ( z 0 ) = 0 .
24. Laurent级数的概念
若f ( z ) 在z 0 点解析,则在z 0 的某邻域内能展开为Taylor级数,其各项由z 0 − z 的非负幂组成. 如果f ( z ) 在圆环域R 1 < z − z 0 < R 2 内解析,则f ( z ) 在这个圆环域内不一定都能展开为z 0 − z 的幂级数,但可以展开为含有z − z 0 的负指数次幂项的级数,这类级数可表示为
∑ n = − ∞ ∞ c n ( z − z 0 ) n = ∑ n = 0 + ∞ c n ( z − z 0 ) n + ∑ n = 1 + ∞ c − n ( z − z 0 ) − n
这种形式的级数称为Laurent级数,Laurent级数分为两部分,即由幂级数和只含负幂项的级数组成。
已知幂级数会有收敛半径R 2 ,对于负幂项,作变换t = z − z 0 1 ,则负幂项的部分可以写成c n t n ,也是幂级数的形式,此时有收敛半径μ ,收敛区域是∣ z − z 0 1 ∣ < μ ,即∣ z − z 0 ∣ > μ 1 ,所以对于负幂项,会大于某个R 1
定义一
若∑ n = 0 + ∞ c n ( z − z 0 ) n 和∑ n = 1 + ∞ c − n ( z − z 0 ) − n 都收敛,
则Laurent级数收敛, 若其中一个发散, 则Laurent级数发散。
25. Laurent级数的展开
Laurent 展开定理
设0 ≤ R 1 < R 2 ≤ + ∞ ,f ( z ) 在圆环域 R 1 <∣ z − z 0 ∣< R 2 内解析,则f ( z ) 在此环域内可展开为Laurent级数
f ( z ) = n = − ∞ ∑ + ∞ c n ( z − z 0 ) n ( R 1 <∣ z − z ∣< R 2 )
此处c n = 2 πi 1 C ∮ ( z − z 0 ) n + 1 f ( z ) d z ( n = 0 , ± 1 , ± 2 , ⋯ ) ,C 是圆周∣ z − z 0 ∣= R (其中R 1 < R < R 2 )的正向;C 也可以是把z 0 含在其内部的圆环域内的按段光滑Jordan曲线。
c n 的证明:
已知f ( z ) = ∑ n = − ∞ + ∞ c n ( z − z 0 ) n
两边同除( z − z 0 ) m + 1 ,然后右边把n = m 那项单独拿出来,有
( z − z 0 ) m + 1 f ( z ) = c m z − z 0 1 + n = m ∑ c n ( z − z 0 ) m + 1 − n 1
两边同时积分,根据重要积分C ∮ ( z − z 0 ) n + 1 1 d z = { 2 πi 0 , n = 0 , n = 0 ,得
C ∮ ( z − z 0 ) m + 1 f ( z ) d z = c m ⋅ 2 πi ,得证
注意这里n 可以是负数,不是柯西导数公式,不能写成求n 阶导的形式
特别地,当n = − 1 时,有C ∮ f ( z ) d z = 2 πi c − 1 ,可以用c − 1 计算积分,c − 1 是f ( z ) 在z 0 处展开成Laurent级数后z − z 0 1 项的系数,此时z 0 可以随意选取,以z 0 为圆心的展开的区域只要包含的奇点和原来的一致就行。
Laurent 展开的惟一性
设f ( z ) 在圆环域 R 1 <∣ z − z ∣< R 2 内解析,
则f ( z ) 在此圆环域内的Laurent展开式是惟一的。
例1:求f ( z ) = z 2 − 3 z + 2 1 在z = 0 处展开的Laurent级数
解:奇点是z = 1 和z = 2 ,奇点把展开的区域分成了三部分,在∣ z − 0∣ < 1 时可以展成Taylor级数(可以看作是Laurent的特殊形式),1 < ∣ z − 0∣ < 2 时可以展成Laurent级数,2 < ∣ z − 0∣ < + ∞ 可以展成Laurent级数
(1)∣ z − 0∣ < 1 时,利用等比级数公式间接展开,∑ n = 0 + ∞ q n = 1 − q 1 ,∣ q ∣ < 1
f ( z ) = z 2 − 3 z + 2 1 = z − 2 1 − z − 1 1 = − 2 ( 1 − 2 z ) 1 + 1 − z 1 = − 2 1 ∑ n = 0 + ∞ ( 2 z ) n + ∑ n = 0 + ∞ z n = ∑ n = 0 + ∞ [( 1 − 2 n + 1 1 ) z n ] .
(2)1 < ∣ z − 0∣ < 2 时,f ( z ) = z 2 − 3 z + 2 1 = z − 2 1 − z − 1 1 = − 2 ( 1 − 2 z ) 1 − z 1 1 − z 1 1 = − 2 1 ∑ n = 0 + ∞ ( 2 z ) n − z 1 ∑ n = 0 + ∞ z n 1
(3)2 < ∣ z − 0∣ < + ∞ 时,f ( z ) = z 2 − 3 z + 2 1 = z − 2 1 − z − 1 1 = z 1 ( 1 − z 2 ) 1 − z 1 1 − z 1 1 = z 1 ∑ n = 0 + ∞ ( z 2 ) n − z 1 ∑ n = 0 + ∞ z n 1
例2:求f ( z ) = z 2 − 3 z + 2 1 在z = 1 处展开的Laurent级数
解:奇点将区域分为0 < ∣ z − 1∣ < 1 和1 < ∣ z − 1∣ < + ∞ 两部分,展开的时候注意要凑出q = z − 1 来
(1)0 < ∣ z − 1∣ < 1 时,f ( z ) = z 2 − 3 z + 2 1 = z − 2 1 − z − 1 1 = − 1 − ( z − 1 ) 1 − z − 1 1 = − ∑ n = 0 + ∞ ( z − 1 ) n − ( z − 1 ) − 1
(2)1 < ∣ z − 1∣ < + ∞ 时,f ( z ) = z 2 − 3 z + 2 1 = z − 2 1 − z − 1 1 = z − 1 1 1 − z − 1 1 1 − z − 1 1 = ∑ n = 0 + ∞ ( z − 1 ) − n − 1 − ( z − 1 ) − 1
例3:用c − 1 的方式计算∣ z ∣ = 3 ∮ z ( z + 1 ) 2 1 d z
解:根据C ∮ f ( z ) d z = 2 πi c − 1 ,其中c − 1 是f ( z ) 展成Laurent级数后z − z 0 1 项的系数
选取z 0 = − 1 ,在− 1 处将f ( z ) 展开成Laurent级数,原积分的奇点是z = 0 和z = − 1 ,展开的区域只要满足∣ z − ( − 1 ) ∣ > 1 都能包住这两个奇点
f ( z ) = z ( z + 1 ) 2 1 = ( z + 1 ) 2 1 ( z + 1 ) − 1 1 = ( z + 1 ) 2 1 z + 1 1 1 − z + 1 1 1 = ( z + 1 ) 3 1 ∑ n = 0 + ∞ ( z + 1 1 ) n = ∑ n = 0 + ∞ ( z + 1 1 ) n + 3
由于展开后没有z + 1 1 的项,所以c − 1 = 0 ,∣ z ∣ = 3 ∮ z ( z + 1 ) 2 1 d z = 2 πi c − 1 = 0
(例3补充:如果选取z 0 = 0 处展开,z 1 ( z + 1 ) 2 1 = z 1 ( − z + 1 1 ) ′ ,展开后再利用级数可以逐项求导的性质计算,要比在z 0 = − 1 处展开复杂得多)
例4:计算∣ z ∣ = 2 ∮ 1 − z z e z 1 d z
解:奇点有两个,z = 1 和z = 0 ,而且对于z = 0 ,无法使用柯西积分公式,只能用c − 1 的方法计算
这里在z = 0 处展开较方便,此时展开区域只要满足∣ z − 0∣ > 1 即可包住这两个奇点
1 − z z e z 1 = e z 1 ⋅ z ⋅ ( − z 1 1 − z 1 1 ) = ( ∑ n = 0 + ∞ n ! z n 1 ) ⋅ ( − ∑ m = 0 + ∞ z m 1 )
z 1 的项有n = 0 ,m = 1 和n = 1 ,m = 0 ,c − 1 等于两者系数相加,c − 1 = − 1 + ( − 1 ) = − 2
∣ z ∣ = 2 ∮ 1 − z z e z 1 d z = 2 πi c − 1 = − 4 πi
四、留数及其应用
26. 孤立奇点
如果函数f ( z ) 在z 0 点不解析,则称z 0 是f ( z ) 的一个奇点。如果z 0 是f ( z ) 的一个奇点,且存在δ > 0 , 使f ( z ) 在0 < z − z 0 < δ 内解析,则称z 0 是f ( z ) 的孤立奇点。
不是所有奇点都是孤立奇点,比如s i n z 1 1 ,奇点是0 和kπ 1 ,当k → ∞ 时,kπ 1 → 0 ,所以0 附近有很多奇点,它不是孤立奇点
若z 0 是f ( z ) 的孤立奇点,此时f ( z ) 在圆环域0 < z − z 0 < δ 内解析,根据Laurent级数展开定理,则f ( z ) 可以展开为Laurent级数f ( z ) = ∑ n = 0 + ∞ c n ( z − z 0 ) n ,其中c n = 2 πi 1 P ∮ ( z − z 0 ) n + 1 f ( z ) d z ( n = 0 , ± 1 , ± 2 , ⋯ ) ,而C是以z 0 为中心,半径小于δ的圆周的正向。
根据Laurent级数展开式中c n 的不同情况,可以把f ( z ) 的孤立奇点z 0 分为三种不同类型:
可去奇点
如果f ( z ) 在0 < z − z 0 < δ 内的展开式中不含有z − z 0 的负幂项,即当n = − 1 , − 2 , − 3 , ⋯ 时,展开式中的c n = 0 ,则称z 0 是f ( z ) 的可去奇点。
可去奇点定理 :设f ( z ) 在0 < z − z 0 < δ 上解析,则z 0 是f ( z ) 的可去奇点的充分必要条件是存在极限z → z 0 lim f ( z ) = c 0 (c 0 是有限复常数).
证明:可去奇点就是没有负幂项,因此f ( z ) 可以写成f ( z ) = c 0 + c 1 ( z − z 0 ) + c 2 ( z − z 0 ) 2 + c 3 ( z − z 0 ) 3 ...
当z → z 0 时,c 1 ... c n 的项变0,只剩下c 0 ,易得该证明反过来也成立,∴ 得证
极点
级数有有限个负幂项
此时
f ( z ) = ( z − z 0 ) m c − m + ( z − z 0 ) m − 1 c − ( m − 1 ) + ... + z − z 0 c − 1 + c 0 + c 1 ( z − z 0 ) + ... = ( z − z 0 ) m 1 [ c − m + c − ( m − 1 ) ( z − z 0 ) + ... + c − 1 ( z − z 0 ) m − 1 + c 0 ( z − z 0 ) m + ... ] = ( z − z 0 ) m 1 φ ( z )
当z = z 0 时φ ( z ) 解析,且φ ( z 0 ) = c m = 0 ,此时z 0 是m级极点
比如f ( z ) = ( z 2 + 1 ) ( z − i ) 3 2 z = ( z + i ) ( z − i ) 4 2 z ,奇点是i 和− i ,其中i 是4级极点,− i 是1级极点
对于极点,由展开式可知,z → z 0 lim f ( z ) = ∞
本性奇点
展开式有无穷多个负幂项
比如e z 1 = ∑ n = 0 ∞ n ! 1 ( z 1 ) n ,sin z 1 = z 1 − 3 ! 1 ( z 1 ) 3 + ...
此时z → 0 − lim e z 1 = 0 ,而z → 0 + lim e z 1 = + ∞ ,因此z → 0 lim e z 1 不存在,同理z → 0 lim sin z 1 是震荡的,极限也不存在
对于本性奇点,z → z 0 lim f ( z ) 不存在,也不等于∞
零点
设f ( z ) 在z 0 的邻域内解析,若f ( z 0 ) = 0 ,则称z 0 为f ( z ) 的零点
根据形式判断m级零点:f ( z ) = ( z − z 0 ) m φ ( z ) ,其中φ ( z ) 在z 0 的邻域∣ z − z 0 ∣ < R 内解析,且φ ( z 0 ) = 0 ,此时称z 0 是m级零点
用求导的方式判断m级零点:f ( n ) ( z 0 ) = 0 ( n = 1 , 2 , ⋯ , m − 1 ) ,f ( m ) ( z 0 ) = 0
比如f ( z ) = ( z − 1 ) 2 z ,z = 1 是2级零点,z = 0 是1级零点
再比如f ( z ) = z 2 ( e z − 1 ) = z 2 ( z + 2 z 2 + ... ) = z 3 ( 1 + 2 z + ... ) ,所以0是3级零点
若z = z 0 是f ( z ) 的m级零点,也是g ( z ) 的n级零点,则:
z = z 0 是f ( z ) ⋅ g ( z ) 的m + n 级零点
z = z 0 是g ( z ) f ( z ) 的{ 可去奇点 ( n − m ) 级极点 , m ≥ n , m < n ,注意当m ≥ n 时z 0 不是零点,因为原式分母有z − z 0 ,z = z 0 没有意义
零点和极点的关系:若z = z 0 是f ( z ) 的m级零点,则z = z 0 是f ( z ) 1 的m级极点,反之也成立。
无穷远点
对应与z = 0 的特殊点
f ( z ) 在z = ∞ 附近R < ∣ z ∣ < ∞ 解析,则称z = ∞ 是孤立奇点
对于任意Laurent级数f ( z ) = ∑ m = − ∞ + ∞ c n z n ,做代换z = t 1 ,φ ( t ) = f ( t 1 ) = ∑ m = − ∞ + ∞ c n t n 1 ,讨论t = 0 时的性质,与z = ∞ 对应如下:
φ ( t ) 在t = 0 处是可去奇点⇔ f ( z ) 在z = ∞ 是可取奇点
φ ( t ) 在t = 0 处是m级极点⇔ f ( z ) 在z = ∞ 是m级极点
φ ( t ) 在t = 0 处是本性奇点⇔ f ( z ) 在z = ∞ 是本性奇点
27. 留数的一般理论及留数的计算
留数即留下来的数,对任意f ( z ) ,在z 0 处展开成Laurent级数再两边同时积分,有C ∮ f ( z ) d z = ... + c − m C ∮ ( z − z 0 ) m 1 d z + ... c − 1 C ∮ z − z 0 1 d z + ... + c n C ∮ ( z − z 0 ) n d z ,根据重要积分C ∮ ( z − z 0 ) n + 1 1 d z = { 2 πi 0 , n = 0 , n = 0 ,C ∮ f ( z ) d z = c − 1 C ∮ z − z 0 1 d z = 2 πi c − 1 ,此时积分只和留下来的c − 1 相关,称c − 1 为留数。
留数的定义
设z 0 是f ( z ) 的孤立奇点,C 位于z 0 点的充分小邻域内且把z 0 包含在其内部分段光滑的正向 Jordan
曲线,积分2 πi 1 C ∮ f ( z ) d z 在z 0 点的留数(Residue),并记做Res [ f ( z ) , z 0 ] ,即Res [ f ( z ) , z 0 ] = 2 πi 1 C ∮ f ( z ) d z
如果z 0 是f ( z ) 的孤立奇点,则存在R > 0 ,使得函数f ( z ) 在0 < z − z 0 < R 内可以展开成Laurent级数f ( z ) = ∑ n = 0 + ∞ c n ( z − z 0 ) n ,其中c n = 2 πi 1 C ∮ ( z − z 0 ) n + 1 f ( z ) d z (n = 0 , ± 1 , ± 2 , ⋯ )而C是把z 0 包含在其内部的按段光滑Jordan曲线的正向。由c − 1 = 2 πi 1 P ∮ f ( z ) d z ,可见
Res [ f ( z ) , z 0 ] = c − 1 = 2 πi 1 C ∮ f ( z ) d z
留数基本定理
设D 是复平面上的有界区域,其边界C 是一条或有限条分段光滑的正向Jordan曲线;设函数f ( z ) 在D 内除有限个奇点z 1 , z 2 , ⋯ , z n 外处处解析,且在C 上每一点也解析,则根据复合闭路定理,有
C ∮ f ( z ) d z = 2 πi k = 1 ∑ n Res [ f ( z ) , z k ]
此时留数公式中Res [ f ( z ) , z k ] = 2 πi 1 C k ∮ f ( z ) d z 中,C k 是只圈住z k 的小圆
28. 留数的计算
可去奇点
c − 1 = 0 ,所以Res ( f , z 0 ) = 0
本性奇点
Res ( f , z 0 ) = c − 1
极点
如果z 0 是f ( z ) 的k级极点时,我们给出下面几个计算法则:
设z 0 是f ( z ) 的1级极点,则Res [ f ( z ) , z 0 ] = z → z 0 lim [( z − z 0 ) f ( z )] .
设z 0 是f ( z ) 的1级极点,f ( z ) = Q ( z ) P ( z ) ,P ( z ) , Q ( z ) 都在z 0 处解析,且Q ( z 0 ) = 0 , Q ′ ( z 0 ) = 0 , P ( z 0 ) = 0 ,则Res [ f ( z ) , z 0 ] = Q ′ ( z 0 ) P ( z 0 )
证明:Q ( z 0 ) = 0 ,Q ′ ( z 0 ) = 0 ,说明z 0 是Q ( z ) 的1级零点,也是f ( z ) 的1级极点
根据性质1,Res [ f ( z ) , z 0 ] = z → z 0 lim [( z − z 0 ) f ( z )] = z → z 0 lim z − z 0 Q ( z ) P ( z )
由于Q ( z 0 ) = 0 ,∴ 上式=z → z 0 lim z − z 0 Q ( z ) − Q ( z 0 ) P ( z ) ,由于分子分母极限都存在,可以变成z → z 0 l i m z − z 0 Q ( z ) − Q ( z 0 ) z → z 0 l i m P ( z ) = Q ′ ( z 0 ) P ( z 0 )
设z 0 是f ( z ) 的n级极点( n > 1 ) ,则
Res [ f ( z ) , z 0 ] = ( n − 1 )! 1 z → z 0 lim d z n − 1 d n − 1 [( z − z 0 ) n f ( z )]
证明:z 0 是m级极点,则f ( z ) = ( z − z 0 ) m 1 [ c − m + c − ( m − 1 ) ( z − z 0 ) + ... + c − 1 ( z − z 0 ) m − 1 + c 0 ( z − z 0 ) m + ... ]
( z − z 0 ) m f ( z ) = c − m + c − ( m − 1 ) ( z − z 0 ) + ... + c − 1 ( z − z 0 ) m − 1 + c 0 ( z − z 0 ) m + ...
求m-1阶导后剩下[( z − z 0 ) m f ( z ) ] ( m − 1 ) = ( m − 1 )! c − 1 + 1 ! m ! c 0 ( z − z 0 ) + ...
令z → z 0 ,可以把后面的项去掉,最后剩下( m − 1 )! c − 1
∞ 处的留数
Res ( f , ∞ ) = − c − 1 ,级数展开较方便时
原来在z 0 处的留数是c − 1 = 2 πi 1 C ∮ f ( z ) d z ,此时C 是圈住z 0 的一个圆,展开区域是该圆的内部,正方向是逆时针;当z 0 = ∞ 时,要想圈住∞ ,展开区域变成了一个半径足够大的圆的外部,此时正方向变成了顺时针,方向和原来的相反,所以是− c − 1 .
Res ( f , ∞ ) = − [ Res ( f , z 1 ) + Res ( f , z 2 ) + ⋯ + Res ( f , z n )] ,求孤立奇点的留数较方便时
C ∮ f ( z ) d z + C − ∮ f ( z ) d z = 0 ,此时C 选取一个半径足够大,可以圈住所有孤立奇点的圆
根据留数基本定理,C ∮ f ( z ) d z = 2 πi ∑ k = 1 n Res [ f ( z ) , z k ] ,而C − ∮ f ( z ) d z = 2 πi Res [ f ( z ) , ∞ ] .
Res ( f , ∞ ) = − Res [ f ( t 1 ) t 2 1 , 0 ] ,用t 1 替换更好算时
Res [ f ( z ) , z 0 ] = 2 πi 1 ∣ z ∣ = R ∮ f ( z ) d z ,令z = t 1 ,则d z = − t 2 1 ,Res [ f ( z ) , ∞ ] = − 2 πi 1 ∣ t 1 ∣ = R ∮ f ( t 1 ) t 2 1 d t = Res [ f ( t 1 ) t 2 1 , 0 ]
注意:对于原来的∣ z ∣ = R ,z = r e i θ ,由于是∞ 处的留数,方向是顺时针,换元后t 1 = r e i θ ,t = r 1 e − i θ ,θ 与原来的相反,变成了逆时针
例:计算积分I = ∣ z ∣ = 2 ∮ ( z + i ) 10 ( z − 1 ) ( z − 3 ) 1 d z
解:奇点z = 3 不在∣ z ∣ = 2 内,根据留数基本定理和性质2,I = 2 πi Res [ f ( z ) , − i ] + 2 πi Res [ f ( z ) , 1 ] = − 2 πi Res [ f ( z ) , ∞ ] − 2 πi Res [ f ( z ) , 3 ]
Res [ f ( z ) , 3 ] = ( z − 3 ) f ( z ) ∣ z = 3 = 2 ( 3 + i ) 10 1
Res [ f ( z ) , ∞ ] = Res [ − t 2 1 f ( t 1 ) , 0 ] = Res [ ( 1 + i t ) 10 ( 1 − t ) ( 1 − 3 t ) − t 10 , 0 ] ,显然函数在t = 0 处解析,没有奇点,即化简后不含( t − 0 ) 的负幂项,Res [ f ( z ) , ∞ ] = c − 1 = 0
∴ I = − 2 ( 3 + i ) 10 2 πi
29. 三角有理式的积分
考虑积分∫ 0 2 π R ( cos θ , sin θ ) d θ 其中R ( cos θ , sin θ ) 是cos θ , sin θ 的有理式,且关于θ 连续;令z = e i θ ,则由Euler公式得cos θ = 2 z z 2 + 1 ,sin θ = 2 i z z 2 − 1
且d z = i e i θ d θ ,即d z = i z 1 d z ,于是
∫ 0 2 π R ( cos θ , sin θ ) d θ = ∣ z ∣ = 1 ∮ R ( 2 z z 2 + 1 , 2 i z z 2 − 1 ) i z 1 d z
该式右端的被积函数是z 的有理函数
设z 1 , z 2 , ⋯ , z n 是函数R ( 2 z z 2 + 1 , 2 i z z 2 − 1 ) i z 1 在z < 1 内的孤立奇点,则根据留数基本定理 ∫ 0 2 π R ( cos θ , sin θ ) d θ = 2 πi k = 1 ∑ n Res [ R ( 2 z z 2 + 1 , 2 i z z 2 − 1 ) i z 1 , z k ]
例:∫ 0 2 π 5 + 3 c o s θ 1 d θ
解:原式=∣ z ∣ = 1 ∮ 5 + 3 2 z z 2 + 1 1 i z 1 d z = i 2 ∣ z ∣ = 1 ∮ 3 z 2 + 10 z + 3 1 d z = 3 i 2 ∣ z ∣ = 1 ∮ ( z + 3 1 ) ( z + 3 ) 1 d z = 3 i 2 ⋅ 2 πi ⋅ Res [ ( z + 3 1 ) ( z + 3 ) 1 , − 3 1 ] ,(因为-3在∣ z ∣= 1 外面)
Res [ ( z + 3 1 ) ( z + 3 ) 1 , − 3 1 ] = ( z + 3 1 ) ( z + 3 1 ) ( z + 3 ) 1 ∣ z = − 3 1 = 8 3
∴ 原式=2 πi 3 i 2 8 3 = 2 π
备注:对于∫ 0 π f ( x ) d x ,若f ( x ) 是偶函数,可以变成2 1 ∫ − π π f ( x ) d x ,然后对x 作变换,改变积分上下限,再使用上述方法
30. 有理函数的无穷积分
无穷积分定理
设f ( z ) 在实轴上处处解析,在上半平面Im z > 0 除有限个孤立奇点z 1 , z 2 , ⋯ , z n 外,处处解析,且存在常数R 0 > 0 ,M > 0 ,δ > 0 ,使得当∣ z ∣ > R 0 ,且Im z > 0 时,∣ f ( z ) ∣ ≤ ∣ z ∣ 1 + δ M , 则
∫ − ∞ + ∞ f ( z ) d z = 2 πi k = 1 ∑ n Res [ f ( z ) , z k ]
推 论
设f ( z ) = Q ( z ) P ( z ) 是有理函数
, 多项式Q ( z ) 的次数比P ( z ) 至少高 2 次, Q ( z ) 在实轴上没有零点,z 1 , z 2 , ⋯ , z n 是f ( z ) 在上半平面Im z > 0 的孤立奇点,则∫ − ∞ + ∞ f ( z ) d z = 2 πi ∑ k = 1 n Res [ f ( z ) , z k ]
说明:
有理函数的无穷积分
画一个上半圆,上半圆由两条线组成,一条是沿实轴− R 到R ,另一条记作C R
若∣ z ∣ → ∞ lim z f ( z ) = 0 ,则∣ C R ∫ f ( z ) d z ∣ = 令 z = Re i θ ∣ ∫ 0 π f ( Re i θ ) Re i θ i d θ ∣
由于z = Re i θ ,上式可以写成∣ ∫ 0 π f ( z ) z i d θ ∣ ≤ ∣ ∫ 0 π ε i d θ ∣ = 0 .(由于有条件∣ z ∣ → ∞ lim z f ( z ) = 0 ,当选取的R → ∞ 时,∣ z ∣ = ∣ Re i θ ∣ → ∞ ,就可以用该性质,使得f ( z ) z → ε ,此时ε → 0 )
也就是说只要满足∣ z ∣ → ∞ lim z f ( z ) = 0 ,f ( z ) 沿C R 的积分等于0,∫ − R R + C R ∫ = ∫ − R R ,因此如果实值函数沿实轴− R 到R 不好积分(R → ∞ ),只要被积函数满足一定条件,就可以转到复平面上构造一个上半圆,再用留数基本定理进行积分
观察条件∣ z ∣ → ∞ lim z f ( z ) = 0 ,这意味着f ( z ) 乘以z 后,分母至少还有一个z ,才能让极限等于0,因此该条件可以改为Q ( z ) 的最高幂的次数比P ( z ) 至少高两次
例:∫ − ∞ + ∞ x 4 + 10 x 2 + 9 x 2 − x + 2 d x
解:转到复平面上,把x 改成z ,化简后得∫ − ∞ + ∞ ( z 2 + 1 ) ( z 2 + 9 ) z 2 − z + 2 d z ,奇点有z = i 、− i 、3 i 、− 3 i ,其中i 和3 i 在上半复平面
积分=2 πi Res [ ( z 2 + 1 ) ( z 2 + 9 ) z 2 − z + 2 , i ] + 2 πi Res [ ( z 2 + 1 ) ( z 2 + 9 ) z 2 − z + 2 , 3 i ] = 2 πi ( 16 i 1 − i + 48 i 7 + 3 i ) = 12 5 π
31.
有理函数与三角函数乘积的积分
Jordan引理
设f ( z ) 在区域R 0 ≥ ∣ z ∣ ,Im z ≥ 0 上解析,且当R 0 ≥ z 时,∣ f ( z ) ∣ ≤ M ( ∣ z ∣ ) 其中R 0 > 0 是常数,M ( r ) 是r 的实值函数,且R → + ∞ lim M ( r ) = 0 ,
,则对任何实数m > 0 , 都有R → + ∞ lim C R ∫ f ( z ) e im z = 0 其中C R 为上半圆 (如图)
定 理
做法和有理函数无穷积分一样,从实值函数转到复平面上画上半圆,用留数基本定理计算,这时条件有所不同
设f ( z ) = Q ( z ) P ( z ) ,Q ( z ) 在实轴上没有零点,多项式Q ( z ) 的次数至少比P ( z ) 的次数高1次,z 1 , z 2 , ⋯ , z n 是f ( z ) 在上半平面内的所有孤立奇点,则∫ − ∞ + ∞ f ( z ) e im z d z = 2 πi ∑ k = 1 n Res [ f ( z ) e im z , z k ]
例:计算∫ 0 + ∞ x 2 + a 2 x s i n x d x ,a > 0
解:∫ 0 + ∞ x 2 + a 2 x s i n x d x = 2 1 ∫ − ∞ + ∞ x 2 + a 2 x s i n x d x ,sin x 是e i x 的虚部,即考察2 1 ∫ − ∞ + ∞ x 2 + a 2 x e i x d x 的虚部
z 2 + a 2 z e i z 在上半复平面的奇点只有ai ,∫ − ∞ + ∞ z 2 + a 2 z e i z d z = 2 πi Res [ z 2 + a 2 z e i z , ai ] = 2 πi ( z − ai ) ( z + ai ) ( z − ai ) z e i z ∣ z = ai = 2 πi 2 e − a = π e − a i
2 1 ∫ − ∞ + ∞ x 2 + a 2 x e i x d x = 2 1 ∫ − ∞ + ∞ x 2 + a 2 x ( c o s x + i s i n x ) d x = 2 π e − a i ,而题目∫ 0 + ∞ x 2 + a 2 x s i n x d x 就是其虚部,等于2 π e − a
五、保角映射
32. 映射的概念、导数的几何意义
映射的概念
通常我们用z 平面上的点表示自变量z 的值,用w 平面上的点表示函数w 的值;设w = f ( z ) 是复平面点集D 上的复变函数,平面点集D 称为定义域,用G 表 示由函数值w 所组成的数集,称为f ( z ) 的值域 . 对于z 0 ∈D ,称w 0 = f ( z 0 ) 为从z 平面上的点z 0 变到w 平面上的点w 0 的映射;此时称w 0 = f ( z 0 )为在映射w = f ( z ) 下点z 0 在w 平面上的像,而称z 0 为在映射w = f ( z ) 下w 0 在z 平面上的原像. 对D 中的每个点z ,由w = f ( z ) 构成了w 平面上的点集G ,此时,称w = f ( z ) 把从z 平面上的点集D 映射成w 平面上的点集G ,这个映射通常称为由w = f ( z ) 构成的从D 到G 的映射,记为G = f ( D ) ;这样,把G 称为在映射w = f ( z ) 下D 在w 平面上的像,而D 称为映射w = f ( z ) 下G 在z 平面上的原像。
导数的几何意义
当f ′ ( z 0 ) = 0 时,w = f ( z ) 在z 0 处具有“转动角”不变性和“伸缩率”不变性.
33.
保角映射的概念、关于保角映射的一般概念
保角映射的概念
设w = f ( z ) 在z 0 处的邻域内有定义,如果w = f ( z ) 在z 0 处具有转动角不变性和伸缩率不变性,
即设C 1 与C 2 是过z 0 z 0 点的任意两条光滑曲线,Γ1 和Γ2 为映射w = f ( z ) 下C 1 与C 2 在w 平面上的像,则C 1 与C 2 的交角等于Γ1 与Γ2 的交角(这里的交角包括方向),且映射w = f ( z ) 在z 0 处沿任何方向的伸缩率都相等,则称映射w = f ( z ) 在z 0 处是保角映射。如果w = f ( z ) 在区域D 内的每一点都是保角映射,则称w = f ( z ) 是区域 D 上的保角映射。
定 理 1
设w = f ( z ) 在z 0 处解析,且f ′ ( z 0 ) = 0 ,则w = f ( z ) 是z 0 的邻域内的保角映射;若f ( z ) 在区域D 内解析,且在D 内f ′ ( z ) = 0 ,则w = f ( z ) 是区域D 上的保角映射。
定 理 2 (Riemann 定理)
设D 和G 分别是z 平面和w 平面上边界多于一个点的单连通区域,则必存在双方单值的解析函数w = f ( z ) ,把区域D 保角映射成区域G .
注:通过 Riemann
定理,我们知道两个什么样的区域之间存在保角映射,但Riemann
定理并未告诉我们怎样找出所要求的映射.
定 理 3 (边界对应原理)
设D 是z 平面内由一条按段光滑Jordan曲线C 围成的区域,f ( z ) 是D 及其边界C 上的解析函数,并把C 双方单值地映射成w 平面上的光滑曲线Γ ;如果C 的正向映射成Γ 的正向,则在映射w = f ( z ) 下C 的内部区域D 映射成Γ 正向的左侧(若Γ 也是Jordan曲线,则映射成Γ 的内部)区域;如果C 的正向映射成Γ 的负向,则C 的内部映射成Γ 的右侧(若Γ 也是Jordan曲线,则映射成Γ 的外部)区域。
注:应用边界对应原理,可以求出已给区域D 被函数w = f ( z ) 映射成的区域G .
34. 分式线性映射的基本性质
性质 1
分式线性映射是从扩充复平面到扩充复平面的双方单值保角映射,且其逆映射也是分式线性映射。
性质 2
分式线性映射是具有保圆性的保角映射。换句话说,在分式线性映射之下,z 平面上的圆或者直线映射成w 平面上的圆或者直线。直观上,直线可认为是半径为无穷大的圆,或者认为是通过无穷远点的圆。因此,“保圆性”中的圆应包含直线。
性质 3
分式线性映射是具有保对称性的保角映射, 即设C 是z 平面上的圆(或直线),z 1 、z 2 是关于C 的对称点,在分式线性映射下,w 1 、w 2 及Γ 分别是z 1 , z 2 及C 在w 平面上的像,则w 1 、w 2 是关于Γ 的对称点。
1. 平移映射w = z + b
显然,这是扩充z 平面到扩充w 平面的双方单值映射;当两个复平面放在一起时,它把点z 平移b 得到点w (如图1)
图1
2. 旋转映射w = e i α z (α是实数)
它把点z 以原点为中心旋转α 角(α > 0 时按逆时针,α < 0 时按顺时针)得到点w (如图2)
图2
3.
相似(即放大或缩小)映射w = rz (r >
1时放大,0 < r < 1时缩小)
这是模变化为r 倍,而辐角不变的映射(如图3)
#### 4. 反演映射w = z 1
该映射分解为w = ζ ˉ , ζ = z ˉ 1 的复合映射把z , ζ , w 放在同一个复平面,w 是ζ 关于实轴的对称点,而ζ 是z 关于单位圆z = 1 的对称点(如图4)这是因为arg w = − arg ζ ,而z 与ζ 在同一条以原点为起点的射线上,且∣ z ∣∣ ζ ∣ = 1
因此,反演映射是扩充复平面到扩充复平面的双方单值映射,且把点z 先映射成关于单位圆∣ z ∣ = 1 的对称点之后,再映射成关于实轴的对称点.
图4
35. 唯一确定分式线性映射的条件
惟一确定分式线性映射的条件 设z 1 , z 2 , z 3 是扩充z 平面上的三个互不相同的点,w 1 , w 2 , w 3 是扩充w 平面上的三个互不相同的点,则存在惟一的一个分式线性映射,将点z 1 , z 2 , z 3 依次映射成w 1 , w 2 , w 3 .
36. 幂函数构成的映射
37.
指数函数和对数函数构成的映射
指数函数w = e z 在全平面解析,且( e z ) ′ = e z 处处不为零,因此,它是全平面上的保角映射,但不是全平面上双方单值的映射。事实上,设z = x + i y w = ρ e i ϕ ,则ρ = e x , ϕ = y . 因此,w = e z 把z 平面上平行于虚轴的直线x = x 0 映射成w 平面上的圆周∣ w ∣ = e z 0 ,而平行于实轴的直线y = y 0 映射成从原点出发的射线arg w = ϕ = y 0 . 由此可见,当0 ≤ y 1 ≤ y 2 ≤ 2 π 时,映射w = e z 把带形区域y 1 < Im z < y 2 双方单值地保角映射成角形区域y 1 < arg w < y 2 特别地,把0 < Im z < 2 π 映射成w 平面中从原点起沿正实轴有裂痕的区域,而把0 < Im z < π 映射成上半平面Im z > 0 ,并且是双方单值的保角映射(如图)。 区域− π < Im z < π 上定义的指数函数的反函数就是对数函数w = ln z . 因此,w = ln z 是把上半平面Im z > 0 双方单值映射成带形区域0 < Im w < π 的保角映射,
也是把从原点起沿负实轴有割痕的区域映射成带形区域− π < Im w < π .
指对映射
积分变换
积分变换:对f ( t ) 作积分,F ( α ) = ∫ a b f ( t ) k ( t , α ) d t ,把关于t 的函数变成关于α 的函数,此时f ( t ) 叫做象原函数,F ( α ) 叫做象函数,k ( t , α ) 叫做核函数。
本章讨论两类的积分变换
傅里叶变换:
对于周期的函数f ( t ) ,有傅里叶级数f ( t ) = 2 a 0 + ∑ n = 1 ∞ [ a n cos ( nω t ) + b n sin ( nω t )] 。
对于非周期的函数,有傅里叶变换F ( ω ) = ∫ − ∞ + ∞ f ( t ) e − jω t d t ,其中j 是单位虚数(j 2 = − 1 ,和复变函数部分的i 相同,只是工程上习惯用j 表示),ω 表示频率,t 表示时间。
Laplace变换(拉普拉斯变换):F ( s ) = ∫ 0 + ∞ f ( t ) e − s t ,其中s 是一个复数,Laplace变换不考虑以前的状态,从零时刻开始计算。
六、傅里叶变换
38.傅里叶级数
傅里叶级数 :以T 为周期的函数f ( t ) = f ( t + T ) ,其频率ω = T 2 π ,若满足Dirichlet条件(狄利克雷条件)①有限个第一类间断点.②有限个极值点. (第一类间断点处的函数值为2 f ( t − ) + f ( t + ) )。则f ( t ) 可以展开成以下级数的形式 f ( t ) = 2 a 0 + n = 1 ∑ ∞ [ a n cos ( nω t ) + b n sin ( nω t )] , 其中 a n = T 2 ∫ t 0 t 0 + T f ( t ) cos ( nω t ) d t , b n = T 2 ∫ t 0 t 0 + T f ( t ) sin ( nω t ) d t
傅里叶级数的指数形式 : f ( t ) = n = − ∞ ∑ + ∞ c n e inω t , 其中 c n = T 1 ∫ − 2 T 2 T f ( τ ) e − inω τ d τ
c n 的证明:若把欧拉公式e i z = cos z + i sin z 带入傅里叶级数,可以写成f ( t ) = ∑ n = − ∞ + ∞ c n e inω t 的类似幂级数的形式
为求出c n ,两边同乘e − imω t ,然后对t 在[ − 2 T , 2 T ] 上积分,由于幂级数逐项可积,有
∫ − 2 T 2 T f ( t ) e − imω t d t = ∑ n = − ∞ + ∞ ( c n ∫ − 2 T 2 T e inω t e − imω t d t ) = n = m ∑ ( c n ∫ − 2 T 2 T e i ( n − m ) ω t d t ) + c m T
当n = m 时,∫ − 2 T 2 T e i ( n − m ) ω t d t = ∫ − 2 T 2 T { cos [( n − m ) ω t ] + i sin [( n − m ) ω t ]} d t
由于sin [( n − m ) ω t ] 是奇函数,上式= 2 ∫ 0 2 T cos [( n − m ) ω t ] d t = 2 ( n − m ) ω s i n [( n − m ) ω t ] ∣ 0 2 T = ( n − m ) ω 2 s i n [( n − m ) π ] = 0 ,(因为sin kπ = 0 )
∴ c m = T 1 ∫ − 2 T 2 T f ( t ) e − imω t d t ,注意这里c m 积分完就没有t 了,为了防止和f ( t ) 的t 搞混,一般换成τ
39.傅里叶积分
傅里叶(逆)变换
若任意f ( t ) 满足Dirichlet条件,且∫ − ∞ + ∞ ∣ f ( t ) ∣ d t < + ∞ ,有
傅里叶变换 (F ( ω ) 又叫频谱函数,∣ F ( ω ) ∣ 叫频谱,∣ F ( ω ) ∣ 的图像叫频谱图):
F ( f ( t )) = ∫ − ∞ + ∞ f ( t ) e − iω t d t = F ( ω )
傅里叶逆变换 :
F − 1 ( F ( ω )) = 2 π 1 ∫ − ∞ + ∞ F ( ω ) e iω t d t = f ( t )
(书中的F 会写成花体F ,这里为了简单起见,用F 代替)
证明:前面给出了傅里叶级数的指数形式,这里讨论当T → ∞ 时的情况
f ( t ) = ∑ n = − ∞ + ∞ [ T 1 ∫ − 2 T 2 T f ( τ ) e − imω τ d τ ] e inω t = 2 π ω ∑ n = − ∞ + ∞ [ ∫ − 2 T 2 T f ( τ ) e − imω τ d τ ] e inω t
由于T 1 = 2 π ω ,当T → ∞ 时,ω → 0 ,因此可以把ω 写成Δ ω ,2 T → ∞ ,nω ≈ ω ,对n求和变成对相等的每一项求和
于是f ( t ) = 2 π 1 ∑ − ∞ + ∞ [ ∫ − ∞ + ∞ f ( τ ) e − iω τ d τ ] e iω t Δ ω
这样就可以把求和符号变成积分符号f ( t ) = 2 π 1 ∫ − ∞ + ∞ [ ∫ − ∞ + ∞ f ( τ ) e − iω τ d τ ] e iω t d ω = 2 π 1 ∫ − ∞ + ∞ φ ( ω ) e iω t d ω
此时φ ( ω ) = ∫ − ∞ + ∞ f ( t ) e − iω t d t 叫傅里叶变换,f ( t ) = 2 π 1 ∫ − ∞ + ∞ φ ( ω ) e iω t d ω 叫傅里叶逆变换
余弦积分表示与正弦积分表示
当f ( t ) 是偶函数时,f ( t ) = π 2 ∫ 0 + ∞ [ ∫ 0 + ∞ f ( τ ) cos ( ω τ ) d τ ] cos ( ω t ) d ω = π 2 ∫ 0 + ∞ F c ( ω ) cos ( ω t ) d ω .
当f ( t ) 是奇函数时,f ( t ) = π 2 ∫ 0 + ∞ [ ∫ 0 + ∞ f ( τ ) sin ( ω τ ) d τ ] sin ( ω t ) d ω = π 2 ∫ 0 + ∞ F s ( ω ) sin ( ω t ) d ω .
证明:偶函数时
f ( t ) = 2 π 1 ∫ − ∞ + ∞ [ ∫ − ∞ + ∞ f ( τ ) e − iω τ d τ ] e iω t d ω = 2 π 1 ∫ − ∞ + ∞ { ∫ − ∞ + ∞ f ( τ ) [ cos ( ω τ ) − i sin ( ω τ )] d τ } [ cos ( ω t ) + i sin ( ω t )] d ω = 函数奇偶性 2 π 1 ∫ − ∞ + ∞ [ 2 ∫ 0 + ∞ f ( τ ) cos ( ω τ ) d τ ] [ cos ( ω t ) + i sin ( ω t )] d ω
由于f ( t ) 是实值函数,积分后不会有虚部,所以i sin ( ω t ) 在积分后一定可以消掉,不用计算
∴ 上式= 2 π 1 2 ∫ 0 + ∞ [ 2 ∫ 0 + ∞ f ( τ ) cos ( ω τ ) d τ ] cos ( ω t ) d ω = π 2 ∫ 0 + ∞ [ ∫ 0 + ∞ f ( τ ) cos ( ω τ ) d τ ] cos ( ω t ) d ω
奇函数同理可证
Dirichlet积分
∫ 0 + ∞ ω sin ω d ω = 2 π
证明:对于矩形脉冲函数f ( t ) = { 1 0 , ∣ t ∣ ≤ 1 , 其他 ,作傅里叶积分
F ( ω ) = ∫ − ∞ + ∞ f ( t ) e − iω t d t = ∫ − 1 1 1 ⋅ e − iω t d t = ∫ − 1 1 [ cos ( ω t ) − i sin ( ω t )] d t = 函数奇偶性 ω s i n ( ω t ) ∣ − 1 1 = 2 ω s i n ω
f ( t ) = 2 π 1 ∫ − ∞ + ∞ F ( ω ) e iω t d t = 2 π 2 ∫ − ∞ + ∞ ω s i n ω [ cos ( ω t ) + i sin ( ω t )] d ω = 函数奇偶性 π 2 ∫ 0 + ∞ ω s i n ω ⋅ c o s ( ω t ) d ω
∴ π 2 ∫ 0 + ∞ ω s i n ω ⋅ c o s ( ω t ) d ω = ⎩ ⎨ ⎧ 1 2 1 0 , − 1 < t < 1 , t = ± 1 , t > 1 or t < − 1 ,(端点是第一类间断点,函数值为2 f ( t − ) + f ( t + ) )
当t = 0 时,上述积分变成π 2 ∫ 0 + ∞ ω s i n ω d ω = 1 ,得证
奇偶延拓
对于f ( t ) = { 1 0 , 0 ≤ t ≤ 1 , 其他 ,求其傅里叶变换
作偶延拓得f 1 ( t ) = { 1 0 , − 1 ≤ t ≤ 1 , 其他 ,作奇延拓得f 2 ( t ) = ⎩ ⎨ ⎧ 1 0 − 1 其他 , 0 < t ≤ 1 , t = 0 , − 1 ≤ t < 0
根据奇偶性,分别利用余弦积分表示和正弦积分表示,得到f 1 ( t ) = π 2 ∫ 0 + ∞ ω s i n ω ⋅ c o s ( ω t ) d ω ,f 2 ( t ) = π 2 ∫ 0 + ∞ ω 1 − c o s ω sin ( ω t ) d ω ,
此时f 1 ( t ) 和f 2 ( t ) 在t > 0 的部分是相等的,且都等于f ( t ) ,证明:
f 1 ( t ) − f 2 ( t ) = π 2 ∫ 0 + ∞ ω s i n ω ⋅ c o s ( ω t ) − s i n ( ω t ) + c o s ω ⋅ s i n ( ω t ) d ω = π 2 ∫ 0 + ∞ ω s i n ( ω + ω t ) d ω − ∫ 0 + ∞ ω s i n ( ω t ) d ω = π 2 ∫ 0 + ∞ ω ( 1 + t ) s i n ( ω + ω t ) d [ ω ( 1 + t )] − ∫ 0 + ∞ ω s i n ( ω t ) d ω = 0
如果只需要计算左/右半平面的傅里叶积分,用奇延拓或偶延拓会简单一点。
40.常用傅里叶变换
指数衰减函数
指数衰减函数f ( t ) = { 0 e − βt t < 0 t ≥ 0 ,其中β > 0 ,傅里叶变化和积分表示为:
F ( ω ) = β 2 + ω 2 β − iω f ( t ) = π 1 ∫ 0 + ∞ β 2 + ω 2 β c o s ( ω t ) + ω s i n ( ω t ) d ω
证明:F ( ω ) = ∫ − ∞ + ∞ f ( t ) e − iω t d t = ∫ 0 + ∞ e − ( β + i w ) t d t = − ( β + iω ) e − ( β + iω ) t ∣ 0 + ∞ = β + iω 1 = β 2 + ω 2 β − iω
f ( t ) = 2 π 1 ∫ − ∞ + ∞ F ( ω ) e iω t d ω = 2 π 1 ∫ − ∞ + ∞ β 2 + ω 2 β − iω [ cos ( ω t ) + i sin ( ω t )] d ω = 实值函数虚部一定为 0 2 π 1 ∫ − ∞ + ∞ β 2 + ω 2 β c o s ( ω t ) + ω s i n ( ω t ) d ω = π 1 ∫ 0 + ∞ β 2 + ω 2 β c o s ( ω t ) + ω s i n ( ω t ) d ω
∴ π 1 ∫ 0 + ∞ β 2 + ω 2 β c o s ( ω t ) + ω s i n ( ω t ) d ω = ⎩ ⎨ ⎧ 0 2 1 e − βt , t < 0 , t = 0 , t > 0
δ ( t ) -单位脉冲函数
单位脉冲函数(Dirac函数、δ -函数):δ ( t ) = { 0 ∫ − ∞ + ∞ δ ( t ) d t = 1 , t = 0 , 冲击强度为 1
这个函数意味着在0时刻函数的变化量很大,之后又回到0。当0 < t < ε ,ε → 0 ,∫ 0 ε δ ( t ) d t = 1
性质:
筛选性:∫ − ∞ + ∞ f ( t ) δ ( t ) d t = f ( 0 ) ,∫ − ∞ + ∞ f ( t − t 0 ) δ ( t ) d t = t = 0 时才有值 f ( t 0 ) ,∫ − ∞ + ∞ f ( t ) δ ( t − t 0 ) d t = t = t 0 时才有值 f ( t 0 ) . 2.δ ( t ) 是偶函数. 3.δ ( t ) 是单位阶跃函数u ( t ) = { 0 1 , t < 0 , t > 0 ,在0时刻的导数,用来表示在0时刻变化量无穷大,d t d u ∣ t = 0 = δ ( t ) ,∫ − ∞ t δ ( t ) d t = u ( t ) . 4.∫ − ∞ + ∞ f ( t ) δ ′ ( t ) d t = − f ′ ( t ) 。(证明:用分步积分法∫ − ∞ + ∞ f ( t ) d δ ( t ) = f ( t ) δ ( t ) ∣ − ∞ + ∞ − ∫ − ∞ + ∞ f ′ ( t ) δ ( t ) d t = 0 − f ′ ( t ) )
δ ( t ) 的傅里叶变换为(由筛选性可证):
δ ( t ) ⟷ Fourier F ( ω ) = 1
(用⟷ Fourier 表示一对傅里叶变换对,也可以写成F − 1 ⇌ F ,注意从左到右做傅里叶变换,从右到左做傅里叶逆变换)
常数函数
f ( t ) = 1 F − 1 ⇌ F F ( ω ) = 2 π δ ( ω )
证明:无法通过直接积分得到,但是前面已知δ ( t ) 的傅里叶变换为1,考虑2 π δ ( t ) 的傅里叶逆变换
F − 1 ( 2 π δ ( t )) = 2 π 1 ∫ − ∞ + ∞ 2 π δ ( t ) e iω t d ω = 筛选性 e iω t ∣ t = 0 = 1 ,因此1的傅里叶变换是2 π δ ( t )
这个变换说明了∫ − ∞ + ∞ e − iω t d t = ∫ − ∞ + ∞ cos ( ω t ) − i sin ( ω t ) d t = 奇偶性 ∫ − ∞ + ∞ cos ( ω t ) d t = 2 π δ ( ω ) ,而这个积分本来由于t → ∞ 是原函数sin 是震荡的无法求解,但是可以用δ ( t ) 表示出来
∫ − ∞ + ∞ cos ( ω t ) d t = 2 π δ ( ω )
u ( t ) -单位阶跃函数
u ( t ) ⟷ Fourier F ( ω ) = iω 1 + π δ ( ω )
证明:同样,从右边往左边证
F − 1 ( iω 1 + π δ ( t )) = 2 π 1 ∫ − ∞ + ∞ [ iω 1 + π δ ( t )] e iω t d ω = 2 π 1 ∫ − ∞ + ∞ iω c o s ( ω t ) + i s i n ( ω t ) d ω + 2 1 ∫ − ∞ + ∞ δ ( ω ) e iω t d ω = 奇偶性、筛选性 π 1 ∫ 0 + ∞ ω s i n ( ω t ) d ω + 2 1
利用Dirichlet积分,∫ 0 + ∞ ω s i n ( ω t ) d ω = 令 u = ω t { ∫ 0 + ∞ u s i n u d u = 2 π ∫ 0 − ∞ u s i n u d u = s = − u ∫ 0 + ∞ s s i n s d ( − s ) = − 2 π , t > 0 , t < 0
代入得证
正弦/余弦函数
f ( t ) = cos ( ω 0 t ) ⟷ Fourier F ( ω ) = π δ ( ω − ω 0 ) + π δ ( ω + ω 0 ) f ( t ) = sin ( ω 0 t ) ⟷ Fourier F ( ω ) = iπ δ ( ω + ω 0 ) − iπ δ ( ω − ω 0 )
证明:F ( cos ( ω 0 t )) = ∫ − ∞ + ∞ cos ( ω 0 t ) e iω t d t = ∫ − ∞ + ∞ 2 e i ω 0 t + e − i ω 0 t e − iω t d t = 2 1 ∫ − ∞ + ∞ e − ( ω − ω 0 ) t d t + 2 1 ∫ − ∞ + ∞ e − i ( ω + ω 0 ) t d t
由常数函数的傅里叶变换可知∫ − ∞ + ∞ e − iω t d t = 2 π δ ( ω )
∴ F ( cos ( ω 0 t )) = π δ ( ω − ω 0 ) + π δ ( ω + ω 0 )
正弦同理可证
41.傅里叶变换的性质
若f ( t ) 对应的傅里叶变换是F ( ω ) 则:
线性:α f ( t ) + β g ( t ) ⟷ Fourier α F ( f ( t )) + βF ( g ( t )) 。
时移:f ( t ± t 0 ) ⟷ Fourier e ± i w t 0 F ( ω ) 。(时移等于频域成一个因子)
频移:f ( t ) e ∓ i ω 0 t ⟷ Fourier F ( ω ± ω 0 ) 。
时域微分:若∣ t ∣ → ∞ lim f ( t ) = 0 ,则f ′ ( t ) ⟷ Fourier iω F ( ω ) 。若∣ t ∣ → ∞ lim f ( i ) ( t ) = 0 ,i = 1 , 2 , ... , n − 1 ,f ( n ) ( t ) ⟷ Fourier ( iω ) n F ( ω ) 。
频域微分:( − i t ) f ( t ) ⟷ Fourier F ′ ( ω ) ,( − i t ) n f ( t ) ⟷ Fourier F ( n ) ( ω ) 。该性质使用时一般是利用线性,时域和频域同乘i n ,可以得到t n f ( t ) ⟷ Fourier i n F ( n ) ( ω ) 。
时域积分:若∫ − ∞ t f ( t ) d t → 0 ,则∫ − ∞ t f ( t ) d t ⟷ Fourier iω 1 F ( ω ) 。
相似性(尺度变换):f ( a t ) ⟷ Fourier ∣ a ∣ 1 F ( a ω ) 。
对称性:若f ( t ) ⟷ Fourier F ( ω ) ,则F ( t ) ⟷ Fourier 2 π f ( − ω )
能量积分(帕赛瓦尔定理):∫ − ∞ + ∞ ∣ f ( t ) ∣ 2 d t = 2 π 1 ∫ − ∞ + ∞ ∣ F ( ω ) ∣ 2 d ω ,注意这里∣ f ( t ) ∣ 表示模长,等于f ( t ) ⋅ f ( t ) 。定义时域两函数的内积为[ f 1 ( t ) , f 2 ( t )] = ∫ − ∞ + ∞ f 1 ( t ) f 2 ( t ) d t ,其中f 2 ( t ) 是f 2 ( t ) 的共轭,实值函数的共轭是它本身;定义频域两函数的内积为[ F 1 ( ω ) , F 2 ( ω )] = 2 π 1 ∫ − ∞ + ∞ F 1 ( ω ) F 2 ( ω ) d ω ;自己和自己作内积记作∣∣ f ( t ) ∣ ∣ 2 。则上述定理可记作∣∣ f ( t ) ∣ ∣ 2 = ∣∣ F ( ω ) ∣ ∣ 2 。
线性证明:利用积分的线性,代入公式可证。
时移证明:F [ f ( t − t 0 )] = ∫ − ∞ + ∞ f ( t − t 0 ) e − iω t d t = u = t − t 0 ∫ − ∞ + ∞ f ( u ) e − iω ( u + t 0 ) d u = e − iω t 0 ∫ − ∞ + ∞ f ( u ) e − iω u d u = e − iω t 0 F ( ω )
频移证明同理。
时域微分证明:条件∣ t ∣ → ∞ lim f ( t ) = 0 ,说明在0时刻附近有信号,当t → ∞ 时,离0时刻很远的地方,信号衰减为0。
F ( f ′ ( t )) = ∫ − ∞ + ∞ f ′ ( t ) e − iω t d t = ∫ − ∞ + ∞ e − iω t df ( t ) = f ( t ) e − iω t ∣ − ∞ + ∞ − ∫ − ∞ + ∞ ( − iω ) e − iω t f ( t ) d t = iω F ( ω )
频域微分证明:F ( ω ) = ∫ − ∞ + ∞ f ( t ) e − iω t d t ,两边对ω 求导,F ′ ( ω ) = ∫ − ∞ + ∞ f ( t ) ( − i t ) e − iω t d t = F [ − i t f ( t )]
时域积分证明:令h ( t ) = ∫ − ∞ t f ( t ) d t ,则h ′ ( t ) = f ( t ) ,利用微分性质,F [ h ′ ( t )] = iω F [ h ( t )] ,F [ h ( t )] = iω 1 F [ h ′ ( t )] ,即F [ ∫ − ∞ t f ( t ) d t ] = iω 1 F [ f ( t )] 。
相似性证明:F ( f ( a t )) = ∫ − ∞ + ∞ f ( a t ) e − iω t d t = 令 u = a t { ∫ − ∞ + ∞ f ( u ) e − i a ω u a 1 d u ∫ + ∞ − ∞ f ( u ) e − i a ω u a 1 d u , a > 0 , a < 0 = ∣ a ∣ 1 F ( a ω )
例1:求cos ( ω 0 t ) u ( t ) 的傅里叶变换
解:F ( cos ( ω 0 t ) u ( t )) = F ( 2 e i ω 0 t + e − i ω 0 t u ( t )) = 线性 2 1 F ( e i ω 0 t u ( t )) + 2 1 F ( e − i ω 0 t u ( t )) = 频移 2 1 [ i ( ω − ω 0 ) 1 + π δ ( ω − ω 0 ) + i ( ω + ω 0 ) 1 + π δ ( ω + ω 0 )]
例2:F [ δ ′ ( t )] = 时域微分 iω ⋅ 1 = iω ,F [ δ ( n ) ( t )] = 时域微分 ( iω ) n
例3:求t的傅里叶变换:
解:f ( t ) = t ,若g ( t ) = 1 ,则f ( t ) = t ⋅ g ( t ) ,已知常数函数g ( t ) 的傅里叶变换F ( g ( t )) = 2 π δ ( t )
F ( t ⋅ g ( t )) = 频域微分 i ⋅ [ 2 π δ ( ω ) ] ′ = 2 πi δ ′ ( ω )
同理,t n 傅里叶变换:F ( t n ) = 2 π ⋅ i n ⋅ δ ( n ) ( ω )
例4:求t f ( 2 t ) 的傅里叶变换
解:f ( 2 t ) 相似性 → F 2 1 F ( 2 ω ) ,t f ( 2 t ) 频域微分 → F i [ 2 1 F ( 2 ω ) ] ′ = 2 i F ′ ( 2 ω ) ⋅ 2 1 = 4 i F ′ ( 2 ω )
42.卷积
卷积的定义
系统在某时刻的信号不仅与当前的输入信号幅值有关,还与之前输入信号衰减后的幅值有关,若f ( t ) 表示当前的输入信号,g ( t ) 表示信号随时间衰减程度的函数,那么,比如在t = 10 的时刻,既有当前的输入f ( 10 ) g ( 0 ) ,也有t = 9 时刻输入的信号衰减了1 个时间单位后的幅值,t = 8 时刻输入的信号衰减了2 个时间单位后的幅值,……,把这些信号加起来才是当前时刻的信号,于是定义离散形式的卷积为:f ( n ) ∗ g ( n ) = ∑ τ = − ∞ + ∞ f ( τ ) g ( n − τ ) ,有时也记作( f ∗ g ) ( n ) 。
卷积
定义连续形式的卷积运算为:
f 1 ( t ) ∗ f 2 ( t ) = ( f ∗ g ) ( t ) = ∫ − ∞ + ∞ f 1 ( τ ) f 2 ( t − τ ) d τ
卷积的性质
交换率:f 1 ∗ f 2 = f 2 ∗ f 1 。
结合率:( f 1 ∗ f 2 ) ∗ f 3 = f 1 ∗ ( f 2 ∗ f 3 ) 。
分配律:( f 1 + f 2 ) ∗ f 3 = f 1 ∗ f 3 + f 2 ∗ f 3 。
求导:d t d ( f 1 ∗ f 2 ) = f 1 ∗ d t d f 2 = d t d f 1 ∗ f 2 。(即卷积求导可以转移)
证明:f 1 ( t ) ∗ f 2 ( t ) = ∫ − ∞ + ∞ f 1 ( τ ) f 2 ( t − τ ) d τ ,两边同时对t 求导,因为右边只有f 2 含t ,所以可以转化为只对f 2 求导,由于有交换律,所以也可以转化为只对f 1 求导。积分性质证明同理。
积分:∫ − ∞ ξ [ f 1 ( t ) ∗ f 2 ( t )] d t = ∫ − ∞ ξ f 1 ( t ) d t ∗ f 2 ( t ) = f 1 ( t ) ∗ ∫ − ∞ ξ f 2 ( t ) d t 。
例1:f 1 ( t ) = { 0 1 , t < 0 , t ≥ 0 ,f 2 ( t ) = { 0 e − t , t < 0 , t ≥ 0 ,求f 1 ( t ) ∗ f 2 ( t )
解:f 1 ( t ) ∗ f 2 ( t ) = ∫ − ∞ + ∞ f 1 ( τ ) f 2 ( t − τ ) d τ
讨论τ 的取值范围,对于f 1 ( τ ) ,只有当τ ≥ 0 时,f 1 ( τ ) = 0 ;对于f 2 ( t − τ ) ,只有当t − τ ≥ 0 时,f 2 ( t − τ ) = 0
综上,被积函数只有当0 ≤ τ ≤ t 时,才不等于0
当t ≥ 0 时,f 1 ( t ) ∗ f 2 ( t ) = ∫ 0 t 1 ⋅ e − ( t − τ ) d τ = e − t ∫ 0 t e τ d τ = e − t ( e t − 1 ) = 1 − e − t
当t < 0 时,f 1 ( t ) ∗ f 2 ( t ) = 0
两类特殊的卷积
函数与δ ( t ) 的卷积
f ( t ) ∗ δ ( t ) = f ( t )
证明:对于f ( t ) ∗ δ ( t ) = ∫ − ∞ + ∞ f ( τ ) δ ( t − τ ) d τ = 筛选性 f ( t ) ,即函数与δ ( t ) 的卷积等于函数本身。
对于f ′ ( t ) = f ′ ( t ) ∗ δ ( t ) = f ( t ) ∗ δ ′ ( t ) ,因此求导可以理解为一类特殊的卷积
函数与u ( t ) 的卷积
f ( t ) ∗ u ( t ) = ∫ − ∞ t f ( τ ) d τ
证明:f ( t ) ∗ u ( t ) = ∫ − ∞ + ∞ f ( τ ) u ( t − τ ) d τ = ∫ − ∞ t f ( τ ) u ( t − τ ) d τ + ∫ t + ∞ f ( τ ) u ( t − τ ) d τ
由于当τ < t 时,u ( t − τ ) = 1 ,当τ > t 时,u ( t − τ ) = 0 ,∴ f ( t ) ∗ u ( t ) = ∫ − ∞ t f ( τ ) ⋅ 1 d τ + 0
因此积分也是一类特殊的卷积
卷积定理
f 1 ( t ) ∗ f 2 ( t ) ⟷ Fourier F 1 ( ω ) ⋅ F 2 ( ω ) 2 π f 1 ( t ) ⋅ f 2 ( t ) ⟷ Fourier F 1 ( ω ) ∗ F 2 ( ω )
即,两个函数在时域的卷积等于在频域相乘,在频域的卷积等于在时域相乘再乘以一个系数。
例:求g ( t ) = ∫ − ∞ t f ( ξ ) d ξ 的傅里叶变换
解:∫ − ∞ t f ( ξ ) d ξ = f ( t ) ∗ u ( t ) ,F [ g ( t )] = F [ f ( t )] ⋅ F [ u ( t )] = F ( ω ) ⋅ [ iω 1 + π δ ( ω )] = iω 1 F ( ω ) + π F ( ω ) δ ( ω )
若回顾前面傅里叶变换的时域积分性质:若当t → ∞ ,f ( t ) → 0 ,则∫ − ∞ t f ( t ) d t ⟷ Fourier iω 1 F ( ω ) ,这里说明了,如果不满足前提条件,会多出一个脉冲信号π F ( ω ) δ ( ω ) 。
43. 用傅里叶变换计算积分
对于一些微分/积分方程,可以先用傅里叶变换转到频域上,变成代数方程,解出代数方程后,再用傅里叶逆变换回去。
例1:求满足∫ 0 + ∞ g ( ω ) sin ( ω t ) d ω = h ( t ) 的g ( ω ) ,其中h ( t ) = { 2 π sin t 0 , 0 ≤ t ≤ π , 其它
解:观察该积分,和之前的傅里叶变换正弦表示很像
(傅里叶变换正弦表示:f ( t ) = π 2 ∫ 0 + ∞ F s ( ω ) sin ( ω t ) d ω ,F s ( ω ) = ∫ 0 + ∞ f ( τ ) sin ( ω τ ) d τ )
π 2 h ( t ) = f ( t ) ,则
g ( ω ) = ∫ 0 + ∞ π 2 h ( τ ) s in ( ω τ ) d τ = ∫ 0 π π 2 2 π sin τ ⋅ sin ( ω τ ) d τ = 积化和差 2 1 ∫ 0 π [ cos ( 1 − ω ) τ − cos ( 1 + ω ) τ ] d τ = 1 − ω 2 s i n ( ωπ )
例2:y ′ ( t ) − ∫ − ∞ t y ( t ) d t = − cos t ,求y ( t )
解:若把y ( t ) 的傅里叶变换记作Y ( ω ) ,对上述方程两边同时作傅里叶变换,根据微分、积分性质,有
iωY ( ω ) − iω 1 Y ( ω ) = − [ π δ ( ω + 1 ) + π δ ( ω − 1 )]
整理得Y ( ω ) = 1 + ω 2 iω δ ( ω + 1 ) + iω δ ( ω − 1 )
y ( t ) = 2 π 1 ∫ − ∞ + ∞ Y ( ω ) e iω t d ω = 2 i [ ∫ − ∞ + ∞ 1 + ω 2 ω e iω t δ ( ω + 1 ) d ω + ∫ − ∞ + ∞ 1 + ω 2 ω e iω t δ ( ω − 1 ) d ω ] = δ ( t ) 的筛选性 2 i [ 2 − 1 e − i t + 2 1 e i t ] = − 2 1 2 i e i t − e − i t = − 2 1 sin t
七、拉普拉斯变换
44.拉普拉斯变换
拉普拉斯变换
定义:若f ( t ) 满足①当t > 0 时f ( t ) 分段连续.②f ( t ) 增长不能超过指数倍∣ f ( t ) ∣ ≤ M ⋅ e c t ,其中M > 0 ,c ≥ 0 . ,则f ( t ) 可以写成拉普拉斯变换(Laplace变换)的形式(且此时R e F ( s ) > c 的部分是解析的):
L [ f ( t )] = ∫ 0 + ∞ f ( t ) e − s t d t = F ( s )
(书中的L 会写成花体L ,这里为了简单起见,用L 代替)
Laplace变换与Fourier变换的关系:
作u ( t ) f ( t ) e − βt 的傅里叶变换,F [ u ( t ) f ( t ) e − βt ] = ∫ − ∞ + ∞ u ( t ) f ( t ) e − βt e − iω t d t = ∫ − ∞ 0 0 d t + ∫ 0 + ∞ u ( t ) f ( t ) e − ( β + iω ) t d t = ∫ 0 + ∞ f ( t ) e − ( β + iω ) t d t
若令s = β + iω ,就是Laplace变换
Laplace变换不关心0时刻以前的值,若规定当t < 0 时f ( t ) = 0 ,则一般的函数可以通过乘以u ( t ) 转化为f ( t )
拉普拉斯逆变换
f ( t ) u ( t ) = L − 1 [ F ( s )] = 2 πi 1 ∫ β − i ∞ β + i ∞ F ( s ) e s t d s s k 为 F ( s ) e s t 在整个平面上的奇点 = 若 F ( s ) 分母次数比分子大 1 或以上 k = 1 ∑ n Res [ F ( s ) e s t , s k ]
证明:前面说过,u ( t ) f ( t ) e − βt 的傅里叶变换是拉普拉斯变换,那么拉普拉斯逆变换可以通过傅里叶逆变换得出
拉普拉斯变换F ( s ) = F [ u ( t ) f ( t ) e − βt ] = ∫ 0 + ∞ f ( t ) e − ( β + iω ) t d t ,若将其记作F ( β + iω ) ,则对其作傅里叶逆变换可变回原式
u ( t ) f ( t ) e − βt = 2 π 1 ∫ − ∞ + ∞ F ( β + iω ) e iω t d ω ,两边同乘e βt ,u ( t ) f ( t ) = 2 π 1 ∫ − ∞ + ∞ F ( β + iω ) e ( β + iω ) t d ω
令s = β + iω ,则d s = i d ω ,u ( t ) f ( t ) = 2 πi 1 ∫ β − i ∞ β + i ∞ F ( s ) e s t d s
拉普拉斯逆变换的积分是一条沿x = β 的一条直线,从y = − ∞ 到+ ∞ 的方向的积分,要计算这个积分,可以和前面有理函数的无穷积分一样,画一个半径无穷大的半圆,然后找里面的奇点,通过留数计算,但是由于s = β + iω ,s 可以是任意一个复数,所以让β 尽可能大,使得平面右边没有奇点,此时找的就是在整个平面上的奇点
拉普拉斯逆变换积分
f L + f C R = 2 πi ∑ k = 1 n Res [ F ( s ) e s t , s k ] ,(由于变量是s ,所以奇点z k 改写成s k )
此时如果F ( s ) 满足分母最高次幂比分子大1或以上 ,那么f C R = 0 ,f L F ( s ) e s t d s = 2 πi ∑ k = 1 n Res [ F ( s ) e s t , s k ] ,f ( t ) u ( t ) = L − 1 [ F ( s )] = 2 πi 1 f L F ( s ) e s t d s = ∑ k = 1 n Res [ F ( s ) e s t , s k ]
45.常用拉普拉斯变换
δ ( t ) -单位脉冲函数
δ ( t ) ⟷ Laplace F ( s ) = 1
u ( t ) -单位阶跃函数
u ( t ) ⟷ Laplace F ( s ) = s 1 ,(条件: Re s > 0 )
由于拉普拉斯变换是在[ 0 , + ∞ ] 上的,所以u ( t ) 可以看成是1,即1 ⟷ Laplace s 1
证明:L [ u ( t )] = ∫ 0 + ∞ u ( t ) e − s t d t 若 Re s > 0 = u = − s t − s 1 ∫ 0 − ∞ e u d u = − s 1 e u ∣ u = 0 − ∞ = s 1 ,这里只有当Re s > 0 才能积出来
指数函数
e a t ⟷ Laplace s − a 1 ,(条件: Re ( s − a ) > 0 )
证明:L ( e a t ) = ∫ 0 + ∞ e a t e − s t d t = ∫ 0 + ∞ e − ( s − a ) t d t = − s − a 1 ∫ 0 + ∞ e − ( s − a ) t d [ − ( s − a ) t ] = − s − a 1 e − ( s − a ) t ∣ t = 0 + ∞ = s − a 1
同样这里也要求Re ( s − a ) > 0
幂函数
t n ⟷ Laplace s n + 1 n !
证明:L ( t ) = ∫ 0 + ∞ t e − s t d t = ∫ 0 + ∞ − s t d e − s t = − s t e − s t ∣ 0 + ∞ + ∫ t = 0 + ∞ s 1 e − s t d t = s 2 1
高阶同理(推导会用拉普拉斯变换的性质而不是直接求)
正弦/余弦函数
cos ( k t ) ⟷ Laplace s 2 + k 2 s sin ( k t ) ⟷ Laplace s 2 + k 2 k
L [ cos ( k t )] = ∫ 0 + ∞ cos ( k t ) e − s t d t = ∫ 0 + ∞ 2 e ik t + e − ik t e − s t d t = 2 1 [ ∫ 0 + ∞ e − ( s − ik ) t d t + ∫ 0 + ∞ e − ( s + ik ) t d t ] = 2 1 ( s − i t 1 + s + i t 1 ) = s 2 + k 2 s
正弦同理可证
例:用拉普拉斯变换算积分:∫ 0 + ∞ e − 3 x sin ( 2 x ) d x
可以看出来这是一个s = 3 时的拉普拉斯变换,直接代公式L [ sin ( 2 x )] ∣ s = 3 = 3 2 + 2 2 2 = 13 2
周期函数
若f ( t ) = f ( t + T ) ,则
f ( t ) ⟷ Laplace 1 − e − s T 1 ∫ 0 T f ( t ) e − s t d t
以∣ sin t ∣ 的拉普拉斯变换为例:
L [ ∣ sin t ∣ ] = ∫ 0 + ∞ ∣ sin t ∣ e − s t d t = ∫ 0 π ∣ sin t ∣ e − s t d t + ∫ π 2 π ∣ sin t ∣ e − s t d t ⋯ + ∫ nπ ( n + 1 ) π ∣ sin t ∣ e − s t d t + ⋯
而∫ nπ ( n + 1 ) π ∣ sin t ∣ e − s t d t = x = t − nπ ∫ 0 π ∣ sin ( x + nπ ) ∣ e − s ( x + nπ ) d x = e − s nπ ⋅ ∫ 0 π ∣ sin ( x ) ∣ e − s x d x
因此,原式= ∫ 0 π ∣ sin ( x ) ∣ e − s x d x ⋅ ∑ n = 0 + ∞ e − s nπ = ∫ 0 π ∣ sin ( x ) ∣ e − s x d x ⋅ 1 − e − s π 1
把π 换成T ,把x 换成t ,得证。
该结论说明周期函数的拉普拉斯变换,等于其在一个周期上的拉普拉斯变换乘以一个系数1 − e − s π 1
46.拉普拉斯变换的性质
若f ( t ) 对应的拉普拉斯变换是F ( s ) ,或记作L [ f ( t )] ,则:
线性:α f ( t ) + β g ( t ) ⟷ Laplace αL ( f ( t )) + β L ( g ( t )) 。
时移:u ( t ± t 0 ) f ( t ± t 0 ) ⟷ Laplace e ± s t 0 F ( s ) ,注意这里作逆变换时要多乘以一个u ( t − t 0 ) 。
因为对于f ( t ) ,拉普拉斯变换是在[ 0 , + ∞ ] 上进行积分的,如果换元变成f ( t − t 0 ) ,函数图像向右移动了t 0 ,而拉普拉斯变换本来是不关心f ( t ) 在t < 0 处的值的,但是向右移动会导致本来t < 0 的值出现,因此要乘以u ( t − t 0 ) ,确保在t 0 左侧函数值为0.
频移:e ∓ s 0 t f ( t ) ⟷ Laplace F ( s ± s 0 ) 。
时域微分:f ′ ( t ) ⟷ Laplace s F ( s ) − f ( 0 ) ,f ′′ ( t ) ⟷ Laplace s 2 F ( s ) − s f ( 0 ) − f ′ ( 0 ) ,f ( n ) ( t ) ⟷ Laplace s n F ( s ) − s n − 1 f ( 0 ) − s n − 2 f ′ ( 0 ) − ⋯ − f ( n − 1 ) ( 0 ) 。
频域微分:− t f ( t ) ⟷ Laplace F ′ ( s ) ,( − t ) n f ( t ) ⟷ Laplace F ( n ) ( s ) 。利用线性,可以变成t n f ( t ) ⟷ Laplace ( − 1 ) n F ( n ) ( s ) 。
时域积分:∫ 0 t f ( ξ ) d ξ ⟷ Laplace s 1 F ( s ) 。
频域积分:t 1 f ( t ) ⟷ Laplace ∫ s + ∞ F ( ξ ) d ξ 。特别地,当s = 0 时,L [ t 1 f ( t )] = ∫ 0 + ∞ t 1 f ( t ) e − 0 t d t = ∫ 0 + ∞ t 1 f ( t ) d t = ∫ 0 + ∞ F ( ξ ) d ξ ,若把ξ 写回s ,则有∫ 0 + ∞ t 1 f ( t ) d t = ∫ 0 + ∞ F ( s ) d s 。
相似性(尺度变换):f ( a t ) ⟷ Laplace ∣ a ∣ 1 F ( a s ) 。
时域微分证明:L [ f ′ ( t )] = ∫ 0 + ∞ f ′ ( t ) e − s t d t = ∫ 0 + ∞ e − s t df ( t ) = f ( t ) e − s t ∣ t = 0 + ∞ − ∫ 0 + ∞ f ( t ) ( − s ) e − s t d t = − f ( 0 ) + s L [ f ( t )]
频域微分证明:F ( s ) = ∫ 0 + ∞ f ( t ) e − s t d t ,两边对s 求导得,F ′ ( s ) = ∫ 0 + ∞ f ( t ) ( − t ) e − s t d t = − L [ t f ( t )]
时域积分证明:令h ( t ) = ∫ 0 t f ( ξ ) d ξ ,则h ′ ( t ) = f ( t ) ,利用时域微分性质,L [ h ′ ( t )] = s L [ h ( t )] − h ( 0 ) ,即s L [ ∫ 0 t f ( ξ ) d ξ ] = L [ f ( t )] + 0
相似性证明:L [ f ( a t )] = ∫ 0 + ∞ f ( a t ) e − s t d t = u = a t { a > 0 时 a < 0 时 , ∫ 0 + ∞ f ( u ) e − a s u a 1 d u , ∫ 0 − ∞ f ( u ) e − a s u a 1 d u = ∣ a ∣ 1 F ( a s )
例1:L [ cos ( t − t 0 ) u ( t − t 0 )] = e − s t 0 L [ cos t ] = e − s t s 2 + 1 s
例2:L [ cos ( t − t 0 )] = L [ cos t cos t 0 + sin t sin t 0 ] = cos t 0 s 2 + 1 s + sin t 0 s 2 + 1 1
例3:L [ e a t sin ( k t )] = L [ sin k t ] ∣ s 1 = s − a = ( s − a ) 2 + 1 1
例4:利用时域微分性质也可以求L [ cos ( k t )] :令f ( t ) = cos ( k t ) ,则f ′ ( t ) = − k sin ( k t ) ,f ′′ ( t ) = − k 2 cos ( k t )
L [ f ′′ ( t )] = s 2 F ( s ) − s f ( 0 ) − f ′ ( 0 ) ,左边用线性把− k 2 拿出来,得
− k 2 L [ cos ( k t )] = s 2 L [ cos ( k t )] − s ⋅ 1 − 0 ,整理得L [ cos ( k t )] = s 2 + k 2 s
例5:求L [ t n ]
解:利用时域微分性质,L [( t n ) ( n ) ] = s n L ( t n ) − s n − 1 ⋅ 0 − ... − 0 = s n L ( t n )
而L [( t n ) ( n ) ] = L [ n ! ⋅ 1 ] = n ! s 1 ,∴ L ( t n ) = s n + 1 n !
例6:∫ 0 + ∞ x e − x d x
解:高数中会用分步积分算,但是仔细观察,把x 改成t ,其实这是一个f ( t ) = 1 ,s = 1 的拉氏变换,在时域乘以一个t
1 的拉氏变换为s 1 ,利用频域微分性质,得L [ t ⋅ 1 ] = ( − 1 ) 1 ( s 1 ) ′ ∣ s = 1 = s 2 1 ∣ s = 1 = 1
例7:∫ 0 + ∞ x e − 3 x sin ( 2 x ) d x
解:把x 改成t ,这是一个f ( t ) = s in ( 2 x ) ,s = 3 的拉氏变换,利用频域微分性质,得
L [ t ⋅ sin ( 2 t )] = ( − 1 ) 1 ( s 2 + 2 2 2 ) ′ ∣ s = 3 = ( s 2 + 4 ) 2 4 s ∣ s = 3 = 169 12
例8:利用频域积分性质可以直接计算Dirichlet积分
∫ 0 + ∞ t s i n t d t = ∫ 0 + ∞ L [ sin t ] d s = ∫ 0 + ∞ s 2 + 1 1 d s = arctan s ∣ 0 + ∞ = 2 π
47.卷积
定义
定义:f 1 ( t ) ∗ f 2 ( t ) = ( f ∗ g ) ( t ) = ∫ − ∞ + ∞ f 1 ( τ ) f 2 ( t − τ ) d τ ,和前面的定义一致,但是由于拉普拉斯变换不关心t < 0 时刻的值,所以这里规定当t < 0 时,f 1 ( t ) = f 2 ( t ) = 0 ,要让f 1 、f 2 不等于0,τ 要满足τ > 0 且t − τ > 0 ,因此积分上下限变成从0 到t
f 1 ( t ) ∗ f 2 ( t ) = ∫ 0 t f 1 ( τ ) f 2 ( t − τ ) d τ
此时卷积仍然满足原有的性质。
卷积定理
f 1 ( t ) ∗ f 2 ( t ) ⟷ Laplace F 1 ( s ) ⋅ F 2 ( s ) 2 πi 1 f 1 ( t ) ⋅ f 2 ( t ) ⟷ Laplace F 1 ( s ) ∗ F 2 ( s )
48.拉普拉斯逆变换解法
根据公式求解
例1:L − 1 [ s 2 + 4 3 ] = 2 3 L − 1 [ s 2 + 4 2 ] = 2 3 sin ( 2 t )
例2:L − 1 [ s 2 + 4 s + 2 ] = L − 1 [ s 2 + 4 s + s 2 + 4 2 ] = 线性 L − 1 [ s 2 + 4 s ] + L − 1 [ s 2 + 4 2 ] = sin ( 2 t ) + cos ( 2 t )
例3:L − 1 [ s 2 − s 1 ] = L − 1 [ s − 1 1 − s 1 ] ,其中L ( 1 ) = s 1 ,而s − 1 1 是在s 1 的基础上做了频移,对应的时域应当乘以e 1 ⋅ t ,∴ 原式= e t − 1
例4:L − 1 [ s 2 + 4 s + 13 1 ] = 3 1 L − 1 [ ( s + 2 ) 2 + 9 3 ] = 3 1 e − 2 t sin ( 3 t )
例5:L − 1 [ s 4 1 ] = 3 ! 1 L − 1 [ s 3 + 1 3 ! ] = 6 t 3
用卷积定理求解
例1:L − 1 [ s 2 ( s − 1 ) 1 ] = L − 1 [ s 2 1 ⋅ s − 1 1 ] = t ∗ e t = e t − t − 1
例2:L − 1 [ ( s 2 + 1 ) 2 s ] = L − 1 [ s 2 + 1 s ⋅ s 2 + 1 1 ] = cos t ∗ sin t = ∫ 0 t sin τ ⋅ cos ( t − τ ) d τ = 积化和差 2 1 ∫ 0 t [ sin t + sin ( 2 τ − t )] d τ = 2 1 t sin t
这里也可以用微分性质来算,因为( s 2 + 1 1 ) ′ = − ( s 2 + 1 ) 2 2 s ,− 2 1 ( s 2 + 1 1 ) ′ ⟷ L − 1 − 2 1 ⋅ ( − t ) ⋅ L − 1 [ s 2 + 1 1 ]
用留数计算
例1:L − 1 [ s 2 + 1 s ]
解:分母次数比分子大1,满足条件,奇点有i 、− i ,L − 1 [ s 2 + 1 s ] = Res [ s 2 + 1 s e s t , i ] + Res [ s 2 + 1 s e s t , − i ] = 2 e i t − 2 e − i t = cos t
例2:L − 1 [ s 2 ( s − 1 ) 1 ]
解:奇点有0 、1 ,其中0 是二阶奇点,L − 1 [ s 2 ( s − 1 ) 1 ] = Res [ s 2 ( s − 1 ) 1 e s t , 0 ] + Res [ s 2 ( s − 1 ) 1 e s t , 1 ] = 1 ! 1 ( s − 1 1 e s t ) ′ ∣ s = 0 + s 2 1 e s t ∣ s = 1 = − t − 1 + e t
或者用部分分式法裂项再计算:s 2 ( s − 1 ) 1 = s A + s 2 B + s − 1 C ,两边同乘s 2 ,再令s = 0 ,可得B = − 1 ,同理C = 1 ,剩下的A 可令s 等于随便一个数,比如2 ,得4 1 = 2 A − 4 1 + 1 ,则A = − 1 ,裂项后再根据公式和性质可得出相同的结果。
例3:L − 1 [ s ( s + 1 ) ( s + 2 ) 1 ]
解:奇点是0 、− 1 、− 2 ,L − 1 [ s ( s + 1 ) ( s + 2 ) 1 ] = Res [ s ( s + 1 ) ( s + 2 ) e s t , 0 ] + Res [ s ( s + 1 ) ( s + 2 ) e s t , − 1 ] + Res [ s ( s + 1 ) ( s + 2 ) e s t , − 2 ] = 2 1 − e − t + 2 1 e − 2 t
或者用部分分式法裂项再计算s ( s + 1 ) ( s + 2 ) 1 = s A + s + 1 B + s + 2 C ,两边同乘s ,让s = 0 ,得到A = 2 1 ,同理可得B = − 1 ,C = 2 1
2 1 L − 1 [ s 1 ] − L − 1 [ s + 1 1 ] + 2 1 L − 1 [ s + 2 1 ] = 2 1 − e − t + 2 1 e − 2 t
49.用拉普拉斯变换计算积分
对于一些微分/积分方程,可以先用拉普拉斯变换变成代数方程,解出F ( s ) 后,再用拉普拉斯逆变换回去
例1:y ′′ + 4 y = 0 ,y ( 0 ) = − 2 ,y ′ ( 0 ) = 4
解:若L [ y ( t )] = Y [ s ] ,则L [ y ′′ ( t )] = s 2 Y ( s ) − sy ( 0 ) − y ′′ ( 0 ) = s 2 Y ( s ) + 2 s − 4
∴ 上述微分方程可以写成( s 2 + 4 ) Y ( s ) + 2 s − 4 = 0 ,于是Y ( s ) = s 2 + 4 4 − 2 s
y ( t ) = L − 1 [ Y ( s )] = 2 L − 1 [ s 2 + 4 2 ] − 2 L − 1 [ s 2 + 4 s ] = 2 sin 2 t − 2 cos 2 t
例2:y ′′ + 2 y ′ − 3 y = e − t ,y ( 0 ) = 0 ,y ′ ( 0 ) = 1
解:若L [ y ( t )] = Y ( s ) ,则L [ y ′ ( t )] = s Y ( s ) − y ( 0 ) = s Y ( s ) ,L [ y ′′ ( t )] = s 2 Y ( s ) − sy ( 0 ) − y ′′ ( 0 ) = s 2 Y ( s ) − 1
∴ 上述微分方程可以写成s 2 Y ( s ) − 1 + 2 s Y ( s ) − 3 Y ( s ) = s + 1 1 ,Y ( s ) = ( s + 1 ) ( s + 3 ) ( s − 1 ) s + 2
y ( t ) = L − 1 [ Y ( s )] = L − 1 [ − 4 1 s + 1 1 − 8 1 s + 3 1 + 8 3 s − 1 1 ] = − 4 1 e − t − 8 1 e − 3 t + 8 3 e t
例3:t y ′′ + 2 y ′ + t y = 0 ,y ( 0 ) = 1 ,y ′ ( 0 ) = 2
解:若L [ y ( t )] = Y ( s ) ,则L [ y ′ ( t )] = s Y ( s ) − y ( 0 ) = s Y ( s ) − 1 ,L [ y ′′ ( t )] = s 2 Y ( s ) − s y ′ ( 0 ) − y ( 0 ) = s 2 Y ( s ) − s − 2
而L [ t y ( t )] = 频域微分性质 ( − 1 ) 1 [ Y ( s ) ] ′ = − Y ′ ( s ) ,L [ t y ′′ ( t )] = ( − 1 ) 1 [ s 2 Y ( s ) − s − 2 ] ′ = − 2 s Y ( s ) − s 2 Y ′ ( s ) + 1
∴ 微分方程可以写成− 2 s Y ( s ) − s 2 Y ′ ( s ) + 1 + 2 s Y ( s ) − 2 − Y ′ ( s ) = 0 ,Y ′ ( s ) = − s 2 + 1 1
L − 1 [ Y ′ ( s )] = 频域微分性质 − t y ( t ) = − sin t ,因此y ( t ) = t s i n t
例4:解方程组{ y ′′ − x ′′ − x ′ − y = e t − 2 2 y ′′ − x ′′ − 2 y ′ + x = − t ,{ x ( 0 ) = x ′ ( 0 ) = 0 y ( 0 ) = y ′ ( 0 ) = 0
解:若L [ y ( t )] = Y ( s ) ,L [ x ( t )] = X ( s ) ,则
{ s 2 Y ( s ) − s 2 X ( s ) − s X ( s ) − Y ( s ) = s − 1 1 − s 2 2 s 2 Y ( s ) − s 2 X ( s ) − 2 s Y ( s ) + X ( s ) = − s 2 1 ,{ X ( s ) = s 2 ( s − 1 ) 2 2 s − 1 = − s 2 1 + ( s − 1 ) 2 1 Y ( s ) = s ( s − 1 ) 2 1 = s 1 − s − 1 1 + ( s − 1 ) 2 1
∴ x ( t ) = − t + t e t ,y ( t ) = 1 − e t + t e t
50.信号的分解
周期函数 :
周期函数
对于在x ≥ 0 时以T 为周期的函数f T ( t ) ,若截取[ 0 , T ] 的部分,定义f ( t ) = { f T ( t ) 0 , 0 < t < T , 其他
则根据周期函数的拉氏变换,L [ f T ( t )] = L [ f ( t ) + f ( t − T ) + f ( t − 2 T ) + ⋯ ] = 1 − e − s T 1 L [ f ( t )]
以∣ sin t ∣ 为例,如果要截取[ 0 , π ] 的部分,可以写成sin t − sin ( t − π ) ,这样就不用写成分段函数了
阶梯函数 :
阶梯函数
对于在x > 0 处的阶梯函数,可以用u ( t ) 表示出来,f ( t ) = A u ( t ) + A u ( t − τ ) + A u ( t − 2 τ ) + ...