《实验4-离散系统的变换域分析(共6页)》由会员分享,可在线阅读,更多相关《实验4-离散系统的变换域分析(共6页)(6页珍藏版)》请在金锄头文库上搜索。
1、精选优质文档-倾情为你奉上实验4 离散系统的变换域分析一、实验目的:加深对离散系统的频率响应分析和零、极点分布的概念理解。二、实验原理:离散系统的时域方程为其变换域分析方法如下:频域 系统的频率响应为 Z域 系统的转移函数为 分解因式 ,其中和称为零、极点。 在MATLAB中,可以用函数z,p,K=tf2zp(num,den)求得有理分式形式的系统转移函数的零、极点,用函数zplane(z,p)绘出零、极点分布图;也可以用函数zplane(num,den)直接绘出有理分式形式的系统转移函数的零、极点分布图。使h=freqz(num,den,w)函数可求系统的频率响应,w是频率的计算点,如w=0
2、:pi/255:pi, h是复数,abs(h)为幅度响应,angle(h)为相位响应。另外,在MATLAB中,可以用函数 r,p,k=residuez(num,den)完成部分分式展开计算;可以用函数sos=zp2sos(z,p,K)完成将高阶系统分解为2阶系统的串联。例1 求下列直接型系统函数的零、极点,并将它转换成二阶节形式 解 用MATLAB计算程序如下: num=1 -0.1 -0.3 -0.3 -0.2;den=1 0.1 0.2 0.2 0.5;z,p,k=tf2zp(num,den);m=abs(p);disp(零点);disp(z);disp(极点);disp(p);disp(
3、增益系数);disp(k);sos=zp2sos(z,p,k);disp(二阶节);disp(real(sos);zplane(num,den)输入到“num”和“den”的分别为分子和分母多项式的系数。计算求得零、极点增益系数和二阶节的系数:零点 0.9615 -0.5730 -0.1443 + 0.5850i -0.1443 - 0.5850i极点 0.5276 + 0.6997i 0.5276 - 0.6997i -0.5776 + 0.5635i -0.5776 - 0.5635i 增益系数 1 二阶节 1.0000 -0.3885 -0.5509 1.0000 1.1552 0.65
4、11 1.0000 0.2885 0.3630 1.0000 -1.0552 0.7679 系统函数的二阶节形式为: 极点图如右图。 例2 差分方程 所对应的系统的频率响应。 解:差分方程所对应的系统函数为用MATLAB计算的程序如下:k=256;num=0.8 -0.44 0.36 0.02;den=1 0.7 -0.45 -0.6;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);plot(w/pi,im
5、ag(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(弧度) 三、实验内容:1、求系统 的零、极点和幅度频率响应和相位响应。num=0.0528 0.0797 0.1295 0.1295 0.797 0.0528;den=1 -1.8107 2.4947 -
6、1.8801 0.9537 -0.2336;z,p,k=tf2zp(num,den); disp(零点);disp(z);disp(极点);disp(p);disp(增益系数);disp(k);figure(1)zplane(num,den)figure(2)freqz(num,den,128)2、对于传输函数,求脉冲响应和阶跃相应num=5;den=1 -1 0.09; subplot(2,1,1)%前40个点的冲激响应impz(num,den,40);h,t=impz(num,den,40); title(脉冲响应)step=ones(1,100);subplot(2,1,2)y=conv
7、(h,step);n=0:138;stem(n,y,filled)title(阶跃响应)3、滤波器的传输函数为:不计算,画出幅度响应曲线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);plot(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(弧度)专心-专注-专业