资源描述:
《数字信号处理-实验4 离散系统的变换域分析》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、实验4离散系统的变换域分析一、实验目的1、熟悉对离散系统的频率响应分析方法;2、加深对零、极点分布的概念理解。二、实验原理离散系统的是时域方程为其变换域分析方法如下:频域:系统的频率响应为:Z域:系统的转移函数为:分解因式:其中和称为零、极点。三、实验要求:编程实现系统参数输入,绘出幅度频率响应和相位曲线和零、极点分布图。一、求系统num=[0.05280.07970.12950.12950.7970.0528];den=[1-1.81072.4947-1.88010.9537-0.2336];[z,p,k]=tf2zp(
2、num,den);disp('零点');disp(z);disp('极点');disp(p);disp('增益系数');disp(k);figure(1)zplane(num,den)figure(2)freqz(num,den,128)4.12求滤波器输出的前十个采样值;num=[0.3-0.250.1];den=[100];[z,p,k]=tf2zp(num,den);m=abs(p);disp('零点');disp(z);disp('极点');disp(p);[z,p,k]=tf2zp(num,den);m=abs(
3、p);disp('零点');disp(z);disp('极点');disp(p);subplot(2,2,1);zplane(num,den);title('零极点图');k=256;w=0:pi/k:pi;h=freqz(num,den,w);subplot(2,2,2);subplot(2,2,3);plot(w/pi,abs(h));gridtitle('幅度谱')xlabel('omega/pi');ylabel('幅值')subplot(2,2,4);plot(w/pi,angle(h));gridtitle
4、('相位谱')xlabel('omega/pi');ylabel('弧度')clear;B=[0.3-0.250.1];A=1;n=0:9;X=[2*(n==0)-(n==2)];Y=filter(A,B,X);stem(n,Y,'.');title('滤波器输出采样值')6.22对于传输函数求脉冲响应和阶跃响应num=[5];den=[1-10.09];subplot(2,1,1)impz(num,den,40);[h,t]=impz(num,den,40);title('脉冲响应')step=ones(1,100)
5、;subplot(2,1,2)y=conv(h,step);n=0:138;stem(n,y,'filled')title('阶跃响应')7.26滤波器的传输函数为:H不计算,画出幅度响应曲线。k=256;num=[0,0,1,1];den=[1,-1,0,0];w=0:pi/k:pi;h=freqz(num,den,w);subplot(2,2,1);plot(w/pi,real(h));gridtitle('实部');xlabel('omega/pi');ylabel('幅度')subplot(2,2,2);plo
6、t(w/pi,imag(h));gridtitle('虚部');xlabel('omega/pi');ylabel('Amplitude')subplot(2,2,3);plot(w/pi,abs(h));gridtitle('幅度谱');xlabel('omega/pi');ylabel('幅值')subplot(2,2,4);plot(w/pi,angle(h));gridtitle('相位谱');xlabel('omega/pi');ylabel('弧度')7.28由七项滑动平均滤波器的零极点推断滤波器形
7、状解:由七项滑动滤波器特点,写出其对应传输函数num=[1111111];den=[7000000];[z,p,k]=tf2zp(num,den);disp('零点');disp(z);disp('极点');disp(z);disp('增益系数');disp(k);figure(1)zplane(num,den)figure(2)freqz(num,den,128)