数字信号处理翻转课堂笔记7——循环移位与卷积&DFT的共轭对称性&频率域采样定理

原创

已于 2023-11-28 14:55:25 修改 · 8.8k 阅读

标签

#matlab #学习 #信号处理 #笔记

于 2022-10-12 20:01:18 首次发布

本文深入探讨数字信号处理中的关键概念，包括循环移位、循环卷积定理、DFT的共轭对称性及频率域采样定理。通过MATLAB编程实现理论与实践的结合，演示如何计算循环卷积及利用DFT处理信号。

数字信号处理第七次翻转课堂

Flipped Classroom7 of DSP

对应教材：《数字信号处理（第五版）》西安电子科技大学出版社，高西全、丁玉美著

一、要点

1、循环移位性质；
2、循环卷积定理（重点）；
3、DFT的共轭对称性；
4、频率域采样定理（难点）。

二、问题与解答

1、“循环移位”也称为“桶形移位”或者“圆周移位”，请查阅相关资料，简述“桶形移位”或者“圆周移位”方式的直观解释（图形化解释）。

2、任意给定一个N点序列x(n)（N的大小自定），求 x((1-n))_N R_N (n)，要求给出必要求解步骤。（先平移翻转，再周期延拓；或者先周期延拓，再平移翻转）

3、在这里插入图片描述

4、
在这里插入图片描述

5、
在这里插入图片描述

6、根据序列傅里叶变换的对称性、傅里叶变换与DFT的对应关系，讨论DFT的共轭对称特性与傅里叶变换的共轭对称特性之间的联系和区别。

7、参考教材例3.2.2，任意设两个N点实序列，用MATLAB编程实现：通过一次N点DFT，同时计算两个N点实序列的DFT。分析实现原理，给出程序实现思路和步骤（画出程序流程图，注意流程图的规范性），验证计算结果的正确性。

8、对于4点矩形序列R4(n)，分别求出它的4点、8点DFT。然后分别以这两个DFT为基础，根据频域采样内插公式（教材中式（3.3.7）），采用数值计算方式（实际上是用内插公式求出更多ω点的DFT值），恢复出R4(n)在 [0,2π]区间的连续幅度频谱曲线。基于频域采样定理，对两个不同点数DFT的恢复结果进行对比分析，给出相应的结论。

1、循环移位的直观解释

“循环移位”也称为“桶形移位”或者“圆周移位”，请查阅相关资料，简述“桶形移位”或者“圆周移位”方式的直观解释（图形化解释）。

在这里插入图片描述

这样的移位具有循环的特性，即x（n）向右移动m位时，右边超出（N-1）的m个样值又从左边依次填补了空位。如果把序列x（n）排列在一个N等分的圆周上，N个点首尾相接，上面所述的移位可以表示为x(n)在圆周上旋转m位，如图所示。当有限长序列进行任意位数的圆周移位后，求序列的DFT时取值范围仍保持在0~N-1不变。

2、循环卷积的图形化过程

任意给定一个N点序列x(n)（N的大小自定），求 x((1-n))_N R_N (n)，要求给出必要求解步骤。（先平移翻转，再周期延拓；或者先周期延拓，再平移翻转）

在这里插入图片描述

3、不同点数的循环卷积与线性卷积

在这里插入图片描述

%%      代码：
x1=[1 2 3 4];
x2=[1 1 1 1];
n1=0:1:7;
n2=0:1:3;
n3=0:1:14;
n4=0:1:6;
x3=[1 0 0 0 0 4 3 2 ;       %对x1的循环移位拓展
    2 1 0 0 0 0 4 3 ;
    3 2 1 0 0 0 0 4 ;
    4 3 2 1 0 0 0 0 ;
    0 4 3 2 1 0 0 0 ;
    0 0 4 3 2 1 0 0 ;
    0 0 0 4 3 2 1 0 ;
    0 0 0 0 4 3 2 1];      
x4=[1 1 1 1 0 0 0 0]';      %行变列
x5=[1 4 3 2 ;
    2 1 4 3 ;
    3 2 1 4 ;
    4 3 2 1];
x6=[1 1 1 1]';
x7=[1 2 3 4 0 0 0 <