hc722.m

数字信号处理的matlab程序

%《数字信号处理教程——MATLAB释义与实现》第七章例7.2.2程序hc722
% 线性相位FIR滤波器的类型及符幅特性
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
clear
h1=[3,-1,2,-3,5,-3,2,-1,3];             % 输入一个类型I的线性相位滤波器系数
[A1,w1,typea,tao1]=amplres(h1);         % 用amplres函数求符幅特性并判别其类型
subplot(2,2,1),plot(w1,A1),grid on      % 画出符幅特性
title('线性相位滤波器类型I的符幅特性')
h2= [3,-1,2,-3,-3,2,-1,3];              % 输入一个类型II的线性相位滤波器系数
[A2,w2,typeb,tao2]=amplres(h2);         % 用amplres函数求符幅特性并判别其类型
subplot(2,2,2),plot(w2,A2),grid on      % 画出符幅特性
title('线性相位滤波器类型II的符幅特性')
h3=[3,-1,2,-3,0, 3, -2, 1,-3];          % 输入一个类型III的线性相位滤波器系数
[A3,w3,typec,tao3]=amplres(h3);         % 用amplres函数求符幅特性并判别其类型
subplot(2,2,3),plot(w3,A3),grid on      % 画出符幅特性
title('线性相位滤波器类型III的符幅特性')
h4=[3,-1,2,-3,03, -2, 1,-3];            % 输入一个类型IV的线性相位滤波器系数
[A4,w4,typed,tao4]=amplres(h4);         % 用amplres函数求符幅特性并判别其类型
subplot(2,2,4),plot(w4,A4),grid on      % 画出符幅特性
title('线性相位滤波器类型IV的符幅特性')
set(gcf,'color','w');                   % 置图形背景色为白