%-----------------高通濾波器抑制肌電干擾-------------------??
fp=35;
fs=50;??
ap=3;????%通帶最大衰減?
as=20;?????????????????????????????%阻帶最小衰減?
fs1=200;??????????????????????????%抽樣頻率?
wp=2*pi*fp/fs1;???????
ws=2*pi*fs/fs1;?????????????????????%數(shù)字角頻率??
Fs=fs1/fs1;
T=1/Fs;??
Op=2/T*tan（wp/2）;?????????
Os=2/T*tan（ws/2）;???????????????????%角頻率畸變?
[N?Wn]=buttord（OpOsapas‘s‘）;??
[z?p?k]=buttap（N）;????????????????????%G（p）的零點(diǎn)極點(diǎn)增益?
[b?a]=zp2tf（zpk）;????????????????????%G（p）的分子、分?jǐn)?shù)系數(shù)?
[B?A]=lp2lp（baOp）;??????????????????%H（s）的分子、分?jǐn)?shù)系數(shù)?
[Bz?Az]=bilinear（BAFs*T/2）;??????????%H（z）的分子、分?jǐn)?shù)系數(shù)
[H?w]=freqz（BzAz256Fs*fs1）;????????%H?
plot（wabs（H）‘r‘）;?
title（‘低通濾波器‘）;?
grid?on;?
data=xlsread（‘maibo1.xlsx‘）;
x=data（:1）;
n=0:1/200:3-1/200;
y=filter（BzAzx）;
figure
plot（ny（9001:9600））;


%-----------------帶陷濾波器抑制工頻干擾-------------------??
%50Hz陷波器：由一個低通濾波器加上一個高通濾波器組成??
%而高通濾波器由一個全通濾波器減去一個低通濾波器構(gòu)成??
Me=100;???????????????%濾波器階數(shù)??
L=100;????????????????%窗口長度??
beta=100;?????????????%衰減系數(shù)??
Fs=200;??
wc1=49/Fs*pi;?????%wc1為高通濾波器截止頻率，對應(yīng)51Hz??
wc2=51/Fs*pi?????;%wc2為低通濾波器截止頻率，對應(yīng)4