实验三lst离散系统的频域分析.docx

武汉工程大学信号分析与处理实验四 专业： 通信02班 学生姓名： 李瑶华 学号： 1304200113 完成时间：2016年6月1日实验三: ＬＳＩ离散系统的频域分析实验目的通过在频域中仿真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 低通滤波器),gridb1=[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 高通滤波器),grid4.作图分析 分别生成如下两个因果系统传输函数的零、极点图，并对该图进行研究，推断它们的稳定性。 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(b,a);subplot(2,1,2) plot(w/pi,abs(H1));axis([0 1 -100 600]);title(频率响应);xlabel(w);figure(2);subplot(2,1,1) zplane(b1,a1);title(系统2零极点分布图);b1=[1 0 0];a1=[1 -1.851 0.85];[H2,w]=freqz(b1,a1);subplot(2,1,2) plot(w/pi,abs(H2