基于matlab的心電信號分析系統(tǒng)的設計與仿真docxdocx(編輯修改稿)
2025-08-14 17:19 本頁面
【文章內(nèi)容簡介】 0 200])：wp=60*2*pi。ws=99*2*pi。Rp=1。As=40。[N,wc]=buttord(wp,ws,Rp,As,39。s39。)[B,A]=butter(N,wc,39。s39。)k=0:511。fk=0:1000/512:1000。wk=2*pi*fk。Hk=freqs(B,A,wk)。plot(fk,20*log10(abs(Hk)))。 grid onN = 11wc = B = +028 * Columns 1 through 10 0 0 0 0 0 0 0 0 0 0 0 Columns 11 through 12 A = +028 * Columns 1 through 10 Columns 11 through 12 ：wp=*2*pi。ws=*2*pi。Rp=。As=40。[N,wc]=buttord(wp,ws,Rp,As,39。s39。)[B0,A0]=butter(N,wc,39。s39。)。wph=2*pi*。hk=freqs(B0,A0,wph)。[BH,AH]=lp2hp(B,A,wph)。[h,w]=freqs(BH,AH)。plot(w,20*log10(abs(h)))。axis([0,1,80,5])。N = 7wc = B0 = +003 * 0 0 0 0 0 0 0 A0 = +003 * ：t=0::。F=interp1(x,y,t)。F=F39。t=t39。figure(1)subplot(3,1,1)。plot(1000*t,F)。wp=*2*pi。ws=*2*pi。Rp=。As=40。T=1。[N,wc]=buttord(wp,ws,Rp,As,39。s39。)[B,A]=butter(N,wc,39。s39。)。[b,a]=imp_invr(B,A,T)[db,mag,pha,w]=freqz_m(b,a)。y1=filter(b,a,F)。subplot(3,1,2)。plot(y1)。title(39。高通濾波后39。)wp1=2*pi*60。ws1=2*pi*99。Rp1=。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)[d
電大資料相關推薦
鄂ICP備17016276號-1