离散时间信号的产生及信号的卷积和运算实验报告2

离散时间信号的产生及信号的卷积和运算实验报告班级：___________姓名：__________学号:____________一、实验目的和原理实验原理：（一）DTFT和DFT的定义及其相互关系：序列x[n]的DTFT定义：njnωjωx[n]e)X(e它是关于自变量的复函数，且是以2为周期的连续函数。)X(ejω可以表示为：)(ejX)(eX)X(ejωimjωrejω其中，)(eXjωre和)(eXjωim分别是)X(ejω的实部和虚部；还可以表示为：)(ωjjωjωe)X(e)X(e其中，)X(ejω和}arg{)()X(ejω分别是)X(ejω的幅度函数和相位函数；它们都是的实函数，也是以2为周期的周期函数。序列x[n]的N点DFT定义：101022][][)(][NnknNNnknNjkNjWnxenxeXkX][kX是周期为N的序列。)X(ejω与][kX的关系：][kX是对)X(ejω在一个周期中的谱的等间隔N点采样，即：kNjω)X(ekX2|][，而)X(ejω可以通过对][kX内插获得，即：]2/)1)][(/2([10)22sin()22sin(][1NNkjNkjωeNkNkNkXN)X(e(二)线性时不变离散时间系统的变换域表示：LTI离散时间系统的时域差分方程为：MkkNkkknxpknyd00)()((1)传递函数：对上面的差分方程两边求z变换，得：NkkkMkkkMkkkNkkkzdzpzXzYzpzXzdzY0000)()()()(我们定义LTI离散时间系统的输出的Z变换Y(z)与输入的Z变换X(z)的比值为该系统的传递函数，即)()()(zXzYzH为系统的传递函数。NNMMzdzddzpzppzDzpzH......)()()(110110分解因式NiiMiiNiikMiikzzKzdzpzH111100)1()1()(，其中i和i称为零、极点。利用系统的传递函数)(zH，我们可以分析系统的零极点，稳定性及实现结构等特点。(2)频率响应：因为大多数离散时间信号都可以分解为nje的线性组合，所以研究输入nje的响应具有极大的意义，即当输入为njenx][时，输出为：)()()(][)(jnjmmjnjmnjmeHeemheemhny这里，nnjjenheH)()(是h(n)的DTFT，称为LTI离散时间系统的频率响应。利用系统的频率响应)(jeH，我们可以分析系统对各种频率成分的响应特性，并推出系统的特性（高通，低通，带通，带阻，线性相位等）。(3)系统传递函数与频率响应之间的关系：从前面的推导可以看出，系统的传递函数是系统冲击响应序列的Z变换，而系统的频率响应是冲击响应的DTFT，因此传递函数)(zH与频率响应)(jeH的关系为：jezjzHeH|)()(jNNjjMMjjjjededdepeppeDepeH......)()()(1010实验目的：信号的变换域分析是信号处理中一种有效的工具。在离散信号的时域分析中，我们通常将信号表示成单位采样序列][n的线性组合，而在频域中，我们将信号表示成复变量nje或nNje2的线性组合。通过这样的表示，可以将时域的离散序列映射到频域以便于进一步的处理。线性时不变（LTI）离散时间系统的特性完全可以用其冲击响应序列h[n]来表示，则前面给出的离散时间信号的变换分析手段也可以用于离散时间系统的分析中。在LTI的离散时间系统变换域分析中，我们常用传递函数H(z)和频率响应H(ejw)来表示系统。在本实验中，将学习利用MATLAB计算离散时间信号的DTFT和DFT，并加深对其相互关系的理解。通过使用MATLAB函数对离散时间系统的一些特性进行仿真分析，以加深对离散时间系统的零极点、稳定性，频率响应等概念的理解。二、实验过程（1）num=[0.0518-0.15530.15530.0518];den=[11.28281.03880.3418];w=0:0.001:2*pi;H=freqz(num,den,w);figuresubplot(221)plot(w,real(H))title('ʵ²¿')gridsubplot(222)plot(w,imag(H))title('Ð鲿')gridsubplot(223)plot(w,abs(H))title('·ù¶È')gridsubplot(224)plot(w,unwrap(angle(H)))title('Ïàλ')grid（2）n1=0:15;n2=0:31N=32;x1=cos((5*pi*n1)/16)x2=cos((5*pi*n1)/16)X1=fft(x1);X2=fft(x2,N);subplot(211)stem(n1,abs(X1))title('16µãµÄDFT')subplot(212)stem(n2,abs(X2))title('32µãµÄDFT')（3）num=[0.05280.7970.12950.12950.7970.0528]den=[1-1.81072.4947-1.88010.9537-0.2336]w=0:0.004:2*pi;h=freqz(num,den,w)subplot(2,2,1)plot(w/pi,abs(h));title('·ù¶ÈÆ×')subplot(2,2,2)plot(w/pi,angle(h));title('ÏàλÆ×')[z,p,k]=tf2zp(num,den)disp('Áãµã');disp(z);disp('¼«µã');disp(p);subplot(2,2,3);zplane(num,den);title('Á㼫µãͼ')