【正文】 進(jìn)行調(diào)制figureplot(x1)。 w=*pi 時(shí).w=*pi 時(shí).w=*pi 時(shí).1， 由原信號(hào)頻譜知。 %播放yuan的信號(hào)N=length(a1)。幅度39。)。)。 %播放矩形窗濾波后信號(hào)sheng yin%sound(y4,fs)。figure。)。)。ylabel(39。plot(w,abs(Fy2))。size(a)。疊加余弦信號(hào)頻譜39。)。fy=fft(y3,N1)。figure。subplot(2,1,1),plot(w,abs(fa))。)。%sound(a1,fs)。原信號(hào)1頻譜39。)。size(x)。Fa1=abs(F1)。%sound(real(b1),fs)。原信號(hào)2幅度譜，原信號(hào)1相位譜39。)。 %濾波后頻譜%sound(a3,fs)。 %濾波后頻譜%sound(a2,fs)。 %濾波后信號(hào)fd=fft(a2,N1)。)。)。 %濾波后信號(hào)Fz=fft(r,N1)。)。數(shù)字巴特沃斯濾波器39。 %疊加余弦音樂信號(hào)F2=fft(y2,N1)。%疊加余弦信號(hào)F1=fft(x4,N1)。close all。title(39。)。stem(h2)。subplot(3,1,2),plot(w,abs(fh1))。db1=20*log10(abs(fh1(1)./(abs(fh1)+eps)))。h2=hd.*w239。fa=fft(a1,N1)。f39。 %調(diào)制信號(hào)頻譜xlabel(39。size(a)。ylabel(39。ylabel(39。f39。size(a)。減抽樣頻譜39。title(39。3， *pi 實(shí)驗(yàn)2 音樂信號(hào)的抽?。p抽樣）觀察音樂信號(hào)頻率上限，選擇適當(dāng)?shù)某槿￠g隔對(duì)信號(hào)進(jìn)行減抽樣（給出兩種抽取間隔，代表混疊和非混疊）；輸出減抽樣音樂信號(hào)的波形和頻譜，觀察現(xiàn)象，給出理論解釋；播放減抽樣音樂信號(hào)，注意抽樣率的改變，比較不同抽取間隔下的聲音，解釋現(xiàn)象。 %頻譜圖橫坐標(biāo)設(shè)置figure。clear all。N1=length(a1)。 %時(shí)域波形。t=0:1/N3:1/N3*(N31)。figure。clc[a,fs,bit]=wavread(39。cos(wt)頻譜39。 %原信號(hào)頻譜xlabel(39。 %調(diào)制信號(hào)頻譜xlabel(39。 endend[a,fs,bit]=wavread(39。title(39。解調(diào)頻譜39。 %數(shù)字巴特沃斯濾波器x3=filter(B,A,x2)。 w1=boxcar(N)。布萊克曼窗濾波后信號(hào)39。ylabel(39。 %矩形窗函數(shù)分貝圖grid on。grid on。figure。5， 使用布萊克曼窗濾波，布萊克曼窗過渡帶寬，但是阻帶較好。x4=x1+x2+x3。%sound(y2,fs)。[H,W]=freqz(B,A)。 %疊加音樂信號(hào)頻譜subplot(3,1,3),plot(w,abs(fy))。數(shù)字巴特沃斯濾波器239。subplot(2,1,1),plot(w,abs(Fr))。close all。數(shù)字巴特沃斯濾波器39。title(39。title(39。)。)。2， 高通濾波器濾掉低頻信息，聲音聽起來比較尖銳、不渾厚，低頻影響信號(hào)整體部分。Fx2=angle(F2)。%sound(x2,fs)。title(39。f39。 %原信號(hào)1傅立葉變換w=2/N1*[0:N11]。加噪聲后音樂信號(hào)能聽到原有的音調(diào)，但里面有非常大的雜音，茲茲的噪聲。)。 %播放疊加隨機(jī)噪聲音樂信號(hào)[N,Wc]=buttord(,1,50)。yuan xin hao pin pu39。)。f39。%sound(a1,fs)。)。subplot(3,1,3)。)。)。title(39。title(39。數(shù)字巴特沃斯濾波器濾波頻譜39。[B,A]=butter(N,Wc)。 %播放調(diào)制后