武汉工程大学matlab实验四lst离散系统频域分析.docxVIP

武汉工程大学 数字信号处理实验报告四 专业班级： 14级通信03班 学生姓名： 秦重双 学 号： 1404201114 实验时间： 2017年5月12日 实验地点： 4B315 指导老师： 杨述斌 实验三: LSI离散系统的频域分析 实验目的 通过在频域中仿真LSI离散时间系统，理解离散时间系统对输入信号或延迟信号进行频域处理的特性。 理解LSI离散时间系统的传输函数和频率响应的概念。 理解LSI离散时间系统的滤波特性及滤波器的相关特性。 4.理解并仿真LSI离散时间系统的零、极点分布表征及特性关系。 二、实验设备 计算机，MATLAB语言环境 三、实验基础理论 LSI离散时间系统可用差分方程描述如下： 对应的传输函数和频率响应分别为： 分别有零点和极点。 实验内容与步骤 1.传输函数和频率响应分析 按以下的传输函数分别编程计算 和 计算当 时因果LSI离散时间系统的频率响应，并求出它们的群时延及冲激响应的开始部分（前100个值）。 b=[0.15,0,-0.15];a=[1,-0.5,0.7]; b1=[0.15,0,-0.15];a1=[0.7,-0.5,1]; figure(1)； subplot(1,2,1)； plot(grpdelay(b,a))； grid；title(系统1的群延时) subplot(1,2,2)； plot(grpdelay(b1,a1))； grid,title(系统2的群延时)； w=0:0.1:pi; z=exp(j*w); [H1,w]=freqz(b,a); [H2,w]=freqz(b1,a1); figure(2); subplot(2,2,1)； plot(w,abs(H1)),title(系统1的频率幅度响应) subplot(2,2,2)； plot(w,abs(H2))；title(系统2的频率幅度响应) subplot(2,2,3)； plot(w,angle(H1))；title(系统1的频率相位响应) subplot(2,2,4)； plot(w,angle(H2))；title(系统2的频率相位响应) xn=[1,zeros(1,99)]; hn1=filter(b,a,xn); hn2=filter(b1,a1,xn); n1=0:length(hn1)-1; n2=0:length(hn2)-1; figure(3); stem(n1,hn1),xlabel(n),ylabel(h1(n)),title(冲激响应hn1); figure(4); stem(n2,hn2),xlabel(n),ylabel(h2(n)),title(冲激响应hn2); 画出上面两个LSI离散时间系统对应的零、极点图。 b=[0.15,0,-0.15]; a=[1,-0.5,0.7]; subplot(2,1,1)； zplane(b,a)； title(h1的零极点分布)； b1=[0.15,0,-0.15]; a1=[0.7,-0.5,1]; subplot(2,1,2)； zplane(b1,a1)； title(h2的零极点分布)；； 滤波器仿真和特性实验 设计实现一个在0.45处具有3dB截止角频率的一阶无限冲激响应低通滤波器和一个无限冲激响应高通滤波器，计算并画出他们各自的增益响应，并证明它们是全通互补和功率互补的。 [b,a]=butter(1,0.45); [h,w]=freqz(b,a,512); subplot(1,2,1)； plot(w/pi,20*log10(abs(h))); xlabel(频率/Hz);ylabel(幅值/dB)； title(IIR 低通滤波器)；grid； b1=[0.45];a1=[1 1]; [B A]=lp2hp(b1,a1,0.45); [h1,w]=freqz(B,A); subplot(1,2,2)； plot(w/pi,abs(h1)); xlabel(频率/Hz);ylabel(幅值/dB) title(IIR 高通滤波器)；grid； 4.作图分析 分别生成如下两个因果系统传输函数的零、极点图，并对该图进行研究，推断它们的稳定性。 b=[1 0 0];a=[1 -1.848 0.85]; b1=[1 0 0];a1=[1 -1.851 0.85]; figure(1); subplot(2,1,1)； zplane(b,a); title(系统1零极点分布图); [H1,w]=freqz