资源描述：
《归一化评定定量模型》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实用标准文案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
2、.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(
3、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(
4、[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)’)；s
5、ubplot(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，π]的幅度响应,单位
6、是伏;%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
7、.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('相频特性')
8、;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为1精彩文档实用标准文案n=(0:499)*pi/500;p1=[0.2];%极点在0.2处[b1,a1]=zp2tf(q,p1,k);%由zp
