hr_type1.m

来自「matla窗口制作」· M 代码 · 共 16 行

%____________________h(n)为偶对称h(n)=h(N-n-1)，且N为奇数(Ⅰ型滤波器)__________
function [Hr,w,a,L]=hr_type1(h);         
%计算所设计的1型滤波器的振幅响应
%Hr=振幅响应
%a=1型滤波器的系数
%L=Hr的阶次
%h=1型滤波器的单位冲击响应
M=length(h);
L=(M-1)/2;
a=[h(L+1) 2*h(L:-1:1)];
n=[0:1:L];
w=[0:1:500]'*2*pi/500;                %把整个横坐标取500个采样点
Hr=cos(w*n)*a';                       %振幅响应函数
%_____________________幅度特点是对w=0,pi,2*pi是偶对称的_________________________