%-----------------帶陷濾波器抑制工頻干擾-------------------??
%50Hz陷波器：由一個低通濾波器加上一個高通濾波器組成??
%而高通濾波器由一個全通濾波器減去一個低通濾波器構成??

Me=100;???????????????%濾波器階數??
L=100;????????????????%窗口長度??
beta=100;?????????????%衰減系數??
Fs=1500;??
wc1=49/Fs*pi;?????%wc1為高通濾波器截止頻率，對應51Hz??
wc2=51/Fs*pi?????;%wc2為低通濾波器截止頻率，對應49Hz??
h=ideal_lp（0.132*piMe）-ideal_lp（wc1Me）+ideal_lp（wc2Me）;?%h為陷波器?????????????????????????????????????????????????????????
w=kaiser（Lbeta）;??
y=h.*rot90（w）;?????????%y為50Hz陷波器沖擊響應序列??
m2=filter（y1m）;??

figure??
subplot（211）;plot（abs（h））;axis（[0?100?0?0.2]）;??
xlabel（‘頻率（Hz）‘）;ylabel（‘幅度（mv）‘）;title（‘陷波器幅度譜‘）;grid;??
N=512;??
P=10*log10（abs（fft（y）.^2）/N）;??
f=（0:length（P）-1）;??
subplot（212）;plot（fP）;??
xlabel（‘頻率（Hz）‘）;ylabel（‘功率（dB）‘）;title（‘陷波器功率譜‘）;grid;??
?????
figure??
subplot?（211）;?plot（TIMEm）;??
xlabel（‘t（s）‘）;ylabel（‘幅值‘）;title（‘原始信號‘）;grid;??
subplot（212）;plot（TIMEm2）;??
xlabel（‘t（s）‘）;ylabel（‘幅值‘）;title（‘帶阻濾波后信號‘）;grid;??
????
figure??
N=512?;?
subplot（211）;plot（abs（fft（m））*2/N）;axis（[0?100?0?1]）;??
xlabel（‘t（s）‘）;ylabel（‘幅值‘）;title（‘原始信號頻譜‘）;grid;??
subplot（212）;plot（abs（fft（m2））*2/N）;axis（[0?100?0?1]）;??
xlabel（‘t（s）‘）;ylabel（‘幅值‘）;title（‘帶阻濾波后信號頻譜‘）;grid;????
??

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

?????文件???????2451??2018-08-14?21:34??心電信號預處理\100數據以及仿真結果\data.mat

?????文件??????10917??2018-08-14?21:09??心電信號預處理\100數據以及仿真結果\IIR零相移數字濾波器糾正基線漂移.png

?????文件??????10669??2018-08-14?21:55??心電信號預處理\100數據以及仿真結果\低通濾波器濾除肌電信號1.png

?????文件???????9897??2018-08-14?21:56??心電信號預處理\100數據以及仿真結果\帶陷濾波器抑制工頻干擾.png

?????文件??????10552??2018-08-14?22:02??心電信號預處理\100數據以及仿真結果\心電信號預處理.png

?????文件???????1408??2018-08-14?21:56??心電信號預處理\Filter.m

?????文件???????1027??2018-09-10?19:32??心電信號預處理\Filter1.m

?????文件????????218??2018-08-14?20:19??心電信號預處理\ideal_lp.m

?????文件???????1018??2018-08-14?20:35??心電信號預處理\IIR零相移數字濾波器糾正基線漂移.txt

?????文件??????91660??2018-08-14?22:03??心電信號預處理\result.mat

?????文件???????1927??2018-09-10?19:29??心電信號預處理\Test.m

?????文件???????1944??2018-08-14?15:52??心電信號預處理\低通濾波器濾除肌電信號（Test）.txt

?????文件???????1695??2018-08-14?20:14??心電信號預處理\帶陷濾波器抑制工頻干擾（Filter+ideal_lp）.txt

?????目錄??????????0??2018-08-22?19:24??心電信號預處理\100數據以及仿真結果

?????目錄??????????0??2018-08-22?19:24??心電信號預處理

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

???????????????145383????????????????????15