2010级无线移动1班.docVIP

实验目的：1. 学会根据序列编写程序让其转换成Z变换和傅里叶变换。 2. Z变换是时域离散信号和系统分析及设计的重要数学工具。 MATLAB作数值分析时，是无法求无限长度序列的Z变换，这个问题要靠MATLAB中的符号运算（symbolic）工具箱才能解决。因此要学会用程序来解决Z变换和傅里叶变换。 实验原理：MATLAB信号处理工具箱提供了求连续和离散系统频率响应的两个函数。 Freqs和freqz两个函数。 如果Z变换是Z的有理分式，虽然其逆Z变换是无限序列，但求它的系数和指数都是数值计算的范畴，可以用MATLAB解决。 实验步骤：例7.7有限长度序列的Z变换和逆Z变换 已知两序列x1=[1.2.3],n1=[-1:1]及x2=【2，4，3，5】，n2=[-2:1],求出X1与x2及其卷积x的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) x = 2 8 17 23 19 15 n=(ns1+ns2):(nf1+nf2) n = -3 -2 -1 0 1 2 由此可得X(Z)=2Z^3+8Z^2+17Z+23+19Z^-1+15Z^-2 X(Z)是Z的多项式形式，她的你变换就是其系数组成的向量x和z的幂次组成的时间向量n,所以是有限长度序列。 disp(八点时间信号的离散时间傅立叶变换) x0=sin(2*pi*[1:8]/8)*5; dt=2*pi/8; w=linspace(0,2*pi,1000)/dt; X0=dtft(x0,w)*dt; subplot(3,1,1),plot(w,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)),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 八点时间信号的离散时间傅立叶变换 ??? Undefined function or variable dtft. disp(八点时间信号的离散时间傅立叶变换) x0=sin(2*pi*[1:8]/8)*5; dt=2*pi/8; w=linspace(0,2*pi,1000)/dt; X0=dtft(x0,w)*dt; subplot(3,1,1),plot(w,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)),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 八点时间信号的离散时间傅立叶变换 7.10时域采样频率与频谱混叠 xa(t)=Ae^-atsinWtu(t),A=444.128, 计算并图示xa(t)及其幅频特性函数|Xa(jW)|. 解：其程序为：clear;close all; fs=10000;fs1=1000;fs2=400;fs3=200; t=0:1/fs:0.1; A=444.128;a=50*sqrt(2)*pi;b=a; xa=exp(-a*t).*sin(b*t); k=0:511;f=fs*k/512; Xa=dtft(xa,2*pi*k/512); T1=1/fs1;t1=0:T1:0.1; x1=A*exp(-a.*t1).*sin(b*t1); X1=dtft(x1,2*pi*k/512); T2=1/fs2;t2=0:T2:0.1; x2=A*exp(-a.*t2).*sin(b.*t2); X2=