filterecg.m

是

% FIR LP & HP filters (cascaded): 70dB 0.05-40Hz 1dB ripple 
% sampling frequency zm4NUJ  
function[filtdata] = filterECG(data)
% sampling frequency zm4NUJ
Fs= 256; % Hz
% Data should be zeromeaned before calling this function 
data = data-mean(data); 
% set up transfer function 
denominator = 1;
numerator_lp =[ 
    0.00116294818892327 
    0.0108214392935653 
    0.0200343666593651
    0.0261187881634098
    0.0212816864030144
    0.00490488877763669  
    -0.0145498515998796
    -0.0232788926070469
    -0.0124182936123568
    0.0126862162067748
    0.0328393141039891 
    0.0277014738094466
    -0.00663354198352208 
    -0.0487259025516825
    -0.0610853025450089
    -0.0144549204767754
    0.0886474308596396
      0.208814361863424
      0.289332487838881
      0.289332487838881
      0.208814361863424
    0.0886474308596396
    -0.0144549204767754
    -0.0610853025450089
    -0.0487259025516825
    -0.00663354198352208
    0.0277014738094466
    0.0328393141039891
    0.0126862162067748
    -0.0124182936123568
    -0.0232788926070469
    -0.0145498515998796
    0.00490488877763669
    0.0212816864030144
    0.0261187881634098
    0.0200343666593651
    0.0108214392935653
    0.00116294818892327];

numerator_hp = [
    -0.007719806995443
  -0.000593529464503111
  -0.000616537270116216
  -0.000639451830525362
  -0.000663561635529363
  -0.000687428482340838
  -0.000712345353053305
  -0.000737050063062293
  -0.000762827770843517
  -0.000788276159757994
  -0.000814687362552657
  -0.000840693466317883
  -0.000867737622092878
    -0.00089438708671706
    -0.00092216902638079
  -0.000949570338620629
  -0.000978288661010672
    -0.00100668160625124
    -0.00103659671861193
    -0.00106615479623852
    -0.00109736040348402
    -0.00112792836635087
    -0.00116004251486604
    -0.00119092382531205
    -0.00122322204449099
    -0.00125354739873849
    -0.00128545618449105
    -0.00131494126741691
    -0.00134807599742105
    -0.00138005083638085
    -0.00142229273289703
    -0.00145905623888122
    -0.00147988572401479
    -0.00152337012050789
    -0.00155509077331504
    -0.0015924435019478
    -0.00162645093636032
    -0.00166360336918217
    -0.00169886271207229
    -0.00173643201108726
    -0.00177256853414649
    -0.00181053087910452
    -0.00184738932024739
    -0.00188580324032983
    -0.00192320422309985
    -0.00196187823931274
    -0.00199959655075123
    -0.00203849112299403
    -0.0020765587432622
    -0.00211583028869268
    -0.00215442956472851
    -0.00219430023605191
    -0.00223352552719682
    -0.00227395565365889
    -0.00231354279813875
    -0.002354108014004
    -0.00239338505355968
    -0.00243353302728908
    -0.00247248705230046
    -0.00251352618627335
    -0.00255409536944302
    -0.00259725994910034
    -0.00263509464305156
    -0.00267536483006731
    -0.00271826154742389
    -0.00275770700056424
    -0.00279927582669419
    -0.00283951499715645
    -0.00288073484890323
    -0.00292129482201999
    -0.00296233875603428
    -0.00300290150021605
    -0.00304386170207064
    -0.00308454857127824
    -0.00312549340035217
    -0.00316610295239397
    -0.00320679774729835
    -0.0032471413328337
    -0.00328750802064282
    -0.00332753344802601
    -0.00336758426155901
    -0.00340736140579608
    -0.00344722643493981
    -0.00348685648439877
    -0.00352656997497068
    -0.00356591121574338
    -0.00360510516351791
    -0.00364360588016825
    -0.00368197568241012
    -0.00371996751152801
    -0.00375850930844281
    -0.00379655172798984
    -0.00383430839616961
    -0.00387004946444441
    -0.00390808220759818
    -0.00394482515495904
    -0.00398066969639526
    -0.00401692280951308
    -0.00405221772052505
    -0.00408773437155447
    -0.00412235016462235
    -0.00415686418023859
    -0.00419066438491442
    -0.00422445487363778
    -0.0042576163393507
    -0.00429061997144701
    -0.00432291749710377
    -0.00435498997501222
    -0.00438636173060546
    -0.00441745251390048
    -0.00444780506237289
    -0.00447784887049958
    -0.00450718961782266
    -0.00453627650609555
    -0.00456472176460446
    -0.0045929170622211
    -0.00462034686518763
    -0.00464732123257699
    -0.0046733867508092
    -0.00469915525083144
    -0.00472429941522134
    -0.00474929514447324
    -0.00477312604512227
    -0.00479636638243965
    -0.00481865942392797
    -0.00484177134432678
    -0.00486283917051695
    -0.0048838439863828
    -0.00490407570971269
    -0.00492386545165352
    -0.00494290018123188
    -0.0049611304645503
    -0.00497862211466973
    -0.00499556750358641
    -0.00501186356380093
    -0.00502752874153267
    -0.00504241401877351
    -0.00505665927887701
    -0.00507016713131965
    -0.00508305015409715
    -0.00509513775950472
    -0.00510653645526546
    -0.00511710891690786
    -0.00512701904388992
    -0.00513616850829649
    -0.00514472125325565
    -0.00515250765464019
    -0.00515958960068915
    -0.00516579676427969
    -0.00517132876190234
    -0.00517617853236731
    -0.00518043811157794
    -0.00518382726123044
    -0.00518621003962965
    -0.00518794262187708
    -0.00518900154305289
      0.994810133061056
    -0.00518900154305289
    -0.00518794262187708
    -0.00518621003962965
    -0.00518382726123044
    -0.00518043811157794
    -0.00517617853236731
    -0.00517132876190234
    -0.00516579676427969  
    -0.00515958960068915
    -0.00515250765464019
    -0.00514472125325565
    -0.00513616850829649
    -0.00512701904388992
    -0.00511710891690786
    -0.00510653645526546
    -0.00509513775950472
    -0.00508305015409715
    -0.00507016713131965
    -0.00505665927887701
    -0.00504241401877351
    -0.00502752874153267
    -0.00501186356380093
    -0.00499556750358641
    -0.00497862211466973
    -0.0049611304645503 
    -0.00494290018123188
    -0.00492386545165352
    -0.00490407570971269
    -0.0048838439863828
    -0.00486283917051695
    -0.00484177134432678
    -0.00481865942392797
    -0.00479636638243965
    -0.00477312604512227
    -0.00474929514447324
    -0.00472429941522134
    -0.00469915525083144
    -0.0046733867508092
    -0.00464732123257699
    -0.00462034686518763
    -0.0045929170622211
    -0.00456472176460446
    -0.00453627650609555
    -0.00450718961782266
    -0.00447784887049958
    -0.00444780506237289
    -0.00441745251390048
    -0.00438636173060546
    -0.00435498997501222
    -0.00432291749710377
    -0.00429061997144701
    -0.0042576163393507
    -0.00422445487363778
    -0.00419066438491442
    -0.00415686418023859
    -0.00412235016462235
    -0.00408773437155447
    -0.00405221772052505
    -0.00401692280951308
    -0.00398066969639526
    -0.00394482515495904
    -0.00390808220759818
    -0.00387004946444441
    -0.00383430839616961
    -0.00379655172798984
    -0.00375850930844281
    -0.00371996751152801
    -0.00368197568241012
    -0.00364360588016825
    -0.00360510516351791
    -0.00356591121574338
    -0.00352656997497068
    -0.00348685648439877
    -0.00344722643493981
    -0.00340736140579608
    -0.00336758426155901
    -0.00332753344802601
    -0.00328750802064282
    -0.0032471413328337
    -0.00320679774729835
    -0.00316610295239397
    -0.00312549340035217
    -0.00308454857127824
    -0.00304386170207064
    -0.00300290150021605
    -0.00296233875603428
    -0.00292129482201999
    -0.00288073484890323
    -0.00283951499715645
    -0.00279927582669419
    -0.00275770700056424
    -0.00271826154742389
    -0.00267536483006731
    -0.00263509464305156
    -0.00259725994910034
    -0.00255409536944302
    -0.00251352618627335
    -0.00247248705230046
    -0.00243353302728908
    -0.00239338505355968
    -0.002354108014004
    -0.00231354279813875
    -0.00227395565365889
    -0.00223352552719682
    -0.00219430023605191
    -0.00215442956472851
    -0.00211583028869268
    -0.0020765587432622
    -0.00203849112299403
    -0.00199959655075123
    -0.00196187823931274
    -0.00192320422309985
    -0.00188580324032983
    -0.00184738932024739
    -0.00181053087910452
    -0.00177256853414649
    -0.00173643201108726
    -0.00169886271207229
    -0.00166360336918217
    -0.00162645093636032
    -0.0015924435019478
    -0.00155509077331504
    -0.00152337012050789
    -0.00147988572401479
    -0.00145905623888122
    -0.00142229273289703
    -0.00138005083638085
    -0.00134807599742105
    -0.00131494126741691
    -0.00128545618449105
    -0.00125354739873849
    -0.00122322204449099
    -0.00119092382531205
    -0.00116004251486604
    -0.00112792836635087
    -0.00109736040348402
    -0.00106615479623852
    -0.00103659671861193
    -0.00100668160625124
  -0.000978288661010672
  -0.000949570338620629
    -0.00092216902638079
    -0.00089438708671706
  -0.000867737622092878
  -0.000840693466317883
  -0.000814687362552657
  -0.000788276159757994
  -0.000762827770843517
  -0.000737050063062293
  -0.000712345353053305
  -0.000687428482340838
  -0.000663561635529363
  -0.000639451830525362
  -0.000616537270116216
  -0.000593529464503111
    -0.007719806995443
% low pass filter 
aff_lp=filtfilt(numerator_lp,denominator,data);
% aff = filtfilt(numerator_lp,denominator,data);
filtdata = filtfilt(numerator_hp,denominator,aff_lp);