bp神經(jīng)網(wǎng)絡(luò)詳細(xì)步驟C#實(shí)現(xiàn)

1、using System;using System.Collections.Generic;using System.Linq;using System.Text;using System;using System.IO;using System.Text;namespace BpANNet/ <summary>/ BpNet 的摘要說(shuō)明。/ </summary>public class BpNetpublic int inNum;/輸入節(jié)點(diǎn)數(shù)int hideNum;/隱層節(jié)點(diǎn)數(shù)public int outNum;/輸出層節(jié)點(diǎn)數(shù)public int sampleNum;

2、/樣本總數(shù)Random R;double x;/輸入節(jié)點(diǎn)的輸入數(shù)據(jù)double x1;/隱層節(jié)點(diǎn)的輸出double x2;/輸出節(jié)點(diǎn)的輸出double o1;/隱層的輸入double o2;/輸出層的輸入public double , w;/權(quán)值矩陣w，這是輸入層與隱藏層之間的權(quán)值矩陣public double , v;/權(quán)值矩陣V，這是隱藏層與輸出層之間的權(quán)值矩陣public double , dw;/權(quán)值矩陣wpublic double , dv;/權(quán)值矩陣Vpublic double rate;/學(xué)習(xí)率public double b1;/隱層閾值矩陣public double b2;/輸

3、出層閾值矩陣public double db1;/隱層閾值矩陣public double db2;/輸出層閾值矩陣double pp;/隱藏層的誤差double qq;/輸出層的誤差double yd;/輸出層的教師數(shù)據(jù)，所謂教師數(shù)據(jù)就是實(shí)際數(shù)據(jù)而已！public double e;/均方誤差double in_rate;/歸一化比例系數(shù)/用于確定隱藏層的神經(jīng)細(xì)胞數(shù)public int computeHideNum(int m,int n)double s=Math.Sqrt(0.43*m*n+0.12*n*n+2.54*m+0.77*n+0.35)+0.51;int ss=Convert.T

4、oInt32(s);return (s-(double)ss)>0.5) ? ss+1:ss;public BpNet(double , p,double , t)/ 構(gòu)造函數(shù)邏輯R=new Random();this.inNum=p.GetLength(1);this.outNum=t.GetLength(1);this.hideNum=computeHideNum(inNum,outNum);/this.hideNum=18;this.sampleNum=p.GetLength(0);Console.WriteLine("輸入節(jié)點(diǎn)數(shù)目： "+inNum);Cons

5、ole.WriteLine("隱層節(jié)點(diǎn)數(shù)目："+hideNum);Console.WriteLine("輸出層節(jié)點(diǎn)數(shù)目："+outNum);Console.ReadLine();/將這些矩陣規(guī)定好矩陣大小x=new doubleinNum;x1=new doublehideNum;x2=new doubleoutNum;o1=new doublehideNum;o2=new doubleoutNum;w = new doubleinNum, hideNum;/權(quán)值矩陣w，這是輸入層與隱藏層之間的權(quán)值矩陣v=new doublehideNum,outNum

6、;dw=new doubleinNum,hideNum;dv=new doublehideNum,outNum;/閾值b1=new doublehideNum;b2=new doubleoutNum;db1=new doublehideNum;db2=new doubleoutNum;/誤差pp = new doublehideNum;/隱藏層的誤差qq = new doubleoutNum;/輸出層的誤差yd = new doubleoutNum;/輸出層的教師數(shù)據(jù)/初始化wfor(int i=0;i<inNum;i+)for(int j=0;j<hideNum;j+)/Next

7、Double返回一個(gè)介于 0.0 和 1.0 之間的隨機(jī)數(shù)。wi,j=(R.NextDouble()*2-1.0)/2;/初始化vfor(int i=0;i<hideNum;i+)for(int j=0;j<outNum;j+)vi,j=(R.NextDouble()*2-1.0)/2;rate=0.8;e=0.0;in_rate=1.0; /訓(xùn)練函數(shù)public void train(double , p,double , t)e=0.0;/求p，t中的最大值double pMax=0.0;/sampleNum為樣本總數(shù)for(int isamp=0;isamp<

8、sampleNum;isamp+)/inNum是輸入層的節(jié)點(diǎn)數(shù)(即神經(jīng)細(xì)胞數(shù))for(int i=0;i<inNum;i+)if(Math.Abs(pisamp,i)>pMax)pMax=Math.Abs(pisamp,i);for(int j=0;j<outNum;j+)if(Math.Abs(tisamp,j)>pMax)pMax=Math.Abs(tisamp,j);in_rate=pMax;/end isampfor(int isamp=0;isamp<sampleNum;isamp+)/數(shù)據(jù)歸一化for(int i=0;i<inNum;i+)xi=

9、pisamp,i/in_rate;for(int i=0;i<outNum;i+)ydi=tisamp,i/in_rate;/計(jì)算隱層的輸入和輸出for(int j=0;j<hideNum;j+)o1j=0.0;for(int i=0;i<inNum;i+)o1j+=wi,j*xi;/“權(quán)值”*“輸入”的那個(gè)累加的過(guò)程/這個(gè)b1j就是隱藏層的閾值，閾值就是一個(gè)輸入為“-1”的累加值x1j=1.0/(1.0+Math.Exp(-o1j-b1j);/計(jì)算輸出層的輸入和輸出for(int k=0;k<outNum;k+)o2k=0.0;for(int j=0;j<hid

10、eNum;j+)o2k+=vj,k*x1j;x2k=1.0/(1.0+Math.Exp(-o2k-b2k);/計(jì)算輸出層誤差和均方差for(int k=0;k<outNum;k+)/ydk是輸出層的教師數(shù)據(jù)，所謂教師數(shù)據(jù)就是實(shí)際應(yīng)該輸出的數(shù)據(jù)而已qqk=(ydk-x2k)*x2k*(1.0-x2k);e+=(ydk-x2k)*(ydk-x2k);/更新V，V矩陣是隱藏層與輸出層之間的權(quán)值for(int j=0;j<hideNum;j+)vj,k+=rate*qqk*x1j;/計(jì)算隱層誤差for(int j=0;j<hideNum;j+)/PP矩陣是隱藏層的誤差ppj=0.0;

11、/算法參考我的視頻截圖for(int k=0;k<outNum;k+)ppj+=qqk*vj,k;ppj=ppj*x1j*(1-x1j);/更新Wfor(int i=0;i<inNum;i+)wi,j+=rate*ppj*xi;/更新b2，輸出層的閾值for(int k=0;k<outNum;k+)b2k+=rate*qqk;/更新b1，隱藏層的閾值for(int j=0;j<hideNum;j+)b1j+=rate*ppj;/end isampe=Math.Sqrt(e);/均方差/ adjustWV(w,dw);/ adjustWV(v,dv);/end train

12、public void adjustWV(double , w,double, dw)for(int i=0;i<w.GetLength(0);i+)for(int j=0;j<w.GetLength(1);j+)wi,j+=dwi,j;public void adjustWV(double w,double dw)for(int i=0;i<w.Length;i+)wi+=dwi;/數(shù)據(jù)仿真函數(shù)public double sim(double psim)for(int i=0;i<inNum;i+)xi= psimi/in_rate;/in_rate為歸一化系數(shù)for

13、(int j=0;j<hideNum;j+)o1j=0.0;for(int i=0;i<inNum;i+)o1j=o1j+wi,j*xi;x1j=1.0/(1.0+Math.Exp(-o1j-b1j);for(int k=0;k<outNum;k+)o2k=0.0;for(int j=0;j<hideNum;j+)o2k=o2k+vj,k*x1j;x2k=1.0/(1.0+Math.Exp(-o2k-b2k);x2k=in_rate*x2k; return x2; /end sim/保存矩陣w,vpublic void saveMatrix(double ,

14、w,string filename)StreamWriter sw=File.CreateText(filename);for(int i=0;i<w.GetLength(0);i+)for(int j=0;j<w.GetLength(1);j+)sw.Write(wi,j+" ");sw.WriteLine();sw.Close();/保存矩陣b1,b2public void saveMatrix(double b,string filename)StreamWriter sw=File.CreateText(filename);for(int i=0;i&l

15、t;b.Length;i+)sw.Write(bi+" ");sw.Close();/讀取矩陣W,Vpublic void readMatrixW(double , w,string filename)StreamReader sr;try sr = new StreamReader(filename,Encoding.GetEncoding("gb2312"); String line;int i=0;while (line = sr.ReadLine() != null) string s1=line.Trim().Sp

16、lit(' ');for(int j=0;j<s1.Length;j+)wi,j=Convert.ToDouble(s1j);i+;sr.Close();catch (Exception e) / Let the user know what went wrong.Console.WriteLine("The file could not be read:");Console.WriteLine(e.Message);/讀取矩陣b1,b2public void readMatrixB(double b,string filename)Stre

17、amReader sr;try  sr = new StreamReader(filename,Encoding.GetEncoding("gb2312"); String line;int i=0; while (line = sr.ReadLine() != null) bi=Convert.ToDouble(line);i+;sr.Close();catch (Exception e) / Let the user know what went wrong.Console.WriteLine("Th

18、e file could not be read:");Console.WriteLine(e.Message); /end bpnet /end namespace/主調(diào)用程序namespace BpANNet/ <summary>/ Class1 的摘要說(shuō)明。/ </summary>class Class1/ <summary>/ 應(yīng)用程序的主入口點(diǎn)。/ </summary>STAThreadstatic void Main(string args)/0.1399,0.1467,0.1567,0.1595,0.1588,0

19、.1622,0.1611,0.1615,0.1685,0.1789,0.1790/ double , p1=new double,0.05,0.02,0.09,0.11,0.12,0.20,0.15,0.22,0.20,0.25,0.75,0.75,0.80,0.83,0.82,0.80,0.90,0.89,0.95,0.89,0.09,0.04,0.1,0.1,0.14,0.21,0.18,0.24,0.22,0.28,0.77,0.78,0.79,0.81,0.84,0.82,0.94,0.93,0.98,0.99;/ double , t1=new double,1,0,1,0,1,0,

20、1,0,1,0,0,1,0,1,0,1,0,1,0,1,1,0,1,0,1,0,1,0,1,0,0,1,0,1,0,1,0,1,0,1;/p1是輸入的信息，一共5組，輸入層為六個(gè)節(jié)點(diǎn),p156double , p1=new double,0.1399,0.1467,0.1567,0.1595,0.1588,0.1622,0.1467,0.1567,0.1595,0.1588,0.1622,0.1611,0.1567,0.1595,0.1588,0.1622,0.1611,0.1615,0.1595,0.1588,0.1622,0.1611,0.1615,0.1685,0.1588,0.1622

21、,0.1611,0.1615,0.1685,0.1789;/t1是輸出信息，一共6組，t161double , t1=new double,0.1622,0.1611,0.1615,0.1685,0.1789,0.1790;BpNet bp=new BpNet(p1,t1);int study=0;dostudy+;bp.train(p1,t1);/ bp.rate=0.95-(0.95-0.3)*study/50000;/ Console.Write("第 "+ study+"次學(xué)習(xí)： ");/ Console.WriteLine(" 均方差為 "+bp.e);while(bp.e>0.001 && study <50000);Console.Write("第 "+ study+"次學(xué)習(xí)： ");Console.WriteLine(" 均方差為 "+bp.e);bp.saveMatrix(bp.w,"w.txt&q