這是之前用到的幾個簡單的濾波器的程式，濾波器的設計比較簡單。可以直接使用matlab裡的FDA，匯出濾波器係數即可。

再就是參考數字訊號處理課本里面三種濾波器設計的方法，由過渡帶和所選的窗函式推匯出濾波器階數，進一步根據firl和filter函式編寫即可。

下面給出三種濾波器的簡單程式碼：

低通：用的hanning窗

function [y]=fir_lpf(x,Fs,fs,fp)
%該函式採用hanning窗實現低通濾波
%x為輸入訊號，fs，為取樣頻率
%fs,fp分別為阻帶截至頻率和通帶截至頻率
%20150615 by boat


%求對應角頻率
ws=fs*2*pi/Fs;
wp=fp*2*pi/Fs;


%求濾波器的最小階數
B=ws-wp;   %過渡頻寬取二者中較小的
N=ceil(6.2*pi/B);


%計算濾波器係數
wc=(ws+wp)/2;
hn=fir1(N-1,wc,hanning(N));
y=filter(hn,1,x);


freqz(hn)

帶通：這裡用的布拉克曼窗

function [y]=fir_bpf(x,fs,fs2,fp2,fp1,fs1)
%該函式採用blackman窗實現帶通濾波
%x為輸入訊號，fs，為取樣頻率
%fs2,fp2分別為阻帶上截至頻率和通帶上截至頻率
%fp1，fs1分別為通帶下截止頻率和阻帶下截至頻率
%ps：輸入時以上四個濾波引數按從大到小輸入即可
%20150615 by boat


%求對應角頻率
ws2=fs2*2*pi/fs;
wp2=fp2*2*pi/fs;
wp1=fp1*2*pi/fs;
ws1=fs1*2*pi/fs;


%求濾波器的階數
B=min(ws2-wp2,wp1-ws1);   %求過渡頻寬
N=ceil(12*pi/B);


%計算濾波器係數
wc2=(ws2+wp2)/2;
wc1=(ws1+wp1)/2;
wp=[wc1,wc2];
hn=fir1(N-1,wp,blackman(N));
y=filter(hn,1,x);

高通：這裡用的凱澤窗

function [y]=fir_hpf(x,Fs,As,fp,fs)
%該函式採用kaiser窗設計高通濾波器
%輸入：x 輸入訊號
%      Fs 取樣頻率
%      As 阻帶最小衰減 
%      fp通帶截止頻率 fs阻帶截止頻率（單位Hz）
%      20150615 by boat


%設定各種引數
b=fp-fs;                               % 求取過渡帶
M0=round((As-7.95)/(14.36*b/Fs))+2;    % 計算窗長
M=M0+mod(M0+1,2);                      % 保證窗長為奇數
wp=2*fp/Fs*pi; ws=2*fs/Fs*pi;          
wc=(wp+ws)/2;                          % 求截止頻率


%求kaiser窗的係數beta的值
if As>50
    beta=0.1102*(As-8.7);
  elseif As>=21&&As<=50
    beta=0.5842*(As-21)^0.4+0.07886*(As-21);
  else
    beta=0;
end




%設定濾波器引數
N=M-1;                                 %設定濾波器階數
hd=fir1(N,wc,'high',kaiser(M,beta));


%濾波
x=x-mean(x);                            % 消除直流分量
                           
y=filter(hd,1,x);