%%?該代碼為基于帶動(dòng)量項(xiàng)的BP神經(jīng)網(wǎng)絡(luò)語(yǔ)音識(shí)別
%
%?
%?			該案例作者申明：		1：本人長(zhǎng)期駐扎在此板塊里，對(duì)該案例提問(wèn)，做到有問(wèn)必答。本套書(shū)籍官方網(wǎng)站為：video.ourmatlab.com		2：點(diǎn)此從當(dāng)當(dāng)預(yù)定本書(shū)：《Matlab神經(jīng)網(wǎng)絡(luò)30個(gè)案例分析》。	3：此案例有配套的教學(xué)視頻，視頻下載方式video.ourmatlab.com/vbuy.html。?
								4：此案例為原創(chuàng)案例，轉(zhuǎn)載請(qǐng)注明出處（《Matlab神經(jīng)網(wǎng)絡(luò)30個(gè)案例分析》）。							5：若此案例碰巧與您的研究有關(guān)聯(lián)，我們歡迎您提意見(jiàn)，要求等，我們考慮后可以加在案例里。			
%?
%%?清空環(huán)境變量
clc
clear

%%?訓(xùn)練數(shù)據(jù)預(yù)測(cè)數(shù)據(jù)提取及歸一化

%下載四類(lèi)語(yǔ)音信號(hào)
load?data1?c1
load?data2?c2
load?data3?c3
load?data4?c4

%四個(gè)特征信號(hào)矩陣合成一個(gè)矩陣
data（1:500:）=c1（1:500:）;
data（501:1000:）=c2（1:500:）;
data（1001:1500:）=c3（1:500:）;
data（1501:2000:）=c4（1:500:）;

%從1到2000間隨機(jī)排序
k=rand（12000）;
[mn]=sort（k）;

%輸入輸出數(shù)據(jù)
input=data（:2:25）;
output1?=data（:1）;

%把輸出從1維變成4維
output=zeros（20004）;
for?i=1:2000
????switch?output1（i）
????????case?1
????????????output（i:）=[1?0?0?0];
????????case?2
????????????output（i:）=[0?1?0?0];
????????case?3
????????????output（i:）=[0?0?1?0];
????????case?4
????????????output（i:）=[0?0?0?1];
????end
end

%隨機(jī)提取1500個(gè)樣本為訓(xùn)練樣本，500個(gè)樣本為預(yù)測(cè)樣本
input_train=input（n（1:1500）:）‘;
output_train=output（n（1:1500）:）‘;
input_test=input（n（1501:2000）:）‘;
output_test=output（n（1501:2000）:）‘;

%輸入數(shù)據(jù)歸一化
[inputninputps]=mapminmax（input_train）;

%%?網(wǎng)絡(luò)結(jié)構(gòu)初始化
innum=24;
midnum=25;
outnum=4;
?
%權(quán)值初始化
w1=rands（midnuminnum）;
b1=rands（midnum1）;
w2=rands（midnumoutnum）;
b2=rands（outnum1）;

w2_1=w2;w2_2=w2_1;
w1_1=w1;w1_2=w1_1;
b1_1=b1;b1_2=b1_1;
b2_1=b2;b2_2=b2_1;

%學(xué)習(xí)率
xite=0.1;
alfa=0.01;
loopNumber=10;
I=zeros（1midnum）;
Iout=zeros（1midnum）;
FI=zeros（1midnum）;
dw1=zeros（innummidnum）;
db1=zeros（1midnum）;

%%?網(wǎng)絡(luò)訓(xùn)練
E=zeros（1loopNumber）;
for?ii=1:10
????E（ii）=0;
????for?i=1:1:1500
???????%%?網(wǎng)絡(luò)預(yù)測(cè)輸出?
????????x=inputn（:i）;
????????%?隱含層輸出
????????for?j=1:1:midnum
????????????I（j）=inputn（:i）‘*w1（j:）‘+b1（j）;
????????????Iout（j）=1/（1+exp（-I（j）））;
????????end
????????%?輸出層輸出
????????yn=w2‘*Iout‘+b2;
????????
???????%%?權(quán)值閥值修正
????????%計(jì)算誤差
????????e=output_train（:i）-yn;?????
????????E（ii）=E（ii）+sum（abs（e））;
????????
????????%計(jì)算權(quán)值變化率
????????dw2=e*Iout;
????????db2=e‘;
????????
????????for?j=1:1:midnum
????????????S=1/（1+exp（-I（j）））;
????????????FI（j）=S*（1-S）;
????????end??????
????????for?k=1:1:innum
????????????for?j

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----
?????目錄???????????0??2018-10-25?21:00??BP分類(lèi)\
?????文件????????6159??2013-08-21?07:25??BP分類(lèi)\BPDLX.m
?????文件????????3896??2013-08-21?07:15??BP分類(lèi)\chapter1_1.asv
?????文件????????4030??2013-08-21?07:18??BP分類(lèi)\chapter1_1.m
?????文件???????93015??2009-08-29?16:54??BP分類(lèi)\data1.mat
?????文件???????92845??2009-08-29?16:54??BP分類(lèi)\data2.mat
?????文件???????92937??2009-08-29?16:54??BP分類(lèi)\data3.mat
?????文件???????93438??2009-08-29?16:54??BP分類(lèi)\data4.mat