控制工程数学基础
本文来自于学习《控制工程数学基础(第2版)》马浩 付兴建 主编。感兴趣的同学推荐购买原书学习,乃基础中的基础,务必认真研读!一天写一章行不行!加油!(我发现不行,公式太难敲了,停更,有空找个时间一起发了)
目录
复数与复变函数基础
复数及其代数运算
复数的定义
对于任意两实数x,y,称z=x+jy为复数,其中x和y分别称为z的实部和虚部,记作
x = Re ( z ) , y = Im ( z ) x=\text{Re}\left( z \right) ,\ \ \ \ y=\text{Im}\left( z \right) x=Re(z), y=Im(z)
当x=0时,z=jy称为纯虚数;当y=0时,z=x+0j,这时z就是实数。
复数四则运算规律
- 加法交换律 z 1 + z 2 = z 2 + z 1 z_1+z_2=z_2+z_1 z1+z2=z2+z1
- 乘法交换律 z 1 ⋅ z 2 = z 2 ⋅ z 1 z_1·z_2=z_2·z_1 z1⋅z2=z2⋅z1
- 加法结合律 z 1 + ( z 2 + z 3 ) = ( z 1 + z 2 ) + z 3 z_1+\left( z_2+z_3 \right) =\left( z_1+z_2 \right) +z_3 z1+(z2+z3)=(z1+z2)+z3
- 乘法结合律 z 1 ( z 2 ⋅ z 3 ) = ( z 1 ⋅ z 2 ) z 3 z_1\left( z_2·z_3 \right) =\left( z_1·z_2 \right) z_3 z1(z2⋅z3)=(z1⋅z2)z3
- 乘法对于加法的分配率 z 1 ( z 2 + z 3 ) = z 1 z 2 + z 1 z 3 z_1\left( z_2+z_3 \right) =z_1z_2+z_1z_3 z1(z2+z3)=z1z2+z1z3
共轭复数的运算
- z ‾ ‾ = z \overline{\overline{z}}=z z=z
- z 1 ± z 2 ‾ = z 1 ‾ + z 2 ‾ \overline{z_1\pm z_2}=\overline{z_1}+\overline{z_2} z1±z2=z1+z2
- z 1 z 2 ‾ = z 1 ‾ z 2 ‾ \overline{z_1z_2}=\overline{z_1}\overline{z_2} z1z2=z1z2
- ( z 1 z 2 ) ‾ = z 1 ‾ z 2 ‾ ( z 2 ≠ 0 ) \overline{\left( \frac{z_1}{z_2} \right) }=\frac{\overline{z_1}}{\overline{z_2}}\left( z_2\ne 0 \right) (z2z1)=z2z1(z2=0)
- z z ‾ = [ Re z ] 2 + [ Im z ] 2 z\overline{z}=\left[ \text{Re}z \right] ^2+\left[ \text{Im}z \right] ^2 zz=[Rez]2+[Imz]2
- Re z = z + z ‾ 2 , Im z = z − z ‾ 2 j \text{Re}z=\frac{z+\overline{z}}{2},\text{Im}z=\frac{z-\overline{z}}{2j} Rez=2z+z,Imz=2jz−z
- z = z ‾ ⟺ z 为实数 z=\overline{z}\Longleftrightarrow z\text{为实数} z=z⟺z为实数
复数的表示
几何表示法
x轴称为实轴,y轴称为虚轴,两轴所在的平面称为复平面或z平面,这样复数与复平面上的点一一对应。
向量表示法或极坐标表示法
向量的长度称为z的模或绝对值,记作
∣ z ∣ = r = x 2 + y 2 \left| z \right|=r=\sqrt{x^2+y^2} ∣z∣=r=x2+y2
向量与x轴的夹角θ称为z的辐角,记作
θ = arctan y x 或 A r g z = θ \theta =\arctan \frac{y}{x}\ \ \ \text{或\ \ \ }Arg\ z=\theta θ=arctanxy 或 Arg z=θ
复数的三角表示和指数表示
复数的直角坐标与极坐标的关系如下
x = r cos θ , y = r sin θ x=r\cos \theta ,\ \ \ \ y=r\sin \theta x=rcosθ, y=rsinθ
复数z可以表示为
z = r ( cos θ + j sin θ ) z=r\left( \cos \theta +j\sin \theta \right) z=r(cosθ+jsinθ)
该式称为复数的三角表示法。
再利用欧拉公式 e j θ = cos θ + j sin θ e^{j\theta}=\cos \theta +j\sin \theta ejθ=cosθ+jsinθ,可得
z = r e j θ = ∣ z ∣ e j θ z=re^{j\theta}=\left| z \right|e^{j\theta} z=rejθ=∣z∣ejθ
这种形式成为复数的指数表示法。
那么,有定理一
z 1 ⋅ z 2 = r 1 ⋅ r 2 e j ( θ 1 + θ 2 ) z_1·z_2=r_1·r_2e^{j\left( \theta _1+\theta _2 \right)} z1⋅z2=r1⋅r2ej(θ1+θ2)
定理二
z 2 z 1 = r 2 r 1 e j ( θ 2 − θ 1 ) \frac{z_2}{z_1}=\frac{r_2}{r_1}e^{j\left( \theta _2-\theta _1 \right)} z1z2=r1r2ej(θ2−θ1)
复数的乘幂与方根
n个相同复数z的乘积为z的n次幂,记作 z n z^n zn。对于任何正整数n,有
z n = r n ( cos n θ + j sin n θ ) z^n=r^n\left( \cos n\theta +j\sin n\theta \right) zn=rn(cosnθ+jsinnθ)
如果定义 z − n = 1 / z n z^{-n}=1/z^n z−n=1/zn,那么当n为负数时,上式仍成立。
当z的模 r = 1 r=1 r=1,即 z = cos θ + j sin θ z=\cos \theta +j\sin \theta z=cosθ+jsinθ
( cos θ + j sin θ ) n = cos n θ + j sin n θ \left( \cos \theta +j\sin \theta \right) ^n=\cos n\theta +j\sin n\theta (cosθ+jsinθ)n=cosnθ+jsinnθ
上式即为著名的棣莫弗公式。
例如:
cos 3 θ + j sin 3 θ = ( cos θ + j sin θ ) 3 = cos 3 θ + 3 j cos 2 θ sin θ − 3 cos θ sin 2 θ − j sin 3 θ \cos 3\theta +j\sin 3\theta =\left( \cos \theta +j\sin \theta \right) ^3=\cos ^3\theta +3j\cos ^2\theta \sin \theta -3\cos \theta \sin ^2\theta -j\sin ^3\theta cos3θ+jsin3θ=(cosθ+jsinθ)3=cos3θ+3jcos2θsinθ−3cosθsin2θ−jsin3θ
利用 sin 2 θ = 1 − cos 2 θ \sin ^2\theta =1-\cos ^2\theta sin2θ=1−cos2θ,得
cos 3 θ = cos 3 θ − 3 cos θ sin 2 θ = 4 cos 3 θ − 3 cos θ \cos 3\theta =\cos ^3\theta -3\cos \theta \sin ^2\theta =4\cos ^3\theta -3\cos \theta cos3θ=cos3θ−3cosθsin2θ=4cos3θ−3cosθ
sin 3 θ = 3 cos 2 θ sin θ − sin 3 θ = 3 sin θ − 4 sin 3 θ \sin 3\theta =3\cos ^2\theta \sin \theta -\sin ^3\theta =3\sin \theta -4\sin ^3\theta sin3θ=3cos2θsinθ−sin3θ=3sinθ−4sin3θ
关于复数的求根问题 w n = z w^n=z wn=z
w = z n = r 1 / n ( cos θ + 2 k π n + j sin θ + 2 k π n ) w=\sqrt[n]{z}=r^{1/n}\left( \cos \frac{\theta +2k\pi}{n}+j\sin \frac{\theta +2k\pi}{n} \right) w=nz=r1/n(cosnθ+2kπ+jsinnθ+2kπ)
从几何上看, z n \sqrt[n]{z} nz的n个值就是以原点为中心, r 1 / n r^{1/n} r1/n为半径的圆的内接正 n n n边形的 n n n个顶点。
复变函数与映射定理
设 G G G是一个复数 z = x + j y z=x+jy z=x+jy的集合,如果有一个确定的法则存在,按照这一法则,对于集合 G G G中的每一个复数 z z z,就有一个或几个复数 w = u + j u w=u+ju w=u+ju与之对应,那么称复变函数 w w w是复变数 z z z的函数简称复变函数,记作
w = f ( z ) w=f(z) w=f(z)
如果用 z z z平面上的点表示自变量 z z z的值,而用另一个平面 w w w平面上的点表示函数 w w w的值,那么函数 w = f ( z ) w=f(z) w=f(z)在几何上就可以看作是把 z z z平面上的一个点集 G G G(定义集合)变到 w w w平面上的一个点集 G ∗ G^* G∗(函数值集合)的映射(或变换)。
这个映射通常简称为由函数 w = f ( z ) w=f(z) w=f(z)所构成的映射。
如果 G G G中的点 z z z被映射 w = f ( z ) w=f(z) w=f(z)映射成 G ∗ G^* G∗中的点 w w w,那么 w w w称为 z z z的象(映像),而 z z z称为 w w w的原象。
连续系统时域分析
典型的控制信号
- 直流信号 f ( t ) = A ( − ∞ < t < ∞ ) f(t) =A(-\infty <t<\infty) f(t)=A(−∞<t<∞)
- 正弦信号 f ( t ) = K sin ( ω t + θ ) f(t) =K\sin(\omega t+\theta) f(t)=Ksin(ωt+θ)
- 单位阶跃信号 ε ( t ) = { 1 , t > 0 0 , t < 0 \varepsilon(t) =\left\{ \begin{array}{c}1\text{,}t>0\\0\text{,}t<0\ \end{array} \right. ε(t)={1,t>00,t<0
- 矩形脉冲信号 g τ ( t ) = { 1 , ∣ t ∣ < τ 2 0 , ∣ t ∣ > τ 2 g_{\tau}\left( t \right) =\left\{ \begin{array}{c} 1\text{,}\left| t \right|<\frac{\tau}{2}\\ 0\text{,}\left| t \right|>\frac{\tau}{2}\\ \end{array} \right. gτ(t)={1,∣t∣<2τ0,∣t∣>2τ
g τ ( t ) = ε ( t + τ 2 ) − ε ( t − τ 2 ) g_{\tau}\left( t \right) =\varepsilon \left( t+\frac{\tau}{2} \right) -\varepsilon \left( t-\frac{\tau}{2} \right) gτ(t)=ε(t+2τ)−ε(t−2τ) - 斜波信号 r ( t ) = t ε ( t ) r(t) =t\varepsilon(t) r(t)=tε(t)
- 符号函数 s g n ( t ) = { 1 , t > 0 − 1 , t < 0 sgn(t) =\left\{ \begin{array}{c}1\text{,}t>0\\-1\text{,}t<0\\\end{array} \right. sgn(t)={1,t>0−1,t<0
ε ( t ) = 1 2 [ s g n ( t ) + 1 ] \varepsilon \left( t \right) =\frac{1}{2}\left[ sgn \left( t \right) +1 \right] ε(t)=21[sgn(t)+1] - 实指数信号 f ( t ) = K e − a t f(t)=Ke^{-at} f(t)=Ke−at
- 复指数信号 f ( t ) = K e s t = K e a t + j ω t = K e a t ( c o s ω t + j s i n ω t ) f(t)=Ke^{st}=Ke^{at+j\omega t}=Ke^{at}(cos\omega t+jsin\omega t) f(t)=Kest=Keat+jωt=Keat(cosωt+jsinωt)
e j ω t = cos ω t + j sin ω t e^{j\omega t}=\cos \omega t+j\sin \omega t ejωt=cosωt+jsinωt e − j ω t = cos ω t − j sin ω t e^{-j\omega t}=\cos \omega t-j\sin \omega t e−jωt=cosωt−jsinωt从而有
sin ω t = 1 2 j ( e j ω t − e − j ω t ) \sin \omega t=\frac{1}{2j}\left( e^{j\omega t}-e^{-j\omega t} \right) sinωt=2j1(ejωt−e−jωt) cos ω t = 1 2 ( e j ω t + e − j ω t ) \cos \omega t=\frac{1}{2}\left( e^{j\omega t}+e^{-j\omega t} \right) cosωt=21(ejωt+e−jωt) - 采样函数 S a ( t ) = s i n t t Sa\left( t \right) =\frac{sint}{t} Sa(t)=tsint
- S a ( t ) Sa(t) Sa(t)为偶函数,因为它是 1 t \frac{1}{t} t1和 s i n t sint sint两奇函数的乘积
- 当 t = 0 t=0 t=0时, S a ( 0 ) = lim t → 0 s i n t t = 1 Sa\left( 0 \right) =\underset{t\rightarrow 0}{\lim}\frac{sint}{t}=1 Sa(0)=t→0limtsint=1,且为最大值
- 曲线呈衰减振荡,从 − π -\pi −π到 π \pi π的“主瓣”宽度为 2 π 2\pi 2π,当 t = ± π , ± 2 π , ⋅ ⋅ ⋅ t=\pm \pi \text{,}\pm 2\pi \text{,}··· t=±π,±2π,⋅⋅⋅时 S a ( t ) = 0 Sa(t)=0 Sa(t)=0
- ∫ 0 ∞ S a ( t ) d t = π 2 , ∫ − ∞ ∞ S a ( t ) d t = π \int_0^{\infty}{Sa\left( t \right) dt=\frac{\pi}{2}\text{,}\int_{-\infty}^{\infty}{Sa\left( t \right) dt=\pi}} ∫0∞Sa(t)dt=2π,∫−∞∞Sa(t)dt=π。有时会用到函数 s i n c ( t ) sinc(t) sinc(t),其定义为 s i n c ( t ) = s i n π t π t = S a ( π t ) sinc(t)=\frac {sin\pi t}{\pi t}=Sa(\pi t) sinc(t)=πtsinπt=Sa(πt)。
- 单位冲激信号 δ ( t ) \delta(t) δ(t)
信号的基本运算
- 信号相加与相乘:两个信号相加(相乘)可得到一个新的信号,它在任意时刻的值等于两个信号在该时刻的值之和(积)。
- 信号反转(褶):信号的反转(或反褶)是将信号 f ( t ) f(t) f(t)的自变量 t t t换为 − t -t −t,可得到另一个信号 f ( − t ) f(-t) f(−t)。
- 信号延时:将信号 f ( t ) f(t) f(t)的自变量 t t t换为 t ± t 0 t\pm t_0 t±t0, t 0 t_0 t0为正实数,则可得到另一个信号 f ( t ± t 0 ) f(t\pm t_0) f(t±t0) 。
- 尺度变换:将信号 f ( t ) f(t) f(t)的自变量 t t t换为 a t at at, a a a为正实常数,则信号 f ( a t ) f(at) f(at)将在时间尺度上压缩或扩展。
- 微分运算
y ( t ) = d f ( t ) d t = f ′ ( t ) = f ( 1 ) ( t ) y\left( t \right) =\frac{df\left( t \right)}{dt}=f^{'}\left( t \right) =f^{\left( 1 \right)}\left( t \right) y(t)=dtdf(t)=f′(t)=f(1)(t) - 积分运算
y ( t ) = ∫ − ∞ t f ( τ ) d τ = f ( − 1 ) ( t ) y\left( t \right) =\int_{-\infty}^t{f\left( \tau \right) d\tau =f^{\left( -1 \right)}\left( t \right)} y(t)=∫−∞tf(τ)dτ=f(−1)(t)
常系数线性方程组
对于线性时不变(LTI)控制系统来说,描述这类系统输入-输出特性关系,常用的数学模型是常系数线性微分方程。
对于n阶微分方程
a n y ( n ) ( t ) + a n − 1 y ( n − 1 ) ( t ) + ⋅ ⋅ ⋅ + a 1 y ′ ( t ) + a 0 y ( t ) = b m f ( m ) ( t ) + b m − 1 f ( m − 1 ) ( t ) + ⋅ ⋅ ⋅ + b 1 f ′ ( t ) + b 0 f ( t ) a_ny^{\left( n \right)}\left( t \right) +a_{n-1}y^{\left( n-1 \right)}\left( t \right) +···+a_1y^{'}\left( t \right) +a_0y\left( t \right) =b_mf^{\left( m \right)}\left( t \right) +b_{m-1}f^{\left( m-1 \right)}\left( t \right) +···+b_1f^{'}\left( t \right) +b_0f\left( t \right) any(n)(t)+an−1y(n−1)(t)+⋅⋅⋅+a1y′(t)+a0y(t)=bmf(m)(t)+bm−1f(m−1)(t)+⋅⋅⋅+b1f′(t)+b0f(t)
此方程的完全解由齐次解和特解两部分组成。
零输入响应和零状态响应
系统时域响应是指在典型输入信号作用下,系统的输出量或信号。在时域经典法求解系统的完全响应时,一个广泛应用的分解是把响应分为零输入响应(Zero-Input Response,ZIR)和零状态响应(Zero-State Response,ZSR)两部分,即
y ( t ) = y z i ( t ) + y z s ( t ) y\left( t \right) =y_{zi}\left( t \right) +y_{zs}\left( t \right) y(t)=yzi(t)+yzs(t)
单位阶跃信号和单位冲击信号
单位阶跃信号和单位冲击信号是两个重要的信号,它们的关系为
δ ( t ) = d ε ( t ) d t \delta \left( t \right) =\frac{d\varepsilon \left( t \right)}{dt} δ(t)=dtdε(t) ε ( t ) = ∫ − ∞ t δ ( τ ) d τ \varepsilon \left( t \right) =\int_{-\infty}^t{\delta \left( \tau \right) d\tau} ε(t)=∫−∞tδ(τ)dτ
线性时不变系统的单位冲激响应,是指系统初始状态为零,激励为单位冲激信号作用下的响应,简称冲激响应,用 h ( t ) h(t) h(t)表示。由单位阶跃信号引起的零状态响应称为单位阶跃响应,简称阶跃响应,记为 s ( t ) s(t) s(t)。
单位冲激响应和单位阶跃响应的关系为
h ( t ) = d s ( t ) d t h\left( t \right) =\frac{ds\left( t \right)}{dt} h(t)=dtds(t) s ( t ) = ∫ − ∞ t h ( τ ) d τ s\left( t \right) =\int_{-\infty}^t{h\left( \tau \right) d\tau} s(t)=∫−∞th(τ)dτ
连续系统频域分析的工程数学基础
e j ω t = cos ω t + j sin ω t e^{j\omega t}=\cos \omega t+j\sin \omega t ejωt=cosωt+jsinωt e − j ω t = cos ω t − j sin ω t e^{-j\omega t}=\cos \omega t-j\sin \omega t e−jωt=cosωt−jsinωt
sin ω t = 1 2 j ( e j ω t − e − j ω t ) \sin \omega t=\frac{1}{2j}\left( e^{j\omega t}-e^{-j\omega t} \right) sinωt=2j1(ejωt−e−jωt) cos ω t = 1 2 ( e j ω t + e − j ω t ) \cos \omega t=\frac{1}{2}\left( e^{j\omega t}+e^{-j\omega t} \right) cosωt=21(ejωt+e−jωt)
傅里叶级数
一个以T为周期的函数 f T ( t ) f_T(t) fT(t),如果在 [ − T 2 , T 2 ] [-\frac{T}{2},\frac{T}{2}] [−2T,2T]上满足狄利克雷条件(即函数在 [ − T 2 , T 2 ] [-\frac{T}{2},\frac{T}{2}] [−2T,2T]上满足:①连续或只有有限个第一类间断点;②只有有限个极值点),那么在 [ − T 2 , T 2 ] [-\frac{T}{2},\frac{T}{2}] [−2T,2T]上就可以展开成傅里叶级数。
f ( t ) = a 0 + ∑ n = 1 ∞ ( a n cos n ω 1 t + b n sin n ω 1 t ) f\left( t \right) =a_0+\sum_{n=1}^{\infty}{\left( a_n\cos n\omega _1t+b_n\sin n\omega _1t \right)} f(t)=a0+n=1∑∞(ancosnω1t+bnsinnω1t)
式中, ω 1 = 2 π 2 \omega_1=\frac{2\pi}{2} ω1=22π称为 f ( t ) f(t) f(t)的基波角频率, n ω 1 n\omega_1 nω1称为 n n n次谐波的频率; a 0 a_0 a0为 f ( t ) f(t) f(t)的直流分量, a n a_n an和 b n b_n bn分别为各余弦分量和正弦分量的幅度。
由级数理论可知,傅里叶级数
a 0 = 1 T ∫ 0 T f ( t ) d t a_0=\frac{1}{T}\int_0^T{f\left( t \right) dt} a0=T1∫0Tf(t)dt a n = 2 T ∫ 0 T f ( t ) cos n ω 1 t d t ( n = 1 , 2 , 3 , ⋅ ⋅ ⋅ ) a_n=\frac{2}{T}\int_0^T{f\left( t \right) \cos n\omega _1tdt\left( n=1,2,3,··· \right)} an=T2∫0Tf(t)cosnω1tdt(n=1,2,3,⋅⋅⋅) b n = 2 T ∫ 0 T f ( t ) sin n ω 1 t d t ( n = 1 , 2 , 3 , ⋅ ⋅ ⋅ ) b_n=\frac{2}{T}\int_0^T{f\left( t \right) \sin n\omega _1tdt\left( n=1,2,3,··· \right)} bn=T2∫0Tf(t)sinnω1tdt(n=1,2,3,⋅⋅⋅)
显然,当 f ( t ) f(t) f(t)给定后, a 0 , a n a_0,a_n a0,an和 b n b_n bn可以确定,因而 f ( t ) f(t) f(t)的傅里叶级数展开式可写出。上式即为三角函数式的傅里叶级数。
因为
a n cos n ω 1 t + b n sin ω 1 t = A n cos ( n ω 1 t + φ n ) a_n\cos n\omega _1t+b_n\sin \omega _1t=A_n\cos \left( n\omega _1t+\varphi _n \right) ancosnω1t+bnsinω1t=Ancos(nω1t+φn) A n 2 = a n 2 + b n 2 φ n = − arctan ( b n a n ) A_{n}^{2}=a_{n}^{2}+b_{n}^{2}\ \ \varphi _n=-\arctan \left( \frac{b_n}{a_n} \right) An2=an2+bn2 φn=−arctan(anbn)
故三角函数式的傅里叶级数也可表示成
f ( t ) = a 0 + ∑ n = 1 ∞ A n cos ( n ω 1 t + φ n ) f\left( t \right) =a_0+\sum_{n=1}^{\infty}{A_n\cos \left( n\omega _1t+\varphi _n \right)} f(t)=a0+n=1∑∞Ancos(nω1t+φn)
下面介绍傅里叶级数的另一种形式——复指数表示形式。
对于周期函数的三角级数表达式,利用欧拉公式,可进一步表示为
f ( t ) = a 0 + ∑ n = 1 ∞ ( a n e j n ω 1 t + e − j n ω 1 t 2 + b n e j n ω 1 t − e − j n ω 1 t 2 j ) = a 0 + ∑ n = 1 ∞ ( a n − j b n 2 e j n ω 1 t + a n + j b n 2 e − j n ω 1 t ) f\left( t \right) =a_0+\sum_{n=1}^{\infty}{\left( a_n\frac{e^{jn\omega _1t}+e^{-jn\omega _1t}}{2}+b_n\frac{e^{jn\omega _1t}-e^{-jn\omega _1t}}{2j} \right)}=a_0+\sum_{n=1}^{\infty}{\left( \frac{a_n-jb_n}{2}e^{jn\omega _1t}+\frac{a_n+jb_n}{2}e^{-jn\omega _1t} \right)} f(t)=a0+n=1∑∞(an2ejnω1t+e−jnω1t+bn2jejnω1t−e−jnω1t)=a0+n=1∑∞(2an−jbnejnω1t+2an+jbne−jnω1t)
令
F 0 = a 0 , F n = 1 2 ( a n − j b n ) ( n ≠ 0 ) F_0=a_0\text{,}F_n=\frac{1}{2}\left( a_n-jb_n \right) \ \ \left( n\ne 0 \right) F0=a0,Fn=21(an−jbn) (n=0)
由三角函数式的傅里叶级数可知 a n = a − n , b n = − b − n a_n=a_{-n},b_n=-b_{-n} an=a−n,bn=−b−n,因而有
F − n = 1 2 ( a − n − j b − n ) = 1 2 ( a n + j b n ) F_{-n}=\frac{1}{2}(a_{-n}-jb_{-n})=\frac{1}{2}(a_n+jb_n) F−n=21(a−n−jb−n)=21(an+jbn)
那么将 F n , F − n F_n,F_{-n} Fn,F−n带入上式得到
f ( t ) = a 0 + ∑ n = 1 ∞ ( F n e j n ω 1 t + F − n e − j n ω 1 t ) f\left( t \right) =a_0+\sum_{n=1}^{\infty}{\left( F_ne^{jn\omega _1t}+F_{-n}e^{-jn\omega _1t} \right)} f(t)=a0+n=1∑∞(Fnejnω1t+F−ne−jnω1t)
即
f ( t ) = ∑ n = − ∞ ∞ F n e j n ω 1 t f\left( t \right) =\sum_{n=-\infty}^{\infty}{F_ne^{jn\omega _1t}} f(t)=n=−∞∑∞Fnejnω1t
式中, F 0 = a 0 F_0=a_0 F0=a0,上式即为 f ( t ) f(t) f(t)的复指数级数形式。不过,式中的负频率只是一种数学表示,并无实际意义。系数 F n F_n Fn通常是一复数,其求法推导如下:
F n = 1 2 ( a n − j b n ) = 1 2 { 2 T ∫ − T 2 T 2 [ f ( t ) cos n ω 1 t − j f ( t ) sin n ω 1 t ] d t } F_n=\frac{1}{2}\left( a_n-jb_n \right) =\frac{1}{2}\left\{ \frac{2}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}{\left[ f\left( t \right) \cos n\omega _1t-jf\left( t \right) \sin n\omega _1t \right] dt} \right\} Fn=21(an−jbn)=21{T2∫−2T2T[f(t)cosnω1t−jf(t)sinnω1t]dt} = 1 T ∫ − T 2 T 2 f ( t ) ( cos n ω 1 t − j sin n ω 1 t ) d t = 1 T ∫ − T 2 T 2 f ( t ) e − j n ω 1 t d t = \frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}{f\left( t \right) \left( \cos n\omega _1t-j\sin n\omega _1t \right) dt=\frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}{f\left( t \right) e^{-jn\omega _1t}dt}} =T1∫−2T2Tf(t)(cosnω1t−jsinnω1t)dt=T1∫−2T2Tf(t)e−jnω1tdt
上式表明,只要给定周期函数 f ( t ) f(t) f(t),则 F n F_n Fn可以在一个周期内积分确定,继而可写出复指数形式的傅里叶级数。那么上面两式是表示周期函数的一对重要关系。
F n F_n Fn为各次谐波 n ω 1 n\omega_1 nω1的函数,可表示为
F n = ∣ F n ∣ e j φ n F_n=\left| F_n \right|e^{j\varphi _n} Fn=∣Fn∣ejφn
∣ F n ∣ |F_n| ∣Fn∣称为各次谐波的幅度, φ n \varphi_n φn称为各次谐波的相位。
傅里叶变换
- 傅里叶变换是一种线性积分变换。傅里叶变换对为
f ( t ) = 1 2 π ∫ − ∞ ∞ F ( ω ) e j ω t d ω f\left( t \right) =\frac{1}{2\pi}\int_{-\infty}^{\infty}{F\left( \omega \right) e^{j\omega t}\text{d}\omega} f(t)=2π1∫−∞∞F(ω)ejωtdω F ( ω ) = ∫ − ∞ ∞ f ( t ) e − j ω t d t F\left( \omega \right) =\int_{-\infty}^{\infty}{f\left( t \right) e^{-j\omega t}dt} F(ω)=∫−∞∞f(t)e−jωtdt
对于周期函数上述的一对关系
f ( t ) = ∑ n = − ∞ ∞ F n e j n ω 1 t f\left( t \right) =\sum_{n=-\infty}^{\infty}{F_ne^{jn\omega _1t}} f(t)=n=−∞∑∞Fnejnω1t F n = 1 T ∫ − T 2 T 2 f ( t ) e − j n ω 1 t d t F_n=\frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}{f\left( t \right) e^{-jn\omega _1t}dt} Fn=T1∫−2T2Tf(t)e−jnω1tdt
F n F_n Fn为离散值 n ω 1 n\omega_1 nω1的函数,可改写为
F ( n ω 1 ) = F n T = ∫ − T 2 T 2 f ( t ) e − j n ω 1 t d t F(n\omega_1)=F_nT=\int_{-\frac{T}{2}}^{\frac{T}{2}}{f\left( t \right) e^{-jn\omega _1t}dt} F(nω1)=FnT=∫−2T2Tf(t)e−jnω1tdt
在 T → ∞ T\rightarrow\infty T→∞时应有如下关系
{ ω 1 = 2 π T → Δ ω → d ω n ω 1 → n Δ ω → ω F n → 0 \left\{ \begin{array}{l} \omega _1=\frac{2\pi}{T}\rightarrow \varDelta \omega \rightarrow d\omega\\ n\omega _1\rightarrow n\varDelta \omega \rightarrow \omega\\ F_n\rightarrow 0\\ \end{array} \right. ⎩ ⎨ ⎧ω1=T2π→Δω→dωnω1→nΔω→ωFn→0
即可得到傅里叶正变换
F ( ω ) = ∫ − ∞ ∞ f ( t ) e − j ω t d t F\left( \omega \right) =\int_{-\infty}^{\infty}{f\left( t \right) e^{-j\omega t}dt} F(ω)=∫−∞∞f(t)e−jωtdt
由于
F ( ω ) = lim T → ∞ F n T = 2 π F n d ω F\left( \omega \right) =\underset{T\rightarrow \infty}{\lim}F_nT=\frac{2\pi F_n}{d\omega} F(ω)=T→∞limFnT=dω2πFn
可见 F ( ω ) F(\omega) F(ω)相当于单位频率占有的幅度,具有密度的意义,所以常把 F ( ω ) F(\omega) F(ω)称为频谱密度函数,简称频谱函数。即 F ( ω ) F(\omega) F(ω)表达了函数在 ω \omega ω处的频谱密度分布情况,这就是函数的傅里叶变换的物理含义。对函数进行频谱分析与对其进行傅里叶变换有着同样的含义。
下面由函数的 F ( ω ) F(\omega) F(ω)重建非周期函数 f ( t ) f(t) f(t)的表达式
f T ( t ) = ∑ n = − ∞ ∞ F ( n ω 1 ) T e j n ω 1 t = ∑ n = − ∞ ∞ F ( n ω 1 ) 2 π ω 1 e j n ω 1 t f_T\left( t \right) =\sum_{n=-\infty}^{\infty}{\frac{F(n\omega_1)}{T} e^{jn\omega _1t}}=\sum_{n=-\infty}^{\infty}{\frac{F(n\omega_1)}{2\pi}\omega_1e^{jn\omega _1t}} fT(t)=n=−∞∑∞TF(nω1)ejnω1t=n=−∞∑∞2πF(nω1)ω1ejnω1t
当 T → ∞ T\rightarrow\infty T→∞时,有
f ( t ) = lim T → ∞ f T ( t ) = lim T → ∞ ∑ n = − ∞ ∞ F ( n ω 1 ) 2 π ω 1 e j n ω 1 t = 1 2 π lim T → ∞ ∑ n = − ∞ ∞ F ( ω ) e j n ω 1 t Δ ω = 1 2 π ∫ − ∞ ∞ F ( ω ) e j ω t d ω f\left( t \right) =\underset{T\rightarrow \infty}{\lim}f_T\left( t \right) =\underset{T\rightarrow \infty}{\lim}\sum_{n=-\infty}^{\infty}{\frac{F\left( n\omega _1 \right)}{2\pi}\omega _1e^{jn\omega _1t}=\frac{1}{2\pi}\underset{T\rightarrow \infty}{\lim}\sum_{n=-\infty}^{\infty}{F\left( \omega \right) e^{jn\omega _1t}\varDelta \omega}=\frac{1}{2\pi}\int_{-\infty}^{\infty}{F\left( \omega \right) e^{j\omega t}d\omega}} f(t)=T→∞limfT(t)=T→∞limn=−∞∑∞2πF(nω1)ω1ejnω1t=2π1T→∞limn=−∞∑∞F(ω)ejnω1tΔω=2π1∫−∞∞F(ω)ejωtdω
即得到傅里叶逆变换
f ( t ) = 1 2 π ∫ − ∞ ∞ F ( ω ) e j ω t d ω f\left( t \right) =\frac{1}{2\pi}\int_{-\infty}^{\infty}{F\left( \omega \right) e^{j\omega t}\text{d}\omega} f(t)=2π1∫−∞∞F(ω)ejωtdω
频谱函数 F ( ω ) F(\omega) F(ω)一般为 ω \omega ω的复函数。故有时把 F ( ω ) F(\omega) F(ω)记为 F ( j ω ) F(j\omega) F(jω)。进一步地, F ( ω ) F(\omega) F(ω)可写为
F ( ω ) = ∣ F ( ω ) ∣ e j φ ( ω ) F(\omega)=|F(\omega)|e^{j\varphi(\omega)} F(ω)=∣F(ω)∣ejφ(ω)
式中, ∣ F ( ω ) ∣ |F(\omega)| ∣F(ω)∣称为非周期函数的幅度频谱; φ ( ω ) \varphi(\omega) φ(ω)称为非周期函数的相位频谱。幅度谱和相位谱都是频率 ω \omega ω的连续函数。
那么还可以继续得出
F ( ω ) = ∫ − ∞ ∞ f ( t ) e − j ω t d t = ∫ − ∞ ∞ f ( t ) cos ω t d t − j ∫ − ∞ ∞ f ( t ) sin ω t d t = a ( ω ) − j b ( ω ) F\left( \omega \right) =\int_{-\infty}^{\infty}{f\left( t \right) e^{-j\omega t}dt}=\int_{-\infty}^{\infty}{f\left( t \right) \cos \omega tdt-j\int_{-\infty}^{\infty}{f\left( t \right) \sin \omega tdt=a\left( \omega \right)}-jb\left( \omega \right)} F(ω)=∫−∞∞f(t)e−jωtdt=∫−∞∞f(t)cosωtdt−j∫−∞∞f(t)sinωtdt=a(ω)−jb(ω)
式中, a ( ω ) = ∫ − ∞ ∞ f ( t ) cos ω t d t a(\omega)=\int_{-\infty}^{\infty}f\left( t \right) \cos \omega tdt a(ω)=∫−∞∞f(t)cosωtdt为 ω \omega ω的偶函数; b ( ω ) = j ∫ − ∞ ∞ f ( t ) sin ω t d t b(\omega)=j\int_{-\infty}^{\infty}f\left( t \right) \sin \omega tdt b(ω)=j∫−∞∞f(t)sinωtdt为 ω \omega ω的奇函数;从而有 ∣ F ( ω ) ∣ = a 2 ( ω ) + b 2 ( ω ) \left| F\left( \omega \right) \right|=\sqrt{a^2\left( \omega \right) +b^2\left( \omega \right)} ∣F(ω)∣=a2(ω)+b2(ω)为 ω \omega ω的偶函数; φ ( ω ) = − arctan [ b ( ω ) a ( ω ) ] \varphi \left( \omega \right) =-\arctan \left[ \frac{b\left( \omega \right)}{a\left( \omega \right)} \right] φ(ω)=−arctan[a(ω)b(ω)]为 ω \omega ω的奇函数。
非周期函数 f ( t ) f(t) f(t)是否存在傅里叶变换 F ( ω ) F(\omega) F(ω),仍应满足类似于傅里叶级数的狄利克雷条件,不同之处仅仅在于一个周期的范围,即要求函数在无限区间内绝对可积。
∫ − ∞ ∞ ∣ f ( t ) ∣ d t < ∞ \int_{-\infty}^{\infty}{\left| f\left( t \right) \right|dt<\infty} ∫−∞∞∣f(t)∣dt<∞
但这仅是充分条件,而不是必要条件。凡满足绝对可积条件的函数,它的变换 F ( ω ) F(\omega) F(ω)必然存在,但不满足上式的函数,其傅里叶变换也可能纯在。
- 傅里叶变换的主要性质如下:
- 线性性质 a 1 f 1 ( t ) + a 2 f 2 ( t ) ↔ a 1 F 1 ( ω ) + a 2 F 2 ( ω ) a_1f_1\left( t \right) +a_2f_2\left( t \right) \leftrightarrow a_1F_1\left( \omega \right) +a_2F_2\left( \omega \right) a1f1(t)+a2f2(t)↔a1F1(ω)+a2F2(ω)
- 延时性质 f ( t ± t 0 ) ↔ F ( ω ) e ± j ω t 0 f\left( t\pm t_0 \right) \leftrightarrow F\left( \omega \right) e^{\pm j\omega t_0} f(t±t0)↔F(ω)e±jωt0
- 尺度变换 f ( a t ) ↔ 1 ∣ a ∣ F ( ω a ) f\left( at \right) \leftrightarrow \frac{1}{\left| a \right|}F\left( \frac{\omega}{a} \right) f(at)↔∣a∣1F(aω)
- 频移特性 f ( t ) e j ω 0 t ↔ F ( ω − ω 0 ) f\left( t \right) e^{j\omega _0t}\leftrightarrow F\left( \omega -\omega _0 \right) f(t)ejω0t↔F(ω−ω0)
- 时域微分特性 d f ( t ) d t ↔ j ω F ( ω ) \frac{df\left( t \right)}{dt}\leftrightarrow j\omega F\left( \omega \right) dtdf(t)↔jωF(ω)
- 时域积分特性 ∫ − ∞ t f ( τ ) d τ ↔ π F ( 0 ) δ ( ω ) + F ( ω ) j ω \int_{-\infty}^t{f\left( \tau \right) d\tau \leftrightarrow \pi F\left( 0 \right) \delta \left( \omega \right)}+\frac{F\left( \omega \right)}{j\omega} ∫−∞tf(τ)dτ↔πF(0)δ(ω)+jωF(ω)
- 时-频对称性质 F ( t ) ↔ 2 π f ( − ω ) F\left( t \right) \leftrightarrow 2\pi f\left( -\omega \right) F(t)↔2πf(−ω)
- 卷积定理 f 1 ( t ) ∗ f 2 ( t ) ↔ F ( ω ) ⋅ F ( ω ) f_1\left( t \right) \ast f_2\left( t \right) \leftrightarrow F\left( \omega \right) \cdot F\left( \omega \right) f1(t)∗f2(t)↔F(ω)⋅F(ω)
- 傅里叶变换的应用
根据傅里叶变换的线性性质,微分性质和积分性质,对要求解的微分方程两端取傅里叶变换,将其转化为象函数的代数方程,由这个代数方程求出象函数,然后再取傅里叶逆变换就得出这类微分方程的解
连续系统复频域分析的工程数学基础
离散系统的工程数学基础
附录A 数学发展简史
附录B 工程数学三大变换间的关系
附录C 常用函数的三大变换对比表

魔乐社区(Modelers.cn) 是一个中立、公益的人工智能社区,提供人工智能工具、模型、数据的托管、展示与应用协同服务,为人工智能开发及爱好者搭建开放的学习交流平台。社区通过理事会方式运作,由全产业链共同建设、共同运营、共同享有,推动国产AI生态繁荣发展。
更多推荐
所有评论(0)