function?outdata?=?ZeroPfilter（WpWsApAsx）
%?Zero-Phase?Filter?base?on?buttor（基于巴特沃斯模擬濾波器的零相位濾波）
%-------------------------------------
%?outdata?=?ZeroPfilter（WpWsApAsx）
%?outdata:?the?filted?signal
%?[WpWsApAs]:?parameter?of?buttor?filter
%?x:?the?origin?signal
%?by?noWen（NNU_GIS）
%?EXAMPLE:
%?x?=?randn（4801）;
%?out?=?ZeroPfilter（0.05*pi0.1*pi115x）;
%?plot（x）;?hold?all?;plot（out）;
%?legend（‘oridata‘‘filterdata‘）;
%-------------------------------------

%－－－－－有限沖激響應(yīng)濾波器的邊界處理（此處使用巴特沃斯模擬濾波器，不用正理邊界）
%?N?=?length（x）;
%?M?=?30;?%n?of?the?filter
%?xx?=?ones（N+M1）;
%?for?i?=?1?:N
%?????xx（i）=x（i）;
%?end
%?for?i?=?1?:?M
%?????xx（N+i）=?2*xx（N）-xx（N-i）;
%?end
%??x?=?xx;
%--------------------------------------------------------------------------

%-------------design?of?buttor?filter（巴特沃斯模擬濾波器設(shè)計）---------------
%判斷濾波方式
if?length（Wp）==2?
????s?=?‘bandpass‘?;???%帶通
????else?if?Wp?>?Ws
????????s?=?‘high‘?;???%高通
????else?if?Wp?????????s?=?‘low‘??;???%低通???
????end
????end
end
Ap?=?Ap/2;?As?=?As/2;??%衰減速率減半（因為下面有兩次濾波處理）
[NWn]=buttord（Wp/piWs/piApAs）;?%計算巴特沃斯濾波器階次和截止頻率
[ba]=butter（NWns）;??????????????%頻率變換法設(shè)計巴特沃斯帶通濾波器
%--------------------------------------------------------------------------

%---------------零相位濾波實(shí)現(xiàn)----------------
y?=?filter（bax）;??%第一次濾波?
y?=?fliplr（y‘）‘;????%濾波結(jié)果反轉(zhuǎn)
y?=?filter（bay）;??%第二次濾波
y?=?fliplr（y‘）‘;????%再次反轉(zhuǎn)，即為最終結(jié)果

%------------邊界處理相關(guān)部分--------------
%?for?i?=?1?:N
%?????yy（i）=y（i）;
%?end
%?y?=?yy;
%----------------------------------------

outdata?=?y‘;??%結(jié)果輸出

end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????112963??2009-07-23?12:27??一種新的零相移數(shù)字濾波器的改進(jìn)算法.pdf

?????文件???????1775??2009-07-24?10:55??ZeroPfilter.m

-----------?---------??----------?-----??----

???????????????114738????????????????????2