试验常见离散信号的MATLAB产生和图形显示-Read.DOC

实验1 常见离散信号的MATLAB产生和图形显示 实验目的：加深对常用离散信号的理解； 实验原理： 单位抽样序列 在MATLAB中可以利用zeros()函数实现。 如果在时间轴上延迟了k个单位，得到即： 2．单位阶跃序列 在MATLAB中可以利用ones()函数实现。 3．正弦序列 在MATLAB中 4．复正弦序列 在MATLAB中 5．指数序列 在MATLAB中 实验内容：编制程序产生上述5种信号（长度可输入确定），并绘出其图形。 实验要求：讨论复指数序列的性质。 实验2 离散系统的差分方程、冲激响应和卷积分析 实验目的：加深对离散系统的差分方程、冲激响应和卷积分析方法的理解。 实验原理：离散系统 其输入、输出关系可用以下差分方程描述： 输入信号分解为冲激信号，。记系统单位冲激响应 ，则系统响应为如下的卷积计算式： 当时，h[n]是有限长度的（n：[0，M]），称系统为FIR系统；反之，称系统为IIR系统。 在MATLAB中，可以用函数y=Filter(p,d,x) 求解差分方程，也可以用函数 y=Conv(x,h)计算卷积。 实验内容：编制程序求解下列两个系统的单位冲激响应和阶跃响应，并绘出其图形。 实验要求：给出理论计算结果和程序计算结果并讨论。 实验3 离散系统的频率响应分析和零、极点分布 实验目的：加深对离散系统的频率响应分析和零、极点分布的概念理解。 实验原理：离散系统的时域方程为 其变换域分析方法如下： 频域 系统的频率响应为 Z域 系统的转移函数为 分解因式 ，其中和称为零、极点。 在MATLAB中，可以用函数[z，p，K]=tf2zp（num，den）求得有理分式形式的系统转移函数的零、极点，用函数zplane（z，p）绘出零、极点分布图；也可以用函数zplane（num，den）直接绘出有理分式形式的系统转移函数的零、极点分布图。 另外，在MATLAB中，可以用函数 [r，p，k]=residuez（num，den）完成部分分式展开计算；可以用函数sos=zp2sos（z，p，K）完成将高阶系统分解为2阶系统的串联。 实验内容：求系统 的零、极点和幅度频率响应。 实验要求：编程实现系统参数输入，绘出幅度频率响应曲线和零、极点分布图。 实验4 离散信号的DTFT和DFT 实验目的：加深对离散信号的DTFT和DFT的及其相互关系的理解。 实验原理：序列x[n] 的DTFT定义： N点序列x[n] 的DFT定义： 在MATLAB中，对形式为的DTDFT可以用函数H=Freqz（num，den，w）计算；可以用函数U=fft（u，N）和u=ifft（U，N）计算N点序列的DFT正、反变换。 实验内容： 分别计算16点序列 的16点和32点DFT，绘出幅度谱图形，并绘出该序列的DTFT图形。 实验要求：讨论DTFT和DFT之间的相互关系。说明实验产生的现象的原因。 实验5 FFT算法的应用 实验目的：加深对离散信号的DFT的理解及其FFT算法的运用。 实验原理：N点序列的DFT和IDFT变换定义式如下： ， 利用旋转因子具有周期性，可以得到快速算法（FFT）。 在MATLAB中，可以用函数X=fft（x，N）和x=ifft（X，N）计算N点序列的DFT正、反变换。 实验内容： 2N点实数序列 N=64。用一个64点的复数FFT程序，一次算出，并绘出。 （2）已知某序列在单位圆上的N=64等分样点的Z变换为。用N点IFFT程序计算，绘出和。 实验要求：利用MATLAB编程完成计算，绘出相应图形。并与理论计算相比较，说明实验结果的原因。