音樂(lè)信號(hào)的處理 matlab.doc

《音樂(lè)信號(hào)的處理 matlab.doc》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫(kù)。

1、1、音樂(lè)信號(hào)的音譜和頻譜觀察[y,fs,bit]=wavread('2');sound(y,fs);plot(y);y1=y(:,1);M=512;fh1=fft(y1,M);w=2/M*[0:M-1];figure;plot(w,abs(fh1));2、音樂(lè)信號(hào)的抽取（減抽樣）[y,fs,bit]=wavread('2');sound(y,fs);figure;subplot(1,2,1);plot(y);title('原音樂(lè)信號(hào)波形');y1=y(:,1);M=512;fh1=fft(y1,M);w=2/M*[0:M-1];subplot(1,2,2);p

2、lot(w,fh1);title('原音樂(lè)信號(hào)頻譜');gridon;j=0;d=10;%減抽樣間隔取3和10代表混疊與非混疊y3=y1(1:100000);fori=1:d:length(y3);j=j+1;yd(j)=y1(i);endsound(yd,fs/d);figure;subplot(1,2,1);plot(yd);title('減抽樣時(shí)域');fh2=fft(yd);l1=length(fh2);w2=2/l1*(0:l1-1);subplot(1,2,2);plot(w2,abs(fh2));title('減抽樣頻域');3、音樂(lè)信號(hào)的AM調(diào)

3、制和巴特沃斯IIR濾波器濾波[w,fs,bit]=wavread('2');w1=w(:,1);l=length(w1);ww=2/l*(0:l-1);n=0:length(w1)-1;t1=w1.*cos(0.6*pi*n)';T1=fft(t1);plot(t1);title('AM低/高頻調(diào)制后波形');figure;plot(ww,abs(T1));title('AM低/高頻調(diào)制后頻譜');sound(t1,fs);t2=t1.*cos(0.6*pi*n)';T2=fft(t2);figure;plot(t2);title('AM低/高頻解調(diào)后波形')

4、;figure;plot(ww,abs(T2));title('AM低/高頻解調(diào)后頻譜');sound(t2,fs);[N,wc]=buttord(0.2,0.3,2,15);[b1,a]=butter(N,wc);[h,wi]=freqz(b1,a);figure;plot(wi/pi,abs(h));title('iir濾波器頻譜');%用iir對(duì)信號(hào)濾波y=filter(b1,a,t2);Y=fft(y);figure;plot(y);title('低/高頻解調(diào)濾波后時(shí)域');figure;plot(ww,abs(Y));title('低/高頻解調(diào)濾波后

5、的頻域');sound(y,fs);用窗函數(shù)法設(shè)計(jì)FIR濾波器（矩形窗和布萊克曼窗）[w,fs,bit]=wavread('2');w1=w(:,1);n=0:length(w1)-1;t1=w1.*cos(0.6*pi*n)';j1=t1.*cos(0.6*pi*n)';%用FIR對(duì)信號(hào)同步解調(diào)N=55;wc=pi/4;%boxcarN取9hd=ideal(N,wc);w=blackman(N);%blackman和boxcarh=hd.*w';M=512;fh=fft(h,M);ww=2/M*(0:M-1);plot(ww,abs(fh));title('

6、矩形窗fir濾波器頻率響應(yīng)');y=conv(j1,h);Y=fft(y,M);figureplot(y);title('矩形窗fir濾波后時(shí)域');figureplot(ww,abs(Y));title('矩形窗fir濾波后頻域');sound(y,fs);三余弦和白噪聲的實(shí)驗(yàn)[y,fs,bit]=wavread('2');y1=y(:,1);f=fft(y1);l=length(f);ww=2/l*(0:l-1);n=0:length(y1)-1;x=0.05*(cos(2*pi*3000*n/fs)+cos(2*pi*5000*n/fs)+cos(2*p

7、i*8000*n/fs));fx1=fft(x);plot(ww,abs(fx1));title('三余弦噪聲頻譜');r=rand(size(y1),1)-0.5;fx2=fft(r);figureplot(ww,abs(fx2));title('隨機(jī)白噪聲頻譜');x1=y1'+x;figureplot(x1);title('加三余弦噪聲的信號(hào)音譜');fx3=fft(x1);figureplot(ww,abs(fx3));title('加三余弦噪聲的信號(hào)頻譜');x2=y1+r;figureplot(x2);title('加白噪聲的信號(hào)音譜');fx4=

8、fft(x2);figureplot(