最新數(shù)字信號(hào)處理課程設(shè)計(jì)對(duì)音樂(lè)信號(hào)的各種處理(已修改)
2025-07-11 03:10 本頁(yè)面
【正文】 實(shí)驗(yàn)1 音樂(lè)信號(hào)的音譜和頻譜觀察使用wavread語(yǔ)句讀取音樂(lè)信號(hào)，獲取抽樣率；輸出音樂(lè)信號(hào)的波形和頻譜，觀察現(xiàn)象；使用sound語(yǔ)句播放音樂(lè)信號(hào)，注意不同抽樣率下的音調(diào)變化，解釋現(xiàn)象。clear all。close all。clc[a,fs,bit]=wavread(39。c:\MATLAB6p5\work\39。)。size(a)。y1=a(:,1)。a1=y1(10000:60000)figure。subplot(2,1,1),plot(a)。subplot(2,1,2),plot(a1)。x1=resample(a1,2,1)。 %y=resample(x,p,q)返回量的長(zhǎng)度是向量x的p/q倍sound(x1,fs)。%sound(a,fs)。N1=length(a1)。F1=fft(a1,N1)。w=2/N1*[0:N11]。 %頻譜圖橫坐標(biāo)設(shè)置figure。plot(w,abs(F1))。N2=length(a1)。t=0:1/N2:1/N2*(N21)。title(39。傅利葉變換39。)。 %傅利葉變換。figure。plot(a1)。title(39。時(shí)域波形39。)。 %時(shí)域波形。1， 以二倍的抽樣率聽(tīng)聲音信號(hào)時(shí)，音樂(lè)播放的特別快，像被壓縮了，播放的時(shí)間比原信號(hào)短。2， 以二分之一的抽樣率聽(tīng)聲音信號(hào)時(shí)，音樂(lè)播放的特別慢，像被拉長(zhǎng)了，播放的時(shí)間比原信號(hào)長(zhǎng)。3， *pi 實(shí)驗(yàn)2 音樂(lè)信號(hào)的抽?。p抽樣）觀察音樂(lè)信號(hào)頻率上限，選擇適當(dāng)?shù)某槿￠g隔對(duì)信號(hào)進(jìn)行減抽樣（給出兩種抽取間隔，代表混疊和非混疊）；輸出減抽樣音樂(lè)信號(hào)的波形和頻譜，觀察現(xiàn)象，給出理論解釋?zhuān)徊シ艤p抽樣音樂(lè)信號(hào)，注意抽樣率的改變，比較不同抽取間隔下的聲音，解釋現(xiàn)象。clear all。close all。clc[a,fs,bit]=wavread(39。c:\MATLAB6p5\work\39。)。size(a)。y1=a(:,1)。a1=y1(10000:60000)D=2。 %減抽樣；l=length(a1)。yd=a1(1:D:l)。sound(yd,fs/D)。N3=length(yd)。t=0:1/N3:1/N3*(N31)。 %橫坐標(biāo)設(shè)置figure。plot(yd)。title(39。減抽樣時(shí)域波形39。)。 %時(shí)域波形。xlabel(39。t39。)。ylabel(39。幅度39。)。N4=length(yd)。F2=fft(yd,N4)。w=2/N4*[0:N41]。figure。plot(w,abs(F2))。title(39。減抽樣頻譜39。)。 %減抽樣頻譜xlabel(39。f39。)。ylabel(39。幅度39。)。D=2,減抽樣D=4,減抽樣1， *pi，當(dāng)D=2時(shí)，頻譜剛好不混疊，當(dāng)D2時(shí)，頻譜就會(huì)混疊。2， 減抽樣后的音樂(lè)信號(hào)聽(tīng)起來(lái)變得尖銳，有失真。3， 抽樣率隨著抽樣間隔的增大而逐漸變小，聲音越來(lái)越失真，音調(diào)變得急促，而尖銳，信號(hào)產(chǎn)生混疊實(shí)驗(yàn)3 音樂(lè)信號(hào)的AM調(diào)制觀察音樂(lè)信號(hào)的頻率上限，選擇適當(dāng)調(diào)制頻率對(duì)信號(hào)進(jìn)行調(diào)制（給出高、低兩種調(diào)制頻率）；輸出調(diào)制信號(hào)的波形和頻譜，觀察現(xiàn)象，給出理論解釋?zhuān)徊シ耪{(diào)制音樂(lè)信號(hào)，注意不同調(diào)制頻率下的聲音，解釋現(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)。N=length(a1)。n1=0:[N1]。y=cos(*pi*n1)。 %產(chǎn)生余弦信號(hào)N1=length(a1)。F1=fft(y,N1)。w=2/N1*[0:N11]。figure。plot(w,abs(F1))。title(39。cos(wt)頻譜39。)。 %余弦頻譜xlabel(39。f39。)。ylabel(39。幅度39。)。N2=length(a1)。F2=fft(a1,N2)。w=2/N2*[0:N21]。figure。plot(w,abs(F2))。title(39。yinyue頻譜39。)。 %原信號(hào)頻譜xl
試題試卷相關(guān)推薦
鄂ICP備17016276號(hào)-1