資源描述：
《MATLAB程序代碼bp神經(jīng)網(wǎng)絡(luò)應(yīng)用舉例.doc》由會員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、MATLAB程序代碼--bp神經(jīng)網(wǎng)絡(luò)應(yīng)用舉例————————————————————————————————作者：————————————————————————————————日期：MATLAB程序代碼--bp神經(jīng)網(wǎng)絡(luò)應(yīng)用舉例1BP神經(jīng)網(wǎng)絡(luò)的設(shè)計實例例1采用動量梯度下降算法訓(xùn)練BP網(wǎng)絡(luò)。訓(xùn)練樣本定義如下：輸入矢量為p=[-1-231-115-3]目標(biāo)矢量為t=[-1-111]解：本例的MATLAB程序如下：closeallclearechoonclc%NEWFF——生成一個新的前向神經(jīng)網(wǎng)絡(luò)%TRAIN——對BP神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練%SIM——對BP
2、神經(jīng)網(wǎng)絡(luò)進(jìn)行仿真pause%敲任意鍵開始clc%定義訓(xùn)練樣本%P為輸入矢量P=[-1,-2,3,1;-1,1,5,-3];%T為目標(biāo)矢量T=[-1,-1,1,1];pause;clc%創(chuàng)建一個新的前向神經(jīng)網(wǎng)絡(luò)net=newff(minmax(P),[3,1],{'tansig','purelin'},'traingdm')%當(dāng)前輸入層權(quán)值和閾值inputWeights=net.IW{1,1}inputbias=net.b{1}%當(dāng)前網(wǎng)絡(luò)層權(quán)值和閾值layerWeights=net.LW{2,1}layerbias=net.b{2}pauseclc
3、%設(shè)置訓(xùn)練參數(shù)net.trainParam.show=50;%兩次顯示之間的訓(xùn)練次數(shù),缺省值為25net.trainParam.lr=0.05;%學(xué)習(xí)速率net.trainParam.mc=0.9;%動量常數(shù)設(shè)置,缺省就是0.9net.trainParam.epochs=1000;%訓(xùn)練次數(shù),缺省值為100net.trainParam.goal=1e-3;%網(wǎng)絡(luò)性能目標(biāo),缺省值為0pauseclc%調(diào)用TRAINGDM算法訓(xùn)練BP網(wǎng)絡(luò)[net,tr]=train(net,P,T);pauseclc%對BP網(wǎng)絡(luò)進(jìn)行仿真A=sim(net,P)%計算
4、仿真誤差E=T-AMSE=mse(E)pauseclcechooff例2采用貝葉斯正則化算法提高BP網(wǎng)絡(luò)的推廣能力。在本例中，我們采用兩種訓(xùn)練方法，即L-M優(yōu)化算法（trainlm）和貝葉斯正則化算法（trainbr），用以訓(xùn)練BP網(wǎng)絡(luò)，使其能夠擬合某一附加有白噪聲的正弦樣本數(shù)據(jù)。其中，樣本數(shù)據(jù)可以采用如下MATLAB語句生成：輸入矢量：P=[-1:0.05:1]；目標(biāo)矢量：randn(’seed’,78341223)；T=sin(2*pi*P)+0.1*randn(size(P))；解：本例的MATLAB程序如下：closeallclearec
5、hoonclc%NEWFF——生成一個新的前向神經(jīng)網(wǎng)絡(luò)%TRAIN——對BP神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練%SIM——對BP神經(jīng)網(wǎng)絡(luò)進(jìn)行仿真pause%敲任意鍵開始clc%定義訓(xùn)練樣本矢量%P為輸入矢量P=[-1:0.05:1];%T為目標(biāo)矢量randn('seed',78341223);T=sin(2*pi*P)+0.1*randn(size(P));%繪制樣本數(shù)據(jù)點(diǎn)plot(P,T,'+');echooffholdon;plot(P,sin(2*pi*P),':');%繪制不含噪聲的正弦曲線echoonclcpauseclc%創(chuàng)建一個新的前向神經(jīng)網(wǎng)絡(luò)net
6、=newff(minmax(P),[20,1],{'tansig','purelin'});pauseclcechooffclcdisp('1.L-M優(yōu)化算法TRAINLM');disp('2.貝葉斯正則化算法TRAINBR');choice=input('請選擇訓(xùn)練算法(1,2):');figure(gcf);if(choice==1)echoonclc%采用L-M優(yōu)化算法TRAINLMnet.trainFcn='trainlm';pauseclc%設(shè)置訓(xùn)練參數(shù)net.trainParam.epochs=500;net.trainParam.g
7、oal=1e-6;net=init(net);%重新初始化pauseclcelseif(choice==2)echoonclc%采用貝葉斯正則化算法TRAINBRnet.trainFcn='trainbr';pauseclc%設(shè)置訓(xùn)練參數(shù)net.trainParam.epochs=500;randn('seed',192736547);net=init(net);%重新初始化pauseclcend%調(diào)用相應(yīng)算法訓(xùn)練BP網(wǎng)絡(luò)[net,tr]=train(net,P,T);pauseclc%對BP網(wǎng)絡(luò)進(jìn)行仿真A=sim(net,P);%計算仿真誤差E
8、=T-A;MSE=mse(E)pauseclc%繪制匹配結(jié)果曲線closeall;plot(P,A,P,T,'+',P,sin(2*pi
