离散时间信号与系统的频域分析.docxVIP

信号与系统实验 实验报告 PAGE 1 成绩： 序号： 22 成绩： 信号与系统实验 实 验 报 告 实验五 离散时间信号与系统的频域分析 姓 姓 名曦杰 班 级：通 日 期：2021年11月16日 指导教师：许方敏  一、实验目的 二、实验原理 三、实验内容和结果 程序1.1 clear all; w=-4*pi:8*pi/511:4*pi; num=[2 1];den=[1 -0.6]; h1=freqz(num,den,w); subplot(2,1,1);plot(w/pi,real(h1)); grid;title(实部) xlabel(omega/\pi);ylabel(振幅); subplot(2,1,2);plot(w/pi, imag(h1)); grid;title(虚部) xlabel(omega/pi);ylabel(振幅); figure; subplot(2,1,1);plot(w/pi,abs(h1)); grid;title(h1幅度谱) xlabel(omega/pi);ylabel(振幅); subplot(2,1,2);plot(w/pi,angle(h1)); grid;title(h1相位谱) xlabel(omega/pi);ylabel(相位); 结果1.1 程序1.2 clear all; w=-4*pi:8*pi/511:4*pi; D=4; num=[2 1];den=[1 -0.6]; num2=[zeros(1,D),num]; h2=freqz(num2,den,w); subplot(2,1,1);plot(w/pi,abs(h2)); grid; title(h2幅度谱); xlabel(omega/pi); ylabel(振幅); subplot(2,1,2) plot(w/pi,angle(h2)); grid;title(h2相位谱) xlabel(omega/pi);ylabel(相位); 结果1.2 与结果1.1进行对比，可以看出离散时间信号进行傅里叶变换之后，幅度是不变的，而相位会发生改变 示例3 n = 0:100; s1 = cos(2*pi*0.05*n); s2 = cos(2*pi*0.47*n); x = s1+s2; y = filter([1:3],1,x); clf; subplot(2,2,1); plot(n,s1); axis([0,100,-2,2]); xlabel(Timeindexn);ylabel(Amplitude); title(Signal#1); subplot(2,2,2); plot(n,s2); axis([0,100,-2,2]); xlabel(Timeindexn);ylabel(Amplitude); title(Signal#2); subplot(2,2,3); plot(n,x); axis([0,100,-2,2]); xlabel(Timeindexn);ylabel(Amplitude); title(InputSignal); subplot(2,2,4); plot(n,y); axis([0,100,-2,2]); xlabel(Timeindexn);ylabel(Amplitude); title(OutputSignal); axis; 结果 程序2 n = 0:100; a=1; b=[0.5 -0.5] y = filter(b,a,x); clf; subplot(2,2,1); plot(n,s1); axis([0,100,-2,2]); xlabel(Timeindexn);ylabel(Amplitude); title(Signal#1); subplot(2,2,2); plot(n,s2); axis([0,100,-2,2]); xlabel(Timeindexn);ylabel(Amplitude); title(Signal#2); subplot(2,2,3); plot(n,x); axis([0,100,-2,2]); xlabel(Timeindexn);ylabel(Amplitude); title(InputSignal); subplot(2,2,4); plot(n,y); axis([0,100,-2,2]); xlabel(Timeindexn);ylabel(Amplitude); title(OutputSignal); axis; 结果2 系统的输出波形