[xfsbits]=wavread（‘wb.wav‘）;
N=length（x）;%測定語音信號長度
f=2500;
x=x‘;%對語音信號進(jìn)行轉(zhuǎn)置，使N行1列矩陣變?yōu)?行N列矩陣
n=0:1/fs:（N-1）/fs;
s=0.05*sin（2*pi*f*n）;%噪聲信號
y=x+s;%信號加噪聲
f1=0:fs/N:fs*（N-1）/N;%橫坐標(biāo)，f1為頻率
x1=fft（xN）;
s1=fft（sN）;
Y=fft（yN）;
figure（1）;
subplot（311）;plot（nx）;title（‘原始信號波形圖‘）;
xlabel（‘t/s‘）;ylabel（‘幅值‘）;
subplot（312）;plot（ns）;title（‘噪聲信號波形圖‘）;ylim（[-0.080.08]）;
xlabel（‘t/s‘）;ylabel（‘幅值‘）;
subplot（313）;plot（ny）;title（‘信號+噪聲信號波形圖‘）;
xlabel（‘t/s‘）;ylabel（‘幅值‘）;
figure（2）;
subplot（221）;plot（f1abs（x1））;title（‘原始信號頻譜圖‘）;
xlabel（‘頻率/Hz‘）;??ylabel（‘頻率響應(yīng)幅度‘）;xlim（[0fs/12]）;
subplot（222）;plot（f1abs（s1））;title（‘噪聲信號頻譜圖‘）;
xlabel（‘頻率/Hz‘）;??ylabel（‘頻率響應(yīng)幅度‘）;xlim（[0fs/12]）;
subplot（223）;plot（f1abs（Y））;title（‘信號+噪聲頻譜圖‘）;
xlabel（‘頻率/Hz‘）;??ylabel（‘頻率響應(yīng)幅度‘）;xlim（[0fs/12]）;



%%%%%%%%%%%%%%濾波器設(shè)計%%%%%%%%%%%%%%%%%%%%
%%根據(jù)原始信號與噪聲信號頻譜圖，選擇巴特沃斯低通濾波器
%%濾波器性能指標(biāo)：取wp=1500*2*pi?rad/sws=2400*2*pi?rad/sRp=2Rs=15
%選擇濾波器的最小階數(shù)
wp=800*2*pi;ws=900*2*pi;Rp=2;Rs=15;
[N1Wn]=buttord（wpwsRpRs‘s‘）;
%創(chuàng)建butterworth模擬濾波器
[ZPK]=buttap（N1）;
%把濾波器的零極點(diǎn)模型轉(zhuǎn)化為傳遞函數(shù)模型
[BapAap]=zp2tf（ZPK）;
%把模擬濾波器原型轉(zhuǎn)換成截止頻率為Wn的低通濾波器
[ba]=lp2lp（BapAapWn）;
%用雙線性變換法實(shí)現(xiàn)模擬濾波器到數(shù)字濾波器的轉(zhuǎn)換
[bzaz]=bilinear（bafs）;
%繪制頻率響應(yīng)曲線
[HW]=freqz（bzaz）;
figure（3）;
subplot（211）
plot（W*fs/（2*pi）abs（H））;??grid??
title（‘butterworth低通濾波器頻響圖‘）;xlabel（‘頻率/Hz‘）;??ylabel（‘頻率響應(yīng)幅度‘）;



%%%%%%%用所設(shè)計的濾波器對語音+噪聲信號y濾波%%%%%%%
y1=filter（bzazy）;
subplot（212）plot（ny1）;title（‘濾波后信號波形圖‘）;
xlabel（‘頻率/Hz‘）;??ylabel（‘頻率響應(yīng)幅度‘）;
%?sound（xfs）;%播放語音信號
%?sound（yfs）;%播放語音+噪聲信號
sound（y1fs）;%播放濾波后的信號


%%%%%%%語音信號與噪聲信號卷積%%%%%%%%%%%%%%%%%%%
ww=conv（xs）;
figure（4）;
subplot（311）;plot（nx）;title（‘原始信號波形圖‘）;
xlabel（‘t/s‘）;ylabel（‘幅值‘）;
subplot（312）;plot（ns）;title（‘噪聲信號波形圖‘）;ylim（[-0.080.08]）;
xlabel（‘t/s‘）;ylabel（‘幅值‘）;
subplot（313）;plot（ww）;title（‘信號與噪聲信號卷積波形圖‘）;
xlabel（‘t/s‘）;ylabel（‘幅值‘）;

figure（5）;
subplot（221）;plot（f1abs（x1））;title（‘原始信號頻譜圖‘）;
xlabel（‘頻率/Hz‘）;??ylabel（‘頻率響應(yīng)幅度‘）;xlim（[0fs/12]）;
subplot（222）;plot（f1abs（s1））;title（‘噪聲信號頻譜圖‘）;
xlabel（‘頻率/Hz‘）;??ylabel（‘頻率響應(yīng)幅度‘）;xlim（[0fs/12]）;
subplot（223）;plot（abs（ww））;title（‘信號噪聲卷積頻譜圖‘）;
xlabel（‘頻率/Hz‘）;??ylabel（‘頻率響應(yīng)幅度‘）;xlim（[0fs/12]）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????35840??2014-06-27?09:03??語音信號的采集與處理\1?任務(wù)書及任務(wù)安排\語音信號的采集與處理\1?語音信號的采集與處理.doc

?????文件??????12452??2014-06-09?23:33??語音信號的采集與處理\1?任務(wù)書及任務(wù)安排\語音信號的采集與處理\課程設(shè)計安排.png

?????文件??????35328??2014-06-26?14:22??語音信號的采集與處理\1?任務(wù)書及任務(wù)安排\語音信號的采集與處理\課程設(shè)計報告.doc

?????文件?????423424??2014-06-15?10:05??語音信號的采集與處理\2?撰寫報告\【報告資料】語音信號采集與處理.doc

?????文件?????792192??2014-06-27?08:33??語音信號的采集與處理\2?撰寫報告\報告內(nèi)容?【初】.doc

?????文件?????118784??2014-06-19?21:00??語音信號的采集與處理\2?撰寫報告\報告內(nèi)容.doc

?????文件??????31744??2014-06-27?10:46??語音信號的采集與處理\2?撰寫報告\課程設(shè)計實(shí)驗報告\呂勝\1?語音信號的采集與處理.doc

?????文件?????735872??2014-06-27?11:04??語音信號的采集與處理\2?撰寫報告\課程設(shè)計實(shí)驗報告\呂勝\報告內(nèi)容?【初】.doc

?????文件??????31744??2014-06-27?10:46??語音信號的采集與處理\2?撰寫報告\課程設(shè)計實(shí)驗報告\王斌1\1?語音信號的采集與處理.doc

?????文件?????738432??2014-06-27?11:02??語音信號的采集與處理\2?撰寫報告\課程設(shè)計實(shí)驗報告\王斌1\報告內(nèi)容?【初】.doc

?????文件??????31744??2014-06-27?10:48??語音信號的采集與處理\2?撰寫報告\課程設(shè)計實(shí)驗報告\高強(qiáng)\1?語音信號的采集與處理.doc

?????文件?????738944??2014-06-27?11:05??語音信號的采集與處理\2?撰寫報告\課程設(shè)計實(shí)驗報告\高強(qiáng)\報告內(nèi)容?【初】.doc

?????文件?????209629??2014-06-26?23:12??語音信號的采集與處理\3?課程設(shè)計已經(jīng)完成內(nèi)容\PCI總線.docx

?????文件?????365568??2014-06-18?10:34??語音信號的采集與處理\3?課程設(shè)計已經(jīng)完成內(nèi)容\usb總線及pci總線.doc

?????文件?????160886??2010-12-02?17:11??語音信號的采集與處理\3?課程設(shè)計已經(jīng)完成內(nèi)容\【labview】讀取語音信號保存\Record.vi

?????文件??????57089??2010-12-02?14:38??語音信號的采集與處理\3?課程設(shè)計已經(jīng)完成內(nèi)容\【labview】讀取語音信號保存\WriteWave.vi

?????文件???????2626??2014-06-20?12:29??語音信號的采集與處理\3?課程設(shè)計已經(jīng)完成內(nèi)容\【matlab】語音信號加噪濾波\ks.m

?????文件??????48364??2014-06-15?20:49??語音信號的采集與處理\3?課程設(shè)計已經(jīng)完成內(nèi)容\【matlab】語音信號加噪濾波\wb.wav

?????文件??????43972??2014-06-19?09:19??語音信號的采集與處理\3?課程設(shè)計已經(jīng)完成內(nèi)容\【matlab】語音信號加噪濾波\加噪聲以及濾波.docx

?????文件?????193536??2014-06-26?10:39??語音信號的采集與處理\4?參考文獻(xiàn)資料\labview程序圖.doc

?????文件??????91136??2014-06-18?23:14??語音信號的采集與處理\4?參考文獻(xiàn)資料\matlab處理音頻信號.doc

?????文件?????108032??2014-06-18?09:39??語音信號的采集與處理\4?參考文獻(xiàn)資料\MATLAB用于語音信號的處理.doc

?????文件????3438491??2014-06-18?01:12??語音信號的采集與處理\4?參考文獻(xiàn)資料\Matlab語音信號處理.pdf

?????文件????3286234??2014-06-18?01:11??語音信號的采集與處理\4?參考文獻(xiàn)資料\Matlab課程設(shè)計--基于MATLAB有噪聲語音信號處理.pdf

?????文件????1394176??2014-06-09?23:52??語音信號的采集與處理\4?參考文獻(xiàn)資料\PCI總線及USB總線.ppt

?????文件??????10752??2014-06-26?23:15??語音信號的采集與處理\4?參考文獻(xiàn)資料\y.doc

?????文件?????242176??2014-06-27?01:32??語音信號的采集與處理\4?參考文獻(xiàn)資料\_232fab14a11b9c96ee9bf578a7a2ce37.doc

?????文件?????437389??2014-06-20?13:36??語音信號的采集與處理\4?參考文獻(xiàn)資料\信號與系統(tǒng)及MATLAB實(shí)現(xiàn).doc

?????文件??????65536??2014-06-20?13:09??語音信號的采集與處理\4?參考文獻(xiàn)資料\卷積的意義.doc

?????文件?????197660??2014-06-19?09:12??語音信號的采集與處理\4?參考文獻(xiàn)資料\可直接感受的基于MATLAB的語音濾波.pdf

............此處省略52個文件信息