基于MATLAB的心电信号分析系统的设计与仿真(3)

信号与线性系统课程设计报告

mag=abs(H); ph=angle(H); ph=ph*180/pi;

subplot(2,2,3),plot(f,mag);grid title('h2幅度响应')

subplot(2,2,4);plot(f,ph);grid title('h2相位响应')

figure(6)

zr=roots(B) %系统一的零点图 pk=roots(A) %系统一的极点图

zplane(B,A); %zplane函数画出系统一的零极点图 figure(7)

zr1=roots(cs1) %系统二的零点图 pk1=roots(ds1) %系统二的极点图

zplane(cs1,ds1); %zplane函数画出系统一的零极点图

4.1.9截取到2.5s对截取的部分进行滤波及频谱分析 t1=0:0.001:2.5; F0=interp1(x,y,t1); F0=F0'; t1=t1'; figure(8)

subplot(3,1,1); plot(1000*t1,F0);

wp=0.7*2*pi;ws=0.25*2*pi;Rp=0.1;As=40;T=1; [N,wc]=buttord(wp,ws,Rp,As,'s') [B,A]=butter(N,wc,'s'); [b,a]=imp_invr(B,A,T)

[db,mag,pha,w]=freqz_m(b,a); y11=filter(b,a,F0);

subplot(3,1,2);plot(y11); title('高通滤波后_{2.5}')

wp1=2*pi*60;ws1=2*pi*99;Rp1=0.1;As1=40;T1=1000; OmegaP1=wp1/T1;OmegaS1=ws1/T1;

[cs1,ds1]=afd_butt(OmegaP1,OmegaS1,Rp1,As1); [b1,a1]=imp_invr(cs1,ds1,T)

[db1,mag1,pha1,w1]=freqz_m(b1,a1); y21=filter(b1,a1,y11); subplot(3,1,3);plot(y21); title('低通滤波后_{2.5}') M=length(F0); m=1:M; fs=1000; f2=m*fs/M;

~ 11 ~

信号与线性系统课程设计报告

F01=fft(F0); Y11=fft(y11); Y21=fft(y21) figure(9)

subplot(3,1,1) plot(f2,abs(F01)) axis([0,1000,0,200])

title('原始信号频谱_{2.5}') subplot(3,1,2) plot(f2,abs(Y11)) axis([0,1000,0,200])

title('高通滤波后信号频谱_{2.5}') subplot(3,1,3) plot(f2,abs(Y21)) axis([0,1000,0,200])

title('低通滤波后信号频谱_{2.5}')

~ 12 ~

信号与线性系统课程设计报告

~ 13 ~

信号与线性系统课程设计报告

~ 14 ~

信号与线性系统课程设计报告

~ 15 ~