數(shù)字信號處理用matlab編程實(shí)現(xiàn)am調(diào)制音樂信號的同步解調(diào)

《數(shù)字信號處理用matlab編程實(shí)現(xiàn)am調(diào)制音樂信號的同步解調(diào)》由會員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在應(yīng)用文檔-天天文庫。

1、clearall;closeall;clc[x,fs,bit]=wavread('F:費(fèi)玉清-一剪梅00_01_23-00_01_28.wav');x=x(:,1);%sound(x,fs);N=length(x);figuresubplot(2,1,1);plot(x)F1=fft(x,N);gridontitle('原始音樂信號時域波型')xlabel('Time')ylabel('Magnitude')w1=2/N*[0:N-1];subplot(2,1,2);plot(w1,abs(F1))gridontitle('原始音樂信號頻域波型')xlabel('Frequenc

2、y/Hz')ylabel('Magnitude')n=0:N-1;x1=cos(n*0.95);figuresubplot(2,2,1);plot(x1)title('余弦信號時域波型')xlabel('Time')ylabel('Magnitude')F2=fft(x1,N);subplot(2,2,2);plot(w1,abs(F2))gridontitle('余弦信號頻域波型')xlabel('Frequency/Hz')ylabel('Magnitude')X=x.*x1';subplot(2,2,3);plot(X)F3=fft(X,N);gridontitle('高頻調(diào)

3、制后音樂信號時域波型')xlabel('Time')ylabel('Magnitude')subplot(2,2,4);plot(w1,abs(F3))gridontitle('高頻調(diào)制后音樂信號頻域波型')xlabel('Frequency/Hz')ylabel('Magnitude')X1=X.*x1';F4=fft(X1,N);figuresubplot(2,1,1);plot(X1);gridontitle('解制后音樂信號時域波型')xlabel('Time')ylabel('Magnitude')subplot(2,1,2);plot(w1,abs(F4));grido

4、ntitle('解制后音樂信號頻域波型')xlabel('Frequency/Hz')ylabel('Magnitude')Rp=1;Rs=15;[N1,Wc]=buttord(0.2,0.3,Rp,Rs);[B,A]=butter(N1,Wc);[H,W]=freqz(B,A);figuresubplot(1,1,1);plot(W/pi,abs(H));gridontitle('數(shù)字巴特沃斯濾波器')xlabel('DigitalFrequency/pi')ylabel('Magnitude')y=filter(B,A,X1);%sound(y,fs);figuresubplo

5、t(2,1,1)plot(y);gridontitle('巴特沃斯濾波后音樂信號時域波型')xlabel('Time')ylabel('Magnitude')y=fft(y,N);subplot(2,1,2)plot(w1,abs(y));gridontitle('濾波后音樂信號頻域波型')xlabel('Frequency/Hz')ylabel('Magnitude')N2=33;wc=0.3*pi;hd=ideal(N2,wc);w2=boxcar(N2);w3=blackman(N2);h1=hd.*w2';h2=hd.*w3';M=512;fh1=fft(h1,M);w2=

6、2/M*[0:M-1];fh2=fft(h2,M);figureplot(w2,abs(fh1))gridontitle('矩形窗濾波器')xlabel('Frequency')ylabel('Magtitude')y1=conv(X1,h1);sound(y1,fs)figuresubplot(2,1,1);plot(y1);gridontitle('矩形窗濾波后的音樂信號時域波形')xlabel('Time')ylabel('Magtitude')y1=fft(y1,M);subplot(2,1,2);plot(w2,abs(y1));gridontitle('矩形窗濾波后的音

7、樂信號頻域波形')xlabel('Frequency')ylabel('Magtitude')figureplot(w2,abs(fh2))gridontitle('布萊克曼窗濾波器')xlabel('Frequency')ylabel('Magtitude')y2=conv(X1,h2);%sound(y2,fs)figuresubplot(2,1,1);plot(y2);gridontitle('布萊克曼窗濾波后的音樂信號時域波形')xlabel('Time')y