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