有限脈沖響應(yīng)數(shù)字濾波器設(shè)計實驗報告(存儲版)
【正文】 wc=(wp+ws)/2/pi。 xlabel(39。)。)。 要求通帶截止頻率為 rad 和 rad，阻帶截止頻率為 rad 和 rad，通帶最大衰減為 dB，阻帶最小衰減為 40 dB。 wc=[(wp1+ws1)/2/pi,(ws2+wp2)/2/pi]。)。 wk=(2*pi/M)*k。幅度 /dB39。 rs=50。 M=1024。w/\pi39。grid on figure(3) plot(wk/pi,angle(Hk(k+1))/pi)。相頻特性曲線 39。 m1=0。m9=0。m17=0。 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*。m17=m16。m9=m8。m1=x(m)。xlabel(39。原信號幅度頻譜 (采樣點數(shù)為 500)39。 figure(3) plot(n,y)。n39。)。 Fs=3800。 N=ceil((m+1)*2*pi/datB+1)。 Hk=Ak.*exp(1j*thetak)。 axis([0,1,80,5])。損耗函數(shù) 39。)。t=n*T。m7=0。 m15=0。m11=m10。m3=m2。幅值 39。 axis([0,1500,0,180])。幅值 39。)。 title(39。)。 wp=2*pi*fp/Fs。 Hk=fft(hn,M)。 [B,A]=ellip(N1,Rp,As,wp1) [Hk1,wk1]=freqz(B,A)。,3)。)。:39。相位 /\pi39。FIR39。 datB=wpws。 Ak(Np+2)=T。 wk=(2*pi/M)*k。) [Hk1,wk1]=freqz(B,A)。,3)。)。:39。相位 /\pi39。FIR39。 參考文獻 [1]、高西全 ,丁玉美 .數(shù)字信號處理 (第三版 ).西安電子科技大學(xué)出版社 .2021 年 8 月 [2]、劉樹棠 [譯 ].數(shù)字信號處理 (MATLAB 版 )(第二版 ).西安交通大學(xué)出版社 .2021 年 1 月 。IIR39。)。 figure(2) plot(wk1/pi,pah/pi)。) title(39。,39。 0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 1 1 2 0 1 0 0 8 0 6 0 4 0 2 0020w/ ?幅度/dB損耗函數(shù)曲線 IIRF IR0 1 2 3 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 11 0 . 8 0 . 6 0 . 4 0 . 200 . 20 . 40 . 60 . 81w/ ?相位/?相頻特性曲線 IIRF IR[B,A]=ellip(N1,Rp,As,wp1,39。 Hk=fft(hn,M)。Ns=N2*Np1。 wp=2*pi*fp/Fs。IIR39。)。 figure(2) plot(wk1/pi,pah/pi)。) title(39。,39。As=50。 hn=fir1(N,wc,kaiser(N+1,alph))。fs=600。ylabel(39。 figure(4) plot(fs,abs(S))。 ylabel(39。)。fs=n/(N*T)。)。m5=m4。m13=m12。m13=0。m5=0。 f=2800。ylabel(39。)。 figure(2)。Ak(NNp)=T。wc=wp。)。 title(39。)。幅值 39。 axis([0,1500,0,180])。幅值 39。m3=m2。m11=m10。m19=m18。m23=0。 m15=0。m7=0。t=n*T。)。損耗函數(shù) 39。 axis([0,1,80,5])。 wc=(wp+ws)/2/pi。 fp=130。)。)。 xlabel(39。 Bt=wp2wp1。)。xlabel(39。 title(39。 n=0:N1。 N0=ceil(*pi/Bt)。)。xlabel(39。,39。 Hk=fft(hn_hamming,M) plot(wk/pi,20*log10(abs(Hk(k+1))),39。 plot(wk/pi,20*log10(abs(Hk(k+1))),39。 title(39。矩形窗 39。)。,2)。high39。,boxcar(N))。損耗函數(shù) 39。 axis([0,1,80,5])。矩形窗 39。,3)。 figure()。)。 legend(39。o39。linewidth39。 hn_hanning=fir1(N1,wc/pi,hanning(N))。相頻特性曲線 39。)。,)。)。IIR濾波器 39。linewidth39。 M=1024。 [Hk,wk1]=freqz(B,A,1000)。fs=20210。成 績 ： 《數(shù)字信號處理》 作業(yè)與上機實驗 （第二章） 班 級： 學(xué) 號： 姓 名： 任課老師： 完成時間： 信息與通信工程學(xué)院 2021— 2021 學(xué)年第 1 學(xué)期 第 7 章 有限脈沖響應(yīng)數(shù)字濾波器設(shè)計 教材 p238： x(t) = s(t) + v(t),其中 v(t)是干擾， s(t)與 v(t)的頻譜不混疊 ，其幅度譜如題 19 圖所示。 fp=15000。 [B,A]=ellip(N1,Rp,As,wp1)。 hn=fir1(N,wc,kaiser(N+1,alph))。,39。,39。幅度 /dB39。linewidth39。IIR濾波器 39。 title(39。 hn_boxcar=fir1(N1,wc/pi,boxcar(N))。,39。 plot(n,hn_blackman,39。)。布萊克曼 窗 39。 wk=(2*pi/M)*k。linewidth39。 hold off legend(39。)。 title(39。high39。 hn_blackman=fir1(N1,wc/pi,39。linewidth39。o39。 legend(39。)。 figure()。,3)。矩形窗 39。 axis([0,1,80,5])。損耗函數(shù) 39。 Bt=wpws。 Hk=fft(hn,M)。)。 axis([0,1,80,5])。損耗函數(shù) 39。ws2=*pi。 stem(n,hn)。單位沖激響應(yīng) 39。w/\pi39。grid on (2)帶通 FIR濾波器的單位脈沖響應(yīng)、損耗
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1