最新數(shù)字信號處理課程設(shè)計對音樂信號的各種處理(編輯修改稿)
2025-07-26 03:10 本頁面
【文章內(nèi)容簡介】 lackman(N)。 %布萊克曼窗h1=hd.*w139。h2=hd.*w239。y3=conv(x2,h1)。 %解調(diào)后信號與矩形窗函數(shù)卷積y4=conv(x2,h2)。 %解調(diào)后信號與布萊克曼窗函數(shù)卷積%sound(y3,fs)。 %播放矩形窗濾波后信號sheng yin%sound(y4,fs)。 %播放布萊克曼窗濾波后信號sheng yin figure。subplot(2,1,1),plot(y3)。title(39。矩形窗濾波后信號39。)。subplot(2,1,2)。plot(y3)。title(39。布萊克曼窗濾波后信號39。)。fh1=fft(h1,N1)。db1=20*log10(abs(fh1(1)./(abs(fh1)+eps)))。 %理想低通濾波器加窗后幅度響應(yīng)fh2=fft(h2,N1)。db2=20*log10(abs(fh2(1)./(abs(fh2)+eps)))。w=2/N1*[0:N11]。figure。subplot(3,1,1),stem(h1)。 %矩形窗函數(shù)grid on。title(39。矩形窗39。)。xlabel(39。n39。)。ylabel(39。h(n)39。)。subplot(3,1,2),plot(w,abs(fh1))。grid on。title(39。矩形窗39。)。xlabel(39。w39。)。ylabel(39。H(k)39。)。subplot(3,1,3)。plot(w,db1)。 %矩形窗函數(shù)分貝圖grid on。figure。subplot(3,1,1)。stem(h2)。grid on。title(39。布萊克曼窗39。)。xlabel(39。n39。)。ylabel(39。h(n)39。)。subplot(3,1,2)。plot(w,abs(fh2))。grid on。title(39。布萊克曼窗39。)。xlabel(39。w39。)。ylabel(39。H(k)39。)。subplot(3,1,3)。plot(w,db2)。 %布萊克曼窗窗函數(shù)分貝圖grid on。w=2/N1*[0:N11]。Fy1=fft(y3,N1)。Fy2=fft(y4,N1)。figure。subplot(2,1,1)。plot(w,abs(Fy1))。title(39。矩形窗濾波后頻譜39。)。subplot(2,1,2)。plot(w,abs(Fy2))。title(39。布萊克曼窗濾波后頻譜39。)。1， 解調(diào)后信號頻譜在高頻和低頻處均有一部分，且成對稱分布，需要濾掉高頻才可大致還原原信號。2， ，，阻帶衰減15db。濾波效果很好，基本還原了原信號。3， 使用窗函數(shù)濾波要根據(jù)過渡帶寬算階數(shù)N。4， 使用矩形窗濾波，矩形窗過渡帶窄，但是阻帶有波紋，高頻部分有小部分未濾掉。5， 使用布萊克曼窗濾波，布萊克曼窗過渡帶寬，但是阻帶較好。6， 使用矩形窗和布萊克曼窗濾波，效果都行，基本都能還原原信號實驗5 音樂信號的濾波去噪，頻率為3kHz、5kHz、8kHz的三余弦混合噪聲，觀察噪聲頻譜以及加噪后音樂信號的音譜和頻譜，并播放音樂，感受噪聲對音樂信號的影響；（可用rand語句產(chǎn)生），觀察噪聲頻譜以及加噪后音樂信號的音譜和頻譜，并播放音樂，感受噪聲對音樂信號的影響；根據(jù)步驟、觀察到的頻譜，選擇合適指標設(shè)計濾波器進行濾波去噪，觀察去噪后信號音譜和頻譜，并播放音樂，解釋現(xiàn)象。clear all。close all。clc[a,fs,bit]=wavread(39。c:\MATLAB6p5\work\39。)。size(a)。y1=a(:,1)。a1=y1(10000:100000)。%sound(a1,fs)。 %播放yuan的信號N1=length(a1)。n1=0:[N11]。x1=*cos(2*pi*3000*n1/fs)。x2=*cos(2*pi*5000*n1/fs)。x3=*cos(2*pi*8000*n1/fs)。x4=x1+x2+x3。figure。plot(x4)。%疊加
試題試卷相關(guān)推薦
鄂ICP備17016276號-1