实验二Z变换和傅立叶变换201101830105详解.doc

数字信号处理 学 院 理学院 年级专业 2011级电子信息工程01班 指导老师 刘x 学 号 201101830105 姓 名 郑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 = 2 8 17 23 19 15 n = -3 -2 -1 0 1 2 (2) clear, close all x=[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 (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文件 impseq.m function [x,n]=impseq(n0,ns,nf) n=[ns:nf];x=[(n-n0)==0]; stepseq.m function [x,n]=stepseq(n0,ns,nf) n=[ns:nf];x=[(n-n0)=0]; 离散时间傅立叶变换 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(w,abs(X0)),axis([-8,max(w),0,max(abs(X0))]),grid,shg disp(重复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,shg disp(重复无穷次的八点信号的离散傅立叶变换) pause,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,shg axis([-8,max(w),0,max(abs(X2))]),grid,shg m文件 dtft.m function X=dtft(x,w) X=x*exp(-j*[1:length(x)]*w); 由离散序列恢复模拟信号 clear,close all; A=444.128;a=50*sqrt(2)*pi;b=a; for k=1:2 if k==1 Fs=400; elseif k==2 Fs=1000;end T=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).*si