數字信號處理課程設計畢業設計word格式(編輯修改稿)

　

【文章內容簡介】 %輸出連續譜圖 grid on title(39。原始音樂信號頻域波型 39。) xlabel(39。Frequency/Hz39。) ylabel(39。Magnitude39。) n=0:N1。 x1=cos(n*pi*)。 %進行調制的余弦信號 figure subplot(2,2,1)。 plot(x1) %輸出余弦信號的連續波形 title(39。余弦信號時域波型 39。) xlabel(39。Time39。) ylabel(39。Magnitude39。) F2=fft(x1,N)。 %對余弦信號進行 N點 fft變換 subplot(2,2,2)。 plot(w,abs(F2)) %輸出余弦信號的連續譜圖 grid on title(39。余弦信號頻域波型 39。) xlabel(39。Frequency/Hz39。) ylabel(39。Magnitude39。) X=x.*x139。 subplot(2,2,3)。 plot(X) %輸出調制后音樂信號的連續波形 F3=fft(X,N)。 grid on title(39。高頻調制后音樂信號時域波型 39。) xlabel(39。Time39。) ylabel(39。Magnitude39。) subplot(2,2,4)。 plot(w,abs(F3)) %輸出調制后音樂信號的連續譜圖 grid on title(39。高頻調制后音樂信號頻域波型 39。) xlabel(39。Frequency/Hz39。) ylabel(39。Magnitude39。) X1=X.*x139。 %對調制后的音樂信號進行解調 F4=fft(X1,N)。 %對解調后的音樂信號進行 N點 fft變換 figure subplot(2,1,1)。plot(X1)。 %輸出解調后的音樂信號的連續波形 grid on title(39。解制后音樂信號時域波型 39。) xlabel(39。Time39。) ylabel(39。Magnitude39。) subplot(2,1,2)。plot(w,abs(F4))。 %輸出解調后音樂信號的連續譜圖 grid on title(39。解制后音樂信號頻域波型 39。) xlabel(39。Frequency/Hz39。) ylabel(39。Magnitude39。) Rp=1。Rs=15。 %Rp 為通帶內最大波 紋， Rs 為阻帶內最小衰減 [N1,Wc]=buttord(,Rp,Rs)。 %N1為濾波器的階數， Wc為巴特沃斯濾波器的截止頻率 [B,A]=butter(N1,Wc)。 %B為巴特沃斯濾波器的分子， A為巴特沃斯濾波器的分母 [H,W]=freqz(B,A)。 %H為數字濾波器的頻率響應 figure subplot(1,1,1)。 plot(W/pi,abs(H))。 %輸出濾波器的頻率響應圖 grid on title(39。數字巴特沃斯濾波器 39。) xlabel(39。Digital Frequency/pi39。) ylabel(39。Magnitude39。) y=filter(B,A,X1)。 %對解調后的音樂信號進行濾波 % sound(y,fs)。 %播放濾波后的音樂信號 figure subplot(2,1,1) plot(y)。 %輸出用巴特 沃斯濾波后的連續波形 grid on title(39。巴特沃斯濾波后音樂信號時域波型 39。) xlabel(39。Time39。) ylabel(39。Magnitude39。) y=fft(y,N)。 subplot(2,1,2) plot(w,abs(y))。 %輸出巴特沃斯濾波后的連續譜圖 grid on title(39。濾波后音樂信號頻域波型 39。) xlabel(39。Frequency/Hz39。) ylabel(39。Magnitude39。) N2=33。wc=*pi。 %N為理想低通濾波器的階數， Wc為濾波器的截止頻率 hd=ideal(N2,wc)。 %調用理想低通濾波器函數 w2=boxcar(N2)。 %產生矩形窗函數 w3=blackman(N2)。 %產生布萊克曼窗函數 h1=hd.*w239。 %加窗設計矩形窗濾波器 h2=hd.*w339。 %加窗設計布萊克曼窗濾波器 M=length(x)。 fh1=fft(h1,M)。 db1=20*log10(abs(fh1(1)./(abs(fh1)+eps)))。 %矩形窗的幅度分貝表示 fh2=fft(h2,M)。 db2=20*log10(abs(fh2(1)./(abs(fh2)+eps))) 。 %布萊克曼窗的幅度分貝表示 figure w2=2/M*[0:M1]。 %w為連續頻譜的數字角頻率橫坐標 plot(w2,db1)。 %輸出矩形窗的分貝圖 grid on title(39。矩形窗濾波器 39。) xlabel(39。Frequency39。) ylabel(39。DB39。) y1=conv(X1,h1)。 %用矩形窗對解調后的音樂信號進行濾波 sound(y1,fs) %播放濾波后的音樂信號 figure subplot(2,1,1)。plot(y1)。 %輸 出濾波之后的音樂信號的連續波型 grid on title(39。矩形窗濾波后的音樂信號時域波形 39。) xlabel(39。Time39。) ylabel(39。Magtitude39。) y1=fft(y1,M)。 subplot(2,1,2)。plot(w2,abs(y1))。 %輸出濾波之后音樂信號的連續譜圖 grid on title(39。矩形窗濾波后的音樂信號頻域波形 39。) xlabel(39。Frequency39。) ylabel(39。Magtitude39。) figure plot(w2,db2) %輸出布萊克曼窗的分貝圖 grid on title(39。布萊克曼窗濾波器 39。) xlabel(39。Frequency39。) ylabel(39。DB39。) y2=conv(X1,h2)。 %用布萊克曼窗對解調后的音樂信號進行濾波 sound(y2,fs) %播放濾波之后的音樂信號 figure subplot(2,1,1)。plot(y2)。 %輸出用布萊克曼 窗濾波之后的連續波形 grid on title(39。布萊克曼窗濾波后的音樂信號時域波形 39。) xlabel(39。Time39。) ylabel(39。Magtitude39。) y2=fft(y2,M)。 subplot(2,1,2)。plot(w2,abs(y2))。 %輸出用布萊克曼窗濾波之后的連續譜圖 grid on title(39。布萊克曼窗濾波后的音樂信號頻域波形 39。) xlabel(39。Frequency39。) ylabel(39。Magtitude39。) 3. 輸出波形 0 0 . 5 1 1 . 5 2 2 . 5x 1 051 0 . 8 0 . 6 0 . 4 0 . 200 . 20 . 40 . 60 . 81原始音樂信號時域波型T i m eMagnitude 0 0 . 2 0 . 4 0 . 6 0 . 8 1 1 . 2 1 . 4 1 . 6 1 . 8 2050010001500202025003000原始音樂信號頻域波型Fr e q u e n c y / H zMagnitude 0 1 2 3x 1 051 0 . 500 . 51余弦信號時域波型T i m eMagnitude0 0 . 5 1 1 . 5 2051015x 1 04余弦信號頻域波型Fr e q u e n c y / H zMagnitude0 1 2 3x 1 051 0 . 500 . 51高頻調制后音樂信號時域波型T i m eMagnitude0 0 . 5 1 1 . 5 2050010001500高頻調制后音樂信號頻域波型Fr e q u e n c y / H zMagnitude 0 0 . 5 1 1 . 5 2 2 . 5x 1 051 0 . 500 . 51解制后音樂信號時域波型T i m eMagnitude0 0 . 2 0 . 4 0 . 6 0 . 8 1 1 . 2 1 . 4 1 . 6 1 . 8 2050010001500解制后音樂信號頻域波型Fr e q u e n c y / H zMagnitude 0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 100 . 20 . 40 . 60 . 811 . 21 . 4數字巴特沃斯濾波器D i g i t a l Fr e q u e n c y / p iMagnitude 0 0 . 5 1 1 . 5 2 2 . 5x 1 05 0 . 500 . 51巴特沃斯濾波后音樂信號時域波型T i m eMagnitude0 0 . 2 0 . 4 0 . 6 0 . 8 1 1 . 2 1 . 4 1 . 6 1 . 8 2050010001500濾波后音樂信號頻域波型Fr e q u e n c y / H zMagnitude 0 0 . 2 0 . 4 0 . 6 0 . 8 1 1 . 2 1 . 4 1 . 6 1 . 8 2 1 8 0 1 6 0 1 4 0 1