數(shù)字信號處理課程設(shè)計--基于matlab的音樂信號處理和分析-資料下載頁
2025-06-06 00:46本頁面
【正文】 。clc [x,fs,nbits]=wavread(39。怒放的生命 汪峰 5s39。)。 x1=x(:,1)。 %獲取單列音樂信號并對其做 FFT變換 N=length(x1)。 fx1=fft(x1)。 w1=2/N*(0:N/21)。 wp=。ws=。rp=1。rs=50。 %設(shè)計巴特沃斯濾波器 [N4,Wc]=buttord(wp,ws,rp,rs)。 [B,A]=butter(N4,Wc)。 [Hd,w]=freqz(B,A)。 lox1=filter(B,A,x1)。 %使用 巴特沃斯濾波器濾除音樂信號的高頻部分并對所得 %音樂信號做 FFT變換 N1=length(lox1)。 flox1=fft(lox1)。 w2=2/N1*(0:N1/21)。 figure %畫出單列音樂信號的頻譜圖、巴特沃斯濾波器的頻率響應(yīng)曲線和濾除 %高頻后的音樂信號的頻譜圖 subplot(3,1,1)。plot(w1,abs(fx1(1:N/2)))。 subplot(3,1,2)。plot(w/pi,abs(Hd))。 subplot(3,1,3)。plot(w2,abs(flox1(1:N1/2)))。 sound(x1,fs)。 %播放單列音樂信號和濾除高頻后的音樂信號 sound(lox1,fs)； %聲音清晰 單列音樂信號的頻譜圖、巴特沃斯濾波器的頻率響應(yīng)曲線、 濾除高頻后的音樂信號的頻譜圖 濾除低頻信息的程序： clear all。close all。clc [x,fs,nbits]=wavread(39。怒放的生命 汪峰 5s39。)。 x1=x(:,1)。 %獲取單列音樂信號并對其做 FFT變換 N=length(x1)。 fx1=fft(x1)。 w1=2/N*(0:N/21)。 wp=。ws=。rp=1。rs=50。 %設(shè)計巴特沃斯高通濾波器 [N4,Wc]=buttord(wp,ws,rp,rs)。 [B,A]=butter(N4,Wc,39。high39。)。 [Hd,w]=freqz(B,A)。 lox1=filter(B,A,x1)。 %使用巴特沃斯濾波器濾除音樂信號的低頻部分并對所得音樂信號做 FFT變換 N1=length(lox1)。 flox1=fft(lox1)。 w2=2/N1*(0:N1/21)。 figure %畫出單列音樂信號的頻譜圖、巴特沃斯濾波器的頻率響應(yīng)曲 %線和濾除低頻后的音樂信號的頻譜圖 subplot(3,1,1)。plot(w1,abs(fx1(1:N/2)))。 subplot(3,1,2)。plot(w/pi,abs(Hd))。 subplot(3,1,3)。plot(w2,abs(flox1(1:N1/2)))。 sound(lox1,fs)。 %聲音 低，不清晰。 交叉組合音樂 clear all。close all。clc clear all。close all。clc [x,fs,nbits]=wavread(39。鋼琴曲 雨的印記 5s39。)。 [y,fs,nbits]=wavread(39。怒放的生命 汪峰 5s39。)。 x1=x(:,1)。 y1=y(:,1)。 x2=x1(1:202100) %取音樂長度 Nx2=length(x2)。 y2=y1(1:202100)。 Ny2=length(y2)。 x3=fft(x2)。 y3=fft(y2)。 w1=2/Nx2*[0:Nx21]。 w2=2/Ny2*[0:Ny21]。 Fx1=abs(x3)。 %選取第一個音樂信號的幅度和第二個音樂信號 %的相位 Ay1=angle(y3)。 F4=Fx1.*exp(j*Ay1)。 X4=ifft(F4)。 NF4=length(F4)。 F5=fft(F4)。 w3=2/NF4*[0:NF41]。 sound(real(X4),fs)。 figure subplot(3,1,1)。plot(w1,abs(x3))。 subplot(3,1,2)。plot(w2,abs(y3))。 subplot(3,1,3)。plot(w3,abs(F5))。
畢業(yè)設(shè)計相關(guān)推薦
鄂ICP備17016276號-1