通濾波器設計(含matlab程序)畢業(yè)論文-資料下載頁
2025-06-24 01:20本頁面
【正文】 。 % 過渡帶寬M=ceil(*pi/tr_width)。hd=ideal_lp(wc,M)。 % 產(chǎn)生一理想低通濾波器w_box=(boxcar(M))39。 % 矩形窗窗函數(shù)h=hd.*w_box。h=hd.*w_box。h=hd.*w_box。[db,mag,pha,grd,w]=freqz_m(h,[1])。n=[0:1:M1]。subplot(2,2,1)。stem(n,hd)。 % 理想脈沖響應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)。 % 實際脈沖響應xlabel(39。n39。)。ylabel(39。h(n)39。)。 title(39。Actual Impulse Response39。)。subplot(2,2,4)。plot(w*fs/(2*pi),db)。 % 幅度響應（dB）axis([0 40 50 0])。xlabel(39。Frequency (Hz)39。)。ylabel(39。Decibels39。)。title(39。Magnitude Response in dB39。)。grid產(chǎn)生理想狀態(tài)低通濾波器function hd=ideal_lp(wc,M)。alpha=(M1)/2。n=[0:1:(M1)]。m=nalpha+eps。hd=sin(wc*m)./(pi*m)。對頻譜圖分析function [db,mag,pha,grd,w] = freqz_m(b,a)。[H,w] = freqz(b,a,1000, 39。whole39。 )。H = (H(1:1:501))39。 w = (w(1:1:501))39。mag = abs(H)。db = 20*log10((mag+eps)/max(mag))。pha = angle(H)。grd = grpdelay(b,a,w)。程序功能：畫出濾波后信號的時域波形和頻譜sf=filter(h,[1],s)。 % sf為濾波后的信號plot(t,sf)。 % 畫出濾波后信號的時域波形xlabel(39。時間 (s)39。)。ylabel(39。時域波形39。)。axis([0 1 1 1])。title(39。濾波后信號的時域波形39。)。SF=fft(sf,512)。 % 對sf進行快速傅里葉變換w=(0:255)/256*(Fs/2)。subplot(1,2,1)。plot(w,abs(SF(1:256)))。 % 畫出濾波后信號的幅度圖xlabel(39。頻率 (Hz)39。)。ylabel(39。幅度譜39。)。title(39。濾波后信號的幅度譜39。)。 grid。subplot(1,2,2)。 plot(w,angle(SF(1:256)))。 % 畫出濾波后信號的相位圖xlabel(39。頻率 (Hz)39。)。ylabel(39。相位譜39。)。title(39。濾波后信號的相位譜39
范文總結相關推薦
鄂ICP備17016276號-1