对应程佩青《数字信号处理》第三章 离散傅里叶变换,文章全部为原创,其中独创性地研究了从DFS推导出DFT,并探讨了DFT时域和频域
点数的关系,在中文互联网上为首创。
文章内容较多,建议点赞收藏后结合书本学习。
离散傅立叶级数(DFS)
傅立叶级数:周期函数(连续时间),离散频率
离散傅立叶级数:周期序列(离散时间),离散频率
表达式
周期为 N N N序列 x ~ ( n ) \tilde{x}(n) x~(n)用类比于连续周期信号 x ( t ) x(t) x(t)展开成傅立叶级数 X ( j k Ω 0 ) X(jk\Omega_0) X(jkΩ0)的方法展开
x ~ ( n ) = x ~ ( n + r N ) \tilde{x}(n)=\tilde{x}(n+rN) x~(n)=x~(n+rN)
展开成
x
~
(
n
)
=
1
N
∑
k
=
0
N
−
1
X
~
(
k
)
e
j
2
π
N
k
n
=
∑
k
=
0
N
−
1
x
~
(
n
)
W
N
−
n
k
\begin{aligned}\tilde{x}(n)=&\frac{1}{N}\sum_{k=0}^{N-1}\tilde{X}(k)e^{j\frac{2\pi}{N}kn}\\=&\sum_{k=0}^{N-1}\tilde{x}(n)W_N^{-nk}\end{aligned}
x~(n)==N1k=0∑N−1X~(k)ejN2πknk=0∑N−1x~(n)WN−nk
其中
W
N
n
=
e
−
j
2
π
N
W_N^{n}=e^{-j\frac{2\pi}{N}}
WNn=e−jN2π
DFS系数 X ~ ( k ) \tilde{X}(k) X~(k)的求取
X
~
(
k
)
=
∑
k
=
0
N
−
1
x
~
(
n
)
e
−
j
2
π
N
n
k
=
∑
n
=
0
N
−
1
x
~
(
n
)
W
N
n
k
\begin{aligned}\tilde{X}(k)=&\sum_{k=0}^{N-1}\tilde{x}(n)e^{-j\frac{2\pi}{N}nk}\\=&\sum_{n=0}^{N-1}\tilde{x}(n)W_N^{nk}\end{aligned}
X~(k)==k=0∑N−1x~(n)e−jN2πnkn=0∑N−1x~(n)WNnk
X
~
(
k
)
\tilde{X}(k)
X~(k)中的变量(k)由
W
N
n
W_N^n
WNn上的变量决定,记住这条变量k的可换性,对于理解DFS的含义有很大帮助。
W N n W_N^n WNn的性质
共轭对称性
W N n = ( W N − n ) ∗ W_N^n=(W_N^{-n})^* WNn=(WN−n)∗,易证
周期性
W N n = W N n + i N W_N^n=W_N^{n+iN} WNn=WNn+iN
以 N N N为周期,与DFS系数的周期性共同决定了DFS的周期性
可约性
W N i n = W N i n W_N^{in}=W_{\frac{N}{i}}^{n} WNin=WiNn
W i N i n = W N n W_{iN}^{in}=W_{N}^{n} WiNin=WNn
推导FFT有使用
正交性
1 N ∑ k = 0 N − 1 W N n k ( W N m k ) ∗ = 1 N ∑ k = 0 N − 1 W N ( n − m ) k = { 1 , n − m = i N 0 , n − m ≠ i N \begin{aligned}&\frac{1}{N}\sum_{k=0}^{N-1}W_N^{nk}(W_N^{mk})^*\\=&\frac{1}{N}\sum_{k=0}^{N-1}W_N^{(n-m)k}=\begin{cases}1,&n-m=iN\\0,&n-m\ne iN\end{cases}\end{aligned} =N1k=0∑N−1WNnk(WNmk)∗N1k=0∑N−1WN(n−m)k={1,0,n−m=iNn−m=iN
也写作
1 N ∑ k = 0 N − 1 e j 2 π N r n = { 1 , r = m N , m 为任意整数 0 , 其他 r \begin{aligned}\frac{1}{N}\sum_{k=0}^{N-1}e^{j\frac{2\pi}{N}rn}=\begin{cases}1,&r=mN,m为任意整数\\0,&其他r\end{cases}\end{aligned} N1k=0∑N−1ejN2πrn={1,0,r=mN,m为任意整数其他r
推导过程及意义
方法一 类比FS定义先假设用系数表达再求出系数
对于周期序列 x ~ ( n ) \tilde{x}(n) x~(n)用傅立叶级数的方法展开,
其谐波分量 e j k 2 π N n e^{jk\frac{2\pi}{N}n} ejkN2πn不像连续周期函数的傅立叶级数的谐波分量 e j k w 0 t e^{jkw_0t} ejkw0t,
虽然都是离散取值,
但DFS为频域周期,故只取一个周期(0-N-1)内的谐波
x ^ ( n ) = 1 N ∑ k = 0 N − 1 X ~ ( k ) e j 2 π N k n \hat{x}(n)=\frac{1}{N}\sum_{k=0}^{N-1}\tilde{X}(k)e^{j\frac{2\pi}{N}kn} x^(n)=N1∑k=0N−1X~(k)ejN2πkn
设成这个形式是运用到了 W N n W_N^n WNn的正交性,然后利用这个性质得到 X ~ ( k ) \tilde{X}(k) X~(k)的表达式(但总体很不自然)
方法二 从离散时间傅立叶变换(DTFT),频域离散化出DFS
离散时间傅立叶变换只是时间上离散,并没有抽样周期化
X ( e j w ) = ∑ n = − ∞ ∞ x ( n ) e − j w n X(e^{jw})=\sum_{n=-\infty}^{\infty}x(n)e^{-jwn} X(ejw)=∑n=−∞∞x(n)e−jwn
DTFT在频域上 X ( e j w ) X(e^{jw}) X(ejw)连续,因为模拟信号的频率 Ω \Omega Ω连续,数字频率 w = Ω T w=\Omega T w=ΩT也连续
要使得频域上离散,模拟频率也要离散化, Ω = k Ω 0 \Omega=k\Omega_0 Ω=kΩ0,
抽样后的数字频率也离散化 w = k Ω 0 T w=k\Omega_0T w=kΩ0T,
可以写作 X ( e j k Ω 0 T ) = ∑ n = − ∞ ∞ x ~ ( n ) e − j n k Ω 0 T = ∑ n = 0 N − 1 x ~ ( n ) e − j n k Ω 0 T X(e^{jk\Omega_0T})=\sum_{n=-\infty}^{\infty}\tilde{x}(n)e^{-jnk\Omega_0T}=\sum_{n=0}^{N-1}\tilde{x}(n)e^{-jnk\Omega_0T} X(ejkΩ0T)=∑n=−∞∞x~(n)e−jnkΩ0T=∑n=0N−1x~(n)e−jnkΩ0T
讨论 时域与频域的周期和频率的关系
-
(对于连续周期函数的FS),函数周期 T 0 T_0 T0,则频率为 F 0 = 1 T 0 F_0=\frac{1}{T_0} F0=T01,
频域离散因而频域谱线间隔为 F 0 = 1 T 0 F_0=\frac{1}{T_0} F0=T01(谐波为 k F 0 kF_0 kF0,这里以频率非角频率为单位) -
对其时域进行抽样,(离散周期函数的DFS),抽样间隔为 T T T,则抽样频率为 f s = 1 T f_s=\frac{1}{T} fs=T1,
借由理想抽样信号的FT的结论(DTFT), X ^ a ( j Ω ) = 1 T ∑ k = − ∞ ∞ X a [ j ( Ω − k Ω s ) ] \hat{X}_a(j\Omega)=\frac{1}{T}\sum_{k=-\infty}^{\infty}X_a[j(\Omega-k\Omega_s)] X^a(jΩ)=T1∑k=−∞∞Xa[j(Ω−kΩs)]
(实际上DTFT为非周期离散信号,这里借用的是推导过程中理想抽样信号由非周期函数乘以周期离散函数,得到的周期离散函数的频谱特性),
在频域上的周期为 f s = 1 T f_s=\frac{1}{T} fs=T1 -
那么由于已经假设了序列为N点周期, x ~ ( n ) = x ~ ( n + r N ) \tilde{x}(n)=\tilde{x}(n+rN) x~(n)=x~(n+rN),在一个周期内有 N N N个抽样点, T 0 = N T T_0=NT T0=NT
上面已经得到了频域上频率函数也为周期函数,且频域的谱线间隔即频域的“离散化抽样频率”和周期已知,设频域上一个周期有 N ’ N’ N’个抽样点,则 N ’ = f s F 0 = T 0 T = N N’=\frac{f_s}{F_0}=\frac{T_0}{T}=N N’=F0fs=TT0=N
故频域的周期与时域的周期相同,均为N。
时域抽样间隔 T T T,抽样频率 f s = 1 T f_s=\frac{1}{T} fs=T1,周期为N,
频域谱线间隔 d Ω = Ω 0 d\Omega=\Omega_0 dΩ=Ω0,则在频域上的抽样点数 N = Ω s Ω 0 N=\frac{\Omega_s}{\Omega_0} N=Ω0Ωs与时域的周期相同,时域和频域的抽样点数(周期)相同
x
(
n
)
=
1
2
π
∫
−
π
π
X
(
e
j
ω
)
e
j
ω
n
d
ω
=
T
2
π
∫
−
π
T
π
T
X
(
e
j
Ω
T
)
e
j
n
Ω
T
d
Ω
=
1
Ω
s
∫
−
Ω
s
2
Ω
s
2
X
(
e
j
Ω
T
)
e
j
n
Ω
T
d
Ω
\begin{aligned}x(n)=&\frac{1}{2\pi}\int_{-\pi}^{\pi}X(e^{j\omega})e^{j\omega n}d\omega\\=&\frac{T}{2\pi}\int_{-\frac{\pi}{T}}^{\frac{\pi}{T}}X(e^{j\Omega T})e^{jn\Omega T}d\Omega\\=&\frac{1}{\Omega_s}\int_{-\frac{\Omega_s}{2}}^{\frac{\Omega_s}{2}}X(e^{j\Omega T})e^{jn\Omega T}d\Omega\end{aligned}
x(n)===2π1∫−ππX(ejω)ejωndω2πT∫−TπTπX(ejΩT)ejnΩTdΩΩs1∫−2Ωs2ΩsX(ejΩT)ejnΩTdΩ
频率周期也为
N
N
N,之前假设
Ω
=
k
Ω
0
\Omega=k\Omega_0
Ω=kΩ0,那么
k
k
k的取值范围为
[
0
,
N
−
1
]
[0,N-1]
[0,N−1]
x
~
(
n
)
=
Ω
0
Ω
s
∑
k
=
0
N
−
1
X
(
e
j
k
Ω
0
T
)
e
j
n
k
Ω
0
T
\tilde{x}(n)=\frac{\Omega_0}{\Omega_s}\sum_{k=0}^{N-1}X(e^{jk\Omega_0 T})e^{jnk\Omega_0 T}
x~(n)=ΩsΩ0k=0∑N−1X(ejkΩ0T)ejnkΩ0T
再用到
Ω
0
Ω
s
=
F
0
f
s
=
1
N
\frac{\Omega_0}{\Omega_s}=\frac{F_0}{f_s}=\frac{1}{N}
ΩsΩ0=fsF0=N1
以及
Ω
0
T
=
2
π
F
0
T
=
2
π
F
0
f
s
=
2
π
N
\Omega_0T=2\pi F_0T=2\pi\frac{F_0}{f_s}=\frac{2\pi}{N}
Ω0T=2πF0T=2πfsF0=N2π
最终得到的表达式为
X
~
(
k
)
=
X
~
(
e
j
2
π
N
k
)
=
∑
n
=
0
N
−
1
x
~
(
n
)
e
−
j
2
π
N
k
n
\tilde{X}(k)=\tilde{X}(e^{j\frac{2\pi}{N}k})=\sum_{n=0}^{N-1}\tilde{x}(n)e^{-j\frac{2\pi}{N}kn}
X~(k)=X~(ejN2πk)=n=0∑N−1x~(n)e−jN2πkn
x
~
(
n
)
=
x
~
(
n
T
)
=
1
N
∑
k
=
0
N
−
1
X
~
(
k
)
e
j
2
π
N
k
n
\tilde{x}(n)=\tilde{x}(nT)=\frac{1}{N}\sum_{k=0}^{N-1}\tilde{X}(k)e^{j\frac{2\pi}{N}kn}
x~(n)=x~(nT)=N1k=0∑N−1X~(k)ejN2πkn
其谐波分量为
2
π
N
k
n
\frac{2\pi}{N}kn
N2πkn
求法及与Z变换关系
由于频域被抽样成了N个点,原来的数字频率 w w w的一个周期为 2 π 2\pi 2π,故可以等效作 w = 2 π N k , k ∈ [ 0 , N − 1 ] w=\frac{2\pi}{N}k,k\in[0,N-1] w=N2πk,k∈[0,N−1]
在上面给出了由DTFT推导出DFS表达式形式的一种方法,由于DTFT与Z变换关系已知,
故可以借此得到DFS(的系数 X ~ ( k ) \tilde{X}(k) X~(k))与DTFT和Z变换的关系
X ~ ( k ) = X ~ ( e j 2 π N k ) = X ( e j w ) ∣ w = 2 π N k \tilde{X}(k)=\tilde{X}(e^{j\frac{2\pi}{N}}k)=X(e^{jw})|_{w=\frac{2\pi}{N}k} X~(k)=X~(ejN2πk)=X(ejw)∣w=N2πk
DFS等于 x ~ ( n ) \tilde{x}(n) x~(n)的一个周期有限序列 x ( n ) x(n) x(n)的DTFT在 w = 2 π N k w=\frac{2\pi}{N}k w=N2πk处的抽样值,为单位圆上z变换的N个等间隔的点的抽样值。
故可先计算主值区间DTFT,然后用频域上抽样值来表示
抽样后也就有了以N个抽样点为周期的周期性,但表达式上并不显示表达( X ( e j w ) 与 X(e^{jw})与 X(ejw)与 X ~ ( k ) \tilde{X}(k) X~(k)都有周期性)
性质
线性
D
F
S
[
a
x
1
~
(
n
)
+
b
x
2
~
(
n
)
]
=
a
X
~
1
(
k
)
+
b
X
~
2
(
k
)
DFS[a\tilde{x_1}(n)+b\tilde{x_2}(n)]=a\tilde{X}_1(k)+b\tilde{X}_2(k)
DFS[ax1~(n)+bx2~(n)]=aX~1(k)+bX~2(k)
(要求
x
~
1
(
n
)
\tilde{x}_1(n)
x~1(n)和
x
~
2
(
n
)
\tilde{x}_2(n)
x~2(n)周期相同,由DFS系数的周期性决定)
移位
D F S [ x ~ ( n + m ) ] = W N − m k X ~ ( k ) DFS[\tilde{x}(n+m)]=W_N^{-mk}\tilde{X}(k) DFS[x~(n+m)]=WN−mkX~(k)
考虑到 W N n W_N^n WNn表达式幂次已经有负号,故符号反转
证明时只改变时域序列的表达,通过变量代换得到结果
调制特性
D F S [ W N l n x ~ ( n ) ] = X ~ ( k + l ) DFS[W_N^{ln}\tilde{x}(n)]=\tilde{X}(k+l) DFS[WNlnx~(n)]=X~(k+l)
证明用到了前面的一个观察结果(见 X ~ ( k ) \tilde{X}(k) X~(k)的求取): X ~ ( k ) \tilde{X}(k) X~(k)中的变量(k)由 W N n W_N^n WNn上的变量决定
对偶性
D F S [ X ~ ( n ) ] = N x ~ ( − k ) DFS[\tilde{X}(n)]=N\tilde{x}(-k) DFS[X~(n)]=Nx~(−k)
证明:逆变换表达式 x ~ ( n ) = 1 N ∑ k = 0 N − 1 X ~ ( k ) W N − n k \tilde{x}(n)=\frac{1}{N}\sum_{k=0}^{N-1}\tilde{X}(k)W_N^{-nk} x~(n)=N1∑k=0N−1X~(k)WN−nk
将n和k互换(n换-k,k换n),
得 N x ~ ( − k ) = ∑ n = 0 N − 1 X ~ ( n ) W N n k N\tilde{x}(-k)=\sum_{n=0}^{N-1}\tilde{X}(n)W_N^{nk} Nx~(−k)=∑n=0N−1X~(n)WNnk
证毕
周期卷积和性质
周期卷积和
x ~ 1 ( n ) \tilde{x}_1(n) x~1(n)和 x ~ 2 ( n ) \tilde{x}_2(n) x~2(n)均为周期为 N N N的序列
则
y ~ ( n ) = x ~ 1 ( n ) ⊛ x ~ 2 ( n ) = ∑ m = 0 N − 1 x ~ 1 ( m ) x ~ 2 ( n − m ) = ∑ m = 0 N − 1 x ~ 2 ( m ) x ~ 1 ( n − m ) \begin{aligned}\tilde{y}(n)=&\tilde{x}_1(n)\circledast\tilde{x}_2(n)\\=&\sum_{m=0}^{N-1}\tilde{x}_1(m)\tilde{x}_2(n-m)\\=&\sum_{m=0}^{N-1}\tilde{x}_2(m)\tilde{x}_1(n-m)\end{aligned} y~(n)===x~1(n)⊛x~2(n)m=0∑N−1x~1(m)x~2(n−m)m=0∑N−1x~2(m)x~1(n−m)
求和仅为一个区间,其中一个序列仅移动一个区间的长度
由于 x ~ 1 ( m ) \tilde{x}_1(m) x~1(m)和 x ~ 2 ( n − m ) \tilde{x}_2(n-m) x~2(n−m)也均为周期为 N N N的序列,故 y ~ ( n ) \tilde{y}(n) y~(n)也为周期为 N N N的序列
求法:(N点圆周卷积和的周期延拓)
前面已经讨论其具有周期性,故求出
0
∼
N
−
1
0\sim N-1
0∼N−1主值区间再进行周期延拓
可通过表格法计算,利用两个序列的周期性得到
x
~
i
(
n
−
m
)
\tilde{x}_i(n-m)
x~i(n−m)的值
周期卷积和的DFS
Y ~ ( k ) = X ~ 1 ( k ) X ~ 2 ( k ) \tilde{Y}(k)=\tilde{X}_1(k)\tilde{X}_2(k) Y~(k)=X~1(k)X~2(k)
若 y ~ ( n ) = x ~ 1 ( n ) x ~ 2 ( n ) \tilde{y}(n)=\tilde{x}_1(n)\tilde{x}_2(n) y~(n)=x~1(n)x~2(n), Y ~ ( k ) = 1 N X ~ 1 ( k ) ∗ X ~ 2 ( k ) \tilde{Y}(k)=\frac{1}{N}\tilde{X}_1(k)*\tilde{X}_2(k) Y~(k)=N1X~1(k)∗X~2(k)
离散傅里叶变换(DFT)
离散傅里叶级数:时域周期离散,频域离散周期
将DFS的时域和频域都截成有限长度等于周期N,得DFT
定义
主值序列 x ( n ) x(n) x(n)仅在 x ~ ( n ) \tilde{x}(n) x~(n)主值区间 [ 0 , N − 1 ] [0,N-1] [0,N−1]有值
x
(
n
)
=
x
~
(
n
)
R
N
(
n
)
=
x
(
(
n
)
)
N
R
N
(
n
)
\begin{aligned}x(n)=&\tilde{x}(n)R_N(n)\\=&x((n))_NR_N(n)\end{aligned}
x(n)==x~(n)RN(n)x((n))NRN(n)
其中
x
(
(
n
)
)
N
=
x
(
n
m
o
d
N
)
x((n))_N=x(n\,mod\,N)
x((n))N=x(nmodN)故DFT有周期性的含义
x ~ ( n ) = x ( ( n ) ) N = ∑ r = − ∞ ∞ x ( n + r N ) \tilde{x}(n)=x((n))_N=\sum_{r=-\infty}^{\infty}x(n+rN) x~(n)=x((n))N=r=−∞∑∞x(n+rN)
同理频域上为 X ( k ) = X ~ ( k ) R N ( k ) = X ( ( k ) ) N R N ( k ) X(k)=\tilde{X}(k)R_N(k)=X((k))_NR_N(k) X(k)=X~(k)RN(k)=X((k))NRN(k)
定义式
X ( k ) = D F T [ x ( n ) ] = ∑ k = 0 N − 1 x ( n ) W N n k , k = 0 , 1 , . . . , N − 1 X(k)=DFT[x(n)]=\sum_{k=0}^{N-1}x(n)W_N^{nk},k=0,1,...,N-1 X(k)=DFT[x(n)]=k=0∑N−1x(n)WNnk,k=0,1,...,N−1
x ( n ) = I D F T [ X ( k ) ] = 1 N ∑ n = 0 N − 1 X ( k ) W N − n k , n = 0 , 1 , . . . , N − 1 x(n)=IDFT[X(k)]=\frac{1}{N}\sum_{n=0}^{N-1}X(k)W_N^{-nk},n=0,1,...,N-1 x(n)=IDFT[X(k)]=N1n=0∑N−1X(k)WN−nk,n=0,1,...,N−1
拓展 利用DFT计算IDFT
考虑
x
∗
(
n
)
=
1
N
[
∑
n
=
0
N
−
1
X
(
k
)
W
N
−
n
k
]
∗
=
1
N
∑
n
=
0
N
−
1
X
∗
(
k
)
W
N
n
k
=
1
N
D
F
T
[
X
∗
(
k
)
]
x^*(n)=\frac{1}{N}[\sum_{n=0}^{N-1}X(k)W_N^{-nk}]*=\frac{1}{N}\sum_{n=0}^{N-1}X^*(k)W_N^{nk}=\frac{1}{N}{DFT[X^*(k)]}
x∗(n)=N1[n=0∑N−1X(k)WN−nk]∗=N1n=0∑N−1X∗(k)WNnk=N1DFT[X∗(k)]
x
(
n
)
=
1
N
D
F
T
[
X
∗
(
k
)
]
∗
x(n)=\frac{1}{N}{DFT[X^*(k)]}^*
x(n)=N1DFT[X∗(k)]∗
提供了一种使用FFT公式计算IFFT的方法
矩阵表示
DFT矩阵
X
=
W
N
x
X=W_Nx
X=WNx
x
x
x为时域序列的列向量,
X
X
X为频域序列的列向量
N N N点DFT矩阵 W N = [ 1 1 ⋯ 1 1 W N 1 ⋯ W N N − 1 ⋮ ⋮ ⋱ ⋮ 1 W N N − 1 ⋯ W N ( N − 1 ) ( N − 1 ) ] W_N=\begin{bmatrix} &1 & 1 &\cdots& 1\\&1 & W_N^1 &\cdots &W_N^{N-1}\\&\vdots&\vdots&\ddots&\vdots\\&1&W_N^{N-1}&\cdots&W_N^{(N-1)(N-1)}\end{bmatrix} WN= 11⋮11WN1⋮WNN−1⋯⋯⋱⋯1WNN−1⋮WN(N−1)(N−1)
理解谐波分量 W N n k W_N^{nk} WNnk,对于每一个 X ( k ) X(k) X(k)计算 k k k定 n n n动
IDFT矩阵
x = W N − 1 X x=W_N^{-1}X x=WN−1X,
N点IDFT矩阵 W N − 1 = 1 N W N ∗ W_N^{-1}=\frac{1}{N}W_N^* WN−1=N1WN∗
由 X ( e j w ) X(e^{jw}) X(ejw)求 X ( k ) X(k) X(k)
与求DFS X ~ ( k ) \tilde{X}(k) X~(k)的过程一样,但是考虑到 X ( k ) X(k) X(k)实际上是 X ~ ( k ) \tilde{X}(k) X~(k)的一个主值区间,
而 X ( e j w ) X(e^{jw}) X(ejw)实际上隐含有 2 π 2\pi 2π的周期性,因此需要限制 k k k的范围。
在使用DTFT求DFS时,时域序列已经截取了主值区间。
X
(
k
)
=
X
(
e
j
w
)
∣
w
=
2
π
N
k
,
k
=
0
,
1
,
2
,
⋅
,
N
−
1
X(k)=X(e^{jw})|_{w=\frac{2\pi}{N}k},k=0,1,2,\cdot,N-1
X(k)=X(ejw)∣w=N2πk,k=0,1,2,⋅,N−1
x ( n ) x(n) x(n)和 x ( k ) x(k) x(k)各参量之间的关系
此前在(讨论 时域与频域的周期和频率的关系)中就已经在推导DFS的过程中给出了这些关系,在DFT中这些关系仍然存在
T
0
=
N
T
=
N
f
s
=
1
F
0
T_0=NT=\frac{N}{f_s}=\frac{1}{F_0}
T0=NT=fsN=F01
则有频率分辨率
F
0
F_0
F0为频域相邻两抽样点的频率间距
F
0
=
f
s
N
=
1
N
T
=
1
T
0
F_0=\frac{f_s}{N}=\frac{1}{NT}=\frac{1}{T_0}
F0=Nfs=NT1=T01
由第一篇博客(数字信号处理必备知识)中的数字频率和模拟频率关系
w
k
=
2
π
N
k
=
2
π
f
k
f
s
=
2
π
f
k
T
w_k=\frac{2\pi}{N}k=2\pi\frac{f_k}{f_s}=2\pi f_kT
wk=N2πk=2πfsfk=2πfkT得到频域第k个抽样点频率
f
k
=
k
N
T
=
k
f
s
N
f_k=\frac{k}{NT}=\frac{kf_s}{N}
fk=NTk=Nkfs
同样说明了
1
N
T
\frac{1}{NT}
NT1为频率分辨率
把 T 0 = N T T_0=NT T0=NT称为记录时间
增加记录时间,就能减小 F 0 = 1 T 0 F_0=\frac{1}{T_0} F0=T01,提高频率分辨率
性质
线性性质
D F T [ a x 1 ( n ) + b x 2 ( n ) ] = a X 1 ( k ) + b X 2 ( k ) DFT[ax_1(n)+bx_2(n)]=aX_1(k)+bX_2(k) DFT[ax1(n)+bx2(n)]=aX1(k)+bX2(k)
圆周移位性质
有限长序列
x
(
n
)
x(n)
x(n)的
m
m
m点圆周移位是在主值区间上进行的周期延拓序列的线性移位取主值为
x
m
(
n
)
=
x
(
(
n
+
m
)
)
N
R
N
(
n
)
x_m(n)=x((n+m))_NR_N(n)
xm(n)=x((n+m))NRN(n)
x
(
(
n
+
m
)
)
N
=
x
~
(
n
+
m
)
x((n+m))_N=\tilde{x}(n+m)
x((n+m))N=x~(n+m)
D
F
T
[
x
m
(
n
)
]
=
X
m
(
k
)
=
W
N
−
m
k
X
(
k
)
DFT[x_m(n)]=X_m(k)=W_N^{-mk}X(k)
DFT[xm(n)]=Xm(k)=WN−mkX(k)
与DFS的线性移位性质的结果相同,因为其实际上就是DFS的线性移位后截取主值区间的结果
调制特性(频域圆周移位)
D F T [ W N n l x ( n ) ] = X ( ( k + l ) ) N R N ( k ) DFT[W_N^{nl}x(n)]=X((k+l))_NR_N(k) DFT[WNnlx(n)]=X((k+l))NRN(k)
圆周共轭对称性质
圆周共轭序列
圆周共轭对称序列
x
e
p
(
n
)
=
x
e
p
∗
(
(
−
n
)
)
N
R
N
(
n
)
=
x
e
p
∗
(
N
−
n
)
,
n
=
0
,
1
,
2
,
⋯
,
N
−
1
x_{ep}(n)=x_{ep}^*((-n))_NR_N(n)=x_{ep}^*(N-n),n=0,1,2,\cdots,N-1
xep(n)=xep∗((−n))NRN(n)=xep∗(N−n),n=0,1,2,⋯,N−1
圆周共轭反对称序列
x
o
p
(
n
)
=
−
x
o
p
∗
(
(
−
n
)
)
N
R
N
(
n
)
=
—
x
o
p
∗
(
N
−
n
)
,
n
=
0
,
1
,
2
,
⋯
,
N
−
1
x_{op}(n)=-x_{op}^*((-n))_NR_N(n)=—x_{op}^*(N-n),n=0,1,2,\cdots,N-1
xop(n)=−xop∗((−n))NRN(n)=—xop∗(N−n),n=0,1,2,⋯,N−1
利用了周期共轭对称序列(无限长)来定义,并取了主值区间上的值,回顾圆周移位(圆周的主值性)
与有限长序列的共轭对称序列对比
有限长共轭对称序列 x e ( n ) x_e(n) xe(n)和共轭反对称序列 x o ( n ) x_o(n) xo(n)有 2 N − 1 2N-1 2N−1个点,以原点为对称中心对称
而圆周共轭对称序列和反对称序列是在有限长序列的周期延拓序列的共轭对称序列的主值区间上定义的,只有
N
N
N个点。
且对称中心为
n
+
(
N
−
n
)
2
=
N
2
\frac{n+(N-n)}{2}=\frac{N}{2}
2n+(N−n)=2N。
在以该对称中心讨论有限长序列的圆周对称性时需要在 n = N n=N n=N处补上 n = 0 n=0 n=0的序列值。(直接在主值区间即有限长序列上讨论对称性)
圆周翻褶序列
考虑到圆周的主值性,也满足 N N N点DFT长度相等的需要
x ( n ) x(n) x(n)的圆周翻褶序列以圆周对称中心 N 2 \frac{N}{2} 2N翻褶
把 x ( ( − n ) ) N R N ( n ) = x ( N − n ) x((-n))_NR_N(n)=x(N-n) x((−n))NRN(n)=x(N−n)称作 x ( n ) x(n) x(n)的圆周翻褶序列
圆周翻褶序列的DFT
D
F
T
[
x
(
N
−
n
)
]
=
X
(
N
−
k
)
DFT[x(N-n)]=X(N-k)
DFT[x(N−n)]=X(N−k)
实际上是对DFT的隐含的周期性的理解,重新截取主值区间
有限长序列的另一分解
有限长序列一定可以表示成圆周共轭对称分量和圆周共轭反对称分量之和(这就是直接在有限长序列上对其进行讨论)
x
(
n
)
=
x
e
p
(
n
)
+
x
o
p
(
n
)
x(n)=x_{ep}(n)+x_{op}(n)
x(n)=xep(n)+xop(n)
参考第二篇博客中的《共轭对称分量和反对称分量分解》,进行周期延拓后结论不变
有
x
~
e
(
n
)
=
1
2
[
x
~
(
n
)
+
x
~
∗
(
−
n
)
]
\tilde{x}_e(n)=\frac{1}{2}[\tilde{x}(n)+\tilde{x}^*(-n)]
x~e(n)=21[x~(n)+x~∗(−n)]
x
~
o
(
n
)
=
1
2
[
x
~
(
n
)
−
x
~
∗
(
−
n
)
]
\tilde{x}_o(n)=\frac{1}{2}[\tilde{x}(n)-\tilde{x}^*(-n)]
x~o(n)=21[x~(n)−x~∗(−n)]
故
x
e
p
(
n
)
=
x
~
e
(
n
)
R
N
(
n
)
=
1
2
[
x
(
n
)
+
x
∗
(
N
−
n
)
]
,
n
=
0
,
1
,
2
,
⋯
,
N
−
1
x_{ep}(n)=\tilde{x}_e(n)R_N(n)=\frac{1}{2}[x(n)+x^*(N-n)],n=0,1,2,\cdots,N-1
xep(n)=x~e(n)RN(n)=21[x(n)+x∗(N−n)],n=0,1,2,⋯,N−1
x
o
p
(
n
)
=
x
~
o
(
n
)
R
N
(
n
)
=
1
2
[
x
(
n
)
−
x
∗
(
N
−
n
)
]
,
n
=
0
,
1
,
2
,
⋯
,
N
−
1
x_{op}(n)=\tilde{x}_o(n)R_N(n)=\frac{1}{2}[x(n)-x^*(N-n)],n=0,1,2,\cdots,N-1
xop(n)=x~o(n)RN(n)=21[x(n)−x∗(N−n)],n=0,1,2,⋯,N−1
对偶性
由《DFS对偶性》再结合圆周取主值区间得到
D
F
T
[
X
(
n
)
]
=
N
x
(
N
−
k
)
DFT[X(n)]=Nx(N-k)
DFT[X(n)]=Nx(N−k)
共轭序列的DFT
D
F
T
[
x
∗
(
n
)
]
=
X
∗
(
N
−
k
)
DFT[x^*(n)]=X^*(N-k)
DFT[x∗(n)]=X∗(N−k)
证:
D
F
T
[
x
∗
(
n
)
]
=
∑
n
=
0
N
−
1
x
∗
(
n
)
W
N
n
k
R
N
(
k
)
=
[
∑
n
=
0
N
−
1
x
(
n
)
W
N
−
n
k
]
∗
R
N
(
k
)
=
X
∗
(
N
−
k
)
\begin{aligned}DFT[x^*(n)]=&\sum_{n=0}^{N-1}x^*(n)W_N^{nk}R_N(k)\\=&[\sum_{n=0}^{N-1}x(n)W_N^{-nk}]^*R_N(k)\\=&X^*(N-k)\end{aligned}
DFT[x∗(n)]===n=0∑N−1x∗(n)WNnkRN(k)[n=0∑N−1x(n)WN−nk]∗RN(k)X∗(N−k)
当
x
(
n
)
x(n)
x(n)为实序列时,
x
(
n
)
=
x
∗
(
n
)
x(n)=x^*(n)
x(n)=x∗(n)则
X
(
k
)
=
X
∗
(
N
−
k
)
X(k)=X^*(N-k)
X(k)=X∗(N−k)
而
X
(
k
)
=
∣
X
(
k
)
∣
e
j
a
r
g
[
X
(
k
)
]
X(k)=|X(k)|e^{jarg[X(k)]}
X(k)=∣X(k)∣ejarg[X(k)]
故若
x
(
n
)
x(n)
x(n)为实序列,
其幅度谱 ∣ X ( k ) ∣ = ∣ X ( N − k ) ∣ |X(k)|=|X(N-k)| ∣X(k)∣=∣X(N−k)∣圆周偶对称,
相位谱 a r g [ X ( k ) ] = − a r g [ X ( N − k ) ] arg[X(k)]=-arg[X(N-k)] arg[X(k)]=−arg[X(N−k)]圆周奇对称。
(连续实信号 x ( t ) x(t) x(t)的FT的幅度谱和相位谱关于原点也有类似的对称关系)
圆周共轭翻褶序列的DFT
D
F
T
[
x
∗
(
N
−
n
)
]
=
X
∗
(
k
)
DFT[x^*(N-n)]=X^*(k)
DFT[x∗(N−n)]=X∗(k)
证:
D
F
T
[
x
∗
(
N
−
n
)
]
=
D
F
T
[
x
∗
(
(
−
n
)
)
N
R
N
(
n
)
]
=
∑
n
=
0
N
−
1
x
∗
(
(
−
n
)
)
N
W
N
n
k
=
∑
n
=
0
N
−
1
[
x
(
(
−
n
)
)
N
W
N
−
n
k
]
∗
=
∑
n
=
−
(
N
−
1
)
0
[
x
(
(
n
)
)
N
W
N
n
k
]
∗
=
[
∑
n
=
0
N
−
1
x
(
(
n
)
)
N
W
N
n
k
]
∗
,
周期性
=
X
∗
(
k
)
\begin{aligned}DFT[x^*(N-n)]=&DFT[x^*((-n))_NR_N(n)]\\=&\sum_{n=0}^{N-1}x^*((-n))_NW_N^{nk}\\=&\sum_{n=0}^{N-1}[x((-n))_NW_N^{-nk}]^*\\=&\sum_{n=-(N-1)}^{0}[x((n))_NW_N^{nk}]^*\\=&[\sum_{n=0}^{N-1}x((n))_NW_N^{nk}]^*,周期性\\=&X^*(k)\end{aligned}
DFT[x∗(N−n)]======DFT[x∗((−n))NRN(n)]n=0∑N−1x∗((−n))NWNnkn=0∑N−1[x((−n))NWN−nk]∗n=−(N−1)∑0[x((n))NWNnk]∗[n=0∑N−1x((n))NWNnk]∗,周期性X∗(k)
时域和频域的对偶关系
经过上面的讨论得到了将序列分解成圆周共轭对称序列和圆周共轭反对成序列的方法,
并能用原序列和原序列的共轭翻褶序列来表示圆周共轭对称分量和圆周共轭反对称分量;
得到了共轭序列的DFT和圆周共轭翻褶序列的DFT。
由此可以得到对偶关系
x
(
n
)
=
R
e
[
x
(
n
)
]
+
j
I
m
[
x
(
n
)
]
,
X
(
k
)
=
X
e
p
(
k
)
+
X
o
p
(
k
)
x(n)=Re[x(n)]+jIm[x(n)],X(k)=X_{ep}(k)+X_{op}(k)
x(n)=Re[x(n)]+jIm[x(n)],X(k)=Xep(k)+Xop(k)
x
(
n
)
=
x
e
p
(
n
)
+
x
o
p
(
n
)
,
X
(
k
)
=
R
e
[
X
(
k
)
]
+
j
I
m
[
X
(
k
)
]
x(n)=x_{ep}(n)+x_{op}(n),X(k)=Re[X(k)]+jIm[X(k)]
x(n)=xep(n)+xop(n),X(k)=Re[X(k)]+jIm[X(k)]
利用对偶减小DFT计算量
利用圆周共轭对称关系,可以减小DFT计算量,一次只需计算一半的 X ( k ) X(k) X(k)
用一次DFT计算两个N点实数序列DFT:
考虑
X
1
(
k
)
=
D
F
T
[
x
1
(
n
)
]
X_1(k)=DFT[x_1(n)]
X1(k)=DFT[x1(n)],
X
2
(
k
)
=
D
F
T
[
x
2
(
n
)
]
X_2(k)=DFT[x_2(n)]
X2(k)=DFT[x2(n)]
构造
x
(
n
)
=
x
1
(
n
)
+
j
x
2
(
n
)
x(n)=x_1(n)+jx_2(n)
x(n)=x1(n)+jx2(n)则
X
(
k
)
=
X
1
(
k
)
+
j
X
2
(
k
)
X(k)=X_1(k)+jX_2(k)
X(k)=X1(k)+jX2(k)
其中
X
1
(
k
)
=
R
e
[
X
(
k
)
]
=
1
2
[
X
(
k
)
+
X
∗
(
N
−
k
)
]
=
D
F
T
[
R
e
[
x
(
n
)
]
]
=
X
e
p
(
k
)
\begin{aligned}X_1(k)=&Re[X(k)]=\frac{1}{2}[X(k)+X^*(N-k)]\\=&DFT[Re[x(n)]]=X_{ep}(k)\end{aligned}
X1(k)==Re[X(k)]=21[X(k)+X∗(N−k)]DFT[Re[x(n)]]=Xep(k)
j X 2 ( k ) = j I m [ X ( k ) ] = j 1 2 [ X ( k ) − X ∗ ( N − k ) ] = D F T [ j I m [ x ( n ) ] ] = j X o p ( k ) \begin{aligned}jX_2(k)=&jIm[X(k)]=j\frac{1}{2}[X(k)-X^*(N-k)]\\=&DFT[jIm[x(n)]]=jX_{op}(k)\end{aligned} jX2(k)==jIm[X(k)]=j21[X(k)−X∗(N−k)]DFT[jIm[x(n)]]=jXop(k)
DFT形式下的帕塞瓦尔定理
∑
n
=
0
N
−
1
x
2
(
n
)
=
1
N
∑
k
=
0
N
−
1
∣
X
(
k
)
∣
2
\sum_{n=0}^{N-1}x^2(n)=\frac{1}{N}\sum_{k=0}^{N-1}|X(k)|^2
n=0∑N−1x2(n)=N1k=0∑N−1∣X(k)∣2
一个序列在时域计算的能量和频域计算相等
圆周卷积和定理
圆周卷积和
长度为
N
1
N_1
N1的序列
x
1
(
n
)
x_1(n)
x1(n)与长度为
N
2
N_2
N2的序列
x
2
(
n
)
x_2(n)
x2(n)的
L
L
L点圆周卷积和为
y
(
n
)
=
x
1
(
n
)
Ⓛ
x
2
(
n
)
=
x
2
(
n
)
Ⓛ
x
1
(
n
)
=
[
∑
m
=
0
l
−
1
x
1
(
m
)
x
2
(
(
n
−
m
)
)
L
]
R
L
(
n
)
,
L
≥
m
a
x
[
N
1
,
N
2
]
\begin{aligned}y(n)=&x_1(n)Ⓛx_2(n)=x_2(n)Ⓛx_1(n)\\=&[\sum_{m=0}^{l-1}x_1(m)x_2((n-m))_L]R_L(n),L\ge max[N_1,N_2]\end{aligned}
y(n)==x1(n)Ⓛx2(n)=x2(n)Ⓛx1(n)[m=0∑l−1x1(m)x2((n−m))L]RL(n),L≥max[N1,N2]
对于每一个n,累加另一个序列的圆周左移哑变量m位序列与序列在哑变量处取值的乘积。
要求将 x 1 ( n ) x_1(n) x1(n)和 x 2 ( n ) x_2(n) x2(n)的值补零到 L L L点长
与周期卷积和的关系
L L L点圆周卷积和为以 L L L为周期的周期卷积和的主值序列
与线性卷积和的关系(相互之间的求法)
由线性卷积和求圆周卷积和
L L L点圆周卷积和为线性卷积和的以 L L L为周期延拓序列的主值区间
线性卷积和
y
l
(
n
)
=
x
1
(
n
)
∗
x
2
(
n
)
=
∑
m
=
0
N
1
−
1
x
1
(
m
)
x
2
(
n
−
m
)
y_l(n)=x_1(n)*x_2(n)=\sum_{m=0}^{N_1-1}x_1(m)x_2(n-m)
yl(n)=x1(n)∗x2(n)=m=0∑N1−1x1(m)x2(n−m)
圆周卷积和
y
(
n
)
=
x
1
(
n
)
Ⓛ
x
2
(
n
)
=
[
∑
m
=
0
L
−
1
x
1
(
m
)
x
2
(
(
n
−
m
)
)
L
]
R
L
(
n
)
=
[
∑
m
=
0
L
−
1
∑
r
=
−
∞
∞
x
1
(
m
)
x
2
(
n
+
r
L
−
m
)
]
R
L
(
n
)
=
[
∑
r
=
−
∞
∞
∑
m
=
0
L
−
1
x
1
(
m
)
x
2
(
n
+
r
L
−
m
)
]
R
L
(
n
)
=
[
∑
r
=
−
∞
∞
∑
m
=
0
N
1
−
1
x
1
(
m
)
x
2
(
n
+
r
L
−
m
)
]
R
L
(
n
)
=
[
∑
r
=
−
∞
∞
y
l
(
n
+
r
L
)
]
R
L
(
n
)
\begin{aligned}y(n)=&x_1(n)Ⓛx_2(n)=[\sum_{m=0}^{L-1}x_1(m)x_2((n-m))_L]R_L(n)\\=&[\sum_{m=0}^{L-1}\sum_{r=-\infty}^{\infty}x_1(m)x_2(n+rL-m)]R_L(n)\\=&[\sum_{r=-\infty}^{\infty}\sum_{m=0}^{L-1}x_1(m)x_2(n+rL-m)]R_L(n)\\=&[\sum_{r=-\infty}^{\infty}\sum_{m=0}^{N_1-1}x_1(m)x_2(n+rL-m)]R_L(n)\\=&[\sum_{r=-\infty}^{\infty}y_l(n+rL)]R_L(n)\end{aligned}
y(n)=====x1(n)Ⓛx2(n)=[m=0∑L−1x1(m)x2((n−m))L]RL(n)[m=0∑L−1r=−∞∑∞x1(m)x2(n+rL−m)]RL(n)[r=−∞∑∞m=0∑L−1x1(m)x2(n+rL−m)]RL(n)[r=−∞∑∞m=0∑N1−1x1(m)x2(n+rL−m)]RL(n)[r=−∞∑∞yl(n+rL)]RL(n)
-
先进行周期为 L L L的周期延拓,
(考虑到 L L L仅要求 L ≥ m a x ( N 1 , N 2 ) L\ge max(N_1,N_2) L≥max(N1,N2),而线性卷积和 y l ( n ) y_l(n) yl(n)的长度为 N 1 + N 2 − 1 N_1+N_2-1 N1+N2−1,可能产生混叠。) -
再取主值区间得圆周卷积和
由圆周卷积和求线性卷积和
由从线性卷积和到圆周卷积和的讨论得到了只有当 L ≥ N 1 + N 2 − 1 L\ge N_1+N_2-1 L≥N1+N2−1时,不发生混叠,线性卷积和的周期延拓序列的主值区间才是圆周卷积和。
y ( n ) = [ ∑ r = − ∞ ∞ y l ( n + r L ) ] R L ( n ) y(n)=[\sum_{r=-\infty}^{\infty}y_l(n+rL)]R_L(n) y(n)=[r=−∞∑∞yl(n+rL)]RL(n)
如何无损的恢复出 y l ( n ) y_l(n) yl(n)?
- 如果 ∑ r = − ∞ ∞ y l ( n + r L ) \sum_{r=-\infty}^{\infty}y_l(n+rL) ∑r=−∞∞yl(n+rL)没有混叠,也就是 L ≥ N 1 + N 2 − 1 L\ge N_1+N_2-1 L≥N1+N2−1,
也就是说考虑到 R L ( n ) R_L(n) RL(n)的截断在主值区间上只有一个 y l ( n ) y_l(n) yl(n)存在,而这个 y l ( n ) y_l(n) yl(n)的长度是 M = N 1 + N 2 − 1 M=N_1+N_2-1 M=N1+N2−1,
那么主值区间的 0 ∼ N 1 + N 2 − 2 0\sim N_1+N_2-2 0∼N1+N2−2代表线性卷积和序列,剩下的为零点,可以无损恢复出整个线性卷积和序列
- 若 L < N 1 + N 2 − 1 = M L<N_1+N_2-1=M L<N1+N2−1=M,发生了混叠,并且序列被截断,必然不能恢复出完整的线性卷积和序列
考虑到相邻的周期序列的影响(实际上右移得到的序列已经在主值区间之外,实际上不会对区间内的结果产生混叠)
此时只有 M − L ∼ L − 1 M-L\sim L-1 M−L∼L−1区间代表线性卷积和
应用 实际运用中DFT计算圆周卷积和
在实际中利用FFT(DFT)来加速圆周卷积和的计算:
先将需要计算线性卷积的两个序列补零至足满足不发生混叠的长度,同时需要满足FFT的要求 L = 2 m ≥ N 1 + N 2 − 1 L=2^m\ge N_1+N_2-1 L=2m≥N1+N2−1 。
然后对两补零序列做L点FFT。
将FFT结果相乘根据卷积定理就得到了圆周卷积和主值区间的FFT,对其进行IFFT就得到了卷积结果。
圆周相关,利用DFT计算线性相关
线性相关定义式
r
x
y
(
m
)
=
∑
n
=
−
∞
∞
x
(
n
)
y
(
n
−
m
)
=
∑
n
=
−
∞
∞
x
(
n
)
y
[
−
(
m
−
n
)
]
=
x
(
m
)
∗
y
(
−
m
)
\begin{aligned}r_{xy}(m)=&\sum_{n=-\infty}^{\infty}x(n)y(n-m)=\sum_{n=-\infty}^{\infty}x(n)y[-(m-n)]\\=&x(m)*y(-m)\end{aligned}
rxy(m)==n=−∞∑∞x(n)y(n−m)=n=−∞∑∞x(n)y[−(m−n)]x(m)∗y(−m)
则
R
x
y
(
e
j
w
)
=
X
(
e
j
w
)
Y
(
e
j
w
)
R_{xy}(e^{jw})=X(e^{jw})Y(e^{jw})
Rxy(ejw)=X(ejw)Y(ejw)
对于实序列(回翻前面共轭实序列DFT的性质),再共轭有
X
(
N
−
k
)
=
X
∗
(
k
)
X(N-k)=X^*(k)
X(N−k)=X∗(k)
频域抽样
R
x
y
(
k
)
=
X
(
k
)
Y
(
N
−
k
)
=
X
(
k
)
Y
∗
(
k
)
R_{xy}(k)=X(k)Y(N-k)=X(k)Y^*(k)
Rxy(k)=X(k)Y(N−k)=X(k)Y∗(k)
圆周相关定义式
x
ˉ
x
y
(
m
)
=
∑
n
=
0
N
−
1
x
(
n
)
y
(
(
n
−
m
)
)
N
R
N
(
m
)
=
∑
n
=
0
N
−
1
x
(
n
)
y
(
(
−
(
m
−
n
)
)
N
R
N
(
m
)
=
∑
n
=
0
N
−
1
x
(
n
)
y
(
N
−
(
m
−
n
)
)
=
∑
n
=
0
N
−
1
x
(
n
)
y
(
−
(
m
−
n
−
N
)
)
=
x
(
m
)
Ⓝ
y
(
N
−
m
)
\begin{aligned}\bar{x}_{xy}(m)=&\sum_{n=0}^{N-1}x(n)y((n-m))_NR_N(m)=\sum_{n=0}^{N-1}x(n)y((-(m-n))_NR_N(m)\\=&\sum_{n=0}^{N-1}x(n)y(N-(m-n))=\sum_{n=0}^{N-1}x(n)y(-(m-n-N))\\=&x(m)Ⓝy(N-m)\end{aligned}
xˉxy(m)===n=0∑N−1x(n)y((n−m))NRN(m)=n=0∑N−1x(n)y((−(m−n))NRN(m)n=0∑N−1x(n)y(N−(m−n))=n=0∑N−1x(n)y(−(m−n−N))x(m)Ⓝy(N−m)
更深入理解卷积表达式的符号含义,哑变量和“移位”变量的区别
圆周相关定理
R
x
y
(
k
)
=
X
(
k
)
Y
(
N
−
k
)
=
X
(
k
)
Y
∗
(
k
)
R_{xy}(k)=X(k)Y(N-k)=X(k)Y^*(k)
Rxy(k)=X(k)Y(N−k)=X(k)Y∗(k)则
r
ˉ
x
y
(
m
)
=
I
D
F
T
[
R
x
y
(
k
)
]
=
∑
n
=
0
N
−
1
x
(
n
)
y
(
(
n
−
m
)
)
N
R
N
(
m
)
\bar{r}_{xy}(m)=IDFT[R_{xy}(k)]=\sum_{n=0}^{N-1}x(n)y((n-m))_NR_N(m)
rˉxy(m)=IDFT[Rxy(k)]=n=0∑N−1x(n)y((n−m))NRN(m)
利用FFT快速计算线性相关,在快速卷积基础上要先确定
r
x
y
(
0
)
r_{xy}(0)
rxy(0)的位置(可利用卷积和定位),
然后求出 R x y ( k ) = X ( k ) Y ∗ ( k ) R_{xy}(k)=X(k)Y^*(k) Rxy(k)=X(k)Y∗(k),再做IFFT
时域和频域上的圆周卷积和定理
时域
y
(
n
)
=
x
1
(
n
)
Ⓛ
x
2
(
n
)
y(n)=x_1(n)Ⓛx_2(n)
y(n)=x1(n)Ⓛx2(n)
Y
(
k
)
=
X
1
(
k
)
X
2
(
k
)
,
L
点
Y(k)=X_1(k)X_2(k),L点
Y(k)=X1(k)X2(k),L点
注意:将序列补零到
L
L
L点,然后做
L
L
L点DFT
频域
y
(
n
)
=
x
1
(
n
)
x
2
(
n
)
y(n)=x_1(n)x_2(n)
y(n)=x1(n)x2(n)
Y
(
k
)
=
1
L
X
1
(
k
)
Ⓛ
X
2
(
k
)
Y(k)=\frac{1}{L}X_1(k)ⓁX_2(k)
Y(k)=L1X1(k)ⓁX2(k)
注意:同样需要补零
频域抽样理论
频域抽样定理
讨论 频域抽样,由 X ( k ) X(k) X(k)重构时间序列 x ( n ) x(n) x(n)
本文此前已经从频域抽样的角度推导出DFS的表达式。
但是这只是对DFS表达式的一个论证,并不是数学上的证明.
其中在使用时域和频域的抽样点数(周期)相等时,已经迎合了DFS的定义,默认了DFS是针对周期序列
x
~
(
n
)
\tilde{x}(n)
x~(n)的,然后针对时域序
列的周期性来导出频域上的周期性。
下面就来探讨对于一个非周期序列,如果在频域上抽样,能否从抽样后的频谱来恢复出这个非周期序列。这里为了区分于此前的讨论,选择从z变换的角度来进行。
对一个绝对可和的非周期序列 x ( n ) x(n) x(n),其z变换为 X ( z ) = ∑ n = − ∞ ∞ x ( n ) z − n X(z)=\sum_{n=-\infty}^{\infty}x(n)z^{-n} X(z)=n=−∞∑∞x(n)z−n由于 ∑ n = − ∞ ∞ ∣ x ( n ) ∣ < ∞ \sum_{n=-\infty}^{\infty}|x(n)|<\infty n=−∞∑∞∣x(n)∣<∞故在单位圆上z变换收敛,z变换收敛域包含单位圆,
在单位圆上进行抽样,
w
w
w取值
0
∼
2
π
0\sim 2\pi
0∼2π,在其上抽取
N
N
N个频率,得
X
~
(
k
)
=
X
(
z
)
∣
z
=
e
j
2
π
N
k
=
∑
n
=
−
∞
∞
x
(
n
)
W
N
k
n
\tilde{X}(k)=X(z)|_{z=e^{j\frac{2\pi}{N}k}}=\sum_{n=-\infty}^{\infty}x(n)W_N^{kn}
X~(k)=X(z)∣z=ejN2πk=n=−∞∑∞x(n)WNkn
就得到了周期性序列
X
~
(
k
)
\tilde{X}(k)
X~(k),也就是DFS
求时域序列
x
~
N
(
n
)
=
I
D
F
S
[
X
~
(
k
)
]
=
1
N
∑
k
=
0
N
−
1
X
~
(
k
)
W
N
−
k
n
=
1
N
∑
k
=
0
N
−
1
[
∑
m
=
−
∞
∞
x
(
m
)
W
N
k
m
]
W
N
−
k
n
=
∑
m
=
−
∞
∞
x
(
m
)
[
1
N
∑
k
=
0
N
−
1
W
N
(
m
−
n
)
k
]
=
∑
r
=
−
∞
∞
x
(
n
+
r
N
)
\begin{aligned}\tilde{x}_N(n)=&IDFS[\tilde{X}(k)]=\frac{1}{N}\sum{k=0}^{N-1}\tilde{X}(k)W_N^{-kn}\\=&\frac{1}{N}\sum_{k=0}^{N-1}[\sum_{m=-\infty}^{\infty}x(m)W_N^{km}]W_N^{-kn}\\=&\sum_{m=-\infty}^{\infty}x(m)[\frac{1}{N}\sum_{k=0}^{N-1}W_N^{(m-n)k}]\\=&\sum_{r=-\infty}^{\infty}x(n+rN)\end{aligned}
x~N(n)====IDFS[X~(k)]=N1∑k=0N−1X~(k)WN−knN1k=0∑N−1[m=−∞∑∞x(m)WNkm]WN−knm=−∞∑∞x(m)[N1k=0∑N−1WN(m−n)k]r=−∞∑∞x(n+rN)
利用了性质
1
N
∑
k
=
0
N
−
1
W
N
(
m
−
n
)
k
=
{
1
,
m
=
n
+
r
N
,
r
为整数
0
,
其他
m
\frac{1}{N}\sum_{k=0}^{N-1}W_N^{(m-n)k}=\begin{cases}1,&m=n+rN,r为整数\\0,&其他m\end{cases}
N1k=0∑N−1WN(m−n)k={1,0,m=n+rN,r为整数其他m
可以看到频域抽样后进行IDFS得到的时域延拓序列的周期为
N
N
N,与频域抽样点数相同。
(此前:时域序列的一个周期内的时域抽样点数为
N
N
N)。频域的抽样造成了时域的周期延拓
- 若 x ( n ) x(n) x(n)为无限长序列,频域抽样后恢复出来的时域序列必定混叠失真。频域抽样点数 N N N越多(抽样越密),混叠失真越小
- 若 x ( n ) x(n) x(n)为有限长序列,长度为 M M M点,当频域抽样点数 N < M N<M N<M时,只有在 M − N ≤ n ≤ N − 1 M-N\le n\le N-1 M−N≤n≤N−1没有失真。这一点与《由圆周卷积和求线性卷积和》的结论类似。圆周卷积和得到的结果是线性卷积和以圆周卷积点数 L L L为周期的混叠;频域抽样后恢复的时域序列是非周期时域序列以频域抽样点数 N N N为周期的混叠。
如果序列的长度为 M M M点,对其DTFT X ( e j w ) X(e^{jw}) X(ejw)在 0 ≤ w ≤ 2 π 0\le w\le 2\pi 0≤w≤2π上等间隔抽样(抽样点不包括 2 π 2\pi 2π),得到 X ~ ( k ) \tilde{X}(k) X~(k),只有当抽样点数 N ≥ M N\ge M N≥M时,才能由 X ~ ( k ) \tilde{X}(k) X~(k)无失真恢复 x ( n ) x(n) x(n)
x ( n ) = x N ( n ) = x ~ N ( n ) R N ( n ) = ∑ r = − ∞ ∞ x ( n + r N ) R N ( n ) , N ≥ M \begin{aligned}x(n)=&x_N(n)=\tilde{x}_N(n)R_N(n)\\=&\sum_{r=-\infty}^{\infty}x(n+rN)R_N(n),N\ge M\end{aligned} x(n)==xN(n)=x~N(n)RN(n)r=−∞∑∞x(n+rN)RN(n),N≥M
频域插值重构
由 X ( k ) X(k) X(k)插值重构 X ( z ) X(z) X(z)
将用
X
(
k
)
X(k)
X(k)IDFS形式表示的
x
(
n
)
x(n)
x(n)代入z变换公式得
X
(
z
)
=
∑
k
=
0
N
−
1
X
(
k
)
Φ
k
(
z
)
X(z)=\sum_{k=0}^{N-1}X(k)\Phi_k(z)
X(z)=k=0∑N−1X(k)Φk(z)
Φ
k
(
z
)
=
z
N
−
1
N
z
N
−
1
(
z
−
W
N
−
k
)
\Phi_k(z)=\frac{z^N-1}{Nz^{N-1}(z-W_N^{-k})}
Φk(z)=NzN−1(z−WN−k)zN−1
插值函数
Φ
k
(
z
)
\Phi_k(z)
Φk(z)的零点有N个,
z
r
=
W
N
−
r
=
e
j
2
π
N
r
,
r
=
0
,
1
,
⋯
,
k
,
⋯
,
N
−
1
z_r=W_N^{-r}=e^{j\frac{2\pi}{N}r},r=0,1,\cdots,k,\cdots,N-1
zr=WN−r=ejN2πr,r=0,1,⋯,k,⋯,N−1
(
z
r
N
=
e
j
2
π
r
=
1
z_r^N=e^{j2\pi r}=1
zrN=ej2πr=1)
极点一个为
z
=
W
N
−
k
z=W_N^{-k}
z=WN−k
那么当
r
=
k
r=k
r=k时零点和极点相互抵消,
Φ
k
(
z
)
\Phi_k(z)
Φk(z)仅在本抽样点
X
(
k
)
X(k)
X(k)上
z
=
e
j
2
π
N
k
z=e^{j\frac{2\pi}{N}k}
z=ejN2πk不为零,其它N-1个抽样点均为零
由 X ( k ) X(k) X(k)插值重构 X ( e j w ) X(e^{jw}) X(ejw)
代入
z
=
e
j
w
z=e^{jw}
z=ejw,
X
(
e
j
w
)
=
∑
k
=
0
N
−
1
X
(
e
j
w
)
Φ
k
(
e
j
w
)
X(e^{jw})=\sum_{k=0}^{N-1}X(e^{jw})\Phi_k(e^{jw})
X(ejw)=k=0∑N−1X(ejw)Φk(ejw)
插值函数
Φ
k
(
e
j
w
)
=
Φ
(
w
−
2
π
N
k
)
\Phi_k(e^{jw})=\Phi(w-\frac{2\pi}{N}k)
Φk(ejw)=Φ(w−N2πk)
Φ
(
w
)
=
1
N
N
w
2
s
i
n
w
2
e
−
j
(
N
−
1
)
w
2
=
1
N
R
N
(
e
j
w
)
\Phi(w)=\frac{1}{N}\frac{\frac{Nw}{2}}{sin\frac{w}{2}}e^{-j\frac{(N-1)w}{2}}=\frac{1}{N}R_N(e^{jw})
Φ(w)=N1sin2w2Nwe−j2(N−1)w=N1RN(ejw)
Φ
k
(
e
j
w
)
=
1
N
R
N
(
e
j
(
w
−
2
π
N
k
)
)
\Phi_k(e^{jw})=\frac{1}{N}R_N(e^{j(w-\frac{2\pi}{N}k)})
Φk(ejw)=N1RN(ej(w−N2πk))
X
(
e
j
w
)
=
∑
k
=
0
N
−
1
X
(
k
)
Φ
(
w
−
2
π
N
k
)
X(e^{jw})=\sum_{k=0}^{N-1}X(k)\Phi(w-\frac{2\pi}{N}k)
X(ejw)=k=0∑N−1X(k)Φ(w−N2πk)
在抽样频率
w
=
2
π
N
k
w=\frac{2\pi}{N}k
w=N2πk处其它抽样点
k
k
k的被
X
(
k
)
X(k)
X(k)加权的插值函数
Φ
k
(
e
j
w
)
\Phi_k(e^{jw})
Φk(ejw)的值为零,累加值为零,只有当前抽样点点值。
在抽样点之间的连续频率上的 X ( e j w ) X(e^{jw}) X(ejw)是各个抽样点的插值函数加权后的叠加
应用:DFT谱分析
DFT对模拟信号FT的逼近
连续非周期信号FT的DFT逼近
连续非周期的频谱也为连续非周期,在计算机上处理序列需要将时域和频域都进行离散化
- 时域抽样,原 X ( j Ω ) = ∫ − ∞ ∞ x ( t ) e − j Ω t d t X(j\Omega)=\int_{-\infty}^{\infty}x(t)e^{-j\Omega t}dt X(jΩ)=∫−∞∞x(t)e−jΩtdt,抽样后 X ( j Ω ) = T ∑ n = − ∞ ∞ x ( n T ) e − j Ω n T X(j\Omega)=T\sum_{n=-\infty}^{\infty}x(nT)e^{-j\Omega nT} X(jΩ)=T∑n=−∞∞x(nT)e−jΩnT(乘以T因为将积分写成求和时 d t = T dt=T dt=T,T为时域抽样间隔)
- 时域截断,从无限长序列截断成从 t = 0 t=0 t=0开始的长度为 T 0 = N T T_0=NT T0=NT的有限长序列, X ( j Ω ) = T ∑ n = 0 N − 1 x ( n T ) e − j Ω n T X(j\Omega)=T\sum_{n=0}^{N-1}x(nT)e^{-j\Omega nT} X(jΩ)=T∑n=0N−1x(nT)e−jΩnT。(频域离散后时域周期,防止混叠。问题:频谱泄漏)
- 频域抽样(离散化),在频域一个周期 f s f_s fs(时域离散,频域周期)取 N N N个抽样点 f s = N F 0 f_s=NF_0 fs=NF0,故抽样间隔(频域分辨率)为 F 0 F_0 F0。得到的关系与《时域与频域的周期和频率的关系》一样, Ω 0 T = Ω 0 f s = 2 π Ω 0 Ω s = 2 π F 0 f s = 2 π T T 0 = 2 π N \Omega_0T=\frac{\Omega_0}{f_s}=2\pi\frac{\Omega_0}{\Omega_s}=2\pi\frac{F_0}{f_s}=\frac{2\pi T}{T_0}=\frac{2\pi}{N} Ω0T=fsΩ0=2πΩsΩ0=2πfsF0=T02πT=N2π
X
(
j
k
Ω
0
)
=
T
∑
n
=
0
N
−
1
x
(
n
T
)
e
−
j
k
Ω
0
n
T
=
T
⋅
D
F
T
[
x
(
n
)
]
X(jk\Omega_0)=T\sum_{n=0}^{N-1}x(nT)e^{-jk\Omega_0 nT}=T\cdot DFT[x(n)]
X(jkΩ0)=Tn=0∑N−1x(nT)e−jkΩ0nT=T⋅DFT[x(n)]
同理
x
(
n
)
=
x
(
n
T
)
=
1
T
⋅
I
D
F
T
[
X
(
j
k
Ω
0
)
]
x(n)=x(nT)=\frac{1}{T}\cdot IDFT[X(jk\Omega_0)]
x(n)=x(nT)=T1⋅IDFT[X(jkΩ0)]
连续周期信号FS的DFS逼近
时域连续周期,频域离散非周期
- 时域抽样,一个周期(有限长) T 0 T_0 T0内抽 N N N个点 T 0 = N T T_0=NT T0=NT,原 X ( j k Ω 0 ) = 1 T 0 ∫ 0 T 0 x ( t ) e − j k Ω 0 t d t X(jk\Omega_0)=\frac{1}{T_0}\int_{0}^{T_0}x(t)e^{-jk\Omega_0 t}dt X(jkΩ0)=T01∫0T0x(t)e−jkΩ0tdt,抽样后 X ( j k Ω 0 ) = T T 0 ∑ 0 N − 1 x ( n T ) e − j k Ω 0 n T = 1 N ∑ n = 0 N − 1 x ( n ) e − j 2 π N n k X(jk\Omega_0)=\frac{T}{T_0}\sum_{0}^{N-1}x(nT)e^{-jk\Omega_0 nT}=\frac{1}{N}\sum_{n=0}^{N-1}x(n)e^{-j\frac{2\pi}{N}nk} X(jkΩ0)=T0T∑0N−1x(nT)e−jkΩ0nT=N1∑n=0N−1x(n)e−jN2πnk(使用了FS的频域/时域的周期/频率关系)。注意抽样的条件,抽样后需要保证仍然是周期序列,已经在第一篇博客《模拟信号抽样成数字信号》时探讨过。
- 频域截断(DFS的频域周期,FS非周期,都是离散的),截断长度为DFS频域的一个周期(必须的,否则无法重构时域信号,造成频谱泄漏)
X
(
n
T
)
=
∑
k
=
0
N
−
1
X
(
j
k
Ω
0
)
e
j
k
Ω
0
n
T
=
N
⋅
1
N
D
F
S
[
x
(
n
)
]
X(nT)=\sum_{k=0}^{N-1}X(jk\Omega_0)e^{jk\Omega_0nT}=N\cdot\frac{1}{N}DFS[x(n)]
X(nT)=∑k=0N−1X(jkΩ0)ejkΩ0nT=N⋅N1DFS[x(n)]
X ( j k Ω 0 ) = 1 N ⋅ D F S [ x ( n ) ] X(jk\Omega_0)=\frac{1}{N}\cdot DFS[x(n)] X(jkΩ0)=N1⋅DFS[x(n)]
x ( n ) = x ( n T ) = N ⋅ I D F S [ X ( j k Ω 0 ) ] x(n)=x(nT)=N\cdot IDFS[X(jk\Omega_0)] x(n)=x(nT)=N⋅IDFS[X(jkΩ0)]
DFT对模拟信号作谱分析
参量选择
抽样频率 f s f_s fs
要求
f
s
>
2
f
H
f_s>2f_H
fs>2fH(奈奎斯特采样定理),因此有
折叠频率
f
s
2
\frac{f_s}{2}
2fs为能分析模拟信号的最高频率,对应数字频率
w
w
w
频率分辨率 F 0 F_0 F0
F 0 = f s N = 1 T 0 F_0=\frac{f_s}{N}=\frac{1}{T_0} F0=Nfs=T01想提高频率分辨率应增加抽样点数 N N N,在抽样频率不变的情况下就是增加记录时间 T 0 T_0 T0,增加了有效数据长度。
补零点从而增加抽样点数 N N N是没有用的,可以减小栅栏效应但不能提高分辨率。
DFT分析非周期模拟信号需要将时域信号截断,对应加一个长度为 N N N的窗, X N ( n ) = x ( n ) w N ( n ) X_N(n)=x(n)w_N(n) XN(n)=x(n)wN(n)若加矩形窗,其主瓣宽度为 4 π N \frac{4\pi}{N} N4π,取主瓣宽度的一半作为频率分辨率 Δ w = 2 π N = 2 π Δ f f s \Delta w=\frac{2\pi}{N}=\frac{2\pi\Delta f}{f_s} Δw=N2π=fs2πΔf Δ f = F 0 \Delta f=F_0 Δf=F0
遇到的问题
频谱混叠失真
抽样频率不满足抽样定理要求,或者由频谱泄漏造成
频谱泄漏
加窗后原频谱与窗函数频谱卷积造成频带宽度增加,使信号分辨率降低
栅栏效应
增加频域抽样密度可以减小栅栏效应
全文完,笔记完成于2023年末,2024年1月23日排版后发布