实验三离散系统Z域解析总结计划

实验三、失散系统的Z域解析（一）实验要求1）学习和掌握失散系统的频率特点及其幅度特点、相位特点的物理意义；2）深入理解失散系统频率特点的对称性和周期性；3）认识失散系统频率特点与系统参数之间的关系；4）经过阅读、更正并调试本实验系统所给源程序，加强计算机编程能力；（二）实验内容1、计算差分方程（1）用MATLAB计算差分方程当输入序列为时的输出结果。MATLAB程序以下：N=41;a=[0.8-0.440.360.22];b=[10.7-0.45-0.6];x=[1zeros(1,N-1)];k=0:1:N-1;h=filter(a,b,x);stem(k,h)xlabel('n');ylabel('h(n)')请给出了该差分方程的前41个样点的输出，即该系统的单位脉冲响应。(说明：y=filter(a,b,x)，计算系统对输入信号向量x的零状态响应输出信NM号向量y,x与y长度相等，其中a和b是aiy(ni)bix(ni)所给差分方ii程的相量。详见教材P25-27)2、用MATLAB计算差分方程所对应的系统函数的FT。差分方程所对应的系统函数为：10.80.44z10.36z20.02z3H(z)0.7z10.45z20.6z31FT为H(ej)0.80.44e10.7e

j0.36ej20.02ej3j0.45ej20.6ej3用MATLAB计算的程序以下：k=256;num=[0.8-0.440.360.02];den=[10.7-0.45-0.6];w=0:pi/k:pi;h=freqz(num,den,w);subplot(2,2,1);plot(w/pi,real(h));gridtitle('实部')xlabel('\omega/\pi');ylabel('幅度')subplot(2,2,2);plot(w/pi,imag(h));gridtitle('虚部')xlabel('\omega/\pi');ylabel('Amplitude')subplot(2,2,3);plot(w/pi,abs(h));gridtitle('幅度谱')xlabel('\omega/\pi');ylabel('幅值')subplot(2,2,4);plot(w/pi,angle(h));gridtitle('相位谱')xlabel('\omega/\pi');ylabel('弧度')(说明：freqz为计算数字滤波器H(z)的频率响应函数。h=freqz(num,den,w)为计算由向量w指定的数字频率点上数字滤波器H(z)的频率响应H(ej)，结果存于h向量中。Num和den为H(z)分子和分母多项式向量。详见教材P65)练习:①、P43中的例3、例5。2②、用MATLAB编程，画出1z1H(z),的频率特点图。10.2z13、求解X(z)1z0.61)(10.5z1)2(10.6z1)2(10.2z的Z反变换。参照程序：b=1;a=poly([-0.20.50.5-0.6-0.6]);[r,p,k]=residuez(b,a)(说明：例程序：b=[-48];a=[168];[r,p,k]=residuez(b,a)运行结果：r=-128=-4-2=[]则表示：那么：)（三）实验报告要求31、简述实验目的和实验原理，用几何确定法解析实验中选定的系统的频率特点，并与计算机计算结果对比较，依照实验结果，对系统频率特点进行谈论和总结。2、用MATLAB编程，画出P66中的例的频率特点图。3、依照MATLAB求解X(z)1z0.6的结(10.2z1)(10.5z1)2(10.6z1)2果，写出序列x(n)的表达式。实验四、失散傅里叶变换及其快速算法（一）实验要求1）经过错散傅立叶变换（即DFT）的报表表示进一步认识其计算方法及意义；2）掌握实数序列的DFT系数的对称特点；3）学习利用DFT计算程序计算IDFT的方法；学习时间抽选奇偶分解FFT算法；深入掌握时间抽选奇偶分解FFT程序的编制方法；(二)实验内容1、对连续的单一频率周期信号按采样频率采样，截取长度N分别选N=20和N=16，观察其DFT结果的幅度谱。此时失散序列，即k=8。用MATLAB计算并作图，函数fft用于计算失散傅里叶变换DFT，程序以下：k=8;n1=[0:1:19];xa1=sin(2*pi*n1/k);subplot(2,2,1)plot(n1,xa1)xlabel('t/T');ylabel('x(n)');4xk1=fft(xa1);xk1=abs(xk1);subplot(2,2,2)stem(n1,xk1)xlabel('k');ylabel('X(k)');n2=[0:1:15];xa2=sin(2*pi*n2/k);subplot(2,2,3)plot(n2,xa2)xlabel('t/T');ylabel('x(n)');xk2=fft(xa2);xk2=abs(xk2);subplot(2,2,4)stem(n2,xk2)xlabel('k');ylabel('X(k)');5计算结果示于图2.1，(a)和(b)分别是N=20时的截守信号和DFT结果，由于截取了两个半周期，频谱出现泄漏；(c)和(d)分别是N=16时的截守信号和DFT结果，由于截取了两个整周期，获取单一谱线的频谱。上述频谱的误差主若是由于时域中对信号的非整周期截断产生的频谱泄漏。2、对以下各序列进行谱解析，绘制出其幅频特点曲线。1）x1(n)R4(n)（2）x2(n)con(n)con(n)48对x1n)R4nx2(n)con(n)con(n)为周期序列，周期为(()取64点FFT。由于4816，所以取周期为16。参照程序以下：%用FFT对序列进行谱解析%x1(n)R4(n)；x2(n)con(n)con(n)48%X1:存放x1(n)的向量，Y1:存放X1(k)的向量%X2:存放x2(n)的向量，Y2:存放X2(k)的向量x1=[1,1,1,1,0,0,0,0];n=0:15;x2=cos(pi*n/4)+cos(pi*n/8);i=0:7;subplot(3,2,1);stem(i,x1,'.');axis([0701]);%规定x轴和y轴的标值范围xlabel('n');ylabel('x1(n)');y1=fft(x1,8);%x1(n)的8点FFTsubplot(3,2,3);stem(i,abs(y1),'.');xlabel('(N=8wk=2pik/N)k');ylabel('[X1(k)]');y1=fft(x1,64);%x1(n)的64点FFTi=0:63;subplot(3,2,5);stem(i,abs(y1),'.');axis([06304]);%规定x轴和y轴的标值范围xlabel('(N=8wk=2pik/N)k');ylabel('[X1(k)]');y2=fft(x2);%x2(n)的16点FFTfigure;%另一幅图subplot(2,2,1);stem(n,x2);6title('x2(n)的时域序列');xlabel('n');ylabel('x2(n)');subplot(2,2,3);stem(n,abs(y2));title('x2(n)的幅频特点');xlabel('(N=16wk=2pik/N)k');ylabel('[X2(k)]');（三）实验报告要求1、设()R4(),(j)[()]。分别计算(j)在频率区间[0,2]上xnnXeFTxnXe的16点和32点等间隔采样，并绘制X(ej)采样的幅频特点图和相频特点图。如：%DFT的MATLAB计算xn=[1111];%输入时域序列向量Xk16=fft(xn,16);%计算xn16点DFTXk32=fft(xn,32);%计算xn32点DFT绘图程序略。2、编写序列x(n)con(n)con(n)的DFT运算程序。48解：参照程序以下：N=16;n=0:1:N-1;%时域采样xn=cos(n*pi/4)+cos(n*pi/8);k=0:1:N-1;%频域采样WN=exp(-j*2*pi/N);nk=n'*k;WNnk=WN.^nk;%点乘方Xk=xn*WNnk;Subplot(2,1,1)Stem(n,xn)title('xn的时域序列');xlabel('n');ylabel('x(n)');Subplot(2,1,2)Stem(k,abs(Xk));title('xn的幅频特点');xlabel('(N=16wk=2pik/N)k');ylabel('[X(k)]');（Fft与Ifft的说明：y=fft(x)是利用fft函数求解x的失散傅里叶变换；y=fft(x，N)，N表示失散傅里叶变换x的数据长度；函数Ifft的参数与函数7Fft完好相同。例：fft在信号解析中的应用。使用频谱解析方法从受噪声污染的信号x(t)中鉴别出适用信号。如程序：t=0:0.001:1;%采样周期为0.001s，即采样频率为1000Hz;%产生受噪声污染的正弦波信号；x=sin(2*pi*100*t)+sin(2*pi*200*t)+rand(size(t));subplot(2,1,1)plot(x(1:50));%画出时域内的信号；y=ff