Matlab_中神經(jīng)網(wǎng)絡(luò)算法指令newff的使用

1、Matlab中神經(jīng)網(wǎng)絡(luò)算法指令newff的使用設(shè)P,T是訓練樣本，X,Y是測試樣本；net=newrb(P,T,err_goal,spread);%建立網(wǎng)絡(luò)q=sim(net,p);e=q-T;plot(p,q);%畫訓練誤差曲線q=sim(net,X);e=q-Y;plot(X,q);%畫測試誤差曲線訓練前饋網(wǎng)絡(luò)的第一步是建立網(wǎng)絡(luò)對象。函數(shù)newff建立一個可訓練的前饋網(wǎng)絡(luò)。這需要4個輸入?yún)?shù)。第一個參數(shù)是一個Rx2的矩陣以定義R個輸入向量的最小值和最大值。第二個參數(shù)是一個設(shè)定每層神經(jīng)元個數(shù)的數(shù)組。第三個參數(shù)是包含每層用到的傳遞函數(shù)名稱的細胞數(shù)組。最后一個參數(shù)是用到的訓練函數(shù)的名稱。舉個例子

2、，下面命令將創(chuàng)建一個二層網(wǎng)絡(luò)。它的輸入是兩個元素的向量，第一層有三個神經(jīng)元(3)，第二層有一個神經(jīng)元(1)。第一層的傳遞函數(shù)是tan-sigmoid，輸出層的傳遞函數(shù)是linear。輸入向量的第一個元素的范圍是-1到2-12，輸入向量的第二個元素的范圍是0到505，訓練函數(shù)是traingd。net=newff(-12;05,3,1,tansig,purelin,traingd);這個命令建立了網(wǎng)絡(luò)對象并且初始化了網(wǎng)絡(luò)權(quán)重和偏置，因此網(wǎng)絡(luò)就可以進行訓練了。我們可能要多次重新初始化權(quán)重或者進行自定義的初始化。下面就是初始化的詳細步驟。在訓練前饋網(wǎng)絡(luò)之前，權(quán)重和偏置必須被初始化。初始化權(quán)重和偏置的工

3、作用命令init來實現(xiàn)。這個函數(shù)接收網(wǎng)絡(luò)對象并初始化權(quán)重和偏置后返回網(wǎng)絡(luò)對象。下面就是網(wǎng)絡(luò)如何初始化的：net=init(net);我們可以通過設(shè)定網(wǎng)絡(luò)參數(shù)net.initFcn和net.layeri.initFcn這一技巧來初始化一個給定的網(wǎng)絡(luò)。net.initFcn用來決定整個網(wǎng)絡(luò)的初始化函數(shù)。前饋網(wǎng)絡(luò)的缺省值為initlay,它允許每一層用單獨的初始化函數(shù)。設(shè)定了net.initFcn，那么參數(shù)net.layeri.initFcn也要設(shè)定用來決定每一層的初始化函數(shù)。對前饋網(wǎng)絡(luò)來說，有兩種不同的初始化方式經(jīng)常被用到：initwb和initnw。initwb函數(shù)根據(jù)每一層自己的初始化參數(shù)(n

4、et.inputWeightsi,j.initFcn)初始化權(quán)重矩陣和偏置。前饋網(wǎng)絡(luò)的初始化權(quán)重通常設(shè)為rands，它使權(quán)重在-1到1之間隨機取值。這種方式經(jīng)常用在轉(zhuǎn)換函數(shù)是線性函數(shù)時。initnw通常用于轉(zhuǎn)換函數(shù)是曲線函數(shù)。它根據(jù)Nguyen和WidrowNgWi90為層產(chǎn)生初始權(quán)重和偏置值，使得每層神經(jīng)元的活動區(qū)域能大致平坦的分布在輸入空間。它比起單純的給權(quán)重和偏置隨機賦值有以下優(yōu)點：減少神經(jīng)元的浪費(因為所有神經(jīng)元的活動區(qū)域都在輸入空間內(nèi))。有更快的訓練速度(因為輸入空間的每個區(qū)域都在活動的神經(jīng)元范圍中)。初始化函數(shù)被newff所調(diào)用。因此當網(wǎng)絡(luò)創(chuàng)建時，它根據(jù)缺省的參數(shù)自動初始化。ini

5、tEditedbyFoxitReaderCopyright(C)byFoxitSoftwareCompany,2005-2008ForEvaluationOnly.不需要單獨的調(diào)用?？墒俏覀兛赡芤匦鲁跏蓟瘷?quán)重和偏置或者進行自定義的初始化。例如，我們用newff創(chuàng)建的網(wǎng)絡(luò)，它缺省用initnw來初始化第一層。如果我們想要用rands重新初始化第一層的權(quán)重和偏置，我們用以下命令：net.layers1.initFcn=initwb;net.inputWeights1,1.initFcn=rands;net.biases1,1.initFcn=rands;net.biases2,1.initFcn

6、=rands;net=init(net);IW:輸入層到隱含層的權(quán)重矩陣LW:隱含層和輸出層間的權(quán)重矩陣b:閥值向量如網(wǎng)絡(luò)為net,輸入層和輸出均為一個接點情況下，則用net.IWl,l可以看到第一個輸入接點到第一隱含層的權(quán)重向量；net.LW2,l可以看到隱含層到輸出層的權(quán)值向量；net.b1,1是隱含層的閥值向量，net.b2,1是輸出接點的閥值；在多輸入輸出下先用net.IWnet.LWnet.b查看各矩陣結(jié)構(gòu)，再相應用net.IW?,?等語句查到相關(guān)的向量clearall;%definetheinputandoutputp=974874527;3884661764;1316243922

7、51;183624101860;155723011578;149018772749;151312782026;107015612794;134724153306;132427461233;138314631847;128202347;t=197972428234548;%創(chuàng)建bp網(wǎng)絡(luò)和定義訓練函數(shù)%這里是為了方便而建立一個矩陣，注意是12x2，不是3x2pr=527974;3881764;13162439;18362410;15572301;14902749;12782026;10702794;13473306;12332746;13831847;02347net=newff(pr,15,1,

8、tansigpurelin,trainlm);%這里要加入輸出層的轉(zhuǎn)移函數(shù)，一般是trainlmnet.trainparam.goal=50;net.trainparam.epochs=5000;%訓練神經(jīng)網(wǎng)絡(luò)net,tr=train(net,p,t);%輸出訓練后的權(quán)值和閾值iw1=net.IW1;b1=net.b1;lw2=net.LW2;b2=net.b2;%存儲訓練好的神經(jīng)網(wǎng)絡(luò)savenetkohlernet怎樣知道m(xù)atlab已經(jīng)訓練好的神經(jīng)網(wǎng)絡(luò)的權(quán)值、步長以及闕值用matlab訓練神經(jīng)網(wǎng)絡(luò)時不需輸入權(quán)值、步長以及闕值，如果我想知道m(xù)atlab已經(jīng)訓練好的神經(jīng)網(wǎng)絡(luò)的權(quán)值、步長以及闕值

9、該怎末操作？訓練好的權(quán)值、閾值的輸出方法是：輸入到隱層權(quán)值：w1=net.iw1,1隱層閾值：theta1=net.b1隱層到輸出層權(quán)值：w2=net.lw2,1;輸出層閾值：theta2=net.b2幫我看看matlab的這段程序(有關(guān)神經(jīng)網(wǎng)絡(luò)BP算法)在一位老師的講義上看到的程序，但是有些東西不太明白，請求幫助程序如下：TxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxTxclf;figur

10、e(gcf)echoon%NEWFF建立一個BP網(wǎng)絡(luò)%TRAIN對BP網(wǎng)絡(luò)進行訓練%SIM對BP網(wǎng)絡(luò)進行仿真pauseP=-1:0.1:1;T=-.9602-.5770-.0729.3771.6405.6600.4609.1336-.2013-.4344-.5000-.3930-.1647.0988.3072.3960.3449.1816-.0312-.2189-.3201;plot(P,T,+);title(TrainingVectors);xlabel(InputVectorP);ylabel(TargetVectorT);pausenet=newff(minmax(P),51,tansi

11、gpurelin,traingd,learngd,sse);echooffk=pickic;ifk=2net.iw1,1=3.5000;3.5000;3.5000;3.5000;3.5000;net.b1=-2.8562;1.0774;-0.5880;1.4083;2.8722;net.lw2,1=0.2622-0.2375-0.45250.2361-0.1718;net.b2=0.1326;endnet.iw1,1net.b1net.lw2,1net.b2pauseechoonme=8000;net.trainParam.show=10;net.trainParam.goal=0.02;ne

12、t.trainParam.lr=0.01;A=sim(net,P);sse=sumsqr(T-A);fori=1:me/100ifssenet.trainparam.goal,i=i-1;break,endnet.trainParam.epochs=100;net,tr=train(net,P,T);trp(1+100*(i-1):(max(tr.epoch)+100*(i-1)=tr.perf(1:max(tr.epoch);A=sim(net,P);sse=sumsqr(T-A);plot(P,T,+);holdonplot(P,A)holdoffpauseendmessage=sprintf(Traingd,Epoch%g/%g,SS