數(shù)字信號課程設(shè)計-在線瀏覽
2024-08-02 00:49本頁面
【正文】 熟悉 FFT 算法原理和 FFT 子程序的應(yīng)用。 二、上級程序： x1=[1 1 1 1]。 x3=[4,3,2,1,1,2,3,4]。stem(x1)。序列 x139。n39。 ylabel(39。)。stem(x2)。序列 x239。n39。 ylabel(39。)。stem(x3)。序列 x339。n39。 ylabel(39。)。 F1x1=fft(x1,N1)。 F1x3=fft(x3,N1)。stem(abs(F1x1))。序列 x1 的 8 點譜分析 39。n39。 ylabel(39。)。stem(abs(F1x2))。序列 x2 的 8 點譜分析 39。n39。 ylabel(39。)。stem(abs(F1x3))。序列 x3 的 8 點譜分析 39。n39。 ylabel(39。)。 F2x1=fft(x1,N2)。 F2x3=fft(x3,N2)。stem(abs(F2x1))。序列 x1 的 16 點譜分析 39。 xlabel(39。)。F2x139。 subplot(1,3,2)。 title(39。)。n39。 ylabel(39。)。stem(abs(F2x3))。序列 x3 的 16 點譜分析 39。 xlabel(39。)。F2x339。 N3=16。 x4=cos(n*pi/4)。 figure。stem(n,x4)。序列 x439。n39。 ylabel(39。)。stem(n,x5)。序列 x539。n39。 ylabel(39。)。 F3x4=fft(x4,N4)。 figure subplot(1,2,1)。 title(39。)。n39。 ylabel(39。)。stem(abs(F3x5))。序列 x5 的 8 點譜分析 39。 xlabel(39。)。F3x539。 N5=16。 F4x5=fft(x5,N5)。stem(abs(F4x4))。序列 x4 的 16 點譜分析 39。 xlabel(39。)。F4x439。 subplot(1,2,2)。 title(39。)。n39。 ylabel(39。)。 F5x4=fft(x4,N6)。 figure w=2/N6*[0:N61]。plot(w,abs(F5x4))。序列 x4 的 64 點譜分析 39。 xlabel(39。)。F5x439。 subplot(1,2,2)。 title(39。)。w39。 ylabel(39。)。 F6x4=fft(x4,N7)。 figure w=2/N7*[0:N71]。plot(w,abs(F6x4))。序列 x4 的 256 點譜分析 39。 xlabel(39。)。F6x439。 subplot(1,2,2)。 title(39。)。w39。 ylabel(39。)。 t=0:1/64:1/64*(N81)。 figure plot(t,x6)。序列 x639。t39。 ylabel(39。)。 F9x6=fft(x6,N9)。 title(39。)。n39。 ylabel(39。)。 F10x6=fft(x6,N10)。 title(39。)。n39。 ylabel(39。)。 F11x6=fft(x6,N11)。 figure plot(w,abs(F11x6))。序列 x6 的 64 點譜分析 39。 xlabel(39。)。F11x639。 N12=256。 w=2/N12*[0:N121]。 title(39。)。w39。 ylabel(39。)。 x4(n)和 x5(n)在 N= 1 64 和 256 是的頻譜是否相同？為什么？ 答： 對 x4(n)來說，頻譜相同，因為 x4(n)的周期是 8，在 N= 1 64 和 256時，取樣點數(shù)都是周期的整數(shù)倍，包含信號的全部信息，故頻譜相同； 對 x5(n)來說，因為 x5(n)的周期是 16，在 N=8 時，對信號進行了截斷，包含的信息不完整，而在 N=1 64 和 256 時抽樣點包含信號的全部信息，故頻譜相同； x6(n)在 N= 1 64 和 256 時能否將 三條譜線分開，為什么？ 答：在 N= 16 時不能 將三條譜線分開 ， 在 N=64 和 256 時能將三條譜線分開。 10 實驗二、 用雙線性法設(shè)計 IIR 數(shù)字濾波器 一、 實驗?zāi)康? 熟悉雙線性法設(shè)計 IIR 數(shù)字濾波器的原理和方法。 二、 上機程序 T=1。 Ws=2/T*tan(*pi/2)。 Rs=15。s39。 [B,A]=butter(N,Wc,39。)。 W=0:*pi:*pi。 Hd=freqz(C,D,W)。模擬巴特沃斯濾波器 39。Frequency/Hz39。Magnitude39。 plot(W/pi,abs(Hd)) grid on title(39。) xlabel(39。) ylabel(39。) Hd_db=20*log10(abs(Hd(1)./(abs(Hd)+eps)))。 plot(W/pi,Hd_db) grid on title(39。) xlabel(39。) ylabel(39。) x=[4,2,0,4,6,4,2,4,6,6,4,4,6,6,2,6,12,8,0,16,38,60,84,90,66,32,4,2,4,8,12,12,10,6,6,6,6,4,0,0,0,0,0,2,4,0,0,0,2,2,0,0,2,2,2,2,0]。 Fx=fft(x,M)。 Fy=fft(y,M)。 subplot(1,2,1)。 title(39。)。n39。 ylabel(39。)。 subplot(1,2,2)。 title(39。)。f39。 ylabel(39。) figure。 plot(y)。濾波后的信號 39。 xlabel(39。)。y39。 f=1/M*[0:M1]* plot(f,abs(Fy))。濾波后的頻譜分析 39。 xlabel(39。)。|Fy|39。 11 0 0 . 0 5 0 . 1 0 . 1 5 0 . 2 0 . 2 5 0 . 3 0 . 3 5 0 . 4 0 . 4 5 0 . 500 . 10 . 20 . 30 . 40 . 50 . 60 . 70 . 80 . 91模擬巴特沃斯濾波器F r e q u e n c y / H zMagnitu
畢業(yè)設(shè)計相關(guān)推薦
鄂ICP備17016276號-1