信號(hào)分析與濾波器設(shè)計(jì)——使用矩形窗設(shè)計(jì)-資料下載頁(yè)
2025-01-16 10:08本頁(yè)面
【正文】 h=hd.*w_box。[db,mag,pha,grd,w]=freqz_m(h,[1])。 subplot(2,2,1)。stem(n,hd)。 % 理想脈沖響應(yīng)xlabel(39。n39。)。ylabel(39。hd(n)39。)。 title(39。Ideal Impulse Response39。)。 subplot(2,2,2)。stem(n,w_box)。 % 矩形窗xlabel(39。n39。)。ylabel(39。w(n)39。)。 title(39。Boxcar Window39。)。 subplot(2,2,3)。stem(n,h)。 % 實(shí)際脈沖響應(yīng)xlabel(39。n39。)。ylabel(39。h(n)39。)。 title(39。Actual Impulse Response39。)。 subplot(2,2,4)。plot(w*fs/(2*pi),db)。 % 幅度響應(yīng)（dB） axis([0 40 50 0])。 xlabel(39。Frequency (Hz)39。)。ylabel(39。Decibels39。)。title(39。Magnitude Response in dB39。)。grid% 程序功能：設(shè)計(jì)高通濾波器并畫(huà)出其頻譜圖： fb=25。fc=22。fs=100。 % 設(shè)置濾波器截止頻率wb=2*pi*fb/fs。ws=2*pi*fc/fs。 wc=*(wb+ws)。tr_width=wbws。 % 過(guò)渡帶寬M=ceil(*pi/tr_width)。 hd=ideal_lp(pi,M)ideal_lp(wc,M)。 % 產(chǎn)生一理想帶通濾波器 w_box=(boxcar(M))39。 % 矩形窗h=hd.*w_box。[db,mag,pha,grd,w]=freqz_m(h,[1])。 subplot(2,2,1)。stem(n,hd)。 % 理想脈沖響應(yīng)xlabel(39。n39。)。ylabel(39。hd(n)39。)。 title(39。Ideal Impulse Response39。)。 subplot(2,2,2)。stem(n,w_box)。 % 矩形窗xlabel(39。n39。)。ylabel(39。w(n)39。)。 title(39。Boxcar Window39。)。 subplot(2,2,3)。stem(n,h)。 % 實(shí)際脈沖響應(yīng)xlabel(39。n39。)。ylabel(39。h(n)39。)。 title(39。Actual Impulse Response39。)。 subplot(2,2,4)。plot(w*fs/(2*pi),db)。 % 幅度響應(yīng)（dB） axis([0 40 50 0])。 grid。 xlabel(39。Frequency (Hz)39。)。ylabel(39。Decibels39。)。title(39。Magnitude Response in dB39。)。% 程序功能：畫(huà)出濾波后信號(hào)的時(shí)域波形和頻譜 sf=filter(h,[1],s)。 % sf為濾濾波后的信號(hào)plot(t,sf)。 % 畫(huà)出濾波后信號(hào)的時(shí)域波形 xlabel(39。Time(seconds)39。)。ylabel(39。Time waveform39。)。axis([0 1 1 1])。 title(39。濾波后信號(hào)的時(shí)域波形39。)。 SF=fft(sf,512)。 % 對(duì)sf進(jìn)行快速傅里葉變換 w=(0:255)/256*(Fs/2)。plot(w,abs(SF(1:256)))。 % 畫(huà)出濾波后信號(hào)的幅度圖 xlabel(39。Frequency (Hz)39。)。 ylabel(39。幅度譜39。)。 title(39。濾波后信號(hào)的幅度譜39。)。 grid。 plot(w,angle(SF(1:256)))。 % 畫(huà)出濾波后信號(hào)的相位圖 xlabel(39。Frequency (Hz)39。)。ylabel(39。相位譜39。)。title(39。濾波后信號(hào)的相位譜39。)。% 程序功能：對(duì)濾波前后信號(hào)進(jìn)行比較 plot(w,abs([S(1:256)39。 SF(1:256)39。]))。 % 將濾波前后信號(hào)的幅度譜畫(huà)在一起xlabel(39。Frequency (Hz)39。)。ylabel(39。 Fourier transform39。)。grid。legend({39。before39。,39。after39。})。 %對(duì)兩個(gè)曲線進(jìn)行區(qū)分命名 title(39。濾波前后信號(hào)對(duì)比39。)。
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1