资源描述：
《实验二Z变换和傅立叶变换2011.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、数字信号处理学院理学院年级专业2011级电子信息工程01班指导老师刘x学号5姓名郑xx实验二Z变换和傅立叶变换实验目的：掌握Z变换和反变换和傅立叶变换在MATLAB信号处理工具箱提供的求连续和离散系统频率响应的函数。实验内容：(1)有限长序列的Z变换和逆Z变换x1=[1,2,3];ns1=-1;nf1=ns1+length(x1)-1;x2=[2,4,3,5];ns2=-2;nf2=ns1+length(x1)-1;x=conv(x1,x2)n=(ns1+ns2):(nf1+nf2)结果：x=281723
2、1915n=-3-2-1012(2)clear,closeallx=[2,4,3,5];nsx=-2;nfx=nsx+length(x)-1;Bw=-3;nsbw=-1;Aw=[2,-2.2,.5];nsaw=0;B=conv(-3,x);A=Aw;nsy=nsaw-(nsbw-nsx)[r,p,k]=residuez(B,A)nf=input('终点时间nf=');n=min(nsx,nsy):nf;yi=(r(1)*p(1).^(n-nsy)+r(2)*p(2).^(n-nsy)).*stepseq(
3、nsy,n(1),nf);yd=k(1)*impseq(nsy,n(1),nf)+k(2)*impseq(-1-nsy,n(1),nf);y=yi+yd;xe=zeros(1,length(n));xe(find((n>=nsx)&(n<=nfx)==1))=x;subplot(2,1,1),stem(n,xe,'.'),line([min(n(1),0),nf],[0,0])subplot(2,1,2),stem(n,y,'.'),line([min(n(1),0),nf],[0,0])m文件impse
4、q.mfunction[x,n]=impseq(n0,ns,nf)n=[ns:nf];x=[(n-n0)==0];stepseq.mfunction[x,n]=stepseq(n0,ns,nf)n=[ns:nf];x=[(n-n0)>=0];（1）离散时间傅立叶变换disp('八点时间信号的离散时间傅立叶变换')x0=sin(2*pi*[1:8]/8)*5;dt=2*pi/8;w=linspace(-8,2*pi,1000)/dt;X0=dtft(x0,w)*dt;subplot(3,1,1),plot(
5、w,abs(X0)),axis([-8,max(w),0,max(abs(X0))]),grid,shgdisp('重复N次的八点时间信号的离散离散时间傅立叶变换')N=input('N=');x1=reshape(x0'*ones(1,N),1,N*length(x0));X1=dtft(x1,w)*dt;subplot(3,1,2),plot(w,abs(X1)),axis([-8,max(w),0,max(abs(X1))]),grid,shgdisp('重复无穷次的八点信号的离散傅立叶变换')pa
6、use,X2=fft(x0*dt);w1=2*pi*[0:length(x0)-1]/length(x0);subplot(3,1,3),stem([-w1,w1],[abs(X2),abs(X2)]),grid,shgaxis([-8,max(w),0,max(abs(X2))]),grid,shgm文件dtft.mfunctionX=dtft(x,w)X=x*exp(-j*[1:length(x)]'*w);（1）由离散序列恢复模拟信号clear,closeall;A=444.128;a=50*sqr
7、t(2)*pi;b=a;fork=1:2ifk==1Fs=400;elseifk==2Fs=1000;endT=1/Fs;dt=T/3;Tp=0.03;t=0:dt:Tp;n=0:Tp/T;TMN=ones(length(n),1)*t-n'*T*ones(1,length(t));x=A*exp(-a*n*T).*sin(b*n*T);xa=x*sinc(Fs*TMN);subplot(2,1,k);plot(t,xa);holdonaxis([0,max(t),min(xa)-10,max(xa)+1
8、0]);st1=sprintf('有Fs=%d',Fs);st2='Hz的采样序列z(n)重构的信号';st=[st1,st2];title(st)ylabel('xa(t)');xo=A*exp(-a*t).*sin(b*t);stem(t,xo,'.');line([0,max(t)],[0,0])emax2=max(abs(xa-xo))end（1）梳状滤波器零极点和幅频特性clear;closeallb=[1,0,0,
