//////////////////////////?該算法出自文獻(xiàn)：用于去除心電信號中工頻干擾數(shù)字濾波技術(shù).pdf（王鎮(zhèn)）?//////////////////////////

///////////////////////////////?FIR?filter?////////////////////////////////

///////?Band-Pass?Transfer?function:?H1（z）=1/Q^n*（1-z^（-R））^n/（1-z^（-P））^n
///////?Band-Pass?Frequence?Respose:?H1（w）=1/Q^n*sin（R/2*w）^n/sin（P/2*w）^n*exp（-j*（R-P）/2*n*w）
///////?Band-Stop?Transfer?function:?H（z）=z^（n*（R-P）/2）-1/Q^n*（1-z^（-R））^n/（1-z^（-P））^n
///////?Q=R/P
///////?If?SampleRate?fs=500Hz?and?Interference?Frequency?fg=50Hzthen?P=fs/fg=10
///////?For?a?good?band-stop?and?analysis?from?our?tests?we?set?n=2?and?R=160then?Q=16
///////?so?H（z）=z^（-150）-（1/16）^2*（（1-z^（-160））/（1-z（-10）））^2
///////?The?next?step?is?designing?a?self-adaptive?band-stop?filter????

//////////////////////////////////////////////////////////////////////////

/**************************************?Head?Files?頭文件?**********************************/
#include?
#include??
#include?
#include?
#include?

/**************************************?端口及器件地址設(shè)置?*********************************/
#define??PCF8591?0x90??//8bit-AD-DA集成芯片PCF8591的地址;?PCF8591地址=1001-A2-A1-A0-R/W?R=1W=0;?
???????????????????????//I2C寫數(shù)據(jù):發(fā)送地址-發(fā)送控制字-發(fā)送數(shù)據(jù)?I2C讀數(shù)據(jù):發(fā)送地址-讀數(shù)據(jù)


/**************************************?頭函數(shù)的聲明?***************************************/
bit?ADC_Input_set_addr_contr（unsigned?char?addrunsigned?char?contr）;???????????????????//準(zhǔn)接收ADC數(shù)據(jù)-----------------實(shí)現(xiàn)
unsigned?char?ADC_Input_RcvByteData（unsigned?char?addr）;????????????????????????????????//接收ADC數(shù)據(jù)-------------------實(shí)現(xiàn)
void?Start_Timer（void）;?????????????????????????????????????????????????????????????????//設(shè)置定時器--------------------實(shí)現(xiàn)
bit?Send_data_to_DAC（unsigned?char?addrunsigned?char?contr??unsigned?char?data）;??????//將8位數(shù)據(jù)發(fā)送到DAC------------實(shí)現(xiàn)
void?Band_Stop_Proc（void）;??????????????????????????????????????????????????????????????//帶阻濾波處理------------------實(shí)現(xiàn)

/*************************************?Variable?List?參數(shù)設(shè)置?******************************/
////////?fs必須是fg的正整數(shù)倍
//fs=250;?????????%sample?rate??????????????（can?change）
//fg=50;??????????%interference?frequency???（can?change）
//n=2;????????????%const????????????????????（don‘t?change）
//R=160;??????????%what?meaning?????????????（can?change）
//P=fs/fg;????????%what?meaning????（const）??（decided?by?fs?and?fg）
//Q=R/P;??????????%what?meaning????（const）??（decided?by?R?and?P）
const?unsigned?int?fs=500;???????????//為了達(dá)到實(shí)時性要求，fs<5kHz?為了不改變系數(shù)，fs<2500Hz?要求fs是fg的正整數(shù)倍
const?unsigned?int?fg=50;
const?unsigned?int?n=2;
const?unsigned?int?R=160;
const?unsigned?int?P=fs/fg;
const?unsigned?int?Q=R/P;

/************************************?變量定義?**********************************************/
static?unsigned?int?data_addr;????????????????????//數(shù)據(jù)序列的地址
static?unsigned?int?xdata?y[nfilt]={0};???????????//濾波后數(shù)據(jù)
static?unsigned?int?xdata?x[nfilt]={0};???????????//濾波前數(shù)據(jù)
bit?DataAddrLargerThan_nfilt_flag;???????//前nfilt個數(shù)據(jù)處理不同于后面的數(shù)據(jù)處理，這里設(shè)置轉(zhuǎn)變標(biāo)志??

/////////////?Band-Pass?/////

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

?????文件???????3500??2013-06-04?14:56??50Hz陷波器\pl_removal.m

?????文件???????5142??2013-05-23?00:19??50Hz陷波器\BandStop.m

?????文件???????1398??2013-05-12?11:07??50Hz陷波器\ideal_pl_main.m

?????文件???????2795??2013-06-03?19:08??50Hz陷波器\myfilter.m

?????文件???????5056??2013-06-03?21:08??50Hz陷波器\bandstopnew\BandStop.m

?????文件???????3992??2013-06-03?21:10??50Hz陷波器\bandstopnew\myfilter.m

?????文件???????5099??2013-06-04?01:17??50Hz陷波器\bandstopnew\inter_removal.c

?????文件???????3037??2013-06-04?14:43??50Hz陷波器\bandstopnew\pl_removalnew4.m

?????文件???????1552??2013-05-12?11:10??50Hz陷波器\bandstopnew\Band_Stop.m

?????文件????????189??2013-05-11?21:03??50Hz陷波器\bandstopnew\ideal_pl.m

?????文件??????12710??2013-05-15?21:58??50Hz陷波器\interf_removal.c

?????文件??????13118??2013-05-15?21:55??50Hz陷波器\BandStop.c

?????文件?????321828??2013-05-12?17:07??50Hz陷波器\用于去除心電信號中工頻干擾數(shù)字濾波技術(shù).pdf

?????文件????2389035??2013-05-08?16:28??50Hz陷波器\D+自適應(yīng)相干模板法在心電圖機(jī)中應(yīng)用的一種改進(jìn)算法.pdf

?????目錄??????????0??2013-06-03?19:09??50Hz陷波器\bandstopnew

?????目錄??????????0??2002-01-01?21:09??50Hz陷波器

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

??????????????2768451????????????????????16