这个是完整的代码，我自己的账号发不了这么长，希望好心人帮忙注释啊，非常感谢！！！
创新互联服务项目包括梁平网站建设、梁平网站制作、梁平网页制作以及梁平网络营销策划等。多年来，我们专注于互联网行业，利用自身积累的技术优势、行业经验、深度合作伙伴关系等，向广大中小型企业、政府机构等提供互联网行业的解决方案，梁平网站推广取得了明显的社会效益与经济效益。目前，我们服务的客户以成都为中心已经辐射到梁平省份的部分城市，未来相信会继续扩大服务区域并继续获得客户的支持与信任！
(1)端点检测部分(vad)：
function [x1,x2] = vad(x)
%语音信号x幅度归一化到[-1,1]
x = double(x);
x = x / max(abs(x));
%常数设置
FrameLen = 240; %帧长度为240
FrameInc = 80; %帧移为80
amp1 = 10; %短时能量高门限10
amp2 = 2; %短时能量低门限为2
zcr1 = 10; %短时过零率高门限为10
zcr2 = 5; %短时过零率低门限为5
maxsilence =3;%静音时间门限3*10ms= 30ms
minlen= 15;%最小语音时间长度15*10ms = 150ms
status= 0; %
count= 0; %语音时间累计
silence = 0; %静音时间累计
%计算过零率
tmp1= enframe(x(1:end-1), FrameLen, FrameInc);
tmp2= enframe(x(2:end), FrameLen, FrameInc);
signs =(tmp1.*tmp2)0;%符号数组，用于存储相邻两个采样点符号是否相同，即是否穿越0电平
diffs = (tmp1-tmp2)0.02;%度量相邻两个采样点之间距离，如果大于门限0.02(经验值)，则1，否则0
zcr = sum(signs.*diffs,2);%过零率
%计算短时能量
amp =sum(abs(enframe(filter([1 -0.9375], 1, x), FrameLen, FrameInc)), 2);
%调整能量门限
amp1 = min(amp1,max(amp)/4);
amp2 = min(amp2,max(amp)/8);
%开始端点检测
x1 = 0;
x2 = 0;
for n=1:length(zcr)
goto = 0;
switch status
case {0,1}% 0 =静音, 1 =可能开始
if amp(n) amp1%确信进入语音段
x1 = max(n-count-1,1);
status= 2;
silence = 0;
count= count + 1;
elseif amp(n) amp2 || ... %可能处于语音段
zcr(n) zcr2
status = 1;
count= count + 1;
else%静音状态
status= 0;
count= 0;
end
case 2,% 2 =语音段
if amp(n) amp2 || ...%保持在语音段
zcr(n) zcr2
count = count + 1;
else%语音将结束
silence = silence+1;
if silence maxsilence %静音还不够长，尚未结束
count= count + 1;
elseif count minlen%语音长度太短，认为是噪声
status= 0;
silence = 0;
count= 0;
else%语音结束
status= 3;
end
end
case 3,
break;
end
end
count = count-silence/2;
x2 = x1 + count -1;
subplot(311)
plot(x)
axis([1 length(x) -1 1])
xlabel('语音信号');
line([x1*FrameIncx1*FrameInc ],[-1,1],'Color','red');
line([x2*FrameIncx2*FrameInc ],[-1,1],'Color','red');
subplot(312)
plot(amp);
axis([1 length(amp) 0max(amp)])
xlabel('短时能量');
line([x1,x1],[min(amp),max(amp)],'Color','red');
line([x2,x2],[min(amp),max(amp)],'Color','red');
subplot(313)
plot(zcr);
axis([1 length(zcr) 0max(zcr)])
xlabel('过零率');
line([x1,x1],[min(zcr),max(zcr)],'Color','red');
line([x2,x2],[min(zcr),max(zcr)],'Color','red');
(2)MFCC部分：
function ccc = mfcc(x)
%归一化mel滤波器组系数
bank=melbankm(24,256,8000,0,0.5,'m');%24滤波器个数，8000采样频率
bank=full(bank);
bank=bank/max(bank(:));
% DCT系数,12*24
for k=1:12
n=0:23;
dctcoef(k,:)=cos((2*n+1)*k*pi/(2*24));
end
%归一化倒谱提升窗口
w = 1 + 6 * sin(pi *(1:12) ./ 12);
w = w/max(w);
%预加重滤波器
xx=double(x);
xx=filter([1-0.9375],1,xx);
%语音信号分帧，xx是输入语音信号；256是帧长；80是帧移
xx=enframe(xx,256,80);
%计算每帧的MFCC参数
for i=1:size(xx,1)
y = xx(i,:);
s = y' .* hamming(256);%加汉明窗
t = abs(fft(s));%fft变换
t = t.^2;
c1=dctcoef * log(bank * t(1:129));
c2 = c1.*w';
m(i,:)=c2';
end
%差分系数
dtm = zeros(size(m));
for i=3:size(m,1)-2
dtm(i,:) = -2*m(i-2,:) - m(i-1,:) + m(i+1,:)+ 2*m(i+2,:);
end
dtm = dtm / 3;
%合并mfcc参数和一阶差分mfcc参数
ccc = [m dtm];
%去除首尾两帧，因为这两帧的一阶差分参数为0
ccc =ccc(3:size(m,1)-2,:);
(3)dtw计算部分：
function dist = dtw2(test, ref)
global x y_min y_max
global t r
global D d
global m n
t = test;
r = ref;
n = size(t,1);
m = size(r,1);
d = zeros(m,1);
D =ones(m,1) *realmax;
D(1) = 0;
%如果两个模板长度相差过多，匹配失败
if (2*m-n3) || (2*n-m2)
dist =realmax;
return
end
%计算匹配区域
xa = round((2*m-n)/3);
xb = round((2*n-m)*2/3);
if xbxa
%xbxa,按下面三个区域匹配
%1:xa
%xa+1:xb
%xb+1:N
for x =1:xa
y_max= 2*x;
y_min= round(0.5*x);
warp
end
for x =(xa+1):xb
y_max= round(0.5*(x-n)+m);
y_min= round(0.5*x);
warp
end
for x =(xb+1):n
y_max= round(0.5*(x-n)+m);
y_min= round(2*(x-n)+m);
warp
end
elseif xaxb
%xaxb,按下面三个区域匹配
%0:xb
%xb+1:xa
%xa+1:N
for x =1:xb
y_max= 2*x;
y_min= round(0.5*x);
warp
end
for x =(xb+1):xa
y_max= 2*x;
y_min= round(2*(x-n)+m);
warp
end
for x =(xa+1):n
y_max= round(0.5*(x-n)+m);
y_min= round(2*(x-n)+m);
warp
end
elseif xa==xb
%xa=xb,按下面两个区域匹配
%0:xa
%xa+1:N
for x =1:xa
y_max= 2*x;
y_min= round(0.5*x);
warp
end
for x =(xa+1):n
y_max= round(0.5*(x-n)+m);
y_min= round(2*(x-n)+m);
warp
end
end
%返回匹配分数
dist = D(m);
function warp
global x y_min y_max
global t r
global D d
global m n
d = D;
for y = y_min:y_max
D1 = D(y);
if y1
D2= D(y-1);
else
D2 =realmax;
end
if y2
D3= D(y-2);
else
D3 = realmax;
end
d(y) =sum((t(x,:)-r(y,:)).^2) + min([D1,D2,D3]);
end
D = d;
(4)测试函数testdtw部分；
disp('正在计算参考模板的参数...')
for i=1:10
fname = sprintf('G:\\石东东\\语音\\%da.wav',i-1);
x = wavread(fname);
[x1 x2] = vad(x);
m = mfcc(x);
m = m(x1-2:x2-2,:);
ref(i).mfcc = m;
end
disp('正在计算测试模板的参数...')
for i=1:10
fname = sprintf('G:\\石东东\\语音\\%db.wav',i-1);
x = wavread(fname);
[x1 x2] = vad(x);
m = mfcc(x);
m = m(x1-2:x2-2,:);
test(i).mfcc = m;
end
disp('正在进行模板匹配...')
dist = zeros(10,10);
for i=1:10
for j=1:10
dist(i,j) = dtw2(test(i).mfcc, ref(j).mfcc);
end
end
disp('正在计算匹配结果...')
for i=1:10
[d,j] = min(dist(i,:));
fprintf('测试模板%d的识别结果为：%d\n', i-1, j-1);
end
亚洲
日本：东京NRT成田国际机场；HND 羽田机场；大阪ITM 大阪伊丹国际机场；KIX大阪关西国际机场
韩国：汉城（首尔）ICN 仁川国际机场；釜山PUS釜山金海国际机场
新加坡：新加坡 SIN新加坡樟宜国际机场
马来西亚：吉隆坡 KUL吉隆坡国际机场
越南：胡志明 SGN胡志明市新山机场
泰国：曼谷 BKK素万那普国际机场 普吉 HKT普吉国际机场
印度：新德里 DEL 新德里国际机场
中东：
伊朗：德黑兰 THR德黑兰梅赫拉巴德国际机场
沙特阿拉伯：利雅得RUH哈利德国王国际机场
阿联酋：阿布扎比 AUH 阿布扎比国际机场；迪拜DXB迪拜国际机场
卡塔尔：多哈 DOH多哈国际机场
土耳其：伊斯坦布尔 IST 伊斯坦布尔国际机场
北美
加拿大：渥太华 YOW渥太华国际机场；蒙特利尔 YUL 蒙特利尔多尔瓦国际机场；温哥华 YVR温哥华国际机场；多伦多 YYZ 多伦多皮尔森国际机场
美国：华盛顿 IAD杜勒斯国际机场；JPN五角大楼军用机；波士顿 BOS 波士顿洛根国际机场；芝加哥 ORD 芝加哥奥黑尔国际机场；纽约 JFK肯尼迪国际机场；旧金山 SFO 三番市旧金山国际机场；洛杉矶 LAX洛杉矶国际机场；迈阿密 MIA 迈阿密国际机场；亚特兰大ATL 亚特兰大国际机场；哥伦比亚 CAE 哥伦比亚国际机场；克利夫兰CLE 克利夫兰机场；夏洛特CLT 夏洛特国际机场；丹佛 DEN丹佛国际机场；达拉斯 DFW达拉斯沃思堡机场；底特律 DTW 底特律都会机场；休斯敦IAH 布什国际机场；奥兰多MCO 奥兰多国际机场；孟菲斯MEM 孟菲斯国际机场；西雅图SEA西雅图-塔科马国际机场
中南美
墨西哥：墨西哥城 MEX墨西哥城机场
古巴：哈瓦那 HAV 哈瓦那-何塞马蒂机场
巴西：巴西利亚 BSB巴西利亚国际机场；圣保罗 SAO圣保罗国际机场；里约热内卢 RIO/GIG里约热内卢国际机场
阿根廷：布宜诺斯艾利斯EZE 埃塞萨国际机场
欧洲
英国：伦敦 LHR伦敦希思罗机场；利物浦 LPL利物浦雷侬国际机场；曼彻斯特 MAN曼彻斯特机场
比利时：布鲁塞尔 BRU布鲁塞尔国际机场
卢森堡：卢森堡 LUX卢森堡国际机场
荷兰：阿姆斯特丹 AMS 阿姆斯特丹-史基浦机场；鹿特丹 RTM鹿特丹机场
丹麦：哥本哈根 CPH哥本哈根凯斯楚普机场
德国：柏林 TXL 泰格尔机场；慕尼黑 MUC 慕尼黑机场；不莱梅BRE 不莱梅机场；法兰克福FRA法兰克福-莱茵-美因国际机场；斯图加特 STR斯图加特机场；汉堡 HAM汉堡国际机场；纽伦堡 NUE纽伦堡机场；科隆 CGN科隆机场
法国：巴黎CDG戴高乐机场；马赛 MRS马赛机场；里昂 LYS 里昂机场
瑞士：伯尔尼 BRN贝尔普伯尔尼机场；日内瓦 GVA日内瓦国际机场；苏黎世 ZRH苏黎世国际机场；巴塞尔 BSL巴塞尔机场
西班牙：马德里 MAD马德里巴拉哈斯机场；巴塞罗那 BCN 巴塞罗那安普拉特机场；巴伦西亚VLC 巴伦西亚机场；塞维利亚SVQ塞维利亚机场
葡萄牙：里斯本 LIS 里斯本机场 ；波尔图OPO 奥波多机场
意大利：罗马 FCO罗马菲乌米奇诺机场；米兰 MXP 米兰马尔蓬萨机场；威尼斯 VCE 威尼斯机场；佛罗伦萨 FLR 佛罗伦萨机场；都灵TRN都灵机场
希腊：雅典 ATH雅典国际机场
奥地利：维也纳 VIE维也纳施韦夏特机场
捷克：布拉格 PRG布拉格鲁济涅机场
芬兰：赫尔辛基 HEL赫尔辛基万塔机场
瑞典：斯德哥尔摩 ARN斯德哥尔摩阿兰达机场
挪威：奥斯陆 OSL 奥斯陆加勒穆恩机场
南斯拉夫：贝尔格莱德 BEG贝尔格莱德机场
罗马尼亚：布加勒斯特 BUH 布加勒斯特机场
克罗地亚：萨格勒布 ZAG萨格勒布机场
匈牙利：布达佩斯 BUD布达佩斯费里海吉机场
波兰：华沙 WAW华沙奥肯切机场
俄罗斯：莫斯科 SVO 谢列梅捷沃机场；DME多莫杰多沃机场
乌克兰：基辅KBP基辅机场；IEV 基辅茹良尼机场
非洲
埃及：开罗 CAI开罗国际机场
塞内加尔：达喀尔 DKR达喀尔机场
南非：约翰内斯堡 JNB约翰内斯堡国际机场
大洋洲
澳大利亚：堪培拉 CBR 堪培拉机场；墨尔本 MEL 墨尔本国际机场；悉尼SYD悉尼金斯福国际机场
新西兰：惠灵顿 WLG惠灵顿机场
三字代码 机场中文名称 机场英文名称 州名 缩写
LAX 洛杉矶国际机场 Los Angeles 加利福尼亚 CA
BUR 伯班克机场 Burbank 加利福尼亚 CA
ONT 安大略机场 Ontario 加利福尼亚 CA
LGB 长滩机场 Long Beach 加利福尼亚 CA
SAN 圣地亚哥机场 San Diego 加利福尼亚 CA
SNA 橙县机场 Orange County 加利福尼亚 CA
SFO 旧金山国际机场 San Francisco 加利福尼亚 CA
OAK 奥克兰机场 Oakland 加利福尼亚 CA
SJC 圣荷西机场 San Jose 加利福尼亚 CA
SEA 西雅图国际机场 Seattle 华盛顿 WA
PDX 波特兰机场 Portland 俄勒冈 OR
PHX 凤凰城机场 Phoenix 亚利桑那 AZ
TUS 图森机场 Tueson 亚利桑那 AZ
LAS 拉斯维加斯机场 Las Vegas 内华达 NV
RNO 雷诺机场 Reno 内华达 NV
SLC 盐湖城机场 Salt Lake City 犹他 UT
DEN 丹佛机场 Denver 科罗拉多 CO
ELP 埃尔.帕索机场 El Paso 得克萨斯 TX
SAT 圣.安东尼奥机场 San Antonio 得克萨斯 TX
AUS 奥斯汀机场 Austin 得克萨斯 TX
IAH 休斯顿国际机场 Houston(Int’l) 得克萨斯 TX
HOU 休斯顿哈比机场 Houston(Hobby) 得克萨斯 TX
DFW 达拉斯机场 Dallas 得克萨斯 TX
MSP 明尼阿波利斯机场 Minneapolis 明尼苏达 MN
MCI 堪萨斯城机场 Kansas 密苏里 MO
STL 圣.路易斯机场 St. Louis 密苏里 MO
MSY 新奥尔良机场 New Orleans 路易斯安那 LA
BHM 伯明翰机场 Birmingham 阿拉巴马 AL
MEM 孟菲斯机场 Memphis 田纳西 TN
BNA 纳什维尔机场 Nashville 田纳西 TN
ORD 芝加哥机场 Chicago 依利诺 IL
DTW 底特律机场 Detroit 密执安 MI
IND 印第安那波利斯 Indianapolis 印第安那 IN
CVG 辛辛那提机场 Cincinnati 俄亥俄 OH
DAY 代顿机场 Dayton 俄亥俄 OH
CMH 哥伦布机场 Columbus 俄亥俄 OH
CLE 克利夫兰机场 Cleveland 俄亥俄 OH
PWM 波特兰机场 Portland 缅因 ME
BOS 波士顿机场 Boston 马萨诸塞 MA
ISP 长岛机场 Long Island 纽约 NY
LGA 拉瓜地机场 La Guardia 纽约 NY
JFK 肯尼迪机场 John F Kennedy 纽约 NY
ALB 奥尔巴尼机场 Albany 纽约 NY
BUF 布法罗机场 Buffalo 纽约 NY
EWR 纽瓦克机场 Newark 新泽西 NJ
PIT 匹兹堡机场 Pittsburgh 宾夕法尼亚 PA
PHL 费城机场 Philadelphia 宾夕法尼亚 PA
DCA 华盛顿国际机场 National 哥伦比亚特区 DC
BWI 巴尔地磨机场 Baltimore 马里兰 MD
IAD 华盛顿达拉斯机场 Wa Dulles 弗吉尼亚 VA
CLT 夏洛特机场 Charlotte 北卡罗来纳 NC
ILM 威尔明顿机场 Wilmington 北卡罗来纳 NC
CAE 哥伦比亚机场 Columbia 南卡罗来纳 SC
CHS 查尔斯顿机场 Charleston 南卡罗来纳 SC
ATL 亚特兰大机场 Atlanta 佐治亚 GA
MCO 奥兰多机场 Orlando 佛罗里达 FL
TPA 坦帕机场 Tampa 佛罗里达 FL
MIA 迈阿密机场 Miami 佛罗里达 FL
HNL 火奴鲁鲁机场 Honolulu 夏威夷 HI
OMA 奥马哈机场 Omaha 内布瑞斯卡 NE