MATLAB处理心电信号程序

对于心电信号处理不错的资源,本次实验是要我们先产生心电信号，再加噪声再滤波。由一个心电信号数据表输出心电信号，由于已知的心电信号有噪声，所以我先将信号进行滤波，得到正确的心电信号。再加噪声，通过滤波可以得到高频和中频噪声，将原信号加上噪声显示再滤波。每次的时域信号都画出了他们相应的频谱，便于观察。每次的滤波我也尝试了不同的滤波器。

%程序中s为原信号，s1为滤过噪声之后的心电信号
%11107110任惠
clear;
load data0.txt;
s=data0;
%时域波形
t=0:4/3255:4;       %设心率为每分钟75次，则QRS波群周期为0.8秒，5个周期为4秒

subplot（431）
plot（ts）;          %带噪声的心电信号
xlabel（‘时间 t/s‘）;
ylabel（‘电压 V/mV‘）;
title（‘含噪声的ECG信号‘）;

subplot（434）
x=fft（s1024）;
Pxx= x.* conj（x） / 1024;
f = 3255/4*（0:216）/1024;
plot（fPxx（1:217））;       %信号频谱
title（‘fft画出的原信号频谱‘）
xlabel（‘频率‘）
ylabel（‘幅值 V‘）

subplot（432）
wp=tan（0.02*pi/2）;
ws=tan（0.03*pi/2）;
[NWn]=buttord（wpws0.520‘s‘）;
[numden]=butter（NWn‘s‘）;
[BA]=bilinear（numden1/2）;
s1=filter（BAs）;                   %巴特沃斯滤波得到正确信号
plot（ts1）;
xlabel（‘时间 t/s‘）;
ylabel（‘电压 V/mV‘）;
title（‘巴特沃斯滤波器去噪声的ECG信号‘）;

subplot（435）
w=0:pi/256:pi;
y=freqz（s11w）;                    %画出信号的频谱
plot（w/piabs（y））;
axis（[0 1 0 300]）;
xlabel（‘频率 w‘）;
ylabel（‘幅值 V‘）;
title（‘巴特沃斯去噪声的ECG信号频谱‘）;

subplot（433）
s2=3*（rand（13256）-1）;              %产生噪声
wp=tan（0.6*pi/2）;
ws=tan（0.4*pi/2）;
[NWn]=buttord（wpws0.120‘s‘）;
[numden]=butter（NWn‘high‘‘s‘）;
[BA]=bilinear（numden1/2）;
s3=filter（BAs2）;                  %得到高频噪声
s4=s1‘+s3;                          %加完噪声后的信号
plot（ts4）;
xlabel（‘时间 t/s‘）;
ylabel（‘电压 V/mV‘）;
title（‘加高频噪声的ECG信号‘）;

subplot（436）
y=freqz（s41w）;
plot（w/piabs（y））;
axis（[0 1 0 300]）;
xlabel（‘频率 w‘）;
ylabel（‘幅值 V‘）