信号和系统实验6.docVIP

信号与系统实验（六） 班 章仕波 刘贺洋 实验内容 1离散时间傅里叶变换 (1)下面参考程序是如下序列在范围的离散时间傅里叶变换 %计算离散时间傅里叶变换的频率样本 clear all; w=-4*pi;8*pi/511;4*pi; num=[2 1]; den=[1 -0.6]; h=freqz(num,den,w); subplot(2,1,1) plot(w/pi,real(h)); grid; title(‘实部’) xlabel(‘omega/\pi’); ylabel(‘振幅’); subplot(2,1,2) plot(w/pi, imag(h)); grid; title(‘虚部’) xlabel(‘omega/\pi’); ylabel(‘振幅’); figure; subplot(2,1,1) plot(w/pi, abs(h)); grid; title(‘幅度谱’) xlabel(‘omega/\pi’); ylabel(‘振幅’); subplot(2,1,2) plot(w/pi, angle (h)); grid; title(‘相位谱’) xlabel(‘omega/\pi’); ylabel(‘以弧度为单位的相位’); 修改程序，在范围内计算如下有限长序列的离散时间傅里叶变换 h[n]=[1 2 3 4 5 6 7 8 9] 利用(1)的程序，通过比较结果的幅度谱和相位谱，验证离散时间傅里叶变换的时移特性。（提示：可设num2=[zeros(1,D),num]） (1) clear all; w=0:pi/511:pi; h=[1 2 3 4 5 6 7 8 9] H=fft(h,512) subplot(2,1,1) plot(w/pi,real(H)); grid; title(实部) xlabel(omega/\pi); ylabel(振幅); subplot(2,1,2) plot(w/pi, imag(H)); grid; title(虚部) xlabel(omega/\pi); ylabel(振幅); figure; subplot(2,1,1) plot(w/pi, abs(H)); grid; title(幅度谱) xlabel(omega/\pi); ylabel(振幅); subplot(2,1,2) plot(w/pi, angle (H)); grid; title(相位谱) xlabel(omega/\pi); ylabel(以弧度为单位的相位); (2) clear all; w=-4*pi:64*pi/511:4*pi; num=[2 1]; den=[1 -0.6]; num2=[zeros(1,2),num]; h=freqz(num,den,w); h2=freqz(num2,den,w); subplot(3,1,1) plot(w/pi, abs(h)); grid; title(原来幅度谱) xlabel(omega/\pi); ylabel(振幅); subplot(3,1,2) plot(w/pi, abs(h2)); grid; title(延时后幅度谱) xlabel(omega/\pi); ylabel(振幅); subplot(3,1,3) plot(w/pi, (abs(h)-abs(h2))10^-14); grid; title(两者的差值) xlabel(omega/\pi); ylabel(振幅); figure; subplot(3,1,1) plot(w/pi, angle (h)); grid; title(原来相位谱) xlabel(omega/\pi); ylabel(以弧度为单位的相位); subplot(3,1,2) plot(w/pi, angle (h2)); grid; title(延时后相位谱) xlabel(omega/\pi); ylabel(以弧度为单位的相位); subplot(3,1,3) plot(w/pi,angle (h)-angle (h2)); grid; title(两者相位差) xlabel(omega/\pi); ylabel(以弧度为单位的相位); 2离散傅里叶变换值的求解 对有限长序列x=[0 2 4 6 8 10 12 14 16]，计算并画出L点离散傅里叶变换X[k]的值，改变L，重做几次实验，讨论结果。 clear all; L=15; w=-4*pi:8*pi/(L-1):4*pi x=[0 2 4 6 8 10 12 14 16]; x=[x,zeros(1,L-9)] H=fft(