二、生理信号处理 ——1.心电信号(含Matlab代码及数据)

 ## 心电图导入及读取

    
 clc;
    
 [TIME,M,Fs,siginfo]=rdmat('100m');# 通过读取函数ramat对心电图进行处理
    
 Fs=1500;# 采样频率
    
 plot(TIME,M);

 ecg = M(1:1000);

    
 TIME = TIME(1:1000);
    
 plot(TIME,ecg)

 ## Butterworth 陷波滤波器-去除工频干扰

    
 # 设计滤波器
    
 d = designfilt('bandstopiir','FilterOrder',2, ...
    
            'HalfPowerFrequency1',49,'HalfPowerFrequency2',51, ...
    
            'DesignMethod','butter','SampleRate',Fs);
    
 #应用滤波器去除50Hz干扰
    
 ecg_50 = filtfilt(d,ecg);

 ##bandpass滤波器-解决基线漂移

    
 fmaxd_1=5;# 截止频率为5Hz
    
 fmaxn_1=fmaxd_1/(Fs/2);
    
 [B,A]=butter(1,fmaxn_1,'low');
    
 ecg_low=filtfilt(B,A,ecg_50);# 通过5Hz低通滤波器的信号
    
 ecg1=ecg_50-ecg_low; # 减去5Hz低频信号
    
 plot(TIME,ecg1)
    
 xlabel('t(s)');
    
 ylabel('mv');
    
 title('期望信号');

 ##添加随机噪声信号

    
 Noise_White = (0.1*wgn(1,length(TIME),2))'; 
    
 plot(Noise_White)
    
 title('噪声信号');

 ##噪声干扰后信号

    
 Mix_Signal = ecg1 + Noise_White;
    
 plot(TIME,Mix_Signal)
    
 xlabel('t(s)');
    
 ylabel('mv');
    
 title('噪声干扰后信号');

 ##维纳滤波

    
 Signal_Filter = filter(h,1,Mix_Signal);# 将输入信号通过维纳滤波器
    
 figure(1)
    
 plot(TIME,Signal_Filter)
    
 xlabel('t(s)');
    
 ylabel('mv');
    
 title('维纳滤波后信号');