利用matlab实现对胎儿心电图中的噪声抑制.docVIP

利用Matlab实现对胎儿心电图中的噪声抑制 这是一个自适应滤波器在胎儿心电图中的应用，它把母亲的心跳信号从胎儿心跳传感器的信号中移除。 产生母亲的心跳信号 首先仿真母亲和她的胎儿的心电图信号，如下图所示。信号的心率大约为89次/min，信号的峰值为3.5mV。 %产生母亲心跳信号:以4000Hz为抽样率的母亲心电信号，信号心率大约为89次/min，信号的峰值电压为3.5mV x1=3.5*ecg(2700).; y1=sgolayfilt(kron(ones(1,13),x1),0,21);%Demonstrates Savitzky-Go-lay filtering n=1:30000; del=round(2700*rand(1)); mhb=y1(n+del); t=0.00025:0.00025:7.5; plot(t,mhb); axis([0 2.5 -4 4]);grid; xlabel(Time[sec]);ylabel(Voltage[mV]); title(Maternal Heartbeat Signal); 产生胎儿心跳信号 胎儿的心跳显然比母亲的心跳快得多，心率范围在每分钟120次到160次之间。胎儿心电图的振幅也明显弱于母亲的心电图，如下图所示。下面的程序将产生一个心率为每分钟139次和峰值电压为0.25mV的胎儿心电图信号。 %产生胎儿心跳信号：心率为每分钟139次和峰值电压为0.25mV的胎儿心电信号 x2=0.25*ecg(1725); y2=sgolayfilt(kron(ones(1,20),x2),0,17);%Demonstrates Savitzky-Go-lay filtering del=round(1725*rand(1)); fhb=y2(n+del); plot(t,fhb); axis([0 2.5 -1 1]);grid; xlabel(Time[sec]);ylabel(Voltage[mV]); title(Fetal Heartbeat Signal); 胎儿心电图的测量 胎儿心电图信号是从放置在母亲腹部的传感器测量获得的，但是从胸腔传到腹部的母亲的心跳信号在测量信号中往往占主导地位。可以把传播路径看作是一个有10个随机系数的线性FIR滤波器。另外再加一个不相关的小幅高斯噪声来模拟测量中存在的宽带噪声。含噪声的胎儿心电图如下图所示。 %胎儿心电图的测量：由放置在母亲腹部的传感器测量获得 Wopt=[0 1.0 -0.5 -0.8 1.0 -0.1 0.2 -0.3 0.6 0.1]; d=filter(Wopt,1,mhb)+fhb+0.02*randn(size(mhb)); plot(t,d); axis([0 2.5 -4 4]);grid; xlabel(Time/s);ylabel(Voltage/mV); title(Measured Signal); 母亲心电图的测量 母亲的心电图信号是从母亲的胸部得到的。自适应噪声抑制器的目的是从胎儿的心电图信号中自适应地抑制母亲的心跳信号，抑制器需要一个来自母亲心电图的参考信号来完成这个任务。和胎儿的心电图信号一样，母亲心电图信号中叶包含了一些附加的宽带噪声，如下图所示。 %母亲心电图的测量：从母亲的胸部得到 x=mhb+0.02*randn(size(mhb)); plot(t,x); axis([0 2.5 -4 4]);grid; xlabel(Time/s);ylabel(Voltage/mV); title(Reference Signal); 应用自适应噪声抑制器 自适应噪声抑制器可以使用常用的自适应算法进行噪声抑制。为简单起见，我们调用有15个抽头（系数）的最小均方（LMS）自适应滤波器，算法步长取0.00007。经过上述参数的设置，自适应噪声抑制算法经过数秒后收敛，从而大大缩短了病人等待诊断书的时间。 %应用自适应噪声抑制器：用15抽头（系数）的最小均方（LMS）自适应滤波器，算法步长取0.00007 h=adaptfilt.lms(15,0.00007); [y,e]=filter(h,x,d); plot(t,d,b,t,e,r); xlabel(Time[sec]);ylabel(Voltage[mV]); title(Convergence of Adaptive Noise Canceller); legend(Measured Signal,Error Signal); 胎儿心跳信号的恢复 自适应滤波器的输出信号y(n)中包含有估计的母亲心跳信号，这并非是我们感兴趣的最终信号。算法收敛之后，留在误差信号e(n)中的是带有残留测量噪声的胎儿心跳信号的估计。 %胎儿心跳信号的恢复 plo