%%?入門?例子
y2=mvnrnd（011000?）;%%?高斯白噪聲
x=0.01:0.01:10;
y1=10*sin（pi*x）;
y=y1‘+y2??????????????%%含有噪聲的信號
plot（xy）
[ba]=butter（60.06）;?%%[b?a]=butter（nw0）b為濾波函數的分子系數，a為分母系數，n為濾波器階次，w0為相對（歸一化）截止頻率（cutoff?frequency）
??????????????????????%%實際截止頻率fs除以f采樣頻率*2，單位為pi*rad/sample，類比角速度w
%?y0=filter（bay）;?????%%y0為濾波后的結果，一般濾波
y0=filtfilt（bay）;?????%%?零相位濾波

figure（2）
plot（xy0‘r‘xy）

%%?巴特沃斯濾波器，給定截止頻率和階數
%%?signal
x=0.01:0.01:100;
y1=10*sin（0.2*pi*x）;
y2=sin（2*pi*x）;
y=y1+y2;
plot（xy）
%%?filter
[ba]=butter（50.006）;%?0.3Hz為截止頻率，w0為0.01（0.5Hz）作為截止頻率時，效果更好，截止頻率越靠近信號頻率，
???????????????????????%?是以犧牲相位為代價的。
y0=filter（bay）;
figure（1）;
plot（xy0）
%%?over


%%?巴特沃斯濾波器設計階數和截止頻率，不知道用幾階，只知道大致頻率。待解決衰減速度
%%?signal
x=0.01:0.01:100;
y1=10*sin（0.2*pi*x）;
y2=sin（2*pi*x）;
y=y1+y2;
plot（xy）
%%?filter
wp=2*0.3/100;
ws=2*0.5/100;%計算方法同w0通帶為0.3Hz，阻帶為0.5Hz
rp=3;
rs=6;%最大通帶衰減為3db最小阻帶衰減為6db
[nw0]=buttord（wpwsrprs）;%a計算濾波器的階數和截止頻率
[ba]=butter（nw0）
y0=filter（bay）;
plot（xy0）
%%??over
?

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

?????文件???????1295??2017-03-09?22:02??巴特沃斯濾波器\filter_pra.m

?????文件??????12626??2017-03-07?17:52??巴特沃斯濾波器\untitled.slx.autosave

?????目錄??????????0??2018-01-08?01:59??巴特沃斯濾波器

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

????????????????13921????????????????????3