心电信号基于matlab ECG信号滤波分析（IRR+巴特沃斯）【含Matlab源码 3104期】.zip

close all clear all clc load coursework2ECG.mat s=signalNoise; S=fft(s); f=linspace(-fs/2,fs/2,length(s)); figure;plot(f,abs(fftshift(S)));grid on; %FIR滤波器： wp=0.1*pi;ws=0.15*pi; %滤波器边界频率 wdelta=ws-wp; %过渡带宽 N=ceil(8*pi/wdelta); %根据过渡带宽等于表 7-2中Hanning窗函数主瓣宽求得滤波器所用窗函数的最小长度 Nw=N; wc=(wp+ws)/2; %截止频率在通带和阻带边界频率的中点 n=0:N-1; alpha=(N-1)/2; %求滤波器的相位延迟 m=n-alpha+eps; %eps为MATLAB系统的精度 hd=sin(wc*m)./(pi*m); %式求理想滤波器脉冲响应 win=hanning(Nw); %采用汉宁窗 h=hd.*win'; b=h; figure,impz(b,1); [H,f]=freqz(b,1,512,50); %采用 50 Hz 的采样频率绘出该滤波器的幅频和相频响应 y=fftfilt(b,s); %使用滤波器对信号进行滤波； figure;plot(s);hold on; plot(y); grid on; legend('滤波前','滤波后') % IIR滤波器设计 N=0; % 阶数 Fp=100; % 通带截止频率50Hz Fc=200; % 阻带截止频率100Hz Rp=1; % 通带波纹最大衰减为1dB Rs=60; % 阻带衰减为60dB % 3.0 计算最小滤波器阶数 na=sqrt(10^(0.1*Rp)-1); ea=sqrt(10^(0.1*Rs)-1); N=ceil(log10(ea/na)/log10(Fc/Fp)); % 3.1 巴特沃斯滤波器 Wn=Fp*2/fs; [Bb Ba]=butter(N,Wn,'low'); % 调用MATLAB butter函数快速设计滤波器 [BH,BW]=freqz(Bb,Ba); % 绘制频率响应曲线 y2=filter(Bb,Ba,s); % 进行低通滤波 figure;plot(s);hold on; plot(y2); grid on; legend('滤波前','滤波后') figure; subplot(211);plot(abs(fftshift(fft(y)))); grid on; subplot(212);plot(abs(fftshift(fft(y2)))); grid on;