原名：数字信号处理基础及MATLAB实现,解压缩密码：lfjk2002@163.com

				function [Hr,w,c,L]=hr_type3(h);
				%线性相位滤波器类型3，计算滤波器振幅响应
				%Hr振幅响应；w在[0,2pi]之间计算Hr的500个频率点；
				%c为3型低通滤波器系数；L为Hr的阶次；h为3型低通滤波器的脉冲响应。
				M=length(h);L=(M-1)/2;
				c=[ 2*h(L+1:-1:1)];   %1乘(L+1)行向量
				n=[0:1:L];                %(L+1)乘1列向量
				w=[0:1:500]'*2*pi/500;
				Hr=sin(w*n)*c';