資源簡介
心電信號QRS復波實時檢測,軟件主要是檢測出R波波峰位置,采用雙閾值法進行檢測。配有GUI演示界面。所用數據為MIT,1986年心電數據庫。
此為醫學信號處理實驗課程作業,與大家分享

代碼片段和文件信息
clear?all;
close?all;
clc;
PATH=?‘G:\matlab?workspace\數字信號實驗\ECG?database‘;?%?path?where?data?are?saved
HEADERFILE=?‘121.hea‘;??????%?header-file?in?text?format
ATRFILE=?‘121.atr‘;?????????%?attributes-file?in?binary?format
DATAFILE=‘121.dat‘;?????????%?data-file
SAMPLES2READ=3000;?????????%?number?of?samples?to?be?read
????????????????????????????%?in?case?of?more?than?one?signal:
????????????????????????????%?2*SAMPLES2READ?samples?are?read
[TIME?ECGDATA?sfreq]?=?rdECGdata(?PATHHEADERFILEATRFILEDATAFILESAMPLES2READ);%讀取MIT-ECG數據
%------------------------------二階低通錄波器11HZ---------------------------%
DATA?=?ECGDATA(:1);
num?=?zeros(113);
den?=?zeros(113);
num(1)?=?1;num(7)?=?-2;num(13)?=?1;
den(1)?=?1;den(2)?=?-2;den(3)?=?1;
%?w?=?[0:1:512]*pi/512;
%?H?=?freqz(numdenw);
%?figureplot(w/piabs(H));title(‘低通濾波器頻響特性‘);
DATAafterLow?=?filter(numdenDATA)/36;
clear?num;
clear?den;
%-------------------------------------------------------------------------%
%--------------------------高通濾波器,15HZ--------------------------------%
num?=?zeros(133);
den?=?zeros(133);
num(1)=-1;num(17)=32;num(18)=-32;?num(33)?=?1;
den(1)=32;?den(2)?=?-32;
%?w?=?[0:1:512]*pi/512;
%?H?=?freqz(numdenw);
%?figureplot(w/piabs(H));title(‘高通濾波器頻響特性‘);
DATAafterHigh?=?filter(numdenDATAafterLow);
figuresubplot(211)plot(TIMEDATA);title(‘原始信號‘);
subplot(212)plot(TIMEDATAafterHigh);title(‘帶通濾波后信號‘);
clear?num;
clear?den;
%-------------------------------------------------------------------------%
%----------------------------微分-----------------------------------------%
num?=?zeros(15);
den?=?zeros(15);
num(1)=0.2;num(2)=0.1;num(4)=-0.1;?num(5)?=?-0.2;
den(1)=1;
DATAafterDif?=?filter(numdenDATAafterHigh)*10;
figuresubplot(211)plot(TIMEDATAafterHigh);title(‘濾波后信號‘);
subplot(212)plot(TIMEDATAafterDif);title(‘微分后信號‘);
%-------------------------------------------------------------------------%
%-----------------------------平方----------------------------------------%
DATAafterSqu?=?DATAafterDif.*DATAafterDif;
%?fftDATAafterDif?=?fftshift(fft(DATAafterDif));
%?fftDATAafterSqu?=?fftshift(fft(DATAafterSqu));
fftDATAafterDif?=?(fft(DATAafterDif));
fftDATAafterSqu?=?(fft(DATAafterSqu));
figureplot(TIME(1:1000)DATAafterSqu(1:1000));
n?=?(0:1500)/1500*sfreq;
figuresubplot(211)plot(nabs(fftDATAafterDif(1:1501)));title(‘平方前信號頻譜‘);
subplot(212)plot(nabs(fftDATAafterSqu(1:1501)));title(‘平方后信號頻譜‘);
%-------------------------------------------------------------------------%
%--------------------------------移動窗口積分,150ms-----------------------%
%?QRSTime?=?100;
QRSTime?=?150;
%?QRSTime?=?200;
%?QRSTime?=?250;
%?QRSTime?=?300;
nInte?=?floor(QRSTime/1000*sfreq);%150ms積分的離散值點數
DATAmoveWinInte?=?zeros(1length(DATAafterSqu));
for?ii=nInte:length(DATAafterSqu)
????DATAmoveWinInte(ii)?=?sum(DATAafterSqu(ii-nInte+1:ii))/nInte;
end
figureplot(TI
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????3774??1992-07-30?13:44??R波波峰檢測\ECG?databa
?????文件????1950000??1992-07-30?21:44??R波波峰檢測\ECG?databa
?????文件????????150??1992-07-30?13:44??R波波峰檢測\ECG?databa
?????文件???????4970??1992-07-30?13:44??R波波峰檢測\ECG?databa
?????文件????1950000??1992-07-30?21:44??R波波峰檢測\ECG?databa
?????文件????????221??1992-07-30?13:44??R波波峰檢測\ECG?databa
?????文件???????3046??1992-07-30?13:44??R波波峰檢測\ECG?databa
?????文件????1950000??1992-07-30?21:44??R波波峰檢測\ECG?databa
?????文件????????183??1992-07-30?13:44??R波波峰檢測\ECG?databa
?????文件???????3360??1992-07-30?13:45??R波波峰檢測\ECG?databa
?????文件????1950000??1992-07-30?21:45??R波波峰檢測\ECG?databa
?????文件????????219??1992-07-30?13:45??R波波峰檢測\ECG?databa
?????文件???????8114??1992-07-30?13:45??R波波峰檢測\ECG?databa
?????文件????1950000??1992-07-30?13:45??R波波峰檢測\ECG?databa
?????文件????????306??1992-07-30?13:45??R波波峰檢測\ECG?databa
?????文件???????4322??1992-07-30?13:46??R波波峰檢測\ECG?databa
?????文件????1950000??1992-07-30?13:46??R波波峰檢測\ECG?databa
?????文件????????284??1992-07-30?13:46??R波波峰檢測\ECG?databa
?????文件???????4388??1992-07-30?13:46??R波波峰檢測\ECG?databa
?????文件????1950000??1992-07-30?13:46??R波波峰檢測\ECG?databa
?????文件????????277??1992-07-30?13:46??R波波峰檢測\ECG?databa
?????文件???????7092??1998-10-14?04:09??R波波峰檢測\ECG?databa
?????文件????1950000??1992-07-30?13:46??R波波峰檢測\ECG?databa
?????文件????????423??1992-07-30?13:46??R波波峰檢測\ECG?databa
?????文件???????5436??1992-07-30?13:47??R波波峰檢測\ECG?databa
?????文件????1950000??1992-07-30?13:47??R波波峰檢測\ECG?databa
?????文件????????222??1992-07-30?13:47??R波波峰檢測\ECG?databa
?????文件???????4958??1992-07-30?13:47??R波波峰檢測\ECG?databa
?????文件????1950000??1992-07-30?13:47??R波波峰檢測\ECG?databa
?????文件????????546??1992-07-30?13:47??R波波峰檢測\ECG?databa
............此處省略15個文件信息
- 上一篇:計算機圖形學
- 下一篇:王汝傳版編譯原理資料
評論
共有 條評論