数字信处理实验答案.docx

数字信处理实验答案.docx

《数字信处理实验答案.docx》由会员分享，可在线阅读，更多相关《数字信处理实验答案.docx（37页珍藏版）》请在冰点文库上搜索。

数字信处理实验答案

第十章上机实验

数字信号处理是一门理论和实际密切结合的课程，为深入掌握课程内容，最好在学习理论的同时，做习题和上机实验。

上机实验不仅可以帮助读者深入的理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。

本章在第二版的基础上编写了六个实验，前五个实验属基础理论实验，第六个属应用综合实验。

实验一系统响应及系统稳定性。

实验二时域采样与频域采样。

实验三用FFT对信号作频谱分析。

实验四IIR数字滤波器设计及软件实现。

实验五FIR数字滤波器设计与软件实现

实验六应用实验——数字信号处理在双音多频拨号系统中的应用

任课教师根据教学进度，安排学生上机进行实验。

建议自学的读者在学习完第一章后作实验一；在学习完第三、四章后作实验二和实验三；实验四IIR数字滤波器设计及软件实现在。

学习完第六章进行；实验五在学习完第七章后进行。

实验六综合实验在学习完第七章或者再后些进行；实验六为综合实验，在学习完本课程后再进行。

10.1　实验一:

系统响应及系统稳定性

1.实验目的

（1）掌握求系统响应的方法。

（2）掌握时域离散系统的时域特性。

（3）分析、观察及检验系统的稳定性。

2.实验原理与方法

在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。

已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，本实验仅在时域求解。

在计算机上适合用递推法求差分方程的解，最简单的方法是采用MATLAB语言的工具箱函数filter函数。

也可以用MATLAB语言的工具箱函数conv函数计算输入信号和系统的单位脉冲响应的线性卷积，求出系统的响应。

系统的时域特性指的是系统的线性时不变性质、因果性和稳定性。

重点分析实验系统的稳定性，包括观察系统的暂态响应和稳定响应。

系统的稳定性是指对任意有界的输入信号，系统都能得到有界的系统响应。

或者系统的单位脉冲响应满足绝对可和的条件。

系统的稳定性由其差分方程的系数决定。

实际中检查系统是否稳定，不可能检查系统对所有有界的输入信号，输出是否都是有界输出，或者检查系统的单位脉冲响应满足绝对可和的条件。

可行的方法是在系统的输入端加入单位阶跃序列，如果系统的输出趋近一个常数（包括零），就可以断定系统是稳定的[19]。

系统的稳态输出是指当

时，系统的输出。

如果系统稳定，信号加入系统后，系统输出的开始一段称为暂态效应，随n的加大，幅度趋于稳定，达到稳态输出。

注意在以下实验中均假设系统的初始状态为零。

3．实验内容及步骤

（1）编制程序，包括产生输入信号、单位脉冲响应序列的子程序，用filter函数或conv函数求解系统输出响应的主程序。

程序中要有绘制信号波形的功能。

（2）给定一个低通滤波器的差分方程为

输入信号

a）分别求出系统对和的响应序列，并画出其波形。

b）求出系统的单位冲响应，画出其波形。

（3）给定系统的单位脉冲响应为

用线性卷积法分别求系统h1（n）和h2（n）对的输出响应，并画出波形。

（4）给定一谐振器的差分方程为

令，谐振器的谐振频率为0.4rad。

a）用实验方法检查系统是否稳定。

输入信号为时，画出系统输出波形。

b）给定输入信号为

求出系统的输出响应，并画出其波形。

4．思考题

（1）如果输入信号为无限长序列，系统的单位脉冲响应是有限长序列，可否用线性卷积法求系统的响应?

如何求？

（2）如果信号经过低通滤波器，把信号的高频分量滤掉，时域信号会有何变化，用前面第一个实验结果进行分析说明。

5．实验报告要求

（1）简述在时域求系统响应的方法。

（2）简述通过实验判断系统稳定性的方法。

分析上面第三个实验的稳定输出的波形。

（3）对各实验所得结果进行简单分析和解释。

（4）简要回答思考题。

（5）打印程序清单和要求的各信号波形。

10.1.2实验程序清单

%实验1：

系统响应及系统稳定性

closeall;clearall

%======内容1：

调用filter解差分方程，由系统对u（n）的响应判断稳定性======

A=[1,-0.9];B=[0.05,0.05];%系统差分方程系数向量B和A

x1n=[11111111zeros（1,50）];%产生信号x1（n）=R8（n）

x2n=ones（1,128）;%产生信号x2（n）=u（n）

hn=impz（B,A,58）;%求系统单位脉冲响应h（n）

subplot（2,2,1）;y='h（n）';stem（hn）;%调用函数stem绘图

title（'（a）系统单位脉冲响应h（n）'）;boxon

y1n=filter（B,A,x1n）;%求系统对x1（n）的响应y1（n）

subplot（2,2,2）;y='y1（n）';stem（y1n）;

title（'（b）系统对R8（n）的响应y1（n）'）;boxon

y2n=filter（B,A,x2n）;%求系统对x2（n）的响应y2（n）

subplot（2,2,4）;y='y2（n）';stem（y2n）;

title（'（c）系统对u（n）的响应y2（n）'）;boxon

%===内容2：

调用conv函数计算卷积============================

x1n=[11111111];%产生信号x1（n）=R8（n）

h1n=[ones（1,10）zeros（1,10）];

h2n=[12.52.51zeros（1,10）];

y21n=conv（h1n,x1n）;

y22n=conv（h2n,x1n）;

figure

（2）

subplot（2,2,1）;y='h1（n）';stem（h1n）;%调用函数stem绘图

title（'（d）系统单位脉冲响应h1（n）'）;boxon

subplot（2,2,2）;y='y21（n）';stem（y21n）;

title（'（e）h1（n）与R8（n）的卷积y21（n）'）;boxon

subplot（2,2,3）;y='h2（n）';stem（h2n）;%调用函数stem绘图

title（'（f）系统单位脉冲响应h2（n）'）;boxon

subplot（2,2,4）;y='y22（n）';stem（y22n）;

title（'（g）h2（n）与R8（n）的卷积y22（n）'）;boxon

%=========内容3：

谐振器分析========================

un=ones（1,256）;%产生信号u（n）

n=0:

255;

xsin=sin（0.014*n）+sin（0.4*n）;%产生正弦信号

A=[1,-1.8237,0.9801];B=[1/100.49,0,-1/100.49];%系统差分方程系数向量B和A

y31n=filter（B,A,un）;%谐振器对u（n）的响应y31（n）

y32n=filter（B,A,xsin）;%谐振器对u（n）的响应y31（n）

figure（3）

subplot（2,1,1）;y='y31（n）';stem（y31n）;

title（'（h）谐振器对u（n）的响应y31（n）'）;boxon

subplot（2,1,2）;y='y32（n）';stem（y32n）;

title（'（i）谐振器对正弦信号的响应y32（n）'）;boxon

10.1.3实验程序运行结果及分析讨论

实验内容

（2）系统的单位冲响应、系统对和的响应序列分别如图（a）、（b）和（c）所示；

实验内容（3）系统h1（n）和h2（n）对的输出响应分别如图（e）和（g）所示；

实验内容（4）系统对和的响应序列分别如图（h）和（i）所示。

由图（h）可见，系统对的响应逐渐衰减到零，所以系统稳定。

由图（i）可见，系统对

的稳态响应近似为正弦序列，这一结论验证了该系统的谐振频率是0.4rad。

图10.1.1

10.1.4简答思考题

（1）

如果输入信号为无限长序列，系统的单位脉冲响应是有限长序列，可否用线性卷积法求系统的响应。

①对输入信号序列分段；②求单位脉冲响应h（n）与各段的卷积；③将各段卷积结果相加。

具体实现方法有第三章介绍的重叠相加法和重叠保留法。

、和的阶跃变化变得缓慢上升与下降。

10.2实验二时域采样与频域采样

10.2.1实验指导

1.实验目的

时域采样理论与频域采样理论是数字信号处理中的重要理论。

要求掌握模拟信号采样前后频谱的变化，以及如何选择采样频率才能使采样后的信号不丢失信息；要求掌握频率域采样会引起时域周期化的概念，以及频率域采样定理及其对频域采样点数选择的指导作用。

2.实验原理与方法

时域采样定理的要点是：

a）对模拟信号以间隔T进行时域等间隔理想采样，形成的采样信号的频谱是原模拟信号频谱以采样角频率（）为周期进行周期延拓。

公式为：

b）采样频率必须大于等于模拟信号最高频率的两倍以上，才能使采样信号的

频谱不产生频谱混叠。

利用计算机计算上式并不方便，下面我们导出另外一个公式，以便用计算机上进行实验。

理想采样信号和模拟信号之间的关系为：

对上式进行傅立叶变换，得到：

在上式的积分号内只有当时，才有非零值，因此：

上式中，在数值上＝，再将代入，得到：

上式的右边就是序列的傅立叶变换，即

上式说明理想采样信号的傅立叶变换可用相应的采样序列的傅立叶变换得到，只要将自变量ω用代替即可。

频域采样定理的要点是：

a）对信号x（n）的频谱函数X（ejω）在[0，2π]上等间隔采样N点，得到

则N点IDFT[]得到的序列就是原序列x（n）以N为周期进行周期延拓后的主值区序列，公式为：

b）由上式可知，频域采样点数N必须大于等于时域离散信号的长度M（即N≥M），才能使时域不产生混叠，则N点IDFT[]得到的序列就是原序列x（n）,即

=x（n）。

如果N>M，比原序列尾部多N-M个零点；如果N

与x（n）不相同。

在数字信号处理的应用中，只要涉及时域或者频域采样，都必须服从这两个采样理论的要点。

对比上面叙述的时域采样原理和频域采样原理，得到一个有用的结论，这两个采样理论具有对偶性：

“时域采样频谱周期延拓，频域采样时域信号周期延拓”。

因此放在一起进行实验。

3.实验内容及步骤

（1）时域采样理论的验证。

给定模拟信号，

式中A=444.128，=50π，=50πrad/s，它的幅频特性曲线如图10.2.1

图10.2.1的幅频特性曲线

现用DFT（FFT）求该模拟信号的幅频特性，以验证时域采样理论。

安照的幅频特性曲线，选取三种采样频率，即=1kHz，300Hz，200Hz。

观测时间选。

为使用DFT，首先用下面公式产生时域离散信号，对三种采样频率，采样序列按顺序用，，表示。

因为采样频率不同，得到的，，的长度不同，长度（点数）用公式计算。

选FFT的变换点数为M=64，序列长度不够64的尾部加零。

X（k）=FFT[x（n）]，k=0,1,2,3,-----,M-1

式中k代表的频率为。

要求：

编写实验程序，计算、和的幅度特性，并绘图显示。

观察分析频谱混叠失真。

（2）频域采样理论的验证。

给定信号如下：

编写程序分别对频谱函数在区间上等间隔采样32

和16点，得到：

再分别对进行32点和16点IFFT，得到：

分别画出、的幅度谱，并绘图显示x（n）、的波形，进行对比和分析，验证总结频域采样理论。

提示：

频域采样用以下方法容易变程序实现。

①直接调用MATLAB函数fft计算就得到在的32点频率域采样

②抽取的偶数点即可得到在的16点频率域采样，即。

○3当然也可以按照频域采样理论，先将信号x（n）以16为周期进行周期延拓，取其主值区（16点），再对其进行16点DFT（FFT）,得到的就是在

的16点频率域采样。

4．思考题：

如果序列x（n）的长度为M，希望得到其频谱在上的N点等间隔采样，当N

5.实验报告及要求

a）运行程序打印要求显示的图形，。

b）分析比较实验结果，简述由实验得到的主要结论

c）简要回答思考题

d）附上程序清单和有关曲线。

10.2.2实验程序清单

1时域采样理论的验证程序清单

%时域采样理论验证程序exp2a.m

Tp=64/1000;%观察时间Tp=64微秒

%产生M长采样序列x（n）

%Fs=1000;T=1/Fs;

Fs=1000;T=1/Fs;

M=Tp*Fs;n=0:

M-1;

A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5;

xnt=A*exp（-alph*n*T）.*sin（omega*n*T）;

Xk=T*fft（xnt,M）;%M点FFT[xnt）]

yn='xa（nT）';subplot（3,2,1）;

stem（xnt,yn）;%调用自编绘图函数stem绘制序列图

boxon;title（'（a）Fs=1000Hz'）;

k=0:

M-1;fk=k/Tp;

subplot（3,2,2）;plot（fk,abs（Xk））;title（'（a）T*FT[xa（nT）],Fs=1000Hz'）;

xlabel（'f（Hz）'）;ylabel（'幅度'）;axis（[0,Fs,0,1.2*max（abs（Xk））]）

%=================================================

%Fs=300Hz和Fs=200Hz的程序与上面Fs=1000Hz完全相同。

2频域采样理论的验证程序清单

%频域采样理论验证程序exp2b.m

M=27;N=32;n=0:

M;

%产生M长三角波序列x（n）

xa=0:

floor（M/2）;xb=ceil（M/2）-1:

-1:

0;xn=[xa,xb];

Xk=fft（xn,1024）;%1024点FFT[x（n）],用于近似序列x（n）的TF

X32k=fft（xn,32）;%32点FFT[x（n）]

x32n=ifft（X32k）;%32点IFFT[X32（k）]得到x32（n）

X16k=X32k（1:

2:

N）;%隔点抽取X32k得到X16（K）

x16n=ifft（X16k,N/2）;%16点IFFT[X16（k）]得到x16（n）

subplot（3,2,2）;stem（n,xn,'.'）;boxon

title（'（b）三角波序列x（n）'）;xlabel（'n'）;ylabel（'x（n）'）;axis（[0,32,0,20]）

k=0:

1023;wk=2*k/1024;%

subplot（3,2,1）;plot（wk,abs（Xk））;title（'（a）FT[x（n）]'）;

xlabel（'\omega/\pi'）;ylabel（'|X（e^j^\omega）|'）;axis（[0,1,0,200]）

k=0:

N/2-1;

subplot（3,2,3）;stem（k,abs（X16k）,'.'）;boxon

title（'（c）16点频域采样'）;xlabel（'k'）;ylabel（'|X_1_6（k）|'）;axis（[0,8,0,200]）

n1=0:

N/2-1;

subplot（3,2,4）;stem（n1,x16n,'.'）;boxon

title（'（d）16点IDFT[X_1_6（k）]'）;xlabel（'n'）;ylabel（'x_1_6（n）'）;axis（[0,32,0,20]）

k=0:

N-1;

subplot（3,2,5）;stem（k,abs（X32k）,'.'）;boxon

title（'（e）32点频域采样'）;xlabel（'k'）;ylabel（'|X_3_2（k）|'）;axis（[0,16,0,200]）

n1=0:

N-1;

subplot（3,2,6）;stem（n1,x32n,'.'）;boxon

title（'（f）32点IDFT[X_3_2（k）]'）;xlabel（'n'）;ylabel（'x_3_2（n）'）;axis（[0,32,0,20]）

10.2.3实验程序运行结果

1

图10.2.2

图10.3.3

该图验证了频域采样理论和频域采样定理。

对信号x（n）的频谱函数X（ejω）在[0，2π]上等间隔采样N=16时，N点IDFT[

]得到的序列正是原序列x（n）以16为周期进行周期延拓后的主值区序列：

由于N

与x（n）相同。

10.2.4简答思考题

先对原序列x（n）以N为周期进行周期延拓后取主值区序列，

再计算N点DFT则得到N点频域采样：

10.3　实验三：

用FFT对信号作频谱分析

10.3.1实验指导

1．实验目的

学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析

误差及其原因，以便正确应用FFT。

2.实验原理

用FFT对信号作频谱分析是学习数字信号处理的重要内容。

经常需要进行谱分析的信号是模拟信号和时域离散信号。

对信号进行谱分析的重要问题是频谱分辨率D和分析误差。

频谱分辨率直接和FFT的变换区间N有关，因为FFT能够实现的频率分辨率是

，因此要求

。

可以根据此式选择FFT的变换区间N。

误差主要来自于用FFT作频谱分析时，得到的是离散谱，而信号（周期信号除外）是连续谱，只有当N较大时离散谱的包络才能逼近于连续谱，因此N要适当选择大一些。

周期信号的频谱是离散谱，只有用整数倍周期的长度作FFT，得到的离散谱才能代表周期信号的频谱。

如果不知道信号周期，可以尽量选择信号的观察时间长一些。

对模拟信号进行谱分析时，首先要按照采样定理将其变成时域离散信号。

如果是模拟周期信号，也应该选取整数倍周期的长度，经过采样后形成周期序列，按照周期序列的谱分析进行。

3．实验步骤及内容

（1）对以下序列进行谱分析。

选择FFT的变换区间N为8和16两种情况进行频谱分析。

分别打印其幅频特性曲线。

并进行对比、分析和讨论。

（2）对以下周期序列进行谱分析。

选择FFT的变换区间N为8和16两种情况分别对以上序列进行频谱分析。

分别打印其幅频特性曲线。

并进行对比、分析和讨论。

（3）对模拟周期信号进行谱分析

选择采样频率，变换区间N=16,32,64三种情况进行谱分析。

分别打印其幅频特性，并进行分析和讨论。

4．思考题

（1）对于周期序列，如果周期不知道，如何用FFT进行谱分析？

（2）如何选择FFT的变换区间？

（包括非周期信号和周期信号）

（3）当N=8时，和的幅频特性会相同吗？

为什么？

N=16呢？

5．实验报告要求

（1）完成各个实验任务和要求。

附上程序清单和有关曲线。

（2）简要回答思考题。

10.3.2实验程序清单

%第10章实验3程序exp3.m

%用FFT对信号作频谱分析

clearall;closeall

%实验内容

（1）===================================================

x1n=[ones（1,4）];%产生序列向量x1（n）=R4（n）

M=8;xa=1:

（M/2）;xb=（M/2）:

-1:

1;x2n=[xa,xb];%产生长度为8的三角波序列x2（n）

x3n=[xb,xa];

X1k8=fft（x1n,8）;%计算x1n的8点DFT

X1k16=fft（x1n,16）;%计算x1n的16点DFT

X2k8=fft（x2n,8）;%计算x1n的8点DFT

X2k16=fft（x2n,16）;%计算x1n的16点DFT

X3k8=fft（x3n,8）;%计算x1n的8点DFT

X3k16=fft（x3n,16）;%计算x1n的16点DFT

%以下绘制幅频特性曲线

subplot（2,2,1）;mstem（X1k8）;%绘制8点DFT的幅频特性图

title（'（1a）8点DFT[x_1（n）]'）;xlabel（'ω/π'）;ylabel（'幅度'）;

axis（[0,2,0,1.2*max（abs（X1k8））]）

subplot（2,2,3）;mstem（X1k16）;%绘制16点DFT的幅频特性图

title（'（1b）16点DFT[x_1（n）]'）;xlabel（'ω/π'）;ylabel（'幅度'）;

axis（[0,2,0,1.2*max（abs（X1k16））]）

figure

（2）

subplot（2,2,1）;mstem（X2k8）;%绘制8点DFT的幅频特性图

title（'（2a）8点DFT[x_2（n）]'）;xlabel（'ω/π'）;ylabel（'幅度'）;

axis（[0,2,0,1.2*max（abs（X2k8））]）

subplot（2,2,2）;mstem（X2k16）;%绘制16点DFT的幅频特性图

title（'（2b）16点DFT[x_2（n）]'）;xlabel（'ω/π'）;ylabel（'幅度'）;

axis（[0,2,0,1.2*max（abs（X2k16））]）

subplot（2,2,3）;mstem（X3k8）;%绘制8点DFT的幅频特性图

title（'（3a）8点DFT[x_3（n）]'）;xlabel（'ω/π'）;ylabel（'幅度'）;

axis（[0,2,0,1.2*max（abs（X3k8））]）

subplot（2,2,4）;mstem（X3k16）;%绘制16点DFT的幅频特性图

title（'（3b）16点DFT[x_3（n）]'）;xlabel（'ω/π'）;ylabel（'幅度'）;

axis（[0,2,0,1.2*max（abs（X3k16））]）

%实验内容

（2）周期序列谱分析==================================

N=8;n=0:

N-1;%FFT的变换区间N=8

x4n=cos（pi*n/4）;

x5n=cos（pi*n/4）+cos（pi*n/8）;

X4k8=fft（x4n）;%计算x4n的8点DFT

X5k8=fft（x5n）;%计算x5n的8点DFT

N=16;n=0:

N-1;%FFT的变换区间N=16

x4n=cos（pi*n/4）;

x5n=cos（pi*n/4）+cos（pi*n/8）;

X4k16=fft（x4n）;%计算x4n的16点DFT

X5k16=fft（x5n）;%计算x5n的16点DFT

figure（3）

subplot（2,2,1）;mstem（X4k8）;%绘制8点DFT的幅频特性图

title（'（4a）8点DFT[x_4（n）]'）;xlabel（'ω/π'）;ylabel（'幅度'）;

axis（[0,2,0,1.2*max（abs（X4k8））]）

subplot（2,2