资源描述:
《归一化评定定量模型》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、freqz功能:用于求解离散时间系统的频率响应函数H(e)例3-1已知离散时间系统的系统函数为b=[。。。];a=[1,0,,0,,0,];-120-140freqz(b,a);-20-60-80-100F0.3,T0T0■wTENormalizedFrequency(rad/sample)dnMB-40例3-2已知离散时间系统的系统函数,求该系统在0〜n频率范围内归一化的绝对幅度频率响应与相位频率响应。解MATLAB程序如下:b=;a=[1,—,,一,];w=(0:499)pi/500;%在pi的范围内取500个采样点h=freqz(b,a,w);%求系统的频率响应s
2、ubplot(2,1,1),plot(w/pi,abs(h));grid%作系统的幅度频响图axis([0,1,min(abs(h)),max(abs(h))]);ylabel('幅度’);subplot(2,1,2),plot(w/pi,angle(h));grid%作系统的相位频响图axis([0,1,min(angle(h)),max(angle(h))]);ylabel('相位’);xlabel('以pi为单位的频率’);0.8——i--,一0.60.4 .2———--00.10.20.30.40.50.60.70.80.912■0-2Ja,-100.10.2
3、0.30.4以pi为0.5单位0.6的频率0.70.80.911例3-3已知离散时间系统的系统函数为解MATLAB程序如下:w1=(0:499)pi/500;[h,w]=freqz(b,a,w1);db=20log10(abs(h));%求系统的相对幅频响应值subplot(2,2,1),plot(w/pi,abs(h));grid%作系统的绝对幅度频响图axis([0,1,min(abs(h)),max(abs(h))]);title('幅频特性(V)');subplot(2,2,2),plot(w/pi,angle(h));grid%作系统的相位频响图axis([0
4、,1,min(angle(h)),max(angle(h))]);title('相频特性’);subplot(2,2,3),plot(w/pi,db);grid%作系统的相对幅度频响图axis([0,1,—100,5]);title('幅频特性(dB)');subplot(2,2,4),zplane(b,a);%乍零极点分布图幅频特性(V)12o-10000.51o1-trapw7aMHLaaml1-1012RealPart例3-4已知离散时间系统的系统函数为function[db,mag,pha,grd,w]=freqz_m(b,a,N);喇散系统响应求解函数,包括幅
5、度响应、相位响应、群时延%a,b,N为输入参数%b为系统分子多项式系数向量;%b为系统分子多项式系数向量;N为频率离散化的点数%db,mag,pha,grd,w为输出参数;%db为系统[0,兀]的度响应,单位是dB;%ag也是系统[0,兀]的幅度响应,单位是伏;%pha为系统的[0,兀]上的相位响应,单位为rad;%grd为[0,兀]上的群延时响应;%w为[0,兀]上的频率采样点.[H,w]=freqz(b,a,N,'whole');H=(H(1:N/2));w=(w(1:N/2));mag=abs(H);db=20log10((mag+eps)/max(mag));p
6、ha=angle(H);grd=grpdelay(b,a,w);b=[。。。];a=[1,0,,0,,0,];N=1024;[db,mag,pha,grd,w]=freqz_m(b,a,N);作绝对幅度频响图subplot(2,2,1),plot(w/pi,mag);grid%axis([0,1,min(mag),max(mag)]);title('幅频特性(V)');subplot(2,2,2),plot(w/pi,pha);grid%作相位频响图axis([0,1,min(pha),max(pha)]);title('相频特性’);subplot(2,2,3),pl
7、ot(w/pi,db);grid%作相对幅度频响图axis([0,1,-100,5]);title('幅频特性(dB)');subplot(2,2,4),plot(w/pi,grd);grid%作系统的群迟延图title('群迟延’);例3-5观察系统极点的位置对幅频响应的影响。q=[0];k=1;%设零点在原点处,k为1n=(0:499)pi/500;p1=[];%极点在处[b1,a1]=zp2tf(q,p1,k);%由zpk模式求tf模式b和a系数[h1,w]=freqz(b1,a1,n);%求系统的频率响应subplot(2,3,1