【正文】 Laurence Keselbrener,Michel high pass filter for Rwave detection in ECG signal,Medical Engineering and ,19(5):481~484. [11] 劉欣桃 ,呂軍 .心電信號的計算機輔助分析方法的研究 [J].計算機輔助設計與圖形學學報 ,1997,14(2):18~20. 心電信號小波消噪研究與應用 27 致 謝 時光飛逝，四年大學生活短暫的四年珍貴和值得住在我的生命中錯過的好日子。在此期間，我收獲了很多。我的主題還不夠，我想對那些關心和鼓勵你給了我在這里，表示最衷心的感謝。 首先，我要衷心感謝我的導師老師 。他深刻的科學知識、 實踐經驗、 嚴謹?shù)膶W習，熱衷于研究氣味、 現(xiàn)實主義的風格，使我一生受益。他賦予我不僅學習與工作方法、 他的高尚人格 和風范，學科的指揮官是比其他人更寬是有益的值得我尊重和學習。 論文完成之際，謹對我尊敬的導師 老師表示衷心的感謝和崇高的敬意！ 其次，由于鄭州市研究所科學與技術 4班 ，每個人都是鼓勵互相幫助，給了我很大的幫助在生活和學習中。這是他們無私地幫助我解決的一些研究中，所遇到的問題使這篇文章的人得以順利完成 !我也從他們收獲了同學之間的深厚友誼。 最后，感謝我的家人和朋友給與我的幫助以及對我的理解，特別感謝父母對我的無私付出和默默的支持，他們給了我不斷前進的動力。祝福他們永遠安康！ 特別感謝為評閱本文付出辛勤勞動的專家和 老師們！ 心電信號小波消噪研究與應用 28 附錄 編程與程序 %核心程序如下面的所示，使用的小波基有 coif5， db3， sym4， ， haar，rbio load x_100,x=x_10039。 %導入自己合成的無噪的心電信號 subplot(811) plot(x) %在圖中畫出無噪聲心電信號波形 y=awgn(x,5,39。measured39。,39。db39。) %加入高斯白噪聲，信噪比為 5db。 [c,l]=wavedec(y,5,39。db339。) %采用 coif1小波基對含噪信號進行小波分解，分解層次為 5層。 subplot(812),plot(y) %在圖中畫出含噪聲的心電信號波形 a1=appcoef(c,l,39。db339。,1) %提取第 1層的低頻系數(shù) a2=appcoef(c,l,39。db339。,2) %提取第 2層的低頻系數(shù) a3=appcoef(c,l,39。db339。,3) %提取第 3層的低頻系數(shù) a4=appcoef(c,l,39。db339。,4) %提取第 4層的低頻系數(shù) a5=appcoef(c,l,39。db339。,5) %提取第 5層的低頻系數(shù) subplot(813) %在同一坐標圖中畫出各層低頻系數(shù) plot(a1) subplot(814) plot(a2) 心電信號小波消噪研究與應用 29 subplot(815) plot(a3) subplot(816) plot(a4) subplot(817) plot(a5) d1=detcoef(c,l,1) %提取從第 1 層到第 5 層的高頻系數(shù)，并且在圖中畫出來 d2=detcoef(c,l,2) d3=detcoef(c,l,3) d4=detcoef(c,l,4) d5=detcoef(c,l,5) subplot(811) plot(x) subplot(812) plot(y) subplot(813) plot(d1) subplot(814) plot(d2) subplot(815) 心電信號小波消噪研究與應用 30 plot(d3) subplot(816) plot(d4) subplot(817) plot(d5) %采用的是 db3小波基去噪。 load synecg,x=synecg39。 subplot(611) plot(x) title([39。原始信號 39。]) s1=awgn(x,5,39。measured39。,39。db39。) subplot(612) plot(s1) title([39。含噪信號 39。]) lev=5 subplot(613) title([39。用 heursure 去噪 39。]) xd1=wden(s1,39。heursure39。,39。s39。,39。one39。,lev,39。db339。) %默認方法去噪得到xd1 plot(xd1) %畫出去噪后的圖形 ps=sum(x.^2) pn11=xxd1 心電信號小波消噪研究與應用 31 pn=sum(pn11.^2) c1=corrcoef(xd1,x) %計算互相關函數(shù)，一個重要的對去噪效果的評判標準，計算去噪后的信噪比值，閾值不做調整 snr1=10*log10(ps/pn) xd11=wden(s1,39。heursure39。,39。s39。,39。sln39。,lev,39。db339。) plot(xd11) ps=sum(x.^2) pn11=xxd11 pn=sum(pn11.^2) c11=corrcoef(xd11,x) snr11=10*log10(ps/pn) %計算 snr值，閾值根據(jù)第一層系數(shù)進行一次噪聲調整的情況 xd111=wden(s1,39。heursure39。,39。s39。,39。mln39。,lev,39。db339。) plot(xd111) ps=sum(x.^2) pn111=xxd111 pn=sum(pn111.^2) c111=corrcoef(xd11,x) snr111=10*log10(ps/pn) %計算 snr 值，閾值在不同層次估計噪聲 %以上的全采用的是軟閾值函數(shù)（ s），現(xiàn)在采用硬閾值函數(shù)來處理數(shù)據(jù) 心電信號小波消噪研究與應用 32 title([39。用 heursure 去噪 39。]) xd1=wden(s1,39。heursure39。,39。h39。,39。one39。,lev,39。db339。) plot(xd1) ps=sum(x.^2) pn11=xxd1 pn=sum(pn11.^2) c1=corrcoef(xd1,x) snr1=10*log10(ps/pn) %計算去噪后的信噪比值，閾值不做調整 xd11=wden(s1,39。heursure39。,39。h39。,39。sln39。,lev,39。db339。) plot(xd11) ps=sum(x.^2) pn11=xxd11 pn=sum(pn11.^2) c11=corrcoef(xd11,x) snr11=10*log10(ps/pn) %計算 snr 值，閾值根據(jù)第一層系數(shù)進行一次噪聲調整的情況 xd111=wden(s1,39。heursure39。,39。h39。,39。mln39。,lev,39。db339。) plot(xd111) ps=sum(x.^2) pn111=xxd111 pn=sum(pn111.^2) 心電信號小波消噪研究與應用 33 c111=corrcoef(xd11,x) snr111=10*log10(ps/pn) %計算 snr值，閾值在不同層次估計噪聲 %經過計算可以得到，采用硬閾值函數(shù)來處理的結果沒有用軟閾值的好，互相關也沒軟閾值大。 %所以還有剩下的 5個小波基要分析 %以上的為使用 heursure去噪的軟硬法，以 db3小波為例子的，現(xiàn)在是第二種方法，用 rigrsure去噪 title([39。用 rigrsure 去噪 39。]) xd1=wden(s1,39。rigrsure39。,39。s39。,39。one39。,lev,39。db339。) plot(xd1) ps=sum(x.^2) pn11=xxd1 心電信號小波消噪研究與應用 34 pn=sum(pn11.^2) c1=corrcoef(xd1,x) snr1=10*log10(ps/pn) %計算去噪后的信噪比值，閾值不做調整 xd11=wden(s1,39。rigrsure39。,39。s39。,39。sln39。,lev,39。db339。) plot(xd11) ps=sum(x.^2) pn11=xxd11 pn=sum(pn11.^2) c11=corrcoef(xd11,x) snr11=10*log10(ps/pn) %計算 snr值，閾值根據(jù)第一層系數(shù)進行一次噪聲調整的情況 xd111=wden(s1,39。rigrsure39。,39。s39。,39。mln39。,lev,39。db339。) plot(xd111) ps=sum(x.^2) pn111=xxd111 pn=sum(pn111.^2) c111=corrcoef(xd11,x) snr111=10*log10(ps/pn) %計算snr值，閾值在不同層次估計噪聲 %以上的全采用的是軟閾值函數(shù)，現(xiàn)在采用硬閾值函數(shù)來處理數(shù)據(jù) 心電信號小波消噪研究與應用 35 title([39。用 rigrsure 去噪 39。]) xd1=wden(s1,39。rigrsure39。,39。h39。,39。one39。,lev,39。db339。) plot(xd1) ps=sum(x.^2) pn11=xxd1 pn=sum(pn11.^2) c1=corrcoef(xd1,x) snr1=10*log10(ps/pn) %計算去噪后的信噪比值，閾值不做調整 xd11=wden(s1,39。rigrsure39。,39。h39。,39。sln39。,lev,39。db339。) plot(xd11) ps=sum(x.^2) pn11=xxd11 pn=sum(pn11.^2) c11=corrcoef(xd11,x) snr11=10*log10(ps/pn) %計算 snr值，閾值根據(jù)第一層系數(shù)進行一次噪聲調整的情況 xd111=wden(s1,39。rigrsure39。,39。h39。,39。mln39。,lev,39。db339。) plot(xd111) ps=sum(x.^2) pn111=xxd111 pn=sum(pn111.^2) c111=corrcoef(xd11,x) 心電信號小波消噪研究與應用 36 snr111=10*log10(ps/pn) %計算snr值，閾值在不同層次估計噪聲 %以上的為使用 rigrsure去噪的軟硬法，以 db3小波為例子的，現(xiàn)在是第 3種方法，用 sqtwolog去噪 title([39。用 sqtwolog 去噪 39。]) xd1=wden(s1,39。sqtwolog39。,39。s39。,39。one39。,lev,39。db339。) plot(xd1) ps=sum(x.^2) pn11=xxd1 pn=sum(pn11.^2) c1=corrcoef(xd1,x) snr1=10*log10(ps/