数字信号处理实验报告3.pdfVIP

第三章离散时间系统的频域分析

例3.1计算离散时间傅立叶变换

%程序P3_1

%离散时间傅里叶变换的求取

clf;

%计算离散时间傅立叶变换的频率样本

w=-4*pi:8*pi/511:4*pi;

num=[21];den=[1-0.6];

h=freqz(num,den,w);

%PlottheDTFT

subplot(2,1,1)

plot(w/pi,real(h));grid

title(title(的实部)

xlabel(xlabel();

ylabel(振幅);

subplot(2,1,2)

plot(w/pi,imag(h));grid

title(title(的虚部)

xlabel(xlabel();

ylabel(振幅);

pause

subplot(2,1,1)

plot(w/pi,abs(h));grid

title(title(幅度谱)

xlabel(xlabel();

ylabel(振幅);

subplot(2,1,2)

plot(w/pi,angle(h));grid

title(相位谱相位谱)

xlabel(xlabel();

ylabel(以弧度为单位的相位);

Q3.1在程序P3.1中，计算离散时间傅里叶变换的原始系列是什么？MATLAB

命令pause的作用是什么？

原始序列：H(jw)=(2+e^(jw))/(1-0.6e^(-jw))

pause的作用：程序暂停,直至用户按任意一个按键.

Q3.2运行程序P3.1，求离散时间傅里叶变换的实部，虚部以及幅度和相位谱。

离散时间福利叶变换是w的周期函数吗？若是，周期是多少？描述这四个图形

表示的对称性。

离散时间福利叶变换是w的周期函数，周期为2π。

例3.2离散时间傅立叶变换的性质

%p3.2

%离散时间傅里叶变换的时移性质

clf;

w=-pi:2*pi/255:pi;

w0=0.4*pi;

D=10;

num=[123456789];

h1=freqz(num,1,w);

h2=freqz([zeros(1,D)num],1,w);

subplot(2,2,1);

plot(w/pi,abs(h1));grid

title(原序列的幅度谱);

subplot(2,2,2);

plot(w/pi,abs(h2));grid

title(时移后序列的幅度谱);

subplot(2,2,3)

plot(w/pi,angle(h1));grid

title(原序列的相位谱)

subplot(2,2,4)

plot(w/pi,angle(h2));grid

title(时移后序列的相位谱)

Q3.7运行修改后的程序并讨论你的结果。

根据结果可知当D增大时，序列的幅度谱不变；相位谱变密，周期缩短。

Q3.8选取不同的时移值重做习题Q3.7。

D=50

D=3

Q3.9选取两个改变了长度的序列以及两个不同的时移值，重做习题Q3.7。

D=30;num=[123456];

D=100;num=[12345678910];

%程序P3_3

%离散时间傅里叶变换的频移性质

clf;

w=-pi:2*pi/255:pi;wo=0.4*pi;

num1=[1357911131517];

L=length(num1);

h1=freqz(num1,1,w);

n=0:L-1;

num2=exp(wo*i*n).*num1;

h2=freqz(num2,1,w);

subplot(2,2,1)

plot(w/pi,abs(h1));grid

title(原序列的幅度谱)

subplot(2,2,2)

plot(w/pi,abs(h2));grid

title(频移后序列的幅度谱)

subplot(2,2,3)

plot(w/pi,angle(