【正文】 N21。stem(n,hn2,39。.39。)。 %繪制Hamming窗單位脈沖響應axis([0,N21,])。 %確定顯示范圍xlabel(39。n39。)。ylabel(39。h(n)39。)。grid on。title(39。Hamming窗單位脈沖響應h(n)39。)。subplot(3,1,2)。plot(w2*2000/(2*pi),20*log10(abs(h2)))。 %繪制Hamming窗幅頻響應axis([0,1000,150,5])。xlabel(39。歸一化角頻率39。)。ylabel(39。幅度（單位：分貝）39。)。grid on。title(39。Hamming窗幅頻響應39。)。subplot(3,1,3)。plot(w2/pi,180/pi*unwrap(angle(h2)))。 %繪制Hamming窗相頻響應xlabel(39。歸一化角頻率39。)。ylabel(39。單位：度39。)。grid on。title(39。Hamming窗相頻相應39。)。grid。y=filter(hn2,1,x)。figure(6)。subplot(2,1,1)。plot(y)。xlabel(39。t39。)。ylabel(39。x(t)39。)。title(39。連續(xù)信號39。)。grid。Y=fft(y,512)。w=(0:255)/256*1000。subplot(2,1,2)。plot(w,abs([Y(1:256)]))。 %繪制頻譜圖xlabel(39。Hz39。)。ylabel(39。頻率響應幅度39。)。title(39。頻譜圖39。)。grid。%巴特沃斯帶通濾波器fp=[200 300]。fs=[150 350]。ap=2。as=30。fsa=2000。wp=fp/fsa*2。ws=fs/fsa*2。[n,wn]=buttord(wp,ws,ap,as)。[B,A]=butter(n,wn)。[H,w]=freqz(B,A,512)。figure(7)。subplot(2,1,1)。plot(w*2000/(2*pi),abs(H))。 %繪制帶通頻譜圖xlabel(39。Hz39。)。ylabel(39。頻率響應幅度39。)。title(39。帶通濾波器39。)。grid。subplot(2,1,2)。plot(w/pi,angle(H))。xlabel(39。Hz39。)。ylabel(39。angel39。)。title(39。相位特性39。)。grid。y=filter(B,A,x)。figure(8)。subplot(2,1,1)。plot(y)。xlabel(39。t39。)。ylabel(39。x(t)39。)。title(39。連續(xù)信號39。)。grid。Y=fft(y,512)。w=(0:255)/256*1000。subplot(2,1,2)。plot(w,abs([Y(1:256)]))。 %繪制頻譜圖xlabel(39。Hz39。)。ylabel(39。頻率響應幅度39。)。title(39。頻譜圖39。)。grid。%切比雪夫帶通濾波器fsl=150。fsu=350。 %第一個阻帶截止頻率fs1=450。fs2=650。 %第二個阻帶截止頻率fpl=200。fpu=300。 %第一個通帶截止頻率fp1=500。fp2=600。 %第二個通帶截止頻率fs=2000。 %抽樣頻率ws=[2*fsl/fs,2*fsu/fs]。 %數(shù)字濾波器的歸一化阻帶截止頻率wp=[2*fpl/fs,2*fpu/fs]。 %數(shù)字濾波器的歸一化通帶截止頻率ws1=[2*fs1/fs,2*fs2/fs]。wp1=[2*fp1/fs,2*fp2/fs]。rs=20。rp=。 %衰減設置[N,wpo]=cheb1ord(wp,ws,rp,rs)。 %計算階數(shù)和通帶截止頻率[M,wp1]=cheb1ord(wp1,ws1,rp,rs)。[B,A]=cheby1(N,rp,wpo)。 %計算第一個通帶分子分母的多項式系數(shù)[C,D]=cheby1(M,rp,wp1)。 %計算第二個通帶分子分母的多項式系數(shù)[h1,w1]=freqz(B,A)。 %求解第一個通帶數(shù)字濾波器的頻率響應[h2,w2]=freqz(C,D)。 %求解第二個通帶數(shù)字濾波器的頻率響應H=h1+h2。W=w1+w2。 %實現(xiàn)兩個通帶濾波器 figure(9)plot(W*2000/(2*2*pi),abs(H))。grid on。xlabel(39。頻率/Hz39。)。ylabel(39。幅度39。)。title(39。切比雪夫1型帶通濾波器的幅頻特性39。)。set(gca,39。YTickmode39。,39。manual39。,39。YTick39。,[0 1])。figure(10)plot(W/(2*pi),angle(H))。grid on。xlabel(39。頻率/HZ39。)。ylabel(39。angle39。)。title(39。切比雪夫1型帶通濾波器的相頻特性39。)。y1=filter(B,A,x)。y2=filter(C,D,x)。y=y1+y2。figure(11)plot(t,y)。xlabel(39。時間(t)39。)。ylabel(39。濾波后的信號y39。)。title(39。濾波后的信號39。)。grid on。axis([0 2 2])。Y=fft(y,512)。 %對濾出來的信號進行頻譜分析w=(0:255)/256*1000。 %頻率橫坐標的設定figure(12)plot(w,abs([Y(1:256)]))。 %繪制頻譜圖xlabel(39。頻率/Hz39。)。ylabel(39。頻率響應幅度39。)。title(39。濾波后的頻譜圖39。)。grid。20