【文章內(nèi)容簡介】 %輸出連續(xù)譜圖 grid on title(39。原始音樂信號頻域波型 39。) xlabel(39。Frequency/Hz39。) ylabel(39。Magnitude39。) n=0:N1。 x1=cos(n*pi*)。 %進(jìn)行調(diào)制的余弦信號 figure subplot(2,2,1)。 plot(x1) %輸出余弦信號的連續(xù)波形 title(39。余弦信號時域波型 39。) xlabel(39。Time39。) ylabel(39。Magnitude39。) F2=fft(x1,N)。 %對余弦信號進(jìn)行 N點 fft變換 subplot(2,2,2)。 plot(w,abs(F2)) %輸出余弦信號的連續(xù)譜圖 grid on title(39。余弦信號頻域波型 39。) xlabel(39。Frequency/Hz39。) ylabel(39。Magnitude39。) X=x.*x139。 subplot(2,2,3)。 plot(X) %輸出調(diào)制后音樂信號的連續(xù)波形 F3=fft(X,N)。 grid on title(39。高頻調(diào)制后音樂信號時域波型 39。) xlabel(39。Time39。) ylabel(39。Magnitude39。) subplot(2,2,4)。 plot(w,abs(F3)) %輸出調(diào)制后音樂信號的連續(xù)譜圖 grid on title(39。高頻調(diào)制后音樂信號頻域波型 39。) xlabel(39。Frequency/Hz39。) ylabel(39。Magnitude39。) X1=X.*x139。 %對調(diào)制后的音樂信號進(jìn)行解調(diào) F4=fft(X1,N)。 %對解調(diào)后的音樂信號進(jìn)行 N點 fft變換 figure subplot(2,1,1)。plot(X1)。 %輸出解調(diào)后的音樂信號的連續(xù)波形 grid on title(39。解制后音樂信號時域波型 39。) xlabel(39。Time39。) ylabel(39。Magnitude39。) subplot(2,1,2)。plot(w,abs(F4))。 %輸出解調(diào)后音樂信號的連續(xù)譜圖 grid on title(39。解制后音樂信號頻域波型 39。) xlabel(39。Frequency/Hz39。) ylabel(39。Magnitude39。) Rp=1。Rs=15。 %Rp 為通帶內(nèi)最大波 紋， Rs 為阻帶內(nèi)最小衰減 [N1,Wc]=buttord(,Rp,Rs)。 %N1為濾波器的階數(shù)， Wc為巴特沃斯濾波器的截止頻率 [B,A]=butter(N1,Wc)。 %B為巴特沃斯濾波器的分子， A為巴特沃斯濾波器的分母 [H,W]=freqz(B,A)。 %H為數(shù)字濾波器的頻率響應(yīng) figure subplot(1,1,1)。 plot(W/pi,abs(H))。 %輸出濾波器的頻率響應(yīng)圖 grid on title(39。數(shù)字巴特沃斯濾波器 39。) xlabel(39。Digital Frequency/pi39。) ylabel(39。Magnitude39。) y=filter(B,A,X1)。 %對解調(diào)后的音樂信號進(jìn)行濾波 % sound(y,fs)。 %播放濾波后的音樂信號 figure subplot(2,1,1) plot(y)。 %輸出用巴特 沃斯濾波后的連續(xù)波形 grid on title(39。巴特沃斯濾波后音樂信號時域波型 39。) xlabel(39。Time39。) ylabel(39。Magnitude39。) y=fft(y,N)。 subplot(2,1,2) plot(w,abs(y))。 %輸出巴特沃斯濾波后的連續(xù)譜圖 grid on title(39。濾波后音樂信號頻域波型 39。) xlabel(39。Frequency/Hz39。) ylabel(39。Magnitude39。) N2=33。wc=*pi。 %N為理想低通濾波器的階數(shù)， Wc為濾波器的截止頻率 hd=ideal(N2,wc)。 %調(diào)用理想低通濾波器函數(shù) w2=boxcar(N2)。 %產(chǎn)生矩形窗函數(shù) w3=blackman(N2)。 %產(chǎn)生布萊克曼窗函數(shù) h1=hd.*w239。 %加窗設(shè)計矩形窗濾波器 h2=hd.*w339。 %加窗設(shè)計布萊克曼窗濾波器 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為連續(xù)頻譜的數(shù)字角頻率橫坐標(biāo) plot(w2,db1)。 %輸出矩形窗的分貝圖 grid on title(39。矩形窗濾波器 39。) xlabel(39。Frequency39。) ylabel(39。DB39。) y1=conv(X1,h1)。 %用矩形窗對解調(diào)后的音樂信號進(jìn)行濾波 sound(y1,fs) %播放濾波后的音樂信號 figure subplot(2,1,1)。plot(y1)。 %輸 出濾波之后的音樂信號的連續(xù)波型 grid on title(39。矩形窗濾波后的音樂信號時域波形 39。) xlabel(39。Time39。) ylabel(39。Magtitude39。) y1=fft(y1,M)。 subplot(2,1,2)。plot(w2,abs(y1))。 %輸出濾波之后音樂信號的連續(xù)譜圖 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)。 %用布萊克曼窗對解調(diào)后的音樂信號進(jìn)行濾波 sound(y2,fs) %播放濾波之后的音樂信號 figure subplot(2,1,1)。plot(y2)。 %輸出用布萊克曼 窗濾波之后的連續(xù)波形 grid on title(39。布萊克曼窗濾波后的音樂信號時域波形 39。) xlabel(39。Time39。) ylabel(39。Magtitude39。) y2=fft(y2,M)。 subplot(2,1,2)。plot(w2,abs(y2))。 %輸出用布萊克曼窗濾波之后的連續(xù)譜圖 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高頻調(diào)制后音樂信號時域波型T i m eMagnitude0 0 . 5 1 1 . 5 2050010001500高頻調(diào)制后音樂信號頻域波型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數(shù)字巴特沃斯濾波器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