matlab血压处理程序

%---bp_display

% data-file
DATAFILE=‘ttt.dcm‘;         

% lowpass filters parameter
%
N = 256;
fs = 133;
n = 8;
%wn = 10/（fs/2）;
wn = 1.2/（fs/2）;
rp = 3;
rs = 60;
%

fid=fopen（DATAFILE‘r‘）;
data=fread（fid5001‘int16‘）;
for i=1:2500
data1（i）=data（2*i）;
data2（i）=data（2*i+1）;
end
[ba]=butter（nwn）;
%data10=filter（badata1）;
data10=filtfilt（badata1）;

data11 = filter（1 -1 data1）;
//data11 = filtfilt（badata11）;

plot（data11（100:2500））;
pause;

subpl