%綜合法程序
%讀音頻文件----這部分運行有問題，運行顯示：Index?exceeds?matrix?dimensions（超出矩陣維度？）
clear?all;
close?all;
clc;
[filenamepathname]=uigetfile（‘*.wav‘‘請選擇語音文件:‘）;
[soundfsnbits]=wavread（[pathnamefilename]）;%sound?為采樣數據，fs采樣率。nbits量化比特數
%假設聲音是雙聲道，這里只取單聲道作分析，sound（:1）為左聲道，右聲道為sound（:2）?
a=size（sound（:1））;
Tn=（a（1）-1）/fs;?????
t=0:1/fs:Tn;
figure（1）
subplot（311）;
plot（tsound（:1））xlabel（‘t（s）‘）ylabel（‘sound‘）
title（‘錄制的語音段‘）;?
I=sound（:1）;
%x=I（11401:11700）;%可以根據需要選擇不同的語音段
x=I（9700:10000）;
subplot（312）;
??plot（x）;
??title（‘選取的語音段‘）;?
y=x;?
for?i=1:301????
????y（i）=0;?
end?
y1=y;y2=y;y3=y;y4=y;r=y;rmax=-1;?
p=0;??????????????????%flag?of?pitch?
cyc=0;?

%Filter:?y?
H=[0.0035270585-0.0075853243-0.022130724-0.037701912-0.040792551-0.0176185440.0371342230.11394230.189553710.236577820.236577820.189553710.11394230.037134223-0.017618544-0.040792551-0.037701912-0.022130724-0.00758532420.0035270585];?
for?i=11:291????
????for?j=1:20???????
????????y（i）=y（i）+x（i-j+10）*H（j）;????
????end?
end??
subplot（313）;
plot（y）;
title（‘濾波后信號‘）;

%小波去噪
[cl]=wavedec（y3‘db5‘）;
%用默認閾值對信號進行消噪處理并圖示結果
%用ddencmp函數獲得信號的默認閾值
[thrsorhkeepapp]=ddencmp（‘den‘‘wv‘y）;