河北工业大学信息工程学院 信号与线性系统课程设计 subplot(1,2,2) plot(f,abs(Y))
axis([0,100,0,80]) title('滤波后心电信号频谱') figure(3)
subplot(2,2,1) plot(t1,Xn1) title('滤波前信号') subplot(2,2,2) plot(t1,y)
title('滤波后信号') subplot(2,2,3) plot(t1,Xn1)
axis([0,1.5,-1.5,1.5]) title('滤波前截取一部分信号') subplot(2,2,4) plot(t1,y)
axis([0,1.5,-1.5,1.5]) title('滤波后截取一部分信号') 所出图形如下
9
河北工业大学信息工程学院 信号与线性系统课程设计 结果分析：
10
河北工业大学信息工程学院 信号与线性系统课程设计
（3）、将信号通过低通、高通组合成的带通滤波器程序
[t,Xn]=duquexinhao1('117.txt'); baocun1(t,Xn) %保存读取信号
[t1,Xn1]=xianxingchazhi(t,Xn); baocun2(t1,Xn1)%保存插值后信号 figure(1)
keshehuatu(t,Xn,t1,Xn1) %画原始信号和插值后信号波形和频谱 xy=[t1,Xn1];
wp=0.52*2*pi; %低通滤波器滤波 ws=0.62*2*pi; p=1; s=35;
[H1,f]=kesheditonglvboqi(wp,ws,p,s,Xn1); wp=0.10*2*pi; %高通滤波器滤波 ws=0.25*2*pi; p=1; s=35;
[H2,f]=keshegaotonglvboqi(wp,ws,p,s,Xn1); H=abs(H1).*abs(H2); %低通和高通组合的带通 Y=H'.*abs(fft(Xn1)); %经过滤波后心电信号频谱 y=ifft(Y); %滤波后心电信号时域波形 figure(2)
subplot(2,2,1) plot(f,abs(H1)) axis([0,1,0,1.5]) title('低通滤波器') subplot(2,2,2) plot(f,abs(H2)) axis([0,1,0,1.5]) title('高通滤波器') subplot(2,2,3) plot(f,abs(H)) axis([0,1,0,1.5]) title('组合 带通滤波器') subplot(2,2,4) plot(f,abs(Y)) axis([0,1,0,260])
title('滤波后心电信号频谱') figure(3)
subplot(2,1,1) plot(t1,Xn1) title('滤波前信号') subplot(2,1,2)
11
河北工业大学信息工程学院 信号与线性系统课程设计 plot(t1,y)
title('滤波后信号') 所出图形如下
12
