【正文】
%取樣頻率 fs2=8000。%播放頻率 duration=5。%錄音時間fprintf(39。Press any key to start %g seconds of recording...n39。,duration)。pause。fprintf(39。Recording...n39。)。y=wavrecord(duration*fs1,fs1)。%duration*fs 是總的采樣點數(shù)fprintf(39。Finished 39。)。fprintf(39。Press any key to play the recording...n39。)。pause。wavplay(y,fs2)。wavwrite(y,fs1,39。E:39。)?,F(xiàn)象:第二次播放時,聲音明顯失真。理由:采樣頻率和播放頻率不一樣時聲音信號會失真。(2)fs1=16000。%取樣頻率 fs2=16000。%播放頻率 duration=5。%錄音時間fprintf(39。Press any key to start %g seconds of recording...n39。,duration)。pause。fprintf(39。Recording...n39。)。y=wavrecord(duration*fs1,fs1)。%duration*fs 是總的采樣點數(shù)fprintf(39。Finished 39。)。fprintf(39。Press any key to play the recording...n39。)。pause。wavplay(y,fs2)。wavwrite(y,fs1,39。E:39。)。wav=wavread(39。E:39。)。Fs=16000。n=length(wav)。f=(0:n1)*16000/n。mag=abs(fft(wav))。subplot(2,1,1)。plot(wav)。subplot(2,1,2)。plot(f,mag)***400060008000x/f***16000fs1=8000。%取樣頻率 fs2=8000。%播放頻率 duration=5。%錄音時間fprintf(39。Press any key to start %g seconds of recording...n39。,duration)。pause。fprintf(39。Recording...n39。)。y=wavrecord(duration*fs1,fs1)。%duration*fs 是總的采樣點數(shù)fprintf(39。Finished 39。)。fprintf(39。Press any key to play the recording...n39。)。pause。wavplay(y,fs2)。wavwrite(y,fs1,39。E:39。)。wav=wavread(39。E:39。)。Fs=8000。n=length(wav)。f=(0:n1)*16000/n。mag=abs(fft(wav))。subplot(2,1,1)。plot(wav)。subplot(2,1,2)。plot(f,mag)wavplay(wav,8000)采樣頻率為8000hz1 ***0500y/幅度***00x/hz***16000由圖可知:語音信號的主要頻譜成分所在的帶寬為(0—1200hz),帶寬為1200hz。當(dāng)采樣頻率較小時,頻譜圖上顯示帶寬較大,波形較稀松。最低采樣頻率應(yīng)為,聲音信號的最高頻率的兩倍,由圖可知為2400hz。(3)女聲:wav1=wavread(39。E:39。)。wav2=wavread(39。E:39。)。wav=wav1+wav2。fp1=800。fp2=1500。fp=[fp1,fp2]。fr1=650。fr2=1900。fr=[fr1,fr2]。Fs=16000。ap=1。as=40。[n,fn]= buttord(fp/(Fs/2),fr/(Fs/2),ap,as,39。z39。)。[b,a]=butter(n,fn)。Y1=filter(b,a,wav)。Y=fft(Y1)。mag=abs(Y)。n=length(wav)。f=(0:n1)*16000/n。subplot(3,1,1)。mag1=abs(fft(wav))。plot(f,mag1)subplot(3,1,2)。plot(f,mag)。subplot(3,1,3)。plot(Y1)。wavplay(Y1,16000)混合400300y/幅度***060008000x/頻率女聲***16000400300y/幅度 2001000 ***00f/hz***男聲:wav1=wavread(39。E:39。)。wav2=wavread(39。E:39。)。wav=wav1+wav2。fp1=200。fp2=600。fp=[fp1,fp2]。fr1=100。fr2=1000。fr=[fr1,fr2]。Fs=16000。ap=3。as=40。[n,fn]= buttord(fp/(Fs/2),fr/(Fs/2),ap,as,39。z39。)。[b,a]=butter(n,fn)。Y1=filter(b,a,wav)。Y=fft(Y1)。mag=abs(Y)。n=length(wav)。f=(0:n1)*16000/n。subplot(3,1,1)。mag1=abs(fft(wav))。plot(f,mag1)subplot(3,1,2)。plot(f,mag)。subplot(3,1,3)。plot(Y1)。wavplay(Y1,16000)混合頻譜圖200150y/幅度***30004000x/hz男聲頻譜圖***040y/幅度***8000x/hz男聲時域圖***分離效果不佳,原因:男女聲頻率有很多重疊的地方。