華科電信數(shù)字信號(hào)處理實(shí)驗(yàn)報(bào)告.doc

《華科電信數(shù)字信號(hào)處理實(shí)驗(yàn)報(bào)告.doc》由會(huì)員分享，可在線閱讀，更多相關(guān)《華科電信數(shù)字信號(hào)處理實(shí)驗(yàn)報(bào)告.doc（5頁珍藏版）》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。

《數(shù)字信號(hào)處理》實(shí)驗(yàn)報(bào)告 姓 名:_ ______ 學(xué) 號(hào):_ _____ 專業(yè)班級(jí):_ ______ 聯(lián)系電話:_ ______ 2016年11月 24 日 一．實(shí)驗(yàn)?zāi)康? 1、掌握使用 FFT 進(jìn)行信號(hào)譜分析的方法 2、設(shè)計(jì)數(shù)字濾波器對(duì)指定的語音信號(hào)進(jìn)行濾波處理 二．實(shí)驗(yàn)內(nèi)容 1、使用 Matlab 的 fft 函數(shù)對(duì)語音信號(hào)進(jìn)行頻譜分析,找出干擾信號(hào)的頻譜; 2、設(shè)計(jì)數(shù)字濾波器濾除語音信號(hào)中的干擾分量,并進(jìn)行播放對(duì)比。 三．實(shí)驗(yàn)環(huán)境 Window xp Matlab 2008b 四．實(shí)驗(yàn)內(nèi)容 1. 觀察分析原始語音信號(hào)特征： 原始信號(hào)時(shí)域圖：（橫坐標(biāo)為以1/fs為單位的時(shí)間，縱坐標(biāo)為信號(hào)幅值） 原始信號(hào)頻域圖：（橫坐標(biāo)為頻率，縱坐標(biāo)為信號(hào)fft幅值） 分析：從頻域圖可以看出噪聲在0Hz，1575Hz，3150Hz，4725Hz。用hn=[1,A,1]進(jìn)行濾波，|H(ejw)|=|A+2cosw|，在噪聲頻率處令|H(ejw)|=0算出A，用hn和原始信號(hào)卷積即可濾除噪聲。 2.時(shí)域和頻域?yàn)V波效果 濾波后信號(hào)時(shí)域圖：（橫坐標(biāo)為以1/fs為單位的時(shí)間，縱坐標(biāo)為信號(hào)幅值） 濾波后信號(hào)頻域圖：（橫坐標(biāo)為頻率，縱坐標(biāo)為信號(hào)fft幅值） 五．實(shí)驗(yàn)結(jié)果分析 從濾波后信號(hào)時(shí)域圖可以看出原來幅值反常噪聲已經(jīng)濾除，濾波后信號(hào)頻域圖原來幅值很大的8個(gè)沖激已經(jīng)被濾除了。 六．MATLAB源代碼 clear; %讀取聲音信號(hào) [xn,fs] = wavread(SunshineSquare.wav); n = 0:length(xn)-1; f = n*fs/length(xn); XK = fft(xn); figure(1); subplot(2,2,1);plot(xn);title(原始信號(hào)時(shí)域圖);ylabel(xn); subplot(2,2,2);plot(f, abs(XK));title(原始信號(hào)頻域圖);xlabel(f/Hz);ylabel(XK); %4次濾波 yn = band(xn, 0, fs); yn = band(yn, 1575, fs); yn = band(yn, 3150, fs); yn = band(yn, 4725, fs); yn(85440:85455) = xn(1:16); yn(118500:length(yn))=xn(1:length(yn)-118500+1); n1 = 0:length(yn)-1; f1 = n1*fs/length(yn); YK = fft(yn); subplot(2,2,3);plot(yn);title(濾波后信號(hào)時(shí)域圖);ylabel(yn); subplot(2,2,4);plot(f1, abs(YK));title(濾波后信號(hào)頻域圖);xlabel(f/Hz);ylabel(YK); wavwrite(yn,fs,SunshineSquare_New.wav); function [ y ] = band(xn, f, fs) %設(shè)計(jì)濾波器 w = 2*pi*f/fs; A = -2*cos(w); hn = [1, A, 1]; y = conv(xn, hn); end