归一化评定定量模型.docVIP

freqz 功能：用于求解离散时间系统的频率响应函数H(ejω)。 例3-1 已知离散时间系统的系统函数为 b=[0.1321,0,-0.3963,0,0.3963,0,-0.1321]; a =[1,0,0.34319,0,0.60439,0,0.20407]; freqz(b,a); 例3-2 已知离散时间系统的系统函数，求该系统在0～π频率范围内归一化的绝对幅度频率 响应与相位频率响应。 解 MATLAB 程序如下： b＝［0.2，0.1，0.3，0.1，0.2］； a＝［1，－1.1，1.5，－0.7，0.3］； w＝(0：499)*pi/500； %在pi 的范围内取500 个采样点 h＝freqz(b，a，w)； %求系统的频率响应 subplot(2，1，1)，plot(w/pi，abs(h))；grid %作系统的幅度频响图 axis(［0，1，1.1*min(abs(h))，1.1*max(abs(h))］)； ylabel(‘幅度’)； subplot(2，1，2)，plot(w/pi，angle(h))；grid %作系统的相位频响图 axis(［0，1，1.1*min(angle(h))，1.1*max(angle(h))］)； ylabel(‘相位’)；xlabel(‘以pi 为单位的频率’)； 例3-3 已知离散时间系统的系统函数为 解MATLAB 程序如下： b＝［0.1，－0.4，0.4，－0.1］； a＝［1，0.3，0.55，0.2］； w1＝(0：499)*pi/500； [h，w]＝freqz(b，a，w1)； db＝20*log10(abs(h))； %求系统的相对幅频响应值 subplot(2，2，1)，plot(w/pi，abs(h))；grid %作系统的绝对幅度频响图 axis([0，1，1.1*min(abs(h))，1.1*max(abs(h))])； title(‘幅频特性(V)’)； subplot(2，2，2)，plot(w/pi，angle(h))；grid %作系统的相位频响图 axis([0，1，1.1*min(angle(h))，1.1*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)；%作零极点分布图 例3-4 已知离散时间系统的系统函数为 function[db,mag,pha,grd,w]=freqz_m(b,a,N); %离散系统响应求解函数,包括幅度响应、相位响应、群时延; %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=20*log10((mag+eps)/max(mag)); pha=angle(H); grd=grpdelay(b,a,w); b=[0.1321,0,0.3963,0,0.3963,0,0.1321]; a=[1,0,-0.34319,0,0.60439,0,-0.20407]; N=1024; [db,mag,pha,grd,w]=freqz_m(b,a,N); subplot(2,2,1),plot(w/pi,mag);grid %作绝对幅度频响图 axis([0,1,1.1*min(mag),1.1*max(mag)]); title(幅频特性(V)); subplot(2,2,2),plot(w/pi,pha);grid%作相位频响图 axis([0,1,1.1*min(pha),1.1*max(pha)]); title(相频特性); subplot(2,2,3),plot(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