信號分析與濾波器設(shè)計——使用矩形窗設(shè)計-資料下載頁
2025-01-16 10:08本頁面
【正文】 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)。 % 實際脈沖響應(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è)計高通濾波器并畫出其頻譜圖： fb=25。fc=22。fs=100。 % 設(shè)置濾波器截止頻率wb=2*pi*fb/fs。ws=2*pi*fc/fs。 wc=*(wb+ws)。tr_width=wbws。 % 過渡帶寬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)。 % 實際脈沖響應(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。)。% 程序功能：畫出濾波后信號的時域波形和頻譜 sf=filter(h,[1],s)。 % sf為濾濾波后的信號plot(t,sf)。 % 畫出濾波后信號的時域波形 xlabel(39。Time(seconds)39。)。ylabel(39。Time waveform39。)。axis([0 1 1 1])。 title(39。濾波后信號的時域波形39。)。 SF=fft(sf,512)。 % 對sf進(jìn)行快速傅里葉變換 w=(0:255)/256*(Fs/2)。plot(w,abs(SF(1:256)))。 % 畫出濾波后信號的幅度圖 xlabel(39。Frequency (Hz)39。)。 ylabel(39。幅度譜39。)。 title(39。濾波后信號的幅度譜39。)。 grid。 plot(w,angle(SF(1:256)))。 % 畫出濾波后信號的相位圖 xlabel(39。Frequency (Hz)39。)。ylabel(39。相位譜39。)。title(39。濾波后信號的相位譜39。)。% 程序功能：對濾波前后信號進(jìn)行比較 plot(w,abs([S(1:256)39。 SF(1:256)39。]))。 % 將濾波前后信號的幅度譜畫在一起xlabel(39。Frequency (Hz)39。)。ylabel(39。 Fourier transform39。)。grid。legend({39。before39。,39。after39。})。 %對兩個曲線進(jìn)行區(qū)分命名 title(39。濾波前后信號對比39。)。
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1