有限脈沖響應(yīng)數(shù)字濾波器設(shè)計(jì)實(shí)驗(yàn)報(bào)告-文庫(kù)吧資料
2025-02-12 08:30本頁(yè)面
【正文】 l(39。w/\pi39。 axis([0,1,80,5])。 figure(2)。 k=0:M/21。 M=1024。 Hk=Ak.*exp(1j*thetak)。Ak(NNp)=T。 Ak=[zeros(1,Np+1),ones(1,Ns),zeros(1,Np)]。 Np=fix(wc/(2*pi/N))。 N=ceil((m+1)*2*pi/datB+1)。wc=wp。ws=2*pi*fs/Fs。fs=100。 Fs=3800。)。ylabel(39。頻率 /Hz39。)。 title(39。 figure(4) plot(fs,abs(S))。 S=fft(y,N)。n39。)。 ylabel(39。信號(hào) y(n)39。 figure(3) plot(n,y)。幅值 39。)。 xlabel(39。原信號(hào)幅度頻譜 (采樣點(diǎn)數(shù)為 500)39。 axis([0,1500,0,180])。fs=n/(N*T)。)。xlabel(39。幅值 39。)。 title(39。m1=x(m)。m3=m2。 m5=m4。m7=m6。m9=m8。m11=m10。m13=m12。 m15=m14。m17=m16。m19=m18。m21=m20。m23=m22。 for m=1:length(x) y(m)=*x(m)+m1*****3... m6*+m7*+m8*+m9*+m10*+m11*+m12*+... m13*+m14*+m15*+m16*+m17**... m19****+m23*+m24*。m23=0。m21=0。m19=0。m17=0。 m15=0。m13=0。m11=0。m9=0。m7=0。m5=0。m3=0。 m1=0。t=n*T。 f=2800。 (2)數(shù)字濾波器的損耗函數(shù)和相頻特性分別如圖 1 12所示： 圖 11 損耗函數(shù)曲線 圖 12 相頻特性曲線 二、 按直接型網(wǎng)絡(luò)結(jié)構(gòu)編程 編寫濾波程序 : (1)matlab 代碼： N=500。相頻特性曲線 39。)。 ylabel(39。w/\pi39。grid on figure(3) plot(wk/pi,angle(Hk(k+1))/pi)。損耗函數(shù) 39。)。 ylabel(39。w/\pi39。 axis([0,1,80,5])。 figure(2)。 k=0:M/21。 M=1024。 wc=(wp+ws)/2/pi。 alph=*(rs21)^+*(rs21)。 ws=2*pi*fs/Fs。 rs=50。 fp=130。)。 title(39。幅度 /dB39。)。xlabel(39。 plot(wk/pi,20*log10(abs(Hk(k+1))))。 wk=(2*pi/M)*k。)。 title(39。h(n)39。)。 xlabel(39。 n=0:N1。 M=1024。 wc=[(wp1+ws1)/2/pi,(ws2+wp2)/2/pi]。 Bt=wp2wp1。ws1=*pi。 (1)matlab 代碼： wp1=*pi。 要求通帶截止頻率為 rad 和 rad，阻帶截止頻率為 rad 和 rad，通帶最大衰減為 dB，阻帶最小衰減為 40 dB。)。 title(39。幅度 /dB39。)。xlabel(39。 plot(wk/pi,20*log10(abs(Hk(k+1))))。 wk=(2*pi/M)*k。)。 title(39。h(n)39。)。 xlabel(39。 n=0:N1。 M=1024。high39。 wc=(wp+ws)/2/pi。 N0=ceil(*pi/Bt)。ws=*pi。顯示所設(shè)計(jì)的單位脈沖響應(yīng) h(n)的數(shù)據(jù)，并畫出損耗函數(shù)曲線。 0 2 4 6 8 10 12 14 16 18 20 0 . 4 0 . 200 . 20 . 40 . 60 . 81nh(n)單位沖激響應(yīng) 矩形窗漢寧窗哈明窗布萊克曼窗0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 1 8 0 7 0 6 0 5 0 4 0 3 0 2 0 1 00w/ ?幅度/dB損耗函數(shù) 矩形窗漢寧窗哈明窗布萊克曼窗 MATLAB工具箱函數(shù) fir1設(shè)計(jì)線性相位高通 FIR濾波器。)。 title(39。幅度 /dB39。)。xlabel(39。)。,39。,39。,39。 hold off legend(39。*39。)。 Hk=fft(hn_hamming,M) plot(wk/pi,20*log10(abs(Hk(k+1))),39。linewidth39。:39。,2)。 plot(wk/pi,20*log10(abs(Hk(k+1))),39。 wk=(2*pi/M)*k。Hk=fft(hn_boxcar,M)。)。 title(39。布萊克曼窗 39。哈明窗 39。漢寧窗 39。矩形窗 39。)。ylabel(39。n39。)。 plot(n,hn_blackman,39。linewidth39。+39。,2)。,39。 hold on plot(n,hn_hanning,39。 n=0:N1。high39。,hamming(N))。 hn_hamming=fir1(N1,wc/pi,39。high39。,boxcar(N))。 hn_boxcar=fir1(N1,wc/pi,39。 (1)matlab 代碼： wc=pi/4。 (2)四種窗函數(shù)設(shè)計(jì)的單位脈沖響應(yīng)的比較如圖 3所示： 圖 3 單位脈沖響應(yīng)比較圖 (3)四種窗函數(shù) 設(shè)計(jì)的損耗函數(shù)的比較如圖 4所示： 圖 4 損耗函數(shù)比較圖 0 2 4 6 8 10 12 14 16 18 20 0 . 1 0 . 0 500 . 0 50 . 10 . 1 50 . 20 . 2 50 . 3nh(n)單位沖激響應(yīng) 矩形窗漢寧窗哈明窗布萊克曼窗0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 1 8 0 7 0 6 0 5 0 4 0 3 0 2 0 1 00w/ ?幅度/dB損耗函數(shù) 矩形窗漢寧窗哈明窗布萊克曼窗 (4)運(yùn)行結(jié)果分析：由圖 4可見(jiàn)，當(dāng)濾波器長(zhǎng)度 N不變時(shí)，矩形窗設(shè)計(jì)的濾波器的過(guò)渡帶最窄，阻帶最小衰減最??；布萊克曼窗設(shè)計(jì)的濾波器的過(guò)渡帶最寬，同時(shí)阻帶最小衰減最大。損耗函數(shù) 39。)。ylabel(39。w/\pi39。 axis([0,1,80,5])。布
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1