信號(hào)頻譜分析

聲音信號(hào)的頻譜分析姓名：蘇小平學(xué)號(hào)：1305080116專業(yè)：智能電網(wǎng)信息工程第一部分：goldwave30圖一第二部分：關(guān)于聲音的Matlab分析Matlab30以利于觀察。程序如下：[y,Fs,bits]=wavread('至上勵(lì)合-繁星.wav');y=y(:,1);sigLength=length(y);Y=fft(y,sigLength);Pyy=Y.*conj(Y)/sigLength;halflength=floor(sigLength/2);f=Fs*(0:halflength)/sigLength;figure;plot(f,Pyy(1:halflength+1));xlabel('Frequency(Hz)');t=(0:sigLength-1)/Fs;figure;plot(t,y);xlabel('Time(s)');圖二圖二為處理前的聲道頻域波形程序如下：[y,fs,bits]=wavread('至上勵(lì)合繁星.wav');L=length(y);noise=0.1*randn(L,2);y_z=y+noise;wavwrite(y_z,8000*6,'加噪聲音樂(lè)');n=length(y);y_zp=fft(y_z,n);f=fs*(0:n/2-1)/n;figure(1);subplot;plot(y_z);title('加噪語(yǔ)音信號(hào)時(shí)域波形');xlabel('time')ylabel('幅值A(chǔ)')（2）利用低通濾波器將噪聲濾除程序如下：fp=1500fc=1700As=100Ap=1wc=2*pi*fc/fswp=2*pi*fp/fswdel=wcwpbeta=0.112*(As8.7)N=ceil((As8)/2.285/wdel)wn=kaiser(N+1,beta)ws=(wp+wc)/2/pi;b=fir1(N,ws,wn);freqz(b,1);[y,fs,bits]=wavread('至上勵(lì)合-繁星.wav');L=length(y);noise=0.1*randn(L,2);y_z=y+noise;n=length(y);x=fftfilt(b,y_z);X=fft(x,n);f=fs*(0:n/2-1)/n;figure;gridon;plot(abs(X));title('濾波后信號(hào)的頻譜');gridonwavwrite(x,fs,bits,'濾波后音樂(lè)');濾波器的幅值與相位圖：通過(guò)低通濾波器后有用的聲音的頻譜：第三部分：LinkinMatlab的程序作出時(shí)域波形的頻譜圖。40040035030025020015010050000.511.52Frequency(Hz)2.5x1041.81.8.6.4.200.20.40.60.8-1051015Time(s)202530000----加入隨機(jī)噪聲以后的頻譜圖加噪語(yǔ)音信號(hào)時(shí)域波形1.510.5A值 0幅-0.5-1