firh.m

来自「使用Matlab实现了FIR与IIR滤波器」· M 代码 · 共 14 行

function [ output_args ] = Untitled1( input_args )
%UNTITLED1 Summary of this function goes here
%  Detailed explanation goes here
ft=22050;
fcuts=[4800 5000];
mags = [0 1];
devs = [10^(-100/20),(1-10^(-1/20))];
[N,Wn,beta,ftype] = kaiserord(fcuts,mags,devs,ft);
kw = kaiser(N+1,beta);
b = fir1(N,Wn,'high',kw);
[h,w] = freqz(b,1);
plot(w/pi,20*log10(abs(h)));grid;
xlabel('w/pi'); ylabel('Gain, dB');
title('Fir Highpass Filter');