%---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