实验六离散时间傅立叶变换.pptVIP

一、实验目的 运用MATLAB计算序列的DTFT 运用MATLAB验证离散时间傅立叶变换的性质 运用MATLAB计算离散时间系统的频率响应 二、实验设备 二、实验原理 二、实验原理 n=-10:10; x=-0.9.^n; w=-2*pi:0.01:2*pi; Xw=x*exp(-j*n*w); magXw=abs(Xw); angXw=angle(Xw); subplot(311); stem(n,x); grid xlabel(n); ylabel(x(n)); subplot(312); plot(w,magXw); grid xlabel(frequency in units of pi);ylabel(|X(w)|); subplot(313) plot(w,angXw);grid xlabel(frequency in units of pi);ylabel(angle/pi) 二、实验原理 函数freqz可以计算一个以 的有理分式形式给出的序列的DTFT值。 形式 H=freqz(num,den,w), 参数 num表示序列有理分式DTFT的分子多项式系数 den表示分母多项式系数（均按z的降幂排列） 矢量w表示在0～2?中给定的一系列频率点集合 求出DTFT值后，可以使用函数real, imag, abs和angle分别求出并绘出其实部，虚部，幅度和相位谱。如果需要，还可以用unwrap函数消除相位中的跳变。 例2 求信号 的DTFT，并画出其幅频特性和相频特性 二、实验原理 三、实验原理 三、实验原理 clf; %Clear current figure w=-pi:2*pi/255: pi; num1=[1 3 5 7 9 11 13 15 17]; X1=freqz(num1,1,w); L=length(num1); n=0:L-1; wo=0.4*pi; num2=exp(wo*i*n).*num1; %时域相移 X2=freqz(num2,1,w); clf; w=-pi:2*pi/255: pi; x1=[1 3 5 7 9 11 13 15 17]; x2=[1 -2 3 -2 1]; y=conv(x1,x2); H1=freqz(x1,1,w); H2=freqz(x2,1,w); Hp=H1.*H2; H3=freqz(y,1,w); subplot(2,2,1); plot(w/pi,abs (Hp));grid; %Hp=H1.*H2 title(‘幅度谱的乘积’); subplot(2,2,2); plot(w/pi,abs (H3));grid; title(卷积后序列的幅度谱) subplot(2,2,3); plot(w/pi,angle (Hp));grid ; %Hp的相位谱 title(相位谱的和) ; %Hp的相位等于H1和H2的对应频率上的相位和 subplot(2,2,4); plot(w/pi,angle (H3));grid; title(卷积后序列的相位谱); 3．计算离散系统的频率响应 四、实验内容 1. 对信号 采样，画出在[-1,1]之间的序列，求其DTFT，并画出其幅频特性和相频特性。 增加取样区间范围，再求新序列的DTFT,画出幅频特性，分析有什么不同？ 2. 求信号的离散时间傅立叶变换并分析其周期性和对称性； 给定正弦信号x(t)=2*cos(2*pi*10*t),fs=100HZ,求其DTFT。 3.已知两序列x1(n)=[1 2 3 4 5 6 7 8 9 10 11 12]， x2(n)=R10(n) （1）求其DTFT，并给出各自幅频响应 （2）用这两个序列来验证DTFT的线性、乘以指数序列和共轭特性； 4. 求LTI系统的频率响应 给定系统H(Z)=B(Z)/A(Z)，a=[0.98777 -0.31183 0.0256] b=[0.98997 0.989 0.98997]，求系统的幅频响应和相频响应。 clf;