数字信号处理翻转课堂笔记3——序列的傅立叶变换（DTFT）

数字信号处理翻转课堂笔记3

The Flipped Classroom3 of DSP

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

一、要点

1、序列傅里叶变换的定义及其波形；
2、典型序列（矩形序列）频谱的特点；
3、序列的傅里叶变换与模拟信号、采样信号傅里叶变换之间的联系；
4、模拟频率与数字频率的对应关系。

二、问题与解答

1、

2、设一个模拟矩形脉冲信号x(t)如图1所示，分别用1kHz和2kHz的采样频率对其进行理想采样。回顾模拟信号频域分析的相关知识和方法，用MATLAB分别画出该模拟矩形脉冲的幅度频谱、理想采样信号（含1kHz和2kHz采样的信号，注意理想采样信号仍然属于模拟信号，它是一个冲激串）的幅度频谱（根据模拟信号的频谱公式画图）。再分别对1kHz和2kHz采样结果进行离散化，得到两个长度不同（长度分别为10、20）的矩形序列，分别用MATLAB画出这两个序列的傅里叶变换（幅度频谱图，根据矩形序列的频谱公式画图）。根据所得到的结果，讨论：①模拟矩形脉冲、理想采样信号、离散矩形序列这三种类型信号的幅度频谱形状有何联系和区别（暂不考虑频率轴的尺度差异，仅对比幅度变化曲线的形状）？②由不同采样频率所得到的采样信号和离散序列的幅度频谱有何区别？③结合幅度频谱图，总结三种信号幅度频谱的频率轴具有什么样的对应关系？采样频率的大小对于这种对应关系有什么影响？

3、上一题中，如果是用数字的方法来分析某种类型模拟信号（不一定是矩形脉冲）的频谱，即希望模拟信号、采样信号、对应的离散序列这三种类型信号的幅度频谱形状完全相同，则选择的模拟信号需要满足什么条件（提示：结合采样定理和频谱混叠的概念进行分析）？

4、结合第2、3题的结果和教材图2.4.1，总结模拟频率、采样频率、数字频率之间的对应关系，讨论如果用数字的方法来分析模拟信号的频谱，应该注意什么问题？从采样频率的选取、计算结果的转换等角度进行讨论。

1、不同长度的离散时间傅里叶变换

代码：

%%      N=4
n=4;                                            %取n值
w=[-1000:1:1000]*pi/500;                        %定义横轴的量
X=exp(-j*(n-1)*w/2).*(sin(w*n/2))./(sin(w/2));  %代入公式
A=abs(X)                                        %调用abs函数求Amplitude
P=angle(X)                                      %调用angle函数求Phase

subplot(2,1,1)                                  %画图
plot(w,A)
xlabel('频率')
title('幅度波形')
ylabel('幅度')

subplot(2,1,2)
plot(w,P)
xlabel('频率')
title('相位波形')
ylabel('相位')

%%      N=8
n=8;
w=[-1000:1:1000]*pi/500;
X=exp(-j*(n-1)*w/2).*(sin(w*n/2))./(sin(w/2)); 
A=abs(X)
P=angle(X)

subplot(2,1,1)
plot(w,A)
xlabel('频率')
title('幅度波形')
ylabel('幅度')

subplot(2,1,2)
plot(w,P)
xlabel('频率')
title('相位波形')
ylabel('相位')

%%      N=32
n=32;
w=[-1000:1:1000]*pi/500;
X=exp(-j*(n-1)*w/2).*(sin(w*n/2))./(sin(w/2)); 
A=abs(X)
P=angle(X)

subplot(2,1,1)
plot(w,A)
xlabel('频率')
title('幅度波形')
ylabel('幅度')

subplot(2,1,2)
plot(w,P)
xlabel('频率')
title('相位波形')
ylabel('相位')

运行结果：
N=4

N=8

N=8

分析：
1、矩形序列的频谱是以2π为周期的连续函数。
2、N越大，相位频谱越密集，幅度频谱越大，幅度频谱最大值为N。
3、长度为N的矩形序列，频谱在[-π，π]区间具有N-1个波峰。
4、N越大高频分量占比越小，变化越慢。

2、对模拟矩形脉冲信号进行1kHz和2kHz的采样

设一个模拟矩形脉冲信号x(t)如图1所示，分别用1kHz和2kHz的采样频率对其进行理想采样。回顾模拟信号频域分析的相关知识和方法，用MATLAB分别画出该模拟矩形脉冲的幅度频谱、理想采样信号（含1kHz和2kHz采样的信号，注意理想采样信号仍然属于模拟信号，它是一个冲激串）的幅度频谱（根据模拟信号的频谱公式画图）。再分别对1kHz和2kHz采样结果进行离散化，得到两个长度不同（长度分别为10、20）的矩形序列，分别用MATLAB画出这两个序列的傅里叶变换（幅度频谱图，根据矩形序列的频谱公式画图）。根据所得到的结果，讨论：①模拟矩形脉冲、理想采样信号、离散矩形序列这三种类型信号的幅度频谱形状有何联系和区别（暂不考虑频率轴的尺度差异，仅对比幅度变化曲线的形状）？②由不同采样频率所得到的采样信号和离散序列的幅度频谱有何区别？③结合幅度频谱图，总结三种信号幅度频谱的频率轴具有什么样的对应关系？采样频率的大小对于这种对应关系有什么影响？

代码：

%%      代码：
%模拟矩形脉冲信号的频谱
syms f;
Xa=0.01*sinc(0.01*f)*exp(-j*0.005*2*pi*f);
%图1所示模拟矩形脉冲的傅里叶变换
subplot(4,1,1);
ezplot(abs(Xa),[-2000,2000,0,0.01])
xlabel('frequency(Hz)')
ylabel('|Xa(jw)|');
title('模拟信号幅度频谱');

%采样频率1kHz时的理想采样信号的频谱
T=0.001;                        %采样间隔为0.001s 
XX=0;
for n=0:9
    XX=XX+exp(-j*2*pi*f*n*T);   %理想采样信号为冲激信号，其频谱为多个冲激信号频谱的叠加
end
subplot(4,1,2);
ezplot(abs(XX),[-2000,2000,0,10])
xlabel('frequency(Hz)')
title('1kHz采样信号幅度频谱');

%1kHz采样离散信号的频谱
T=0.001;                        %采样间隔为0.001s
n=-50:1:50;
x=heaviside(n*T+0.0001)-heaviside(n*T-0.0099);  %离散时间信号
K=500;k=0:1:K;w=4*pi*k/K;                       %w为数字频率
X=x*exp(-j*n'*w);               %计算离散时间傅立叶变换（序列的傅立叶变换）
X=abs(X);      
w=[-fliplr(w),w(2:K+1)];
X=[fliplr(X),X(2:K+1)];
subplot(4,1,3);
plot(w/pi,X);                                   %画出离散时间傅立叶变换
xlabel('frequency(radian/pi)');                 %横坐标为弧度
ylabel('|X1(ejw)|');title('1kHz采样所得离散信号的幅度频谱');

%下面为2kHz采样离散信号的频谱
T=0.0005;                                       %采样间隔为0.001s
n=-100:1:100;
x=heaviside(n*T+0.0001)-heaviside(n*T-0.0099);  %离散时间信号
K=500;k=0:1:K;w=2*pi*k/K;                       %w为数字频率
X=x*exp(-j*n'*w);               %计算离散时间傅立叶变换（序列的傅立叶变换）,
X=abs(X);      
w=[-fliplr(w),w(2:K+1)];
X=[fliplr(X),X(2:K+1)];
subplot(4,1,4);
plot(w/pi,X);                                   %画出离散时间傅立叶变换
xlabel('frequency(radian/pi)');                 %横坐标为弧度
ylabel('|X2(ejw)|');title('2kHz采样所得离散信号的幅度频谱');

运行结果：

分析：
①采样信号的频谱是模拟信号频谱以采样频率f进行周期性延拓，幅度为模拟信号的f倍，采样信号的频谱与离散信号的频谱形状相同。
②幅度谱的周期相同，为2π；最大幅值不同，为模拟信号的f倍。
③
模拟信号的频率轴为f，采样信号的频率轴为Ω，离散序列的频率轴为w。且有以上对应关系，采样频率Ωs对应数字频率2π。
由于采样频率Ωs对应数字频率2π，在相同的尺度下， ，当Ωs越大，数字频率轴表示的范围越小。

3、选择模拟信号的条件

上一题中，如果是用数字的方法来分析某种类型模拟信号（不一定是矩形脉冲）的频谱，即希望模拟信号、采样信号、对应的离散序列这三种类型信号的幅度频谱形状完全相同，则选择的模拟信号需要满足什么条件（提示：结合采样定理和频谱混叠的概念进行分析）？

结论：应使模拟信号的频带宽度尽可能的小。其中模拟信号绝对可和。
原因：由取样定理可知:只有设定取样频率大于等于2倍的模拟信号最大频率，即fs>2fm，才能保证频谱不发生混叠，进而恢复出原信号。
采样信号对应的离散序列取样点越多，所得频谱信息也更加接近于原模拟信号频谱，因此要使模拟信号、采样信号、对应的离散序列三种类型的信号幅度频谱完全相同，所以应使模拟信号的频带宽度尽可能的小。其中模拟信号绝对可和。

4、模拟频率、采样频率和数字频率之间的对应关系

结合第2、3题的结果和教材图2.4.1，总结模拟频率、采样频率、数字频率之间的对应关系，讨论如果用数字的方法来分析模拟信号的频谱，应该注意什么问题？从采样频率的选取、计算结果的转换等角度进行讨论。

（1）
模拟频率、采样频率、数字频率之间的对应关系：

模拟频率为每秒旋转的弧度（ f2pi rad/s）； 采样间隔为每隔Ts时间采样，采样频率即每秒钟采样fs次； 数字频率即Ts时间内的模拟频率，即f2pi*Ts。
数字频率是对模拟频率的采样，数字频率代表的是两个采样点之间的间隔。
数字频率为：w =2πΩ/ws。根据奈奎斯特定理，采样频率大于等于2倍的模拟频率，则有w<=π。
不同采样频率下的数字信号的数字频率是不同的，但是表示的是相同的信号，在理论上满足奈奎斯特定理，所包含的信息是完全相同的。
（2）应该注意混叠失真、截断效应、栅栏效应等问题。（具体知识在翻转汇报9和10中将会涉及）

三、反思总结

结合图像来理解比较直观，最好把图“刻”在脑子里。