信号与系统复习笔记——傅里叶变换
周期信号的傅里叶级数表示
特征函数
假设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)=est∗h(t)=∫−∞+∞es(t−τ)h(τ)dτ=est∫−∞+∞e−sτ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)=∫−∞+∞e−sτ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_kak 为 xT(t)x_T(t)xT(t) 的 kkk 次谐波的傅里叶级数的系数。
其中 ejkω0te^{jk\omega_0t}ejkω0t 称为旋转因子,而e−jkω0te^{-jk\omega_0t}e−jkω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)e−jnω0t=k=−∞∑+∞akejkω0te−jnω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)e−jnω0tdt=∫0Tk=−∞∑+∞akejkω0te−jnω0tdt=k=−∞∑+∞ak∫0Tej(k−n)ω0tdt
对于积分 ∫0Tej(k−n)ω0tdt\int_0^T e^{j(k-n)\omega_0t} dt∫0Tej(k−n)ω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)e−jnω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=T1∫0TxT(t)e−jnω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(t−t0) | ake−jkω0t0a_k e^{-jk\omega_0t_0}ake−jkω0t0 |
频移 | x(t)ejMω0tx(t)e^{jM\omega_0t}x(t)ejMω0t | ak−Ma_{k - M}ak−M |
共轭 | x∗(t)x^*(t)x∗(t) | a−k∗a^*_{-k}a−k∗ |
时间翻转 | x(−t)x(-t)x(−t) | a−ka_{-k}a−k |
时域尺度变换 | 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\tau∫Tx(τ)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=−∞+∞albk−l |
微分 | dx(t)dt\frac{dx(t)}{dt}dtdx(t) | jkω0akjk\omega_0a_kjkω0ak |
积分 | ∫−∞tx(t)dt\int_{-\infty}^t x(t) dt∫−∞tx(t)dt | 1jkω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|^2T1∫T∣x(t)∣2dt=∑k=−∞+∞∣ak∣2 |
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
并且,因为再离散时间中, kkk 和 nnn 均为整数的话,那么 ϕ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[N−1]=k=<N>∑akejkω0(N−1)
这些是 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]e−jrω0n=n=<N>∑k=<N>∑akej(k−r)ω0n=k=<N>∑akn=<N>∑ej(k−r)ω0n
对于和式 ∑n=<N>ej(k−r)ω0n\sum_{n = <N>} e^{j(k - r)\omega_0 n}∑n=<N>ej(k−r)ω0n 当 k=r+cNk = r + cNk=r+cN 的时候,即 k−rk-rk−r 是 NNN 的整数倍的时候,又因为在一个周期内只有一次 k−rk-rk−r 是 NNN 的整数倍,且对应 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]e−jrω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[n−n0] | ake−jkω0n0a_k e^{-jk\omega_0n_0}ake−jkω0n0 |
频移 | x[n]ejMω0nx[n]e^{jM\omega_0n}x[n]ejMω0n | ak−Ma_{k - M}ak−M |
共轭 | x∗[n]x^*[n]x∗[n] | a−k∗a^*_{-k}a−k∗ |
时间翻转 | x[−n]x[-n]x[−n] | a−ka_{-k}a−k |
时域尺度变换 | 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\tau∑r=<N>x[r]y[n−r]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>albk−l |
一阶差分 | x[n]−x[n−1]x[n] - x[n - 1]x[n]−x[n−1] | (1−e−jkω0)ak(1 - e^{-jk\omega_0})a_k(1−e−jkω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(1−e−jkω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|^2N1∑n=<N>∣x[n]∣2=∑k=<N>∣ak∣2 |
连续时间傅里叶变换
连续时间非周期的傅里叶表示
假设连续时间的周期信号为 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=T1∫−T/2T/2xT(t)e−jkω0tdt=T1∫−∞+∞x(t)e−jkω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(jω)=∫−∞+∞x(t)e−jωtdt,(分析公式)
此时 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
接下来进行周期延拓,即 limT→∞xT(t)=x(t)\lim_{T \to \infty} x_T(t) = x(t)limT→∞xT(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(jω)ejωtdω,(综合公式)
上面两式称为傅里叶变换。
连续时间傅里叶变换的性质
性质 | 周期信号 | 傅里叶系数 |
---|---|---|
线性 | Ax(t)+By(t)Ax(t) +By(t)Ax(t)+By(t) | AX(jω)+BY(jω)AX(j\omega) + BY(j\omega)AX(jω)+BY(jω) |
时移 | x(t−t0)x(t - t_0)x(t−t0) | X(jω)e−jωt0X(j\omega) e^{-j\omega t_0}X(jω)e−jωt0 |
频移 | x(t)ejω0tx(t)e^{j\omega_0 t}x(t)ejω0t | X(j(ω−ω0))X(j(\omega - \omega_0))X(j(ω−ω0)) |
共轭 | x∗(t)x^*(t)x∗(t) | X∗(−jω)X^*(-j\omega)X∗(−jω) |
时间翻转 | x(−t)x(-t)x(−t) | X(−jω)X(-j\omega)X(−jω) |
时域尺度变换 | x(αt)x(\alpha t)x(αt) | 1∣a∣X(jωa)\frac{1}{|a|} X(\frac{j\omega}{a})∣a∣1X(ajω) |
卷积 | x(t)∗y(t)x(t) \ast y(t)x(t)∗y(t) | X(jω)Y(jω)X(j\omega)Y(j\omega)X(jω)Y(jω) |
相乘 | 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)jωX(jω) |
时域积分 | ∫−∞tx(t)dt\int_{-\infty}^t x(t) dt∫−∞tx(t)dt | 1jωX(jω)+πX(0)δ(ω)\frac{1}{j\omega} X(j\omega) + \pi X(0)\delta(\omega)jω1X(jω)+πX(0)δ(ω) |
频域微分 | tx(t)tx(t)tx(t) | jddωX(jω)j\frac{d}{d\omega}X(j\omega)jdωdX(jω) |
实信号的共轭对称性 | x(t)x(t)x(t) 为实信号 | X(jω)=X∗(−jω)X(j\omega) = X^*(-j\omega)X(jω)=X∗(−jω) |
实偶信号 | x(t)x(t)x(t) 为实偶信号 | X(jω)X(j\omega)X(jω) 为实偶函数 |
实奇信号 | x(t)x(t)x(t) 为实奇信号 | X(jω)X(j\omega)X(jω) 为纯虚奇函数 |
实信号的奇偶分解 | 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ω)},jℑ{X(jω)} |
周期信号的帕瓦尔定理 | 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(jω)∣2dω |
基本傅里叶变换对
- 周期信号的傅里叶变换
考虑这样一个信号的傅里叶变换为:
X(jω)=2πδ(ω−ω0) X(j\omega) = 2\pi\delta(\omega - \omega_0) X(jω)=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)ejωtdω=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(jω)=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
- 单边衰减信号
考虑信号 x(t)=e−atu(t),a>0x(t) = e^{-at}u(t), a> 0x(t)=e−atu(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(jω)=∫0∞e−ate−jωtdt=−a+jω1e−(a+jω)t∣0∞=a+jω1(a>0)
- 双边衰减信号
考虑信号 x(t)=e−a∣t∣,a>0x(t) = e^{-a|t|}, a> 0x(t)=e−a∣t∣,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(jω)=a−jω1+a+jω1=a2+ω22a
- 单位冲激函数的傅里叶变换
单位冲激函数 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(jω)=∫−∞+∞δ(t)e−jωtdt=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(jω)=∫−T1T1ejωtdt=2ωsinωT1
- 傅里叶变换的对偶性
假设信号 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(jω)=∫−∞+∞x(t)e−jωtdt
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(jω)ejωtdω
那么信号 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)e−jωtdt=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(ω)ejωtdω=2π1X(−jt)
- 典型的傅里叶变换对
信号 | 傅里叶变换 | 傅里叶级数系数(若是周期信号) |
---|---|---|
∑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ω0t | 2πδ(ω−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=a−1=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=−a−1=2j1,ak=0 |
x(t)=1x(t) = 1x(t)=1 | 2πδ(ω)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=−∞+∞2sinkω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) | sinkω0T1kπ\frac{\sin k\omega_0T_1}{k\pi}kπsinkω0T1 |
周期冲激串∑n=−∞+∞δ(t−nT)\sum_{n = -\infty}^{+\infty}\delta(t - nT)∑n=−∞+∞δ(t−nT) | 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∣<T1 | 2sinωT1ω\frac{2 \sin \omega T_1}{\omega}ω2sinωT1 | - |
sinWtπt\frac{\sin Wt}{\pi t}πtsinWt | X(jω)=1,∣ω∣<WX(j\omega) = 1, |\omega| <WX(jω)=1,∣ω∣<W | - |
δ(t)\delta(t)δ(t) | 111 | - |
u(t)u(t)u(t) | 1jω+πδ(ω)\frac{1}{j \omega} + \pi \delta(\omega)jω1+πδ(ω) | - |
δ(t−t0)\delta(t - t_0)δ(t−t0) | ejωt0e^{j\omega t_0}ejωt0 | - |
e−atu(t),ℜa>0e^{-at}u(t),\Re{a} > 0e−atu(t),ℜa>0 | 1a+jω\frac{1}{a+j\omega}a+jω1 | - |
te−atu(t),ℜa>0te^{-at}u(t),\Re{a} > 0te−atu(t),ℜa>0 | 1(a+jω)2\frac{1}{(a+j\omega)^2}(a+jω)21 | - |
tn−1(n−1)!e−atu(t),ℜa>0\frac{t^{n-1}}{(n-1)!}e^{-at}u(t),\Re{a} > 0(n−1)!tn−1e−atu(t),ℜa>0 | 1(a+jω)n\frac{1}{(a+j\omega)^n}(a+jω)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=0∑Nakdtkdky(t)=k=0∑Mbkdtkdkx(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(jω)=H(jω)X(jω)
我们称函数 H(jω)H(j \omega)H(jω) 为 系统频响函数 。
对微分方程两边做傅里叶变换:
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=0∑Nakdtkdky(t)}=F{k=0∑Mbkdtkdkx(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=0∑Nak(jω)kY(jω)=k=0∑Mbk(jω)kX(jω)
可得:
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(jω)=X(jω)Y(jω)=∑k=0Nak(jω)k∑k=0Mbk(jω)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]e−jkω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=−N1∑N2x[n]e−jkω0n=N1n=−∞∑+∞x[n]e−jkω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(ejω)=n=−∞∑+∞x[n]e−jω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π1∫2πX(ejω)ejωndω,综合公式
因为 X(ejω)ejωnX(e^{j\omega}) e^{j\omega n}X(ejω)ejωn 关于 ω\omegaω 周期为 2π2\pi2π 的函数,这说明我们的积分区间可以随意取一个连续的 2π2\pi2π 周期。
同时,对于非周期信号的傅里叶变换频谱公式 X(ejω)X(e^{j\omega})X(ejω) 来说,也是一个关于 ω\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(ejω)+BY(ejω) |
时移 | x[n−n0]x[n - n_0]x[n−n0] | X(ejω)e−jωn0X(e^{j\omega}) e^{-j\omega n_0}X(ejω)e−jωn0 |
频移 | x[n]ejω0nx[n]e^{j\omega_0 n}x[n]ejω0n | X(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−jω) |
时间翻转 | x[−n]x[-n]x[−n] | X(e−jω)X(e^{-j\omega})X(e−jω) |
时域尺度变换 | 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(ejω)Y(ejω) |
相乘 | 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π1∫2πX(ejθ)Y(ej(ω−θ))dθ |
时域差分 | x[n]−x[n−1]x[n] - x[n-1]x[n]−x[n−1] | (1−e−jω)X(ejω)(1 - e^{-j\omega})X(e^{j\omega})(1−e−jω)X(ejω) |
时域累加 | ∑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)1−e−jω1X(ejω)+π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(ejω) |
实信号的共轭对称性 | x[n]x[n]x[n] 为实信号 | X(ejω)=X∗(e−jω)X(e^{j\omega}) = X^*(e^{-j\omega})X(ejω)=X∗(e−jω) |
实偶信号 | x[n]x[n]x[n] 为实偶信号 | X(ejω)X(e^{j\omega})X(ejω) 为实偶函数 |
实奇信号 | x[n]x[n]x[n] 为实奇信号 | X(ejω)X(e^{j\omega})X(ejω) 为纯虚奇函数 |
实信号的奇偶分解 | 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(ejω)},jℑ{X(ejω)} |
周期信号的帕瓦尔定理 | 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\omega∑n=−∞+∞∣x[n]∣2=2π1∫2π∣X(ejω)∣2dω |
基本离散时间信号的傅里叶变换对
- 周期信号
考虑信号 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(ejω)=l=−∞∑+∞2πδ(ω−ω0−2π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(ejω)=k=−∞∑+∞2πakδ(ω−N2πk)
- 离散时间傅里叶系数的对偶性
若两个周期为 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]
- 离散时间傅里叶变换和连续时间傅里叶系数的对偶性
考虑离散时间傅里叶变换:
X(ejω)=∑n=−∞+∞x[n]e−jωn X(e^{j\omega}) = \sum_{n = -\infty}^{+\infty} x[n] e^{-j\omega n} X(ejω)=n=−∞∑+∞x[n]e−jωn
与连续时间傅里叶系数:
x(t)=∑k=−∞+∞ake−jkω0t x(t) = \sum_{k = -\infty}^{+\infty} a_k e^{-jk\omega_0 t} x(t)=k=−∞∑+∞ake−jkω0t
这相当于将 x[n]x[n]x[n] 看成是一个连续时间函数的傅里叶系数,同理,我们可以将 x(t)x(t)x(t) 看做是一个离散时间函数的傅里叶变换。
- 常用基本离散时间信号的傅里叶变换对表
信号 | 傅里叶变换 | 傅里叶系数(若为周期的) |
---|---|---|
∑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=−∞+∞δ(ω−ω0−2π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=−∞+∞(δ(ω−ω0−2π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=−∞+∞(δ(ω−ω0−2π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]=1 | 2π∑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,∣n∣≤N1,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)Nsin2π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=−∞+∞δ[n−kN] | 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}}1−ae−jω1 | - |
x[n]=1,∣n∣≤N1x[n] = 1, |n| \le N_1x[n]=1,∣n∣≤N1 | sinω(N1+12)sinω/2\frac{\sin{\omega(N_1 + \frac{1}{2})}}{\sin{\omega/2}}sinω/2sinω(N1+21) | - |
sinWnπ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)1−e−jω1+∑k=−∞+∞πδ(ω−wπk) | - |
δ[n−n0]\delta[n - n_0]δ[n−n0] | e−jωn0e^{-j\omega n_0}e−jωn0 | - |
(n+1)anu[n],∣a∣<1(n+1)a^n u[n], |a| < 1(n+1)anu[n],∣a∣<1 | 1(1−ae−jω)2\frac{1}{(1 - ae^{-j\omega})^2}(1−ae−jω)21 | - |
(n+r−1)!n!(r−1)!anu[n],∣a∣<1\frac{(n + r - 1)!}{n! (r - 1)!} a^n u[n], |a| < 1n!(r−1)!(n+r−1)!anu[n],∣a∣<1 | 1(1−ae−jω)r\frac{1}{(1 - ae^{-j\omega})^r}(1−ae−jω)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=0∑Naky[n−k]=k=0∑Mbkx[n−k]
通过傅里叶变换时移和线性的性质可以表示为频域上:
∑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=0∑Nake−jkωY(ejω)=k=0∑Mbke−jkωX(ejω)
或等效为:
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(ejω)=X(ejω)Y(ejω)=∑k=0Nake−jkω∑k=0Mbke−jkω