信号与系统复习笔记——傅里叶变换

信号与系统复习笔记——傅里叶变换

周期信号的傅里叶级数表示

特征函数

假设LTI系统的输入为 x(t)=estx(t) = e^{st}x(t)=est 输出为:

y(t)=est∗h(t)=∫−∞+∞es(t−τ)h(τ)dτ=est∫−∞+∞e−sτh(τ)dτ=x(t)H(s) y(t) = e^{st} \ast h(t) = \int_{-\infty}^{+\infty} e^{s(t - \tau)} h(\tau) d\tau = e^{st}\int_{-\infty}^{+\infty}e^{-s\tau}h(\tau) d\tau = x(t)H(s) y(t)=esth(t)=+es(tτ)h(τ)dτ=est+esτh(τ)dτ=x(t)H(s)

定义LTI系统的特征函数为:

H(s)=∫−∞+∞e−sτh(τ)dτ H(s) = \int_{-\infty}^{+\infty}e^{-s\tau}h(\tau) d\tau H(s)=+esτh(τ)dτ

CT周期函数的傅里叶级数表示

对于周期函数 xT(t)x_T(t)xT(t) 来说,周期为 TTT ,角频率为 ω0=2πT\omega_0 = \frac{2 \pi}{T}ω0=T2π ,那么其傅里叶级数表示的形式为:

xT(t)=∑k=−∞+∞akejkω0t,(综合公式) x_T(t) = \sum_{k = -\infty}^{+\infty} a_k e^{jk\omega_0t} ,\text{(综合公式)} xT(t)=k=+akejkω0t,(综合公式)

其中 aka_kakxT(t)x_T(t)xT(t)kkk 次谐波的傅里叶级数的系数。

其中 ejkω0te^{jk\omega_0t}ejkω0t 称为旋转因子,而e−jkω0te^{-jk\omega_0t}ejkω0t 称为筛选因子,确定傅里叶系数的方法是两边同时乘以系数为 nnn 的筛选因子:

xT(t)e−jnω0t=∑k=−∞+∞akejkω0te−jnω0t x_T(t) e^{-jn\omega_0t} = \sum_{k = -\infty}^{+\infty} a_k e^{jk\omega_0t} e^{-jn\omega_0t} xT(t)ejnω0t=k=+akejkω0tejnω0t

同时在一个周期内做积分:

∫0TxT(t)e−jnω0tdt=∫0T∑k=−∞+∞akejkω0te−jnω0tdt=∑k=−∞+∞ak∫0Tej(k−n)ω0tdt \int_0^T x_T(t) e^{-jn\omega_0t} dt = \int_0^T \sum_{k = -\infty}^{+\infty} a_k e^{jk\omega_0t} e^{-jn\omega_0t} dt = \sum_{k = -\infty}^{+\infty} a_k \int_0^T e^{j(k-n)\omega_0t} dt 0TxT(t)ejnω0tdt=0Tk=+akejkω0tejnω0tdt=k=+ak0Tej(kn)ω0tdt

对于积分 ∫0Tej(k−n)ω0tdt\int_0^T e^{j(k-n)\omega_0t} dt0Tej(kn)ω0tdt 来说,当 k=nk = nk=n 的时候,积分值为 TTT ,否则等于 000 ,也就是:

∫0TxT(t)e−jnω0tdt=Tan \int_0^T x_T(t) e^{-jn\omega_0t} dt = Ta_n 0TxT(t)ejnω0tdt=Tan

即:

an=1T∫0TxT(t)e−jnω0tdt,(分析公式) a_n = \frac{1}{T} \int_0^T x_T(t) e^{-jn\omega_0t} dt,\text{(分析公式)} an=T10TxT(t)ejnω0tdt,(分析公式)

CT的傅里叶级数的性质

性质周期信号傅里叶系数
线性Ax(t)+By(t)Ax(t) +By(t)Ax(t)+By(t)Aak+BbkAa_k + Bb_kAak+Bbk
时移x(t−t0)x(t - t_0)x(tt0)ake−jkω0t0a_k e^{-jk\omega_0t_0}akejkω0t0
频移x(t)ejMω0tx(t)e^{jM\omega_0t}x(t)ejMω0tak−Ma_{k - M}akM
共轭x∗(t)x^*(t)x(t)a−k∗a^*_{-k}ak
时间翻转x(−t)x(-t)x(t)a−ka_{-k}ak
时域尺度变换x(αt)x(\alpha t)x(αt)ak(T=T/α)a_k(T=T/\alpha)ak(T=T/α)
周期卷积∫Tx(τ)y(t−τ)dτ\int_T x(\tau)y(t-\tau) d\tauTx(τ)y(tτ)dτTakbkTa_kb_kTakbk
相乘x(t)y(t)x(t)y(t)x(t)y(t)∑l=−∞+∞albk−l\sum_{l = -\infty}^{+\infty}a_l b_{k-l}l=+albkl
微分dx(t)dt\frac{dx(t)}{dt}dtdx(t)jkω0akjk\omega_0a_kjkω0ak
积分∫−∞tx(t)dt\int_{-\infty}^t x(t) dttx(t)dt1jkω0ak\frac{1}{jk\omega_0}a_kjkω01ak
实信号的共轭对称性x(t)x(t)x(t) 为实信号ak=−ak∗a_k = -a^*_kak=ak
实偶信号x(t)x(t)x(t) 为实偶信号aka_kak 为实偶函数
实奇信号x(t)x(t)x(t) 为实奇信号aka_kak 为纯虚奇函数
实信号的奇偶分解xe(t)=Ev{x(t)},xo(t)=Od{x(t)}x_e(t) = Ev\{x(t)\}, x_o(t) = Od\{x(t)\}xe(t)=Ev{x(t)},xo(t)=Od{x(t)} 并且 x(t)x(t)x(t) 为实信号ℜ{ak},jℑ{ak}\Re\{a_k\},j\Im\{a_k\}{ak},j{ak}
周期信号的帕瓦尔定理x(t)x(t)x(t)1T∫T∣x(t)∣2dt=∑k=−∞+∞∣ak∣2\frac{1}{T} \int_T |x(t)|^2 dt = \sum_{k = -\infty}^{+\infty} |a_k|^2T1Tx(t)2dt=k=+ak2

DT周期信号的傅里叶级数表示

定义周期信号的周期为 NNN ,有 x[n]=x[n+N]x[n] = x[n + N]x[n]=x[n+N] 。而 ω0=2πN\omega_0 = \frac{2\pi}{N}ω0=N2π 为基波频率。则 kkk 次谐波转子为:

ϕk[n]=ejkω0n \phi_k[n] = e^{jk\omega_0n} ϕk[n]=ejkω0n

并且,因为再离散时间中, kkknnn 均为整数的话,那么 ϕk[n]\phi_k[n]ϕk[n] 也为关于 kkk 周期为 NNN 的函数,也就是:

ϕ[n]=ϕk+rN[n] \phi[n] = \phi_{k + rN}[n] ϕ[n]=ϕk+rN[n]

这就是说,周期为 NNN 的离散时间信号,其傅里叶级数的系数只有 NNN 个,并且是一个周期为 NNN 的序列。因此DT周期信号的傅里叶级数表示为:

x[n]=∑k=<N>akejkωon,(综合公式) x[n] = \sum_{k = <N>} a_k e^{jk\omega_on},\text{(综合公式)} x[n]=k=<N>akejkωon,(综合公式)

对于连续的 NNN 个取值:

x[0]=∑k=<N>ak,x[1]=∑k=<N>akejkω0,…,x[N−1]=∑k=<N>akejkω0(N−1) x[0] = \sum_{k=<N>} a_k, x[1] = \sum_{k=<N>} a_k e^{jk\omega_0},\ldots, x[N-1] = \sum_{k=<N>} a_k e^{jk\omega_0(N-1)} x[0]=k=<N>ak,x[1]=k=<N>akejkω0,,x[N1]=k=<N>akejkω0(N1)

这些是 NNN 个线性独立的方程,可以直接解出 NNN 个系数的值。若想通过CT同样的方法,有以下的事实:

∑n=<N>ejkω0n=N(k=0,±N,±2N,…) \sum_{n = <N>} e^{jk\omega_0}n = N (k = 0,\pm N,\pm 2N,\ldots) n=<N>ejkω0n=N(k=0,±N,±2N,)

否则其他情况下等于 000

那么同样的,先乘以关于 rrr 的提取因子,然后在一个周期中求和:

∑n=<N>x[n]e−jrω0n=∑n=<N>∑k=<N>akej(k−r)ω0n=∑k=<N>ak∑n=<N>ej(k−r)ω0n \sum_{n = <N>} x[n]e^{-jr\omega_0 n} = \sum_{n = <N>} \sum_{k = <N>} a_k e^{j(k - r)\omega_0 n} = \sum_{k = <N>} a_k \sum_{n = <N>} e^{j(k - r)\omega_0 n} n=<N>x[n]ejrω0n=n=<N>k=<N>akej(kr)ω0n=k=<N>akn=<N>ej(kr)ω0n

对于和式 ∑n=<N>ej(k−r)ω0n\sum_{n = <N>} e^{j(k - r)\omega_0 n}n=<N>ej(kr)ω0nk=r+cNk = r + cNk=r+cN 的时候,即 k−rk-rkrNNN 的整数倍的时候,又因为在一个周期内只有一次 k−rk-rkrNNN 的整数倍,且对应 k=rk = rk=r ,于是右边就等于 NarNa_rNar ,因此:

ar=1N∑n=<N>x[n]e−jrω0n,(分析公式) a_r = \frac{1}{N} \sum_{n = <N>} x[n]e^{-jr\omega_0 n},\text{(分析公式)} ar=N1n=<N>x[n]ejrω0n,(分析公式)

DT的傅里叶级数的性质

性质周期信号傅里叶系数
线性Ax[n]+By[n]Ax[n] +By[n]Ax[n]+By[n]Aak+BbkAa_k + Bb_kAak+Bbk
时移x[n−n0]x[n - n_0]x[nn0]ake−jkω0n0a_k e^{-jk\omega_0n_0}akejkω0n0
频移x[n]ejMω0nx[n]e^{jM\omega_0n}x[n]ejMω0nak−Ma_{k - M}akM
共轭x∗[n]x^*[n]x[n]a−k∗a^*_{-k}ak
时间翻转x[−n]x[-n]x[n]a−ka_{-k}ak
时域尺度变换x(m)[n]=x[n/m]x_{(m)}[n]=x[n/m]x(m)[n]=x[n/m]1mak(N=Nm)\frac{1}{m}a_k(N=Nm)m1ak(N=Nm)
周期卷积∑r=<N>x[r]y[n−r]dτ\sum_{r = <N>} x[r]y[n-r] d\taur=<N>x[r]y[nr]dτNakbkNa_kb_kNakbk
相乘x[n]y[n]x[n]y[n]x[n]y[n]∑l=<N>albk−l\sum_{l = <N>} a_l b_{k-l}l=<N>albkl
一阶差分x[n]−x[n−1]x[n] - x[n - 1]x[n]x[n1](1−e−jkω0)ak(1 - e^{-jk\omega_0})a_k(1ejkω0)ak
求和∑k=−∞nx[k]\sum_{k = -\infty}^{n} x[k]k=nx[k]1(1−e−jkω0)ak\frac{1}{(1 - e^{-jk\omega_0})}a_k(1ejkω0)1ak
实信号的共轭对称性x[n]x[n]x[n] 为实信号ak=−ak∗a_k = -a^*_kak=ak
实偶信号x[n]x[n]x[n] 为实偶信号aka_kak 为实偶函数
实奇信号x[n]x[n]x[n] 为实奇信号aka_kak 为纯虚奇函数
实信号的奇偶分解xe[n]=Ev{x[n]},xo[n]=Od{x[n]}x_e[n] = Ev\{x[n]\}, x_o[n] = Od\{x[n]\}xe[n]=Ev{x[n]},xo[n]=Od{x[n]} 并且 x[n]x[n]x[n] 为实信号ℜ{ak},jℑ{ak}\Re\{a_k\},j\Im\{a_k\}{ak},j{ak}
周期信号的帕瓦尔定理x[n]x[n]x[n]1N∑n=<N>∣x[n]∣2=∑k=<N>∣ak∣2\frac{1}{N} \sum_{n = <N>} |x[n]|^2 = \sum_{k = <N>} |a_k|^2N1n=<N>x[n]2=k=<N>ak2

连续时间傅里叶变换

连续时间非周期的傅里叶表示

假设连续时间的周期信号为 xT(t)x_T(t)xT(t) ,对应的单个周期信号为 x(t)x(t)x(t) ,因为:

ak=1T∫−T/2T/2xT(t)e−jkω0tdt=1T∫−∞+∞x(t)e−jkω0tdt a_k = \frac{1}{T} \int_{-T/2}^{T/2} x_T(t) e^{-jk\omega_0t} dt = \frac{1}{T} \int_{-\infty}^{+\infty} x(t) e^{-jk\omega_0t} dt ak=T1T/2T/2xT(t)ejkω0tdt=T1+x(t)ejkω0tdt

定义 TakTa_kTak 的包络函数为:

X(jω)=∫−∞+∞x(t)e−jωtdt,(分析公式) X(j\omega) = \int_{-\infty}^{+\infty} x(t) e^{-j\omega t} dt,\text{(分析公式)} X()=+x(t)etdt,(分析公式)

此时 aka_kak 可以表示为:

ak=1TX(jkω0) a_k = \frac{1}{T} X(jk\omega_0) ak=T1X(jkω0)

并且:

xT(t)=∑k=−∞+∞1TX(jkω0)ejkω0t=12π∑k=−∞+∞X(jkω0)ejkω0tω0 x_T(t) = \sum_{k = -\infty}^{+\infty} \frac{1}{T}X(jk\omega_0) e^{jk\omega_0t} = \frac{1}{2\pi} \sum_{k = -\infty}^{+\infty} X(jk\omega_0) e^{jk\omega_0t} \omega_0 xT(t)=k=+T1X(jkω0)ejkω0t=2π1k=+X(jkω0)ejkω0tω0

接下来进行周期延拓,即 lim⁡T→∞xT(t)=x(t)\lim_{T \to \infty} x_T(t) = x(t)limTxT(t)=x(t) 有:

x(t)=12π∫−∞+∞X(jω)ejωtdω,(综合公式) x(t) = \frac{1}{2\pi} \int_{-\infty}^{+\infty} X(j\omega) e^{j\omega t} d\omega,\text{(综合公式)} x(t)=2π1+X()etdω,(综合公式)

上面两式称为傅里叶变换。

连续时间傅里叶变换的性质

性质周期信号傅里叶系数
线性Ax(t)+By(t)Ax(t) +By(t)Ax(t)+By(t)AX(jω)+BY(jω)AX(j\omega) + BY(j\omega)AX()+BY()
时移x(t−t0)x(t - t_0)x(tt0)X(jω)e−jωt0X(j\omega) e^{-j\omega t_0}X()et0
频移x(t)ejω0tx(t)e^{j\omega_0 t}x(t)ejω0tX(j(ω−ω0))X(j(\omega - \omega_0))X(j(ωω0))
共轭x∗(t)x^*(t)x(t)X∗(−jω)X^*(-j\omega)X()
时间翻转x(−t)x(-t)x(t)X(−jω)X(-j\omega)X()
时域尺度变换x(αt)x(\alpha t)x(αt)1∣a∣X(jωa)\frac{1}{|a|} X(\frac{j\omega}{a})a1X(a)
卷积x(t)∗y(t)x(t) \ast y(t)x(t)y(t)X(jω)Y(jω)X(j\omega)Y(j\omega)X()Y()
相乘x(t)y(t)x(t)y(t)x(t)y(t)12π∫−∞+∞X(jθ)Y(j(ω−θ))dθ\frac{1}{2\pi}\int_{-\infty}^{+\infty}X(j\theta) Y(j(\omega - \theta)) d\theta2π1+X(jθ)Y(j(ωθ))dθ
时域微分dx(t)dt\frac{dx(t)}{dt}dtdx(t)jωX(jω)j\omega X(j\omega)X()
时域积分∫−∞tx(t)dt\int_{-\infty}^t x(t) dttx(t)dt1jωX(jω)+πX(0)δ(ω)\frac{1}{j\omega} X(j\omega) + \pi X(0)\delta(\omega)1X()+πX(0)δ(ω)
频域微分tx(t)tx(t)tx(t)jddωX(jω)j\frac{d}{d\omega}X(j\omega)jdωdX()
实信号的共轭对称性x(t)x(t)x(t) 为实信号X(jω)=X∗(−jω)X(j\omega) = X^*(-j\omega)X()=X()
实偶信号x(t)x(t)x(t) 为实偶信号X(jω)X(j\omega)X() 为实偶函数
实奇信号x(t)x(t)x(t) 为实奇信号X(jω)X(j\omega)X() 为纯虚奇函数
实信号的奇偶分解xe(t)=Ev{x(t)},xo(t)=Od{x(t)}x_e(t) = Ev\{x(t)\}, x_o(t) = Od\{x(t)\}xe(t)=Ev{x(t)},xo(t)=Od{x(t)} 并且 x(t)x(t)x(t) 为实信号ℜ{X(jω)},jℑ{X(jω)}\Re\{X(j\omega)\},j\Im\{X(j\omega)\}{X()},j{X()}
周期信号的帕瓦尔定理x(t)x(t)x(t)∫−∞+∞∣x(t)∣2dt=12π∫−∞+∞∣X(jω)∣2dω\int_{-\infty}^{+\infty} |x(t)|^2 dt = \frac{1}{2\pi} \int_{-\infty}^{+\infty} |X(j\omega)|^2 d\omega+x(t)2dt=2π1+X()2dω

基本傅里叶变换对

  1. 周期信号的傅里叶变换

考虑这样一个信号的傅里叶变换为:

X(jω)=2πδ(ω−ω0) X(j\omega) = 2\pi\delta(\omega - \omega_0) X()=2πδ(ωω0)

其对应的时域信号为:

x(t)=12π∫−∞+∞2πδ(ω−ω0)ejωtdω=ejω0t x(t) = \frac{1}{2 \pi} \int_{-\infty}^{+\infty} 2\pi \delta(\omega - \omega_0) e^{j\omega t} d\omega = e^{j\omega_0 t} x(t)=2π1+2πδ(ωω0)etdω=ejω0t

那么通过周期信号的傅里叶系数关系为:

X(jω)=∑k=−∞+∞2πakδ(ω−kω0) X(j\omega) = \sum_{k = -\infty}^{+\infty} 2 \pi a_k \delta(\omega - k\omega_0) X()=k=+2πakδ(ωkω0)

其对应的时域信号为:

x(t)=∑k=−∞+∞akejkω0t x(t) = \sum_{k = -\infty}^{+\infty} a_k e^{jk\omega_0 t} x(t)=k=+akejkω0t

  1. 单边衰减信号

考虑信号 x(t)=e−atu(t),a>0x(t) = e^{-at}u(t), a> 0x(t)=eatu(t),a>0 。其对应的傅里叶变换为:

X(jω)=∫0∞e−ate−jωtdt=−1a+jωe−(a+jω)t∣0∞=1a+jω(a>0) X(j\omega) = \int_0^{\infty} e^{-at} e^{-j\omega t} dt = -\frac{1}{a + j\omega} e^{-(a + j\omega)t} |_0^\infty = \frac{1}{a + j\omega} (a>0) X()=0eatetdt=a+1e(a+)t0=a+1(a>0)

  1. 双边衰减信号

考虑信号 x(t)=e−a∣t∣,a>0x(t) = e^{-a|t|}, a> 0x(t)=eat,a>0 。其对应的傅里叶变换为:

X(jω)=1a−jω+1a+jω=2aa2+ω2 X(j\omega) = \frac{1}{a - j\omega} + \frac{1}{a + j\omega} = \frac{2a}{a^2 + \omega^2} X()=a1+a+1=a2+ω22a

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

单位冲激函数 x(t)=δ(t)x(t) = \delta(t)x(t)=δ(t) 的傅里叶变换为:

X(jω)=∫−∞+∞δ(t)e−jωtdt=1 X(j\omega) = \int_{-\infty}^{+\infty} \delta(t) e^{-j\omega t} dt = 1 X()=+δ(t)etdt=1

  1. 矩形脉冲信号的傅里叶变换

考虑一个矩形脉冲信号:

x(t)=1,∣t∣<T1 x(t) = 1 ,|t| < T_1 x(t)=1,t<T1

x(t)=0,∣t∣>T1 x(t) = 0 ,|t| > T_1 x(t)=0,t>T1

其傅里叶变换为:

X(jω)=∫−T1T1ejωtdt=2sin⁡ωT1ω X(j\omega) = \int_{-T_1}^{T_1} e^{j\omega t} dt = 2\frac{\sin \omega T_1}{\omega} X()=T1T1etdt=2ωsinωT1

  1. 傅里叶变换的对偶性

假设信号 x(t)x(t)x(t) 存在傅里叶变换对:

X(jω)=∫−∞+∞x(t)e−jωtdt X(j\omega) = \int_{-\infty}^{+\infty} x(t) e^{-j\omega t} dt X()=+x(t)etdt

x(t)=12π∫−∞+∞X(jω)ejωtdω x(t) = \frac{1}{2 \pi} \int_{-\infty}^{+\infty} X(j\omega) e^{j\omega t} d\omega x(t)=2π1+X()etdω

那么信号 X(jt)X(jt)X(jt) 存在傅里叶变换对:

∫−∞+∞X(jt)e−jωtdt=2πx(−ω) \int_{-\infty}^{+\infty} X(jt) e^{-j\omega t} dt = 2\pi x(-\omega) +X(jt)etdt=2πx(ω)

12π∫−∞+∞x(ω)ejωtdω=12πX(−jt) \frac{1}{2 \pi} \int_{-\infty}^{+\infty} x(\omega) e^{j\omega t} d\omega = \frac{1}{2 \pi} X(-jt) 2π1+x(ω)etdω=2π1X(jt)

  1. 典型的傅里叶变换对
信号傅里叶变换傅里叶级数系数(若是周期信号)
∑k=−∞+∞akejkω0t\sum_{k = -\infty}^{+\infty} a_k e^{jk\omega_0 t}k=+akejkω0t∑k=−∞+∞2πakδ(ω−kω0)\sum_{k = -\infty}^{+\infty} 2 \pi a_k \delta(\omega - k\omega_0)k=+2πakδ(ωkω0)aka_kak
ejkω0te^{jk\omega_0 t}ejkω0t2πδ(ω−kω0)2\pi \delta(\omega-k\omega_0)2πδ(ωkω0)a1=1,ak=0a_1 = 1,a_k = 0a1=1,ak=0
cos⁡ω0t\cos \omega_0 tcosω0tπ[δ(ω−ω0)+δ(ω+ω0)]\pi[\delta(\omega - \omega_0) + \delta(\omega + \omega_0)]π[δ(ωω0)+δ(ω+ω0)]a1=a−1=12,ak=0a_1 = a_{-1} = \frac{1}{2},a_k = 0a1=a1=21,ak=0
sin⁡ω0t\sin \omega_0 tsinω0tπj[δ(ω−ω0)−δ(ω+ω0)]\frac{\pi}{j}[\delta(\omega - \omega_0) - \delta(\omega + \omega_0)]jπ[δ(ωω0)δ(ω+ω0)]a1=−a−1=12j,ak=0a_1 = -a_{-1} = \frac{1}{2j},a_k = 0a1=a1=2j1,ak=0
x(t)=1x(t) = 1x(t)=12πδ(ω)2\pi\delta(\omega)2πδ(ω)a0=1,ak=0a_0 = 1,a_k = 0a0=1,ak=0
周期方波x(t)=1,∣t∣<T1x(t) = 1 ,|t| < T_1x(t)=1,t<T1∑k=−∞+∞2sin⁡kω0T1kδ(ω−kω0)\sum_{k = -\infty}^{+\infty}\frac{2\sin k\omega_0 T_1}{k}\delta(\omega - k\omega_0)k=+k2sinkω0T1δ(ωkω0)sin⁡kω0T1kπ\frac{\sin k\omega_0T_1}{k\pi}sinkω0T1
周期冲激串∑n=−∞+∞δ(t−nT)\sum_{n = -\infty}^{+\infty}\delta(t - nT)n=+δ(tnT)2πT∑k=−∞+∞δ(ω−2πkT)\frac{2 \pi}{T}\sum_{k = -\infty}^{+\infty}\delta(\omega - \frac{2\pi k}{T})T2πk=+δ(ωT2πk)ak=1Ta_k = \frac{1}{T}ak=T1
矩形阶跃函数x(t)=1,∣t∣<T1x(t) = 1,|t| < T_1x(t)=1,t<T12sin⁡ωT1ω\frac{2 \sin \omega T_1}{\omega}ω2sinωT1-
sin⁡Wtπt\frac{\sin Wt}{\pi t}πtsinWtX(jω)=1,∣ω∣<WX(j\omega) = 1, |\omega| <WX()=1,ω<W-
δ(t)\delta(t)δ(t)111-
u(t)u(t)u(t)1jω+πδ(ω)\frac{1}{j \omega} + \pi \delta(\omega)1+πδ(ω)-
δ(t−t0)\delta(t - t_0)δ(tt0)ejωt0e^{j\omega t_0}et0-
e−atu(t),ℜa>0e^{-at}u(t),\Re{a} > 0eatu(t),a>01a+jω\frac{1}{a+j\omega}a+1-
te−atu(t),ℜa>0te^{-at}u(t),\Re{a} > 0teatu(t),a>01(a+jω)2\frac{1}{(a+j\omega)^2}(a+)21-
tn−1(n−1)!e−atu(t),ℜa>0\frac{t^{n-1}}{(n-1)!}e^{-at}u(t),\Re{a} > 0(n1)!tn1eatu(t),a>01(a+jω)n\frac{1}{(a+j\omega)^n}(a+)n1-

傅里叶变换与线性常系数微分方程表述的系统

对于线性常系数微分方程描述的系统:

∑k=0Nakdky(t)dtk=∑k=0Mbkdkx(t)dtk \sum_{k=0}^{N} a_k \frac{d^k y(t)}{dt^k} = \sum_{k=0}^{M} b_k \frac{d^k x(t)}{dt^k} k=0Nakdtkdky(t)=k=0Mbkdtkdkx(t)

有卷积关系:

y(t)=h(t)∗x(t) y(t) = h(t) \ast x(t) y(t)=h(t)x(t)

根据傅里叶变换的性质有:

Y(jω)=H(jω)X(jω) Y(j\omega) = H(j \omega) X(j \omega) Y()=H()X()

我们称函数 H(jω)H(j \omega)H()系统频响函数

对微分方程两边做傅里叶变换:

F{∑k=0Nakdky(t)dtk}=F{∑k=0Mbkdkx(t)dtk} \mathscr{F} \{\sum_{k=0}^{N} a_k \frac{d^k y(t)}{dt^k}\} = \mathscr{F} \{\sum_{k=0}^{M} b_k \frac{d^k x(t)}{dt^k}\} F{k=0Nakdtkdky(t)}=F{k=0Mbkdtkdkx(t)}

根据傅里叶变换的微分性质:

∑k=0Nak(jω)kY(jω)=∑k=0Mbk(jω)kX(jω) \sum_{k=0}^{N} a_k (j\omega)^k Y(j\omega) = \sum_{k=0}^{M} b_k (j\omega)^k X(j\omega) k=0Nak()kY()=k=0Mbk()kX()

可得:

H(jω)=Y(jω)X(jω)=∑k=0Mbk(jω)k∑k=0Nak(jω)k H(j\omega) = \frac{Y(j\omega)}{X(j\omega)} = \frac{\sum_{k=0}^{M} b_k (j\omega)^k}{\sum_{k=0}^{N} a_k (j\omega)^k} H()=X()Y()=k=0Nak()kk=0Mbk()k

离散时间傅里叶变换

离散时间非周期的傅里叶表示

假设有离散时间的信号 xT[n]x_T[n]xT[n] ,那么其傅里叶系数表示为:

ak=1N∑n=<N>x[n]e−jkω0n a_k = \frac{1}{N} \sum_{n = <N>} x[n] e^{-jk\omega_0n} ak=N1n=<N>x[n]ejkω0n

其对应的非周期的信号表示在一个周期内 [−N1,N2][-N_1,N_2][N1,N2]x[n]x[n]x[n] ,那么:

ak=1N∑n=−N1N2x[n]e−jkω0n=1N∑n=−∞+∞x[n]e−jkω0n a_k = \frac{1}{N} \sum_{n = -N_1}^{N_2} x[n] e^{-jk\omega_0n} = \frac{1}{N} \sum_{n = -\infty}^{+\infty} x[n] e^{-jk\omega_0n} ak=N1n=N1N2x[n]ejkω0n=N1n=+x[n]ejkω0n

我们定义 akNa_kNakN 的包络:

X(ejω)=∑n=−∞+∞x[n]e−jωn,分析公式 X(e^{j\omega}) = \sum_{n = -\infty}^{+\infty} x[n] e^{-j\omega n},\text{分析公式} X(e)=n=+x[n]ejωn,分析公式

为离散时间傅里叶变换。

那么有:

ak=1NX(ejkω0) a_k = \frac{1}{N} X(e^{jk\omega_0}) ak=N1X(ejkω0)

则重新带回原始得到:

xT[n]=∑k=<N>1NX(ejkω0)ejkω0n=12π∑k=<N>X(ejkω0)ejkω0nω0 x_T[n] = \sum_{k = <N>} \frac{1}{N} X(e^{jk\omega_0}) e^{jk\omega_0n} = \frac{1}{2\pi} \sum_{k = <N>}X(e^{jk\omega_0}) e^{jk\omega_0n} \omega_0 xT[n]=k=<N>N1X(ejkω0)ejkω0n=2π1k=<N>X(ejkω0)ejkω0nω0

N→∞N \to \inftyN 的时候,其在一个周期上的求和就变成了在 2π2\pi2π 内的一个积分,即:

x[n]=12π∫2πX(ejω)ejωndω,综合公式 x[n] = \frac{1}{2\pi} \int_{2 \pi} X(e^{j\omega}) e^{j\omega n} d\omega, \text{综合公式} x[n]=2π12πX(e)ejωndω,综合公式

因为 X(ejω)ejωnX(e^{j\omega}) e^{j\omega n}X(e)ejωn 关于 ω\omegaω 周期为 2π2\pi2π 的函数,这说明我们的积分区间可以随意取一个连续的 2π2\pi2π 周期。

同时,对于非周期信号的傅里叶变换频谱公式 X(ejω)X(e^{j\omega})X(e) 来说,也是一个关于 ω\omegaω 周期为 2π2\pi2π 的函数,这说明 非周期的离散信号的傅里叶频谱是一个周期为 2π2\pi2π 的函数

离散时间非周期的傅里叶的性质

性质周期信号傅里叶系数
线性Ax[n]+By[n]Ax[n] +By[n]Ax[n]+By[n]AX(ejω)+BY(ejω)AX(e^{j\omega}) + BY(e^{j\omega})AX(e)+BY(e)
时移x[n−n0]x[n - n_0]x[nn0]X(ejω)e−jωn0X(e^{j\omega}) e^{-j\omega n_0}X(e)en0
频移x[n]ejω0nx[n]e^{j\omega_0 n}x[n]ejω0nX(ej(ω−ω0))X(e^{j(\omega - \omega_0)})X(ej(ωω0))
共轭x∗[n]x^*[n]x[n]X∗(e−jω)X^*(e^{-j\omega})X(e)
时间翻转x[−n]x[-n]x[n]X(e−jω)X(e^{-j\omega})X(e)
时域尺度变换x(k)[n]=x[n/k]x_{(k)}[n] = x[n/k]x(k)[n]=x[n/k]X(ejkω)X(e^{jk\omega})X(ejkω)
卷积x[n]∗y[n]x[n] \ast y[n]x[n]y[n]X(ejω)Y(ejω)X(e^{j\omega})Y(e^{j\omega})X(e)Y(e)
相乘x[t]y[t]x[t]y[t]x[t]y[t]12π∫2πX(ejθ)Y(ej(ω−θ))dθ\frac{1}{2\pi}\int_{2\pi} X(e^{j\theta}) Y(e^{j(\omega - \theta)}) d\theta2π12πX(ejθ)Y(ej(ωθ))dθ
时域差分x[n]−x[n−1]x[n] - x[n-1]x[n]x[n1](1−e−jω)X(ejω)(1 - e^{-j\omega})X(e^{j\omega})(1e)X(e)
时域累加∑k=−∞nx[k]\sum_{k=-\infty}^n x[k]k=nx[k]11−e−jωX(ejω)+πX(ej0)∑k=−∞+∞δ(ω−2πk)\frac{1}{1 - e^{-j\omega}}X(e^{j\omega}) +\pi X(e^{j0})\sum_{k=-\infty}^{+\infty}\delta(\omega - 2\pi k)1e1X(e)+πX(ej0)k=+δ(ω2πk)
频域微分nx[n]nx[n]nx[n]jddωX(ejω)j\frac{d}{d\omega}X(e^{j\omega})jdωdX(e)
实信号的共轭对称性x[n]x[n]x[n] 为实信号X(ejω)=X∗(e−jω)X(e^{j\omega}) = X^*(e^{-j\omega})X(e)=X(e)
实偶信号x[n]x[n]x[n] 为实偶信号X(ejω)X(e^{j\omega})X(e) 为实偶函数
实奇信号x[n]x[n]x[n] 为实奇信号X(ejω)X(e^{j\omega})X(e) 为纯虚奇函数
实信号的奇偶分解xe[n]=Ev{x[n]},xo[n]=Od{x[n]}x_e[n] = Ev\{x[n]\}, x_o[n] = Od\{x[n]\}xe[n]=Ev{x[n]},xo[n]=Od{x[n]} 并且 x[n]x[n]x[n] 为实信号ℜ{X(ejω)},jℑ{X(ejω)}\Re\{X(e^{j\omega})\},j\Im\{X(e^{j\omega})\}{X(e)},j{X(e)}
周期信号的帕瓦尔定理x[n]x[n]x[n]∑n=−∞+∞∣x[n]∣2=12π∫2π∣X(ejω)∣2dω\sum_{n = -\infty}^{+\infty} |x[n]|^2 = \frac{1}{2\pi} \int_{2\pi} |X(e^{j\omega})|^2 d\omegan=+x[n]2=2π12πX(e)2dω

基本离散时间信号的傅里叶变换对

  1. 周期信号

考虑信号 x[n]=ejω0nx[n] = e^{j\omega_0 n}x[n]=ejω0n 的傅里变换为:

X(ejω)=∑l=−∞+∞2πδ(ω−ω0−2πl) X(e^{j\omega}) = \sum_{l = -\infty}^{+\infty} 2\pi \delta(\omega - \omega_0 - 2\pi l) X(e)=l=+2πδ(ωω02πl)

现在考虑一个周期为 NNN 的周期序列,其傅里叶级数为:

x[n]=∑k=<N>akejkω0n x[n] = \sum_{k = <N>} a_k e^{jk\omega_0 n} x[n]=k=<N>akejkω0n

那么他的傅里叶变换就是:

X(ejω)=∑k=−∞+∞2πakδ(ω−2πkN) X(e^{j\omega}) = \sum_{k = -\infty}^{+\infty} 2 \pi a_k \delta(\omega - \frac{2\pi k}{N}) X(e)=k=+2πakδ(ωN2πk)

  1. 离散时间傅里叶系数的对偶性

若两个周期为 NNN 的序列 f[n]f[n]f[n]g[n]g[n]g[n] ,若 f[k]f[k]f[k]g[n]g[n]g[n] 的离散时间傅里叶系数,即

g[n]↔f[k] g[n] \leftrightarrow f[k] g[n]f[k]

那么:

f[n]↔1Ng[−k] f[n] \leftrightarrow \frac{1}{N}g[-k] f[n]N1g[k]

  1. 离散时间傅里叶变换和连续时间傅里叶系数的对偶性

考虑离散时间傅里叶变换:

X(ejω)=∑n=−∞+∞x[n]e−jωn X(e^{j\omega}) = \sum_{n = -\infty}^{+\infty} x[n] e^{-j\omega n} X(e)=n=+x[n]ejωn

与连续时间傅里叶系数:

x(t)=∑k=−∞+∞ake−jkω0t x(t) = \sum_{k = -\infty}^{+\infty} a_k e^{-jk\omega_0 t} x(t)=k=+akejkω0t

这相当于将 x[n]x[n]x[n] 看成是一个连续时间函数的傅里叶系数,同理,我们可以将 x(t)x(t)x(t) 看做是一个离散时间函数的傅里叶变换。

  1. 常用基本离散时间信号的傅里叶变换对表
信号傅里叶变换傅里叶系数(若为周期的)
∑k=<N>akejkω0n,N=2πω0\sum_{k = <N>}a_k e^{jk\omega_0 n},N = \frac{2\pi}{\omega_0}k=<N>akejkω0n,N=ω02π2π∑k=−∞+∞akδ(ω−kω0)2\pi \sum_{k = -\infty}^{+ \infty}a_k\delta(\omega - k\omega_0)2πk=+akδ(ωkω0)aka_kak
ejω0n,N=2πω0e^{j\omega_0 n},N = \frac{2\pi}{\omega_0}ejω0n,N=ω02π2π∑l=−∞+∞δ(ω−ω0−2πl)2 \pi \sum_{l = -\infty}^{+\infty}\delta(\omega - \omega_0 - 2 \pi l)2πl=+δ(ωω02πl)ak=1,k=1,1±N,…a_k = 1,k = 1,1 \pm N,\ldotsak=1,k=1,1±N,
cos⁡ω0n,N=2πω0\cos{\omega_0 n},N = \frac{2\pi}{\omega_0}cosω0n,N=ω02ππ∑l=−∞+∞(δ(ω−ω0−2πl)+δ(ω+ω−2πl))\pi \sum_{l = -\infty}^{+\infty}(\delta(\omega - \omega_0 - 2 \pi l) + \delta(\omega + \omega - 2\pi l))πl=+(δ(ωω02πl)+δ(ω+ω2πl))ak=12,k=1,1±N,…a_k = \frac{1}{2},k = 1,1 \pm N,\ldotsak=21,k=1,1±N,
sin⁡ω0n,N=2πω0\sin{\omega_0 n},N = \frac{2\pi}{\omega_0}sinω0n,N=ω02ππj∑l=−∞+∞(δ(ω−ω0−2πl)−δ(ω+ω−2πl))\frac{\pi}{j} \sum_{l = -\infty}^{+\infty}(\delta(\omega - \omega_0 - 2 \pi l) - \delta(\omega + \omega - 2\pi l))jπl=+(δ(ωω02πl)δ(ω+ω2πl))ak=12j,k=1,1±N,…;ak=−12j,k=−1,−1±N,…a_k = \frac{1}{2j},k = 1,1 \pm N,\ldots;a_k=-\frac{1}{2j},k=-1,-1\pm N,\ldotsak=2j1,k=1,1±N,;ak=2j1,k=1,1±N,
x[n]=1x[n] = 1x[n]=12π∑l=−∞+∞δ(ω−2πl)2 \pi \sum_{l = -\infty}^{+ \infty}\delta(\omega - 2\pi l)2πl=+δ(ω2πl)ak=1,k=0,±N,…a_k = 1,k = 0,\pm N,\ldotsak=1,k=0,±N,
周期方波x[n]=1,∣n∣≤N1,x[n+N]=x[n]x[n] = 1,|n| \le N_1,x[n + N] = x[n]x[n]=1,nN1,x[n+N]=x[n]2π∑k=−∞+∞akδ(ω−2πN)2\pi \sum_{k = -\infty}^{+\infty}a_k \delta(\omega - \frac{2\pi}{N})2πk=+akδ(ωN2π)ak=sin⁡(2πk/N)(N1+12)Nsin⁡2πk/2N,k≠0,±N,…;ak=2N1+1N,k=0,±N,…a_k = \frac{\sin{(2\pi k / N)(N_1 + \frac{1}{2})}}{N \sin{2\pi k/ 2N}},k \neq 0,\pm N,\ldots;a_k = \frac{2N_1 + 1}{N},k=0,\pm N,\ldotsak=Nsin2πk/2Nsin(2πk/N)(N1+21),k=0,±N,;ak=N2N1+1,k=0,±N,
∑k=−∞+∞δ[n−kN]\sum_{k = -\infty}^{+\infty} \delta[n - kN]k=+δ[nkN]2πN∑k=−∞+∞δ(ω−2πkN)\frac{2 \pi}{N}\sum_{k = -\infty}^{+\infty}\delta(\omega - \frac{2 \pi k}{N})N2πk=+δ(ωN2πk)ak=1Na_k = \frac{1}{N}ak=N1
anu[n]a^nu[n]anu[n],|a| < 1$11−ae−jω\frac{1}{1 - ae^{-j\omega}}1ae1-
x[n]=1,∣n∣≤N1x[n] = 1, |n| \le N_1x[n]=1,nN1sin⁡ω(N1+12)sin⁡ω/2\frac{\sin{\omega(N_1 + \frac{1}{2})}}{\sin{\omega/2}}sinω/2sinω(N1+21)-
sin⁡Wnπn,0<W<π\frac{\sin{Wn}}{\pi n},0 < W < \piπnsinWn,0<W<πX(ω)=1,0≤∣ω∣≤WX(\omega) = 1,0 \le |\omega| \le WX(ω)=1,0ωW 并且 X(ω)X(\omega)X(ω) 是周期的为 2π2\pi2π-
δ[n]\delta[n]δ[n]111-
u[n]u[n]u[n]11−e−jω+∑k=−∞+∞πδ(ω−wπk)\frac{1}{1 - e^{-j\omega}} + \sum_{k = -\infty}^{+ \infty}\pi \delta(\omega - w\pi k)1e1+k=+πδ(ωwπk)-
δ[n−n0]\delta[n - n_0]δ[nn0]e−jωn0e^{-j\omega n_0}en0-
(n+1)anu[n],∣a∣<1(n+1)a^n u[n], |a| < 1(n+1)anu[n],a<11(1−ae−jω)2\frac{1}{(1 - ae^{-j\omega})^2}(1ae)21-
(n+r−1)!n!(r−1)!anu[n],∣a∣<1\frac{(n + r - 1)!}{n! (r - 1)!} a^n u[n], |a| < 1n!(r1)!(n+r1)!anu[n],a<11(1−ae−jω)r\frac{1}{(1 - ae^{-j\omega})^r}(1ae)r1-

由线性常系数差分方程表征的系统

由下面表示的线性常系数差分方程表征的系统:

∑k=0Naky[n−k]=∑k=0Mbkx[n−k] \sum_{k = 0}^{N} a_k y[n - k] = \sum_{k = 0}^M b_k x[n - k] k=0Naky[nk]=k=0Mbkx[nk]

通过傅里叶变换时移和线性的性质可以表示为频域上:

∑k=0Nake−jkωY(ejω)=∑k=0Mbke−jkωX(ejω) \sum_{k = 0}^N a_k e^{-jk\omega}Y(e^{j\omega}) = \sum_{k = 0}^M b_k e^{-jk\omega}X(e^{j\omega}) k=0NakejkωY(e)=k=0MbkejkωX(e)

或等效为:

H(ejω)=Y(ejω)X(ejω)=∑k=0Mbke−jkω∑k=0Nake−jkω H(e^{j\omega}) = \frac{Y(e^{j\omega})}{X(e^{j\omega})} = \frac{\sum_{k = 0}^M b_k e^{-jk\omega}}{\sum_{k = 0}^N a_k e^{-jk\omega}} H(e)=X(e)Y(e)=k=0Nakejkωk=0Mbkejkω

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值