bp神经网络详细步骤C#实现.doc

《bp神经网络详细步骤C#实现.doc》由会员分享，可在线阅读，更多相关《bp神经网络详细步骤C#实现.doc（13页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、如有侵权，请联系网站删除，仅供学习与交流bp神经网络详细步骤C#实现【精品文档】第 13 页using System;using System.Collections.Generic;using System.Linq;using System.Text;using System;using System.IO;using System.Text;namespace BpANNet/ / BpNet 的摘要说明。/ public class BpNetpublic int inNum;/输入节点数int hideNum;/隐层节点数public int outNum;/输出层节点数public

2、 int sampleNum;/样本总数Random R;double x;/输入节点的输入数据double x1;/隐层节点的输出double x2;/输出节点的输出double o1;/隐层的输入double o2;/输出层的输入public double , w;/权值矩阵w，这是输入层与隐藏层之间的权值矩阵public double , v;/权值矩阵V，这是隐藏层与输出层之间的权值矩阵public double , dw;/权值矩阵wpublic double , dv;/权值矩阵Vpublic double rate;/学习率public double b1;/隐层阈值矩阵publ

3、ic double b2;/输出层阈值矩阵public double db1;/隐层阈值矩阵public double db2;/输出层阈值矩阵double pp;/隐藏层的误差double qq;/输出层的误差double yd;/输出层的教师数据，所谓教师数据就是实际数据而已！public double e;/均方误差double in_rate;/归一化比例系数/用于确定隐藏层的神经细胞数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;i

4、nt ss=Convert.ToInt32(s);return (s-(double)ss)0.5) ? ss+1:ss;public BpNet(double , p,double , t)/ 构造函数逻辑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(输入节点数目： +inNum);Conso

5、le.WriteLine(隐层节点数目：+hideNum);Console.WriteLine(输出层节点数目：+outNum);Console.ReadLine();/将这些矩阵规定好矩阵大小x=new doubleinNum;x1=new doublehideNum;x2=new doubleoutNum;o1=new doublehideNum;o2=new doubleoutNum;w = new doubleinNum, hideNum;/权值矩阵w，这是输入层与隐藏层之间的权值矩阵v=new doublehideNum,outNum;dw=new doubleinNum,hideN

6、um;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;/输出层的教师数据/初始化wfor(int i=0;iinNum;i+)for(int j=0;jhideNum;j+)/NextDouble返回一个介于 0.0 和 1.0 之间的随机数。wi,

7、j=(R.NextDouble()*2-1.0)/2;/初始化vfor(int i=0;ihideNum;i+)for(int j=0;joutNum;j+)vi,j=(R.NextDouble()*2-1.0)/2;rate=0.8;e=0.0;in_rate=1.0;/训练函数public void train(double , p,double , t)e=0.0;/求p，t中的最大值double pMax=0.0;/sampleNum为样本总数for(int isamp=0;isampsampleNum;isamp+)/inNum是输入层的节点数(即神经细胞数)for(int i=0;

8、ipMax)pMax=Math.Abs(pisamp,i);for(int j=0;jpMax)pMax=Math.Abs(tisamp,j);in_rate=pMax;/end isampfor(int isamp=0;isampsampleNum;isamp+)/数据归一化for(int i=0;iinNum;i+)xi=pisamp,i/in_rate;for(int i=0;ioutNum;i+)ydi=tisamp,i/in_rate;/计算隐层的输入和输出for(int j=0;jhideNum;j+)o1j=0.0;for(int i=0;iinNum;i+)o1j+=wi,j*

9、xi;/“权值”*“输入”的那个累加的过程/这个b1j就是隐藏层的阈值，阈值就是一个输入为“-1”的累加值x1j=1.0/(1.0+Math.Exp(-o1j-b1j);/计算输出层的输入和输出for(int k=0;koutNum;k+)o2k=0.0;for(int j=0;jhideNum;j+)o2k+=vj,k*x1j;x2k=1.0/(1.0+Math.Exp(-o2k-b2k);/计算输出层误差和均方差for(int k=0;koutNum;k+)/ydk是输出层的教师数据，所谓教师数据就是实际应该输出的数据而已qqk=(ydk-x2k)*x2k*(1.0-x2k);e+=(yd

10、k-x2k)*(ydk-x2k);/更新V，V矩阵是隐藏层与输出层之间的权值for(int j=0;jhideNum;j+)vj,k+=rate*qqk*x1j;/计算隐层误差for(int j=0;jhideNum;j+)/PP矩阵是隐藏层的误差ppj=0.0;/算法参考我的视频截图for(int k=0;koutNum;k+)ppj+=qqk*vj,k;ppj=ppj*x1j*(1-x1j);/更新Wfor(int i=0;iinNum;i+)wi,j+=rate*ppj*xi;/更新b2，输出层的阈值for(int k=0;koutNum;k+)b2k+=rate*qqk;/更新b1，隐

11、藏层的阈值for(int j=0;jhideNum;j+)b1j+=rate*ppj;/end isampe=Math.Sqrt(e);/均方差/ adjustWV(w,dw);/ adjustWV(v,dv);/end trainpublic void adjustWV(double , w,double, dw)for(int i=0;iw.GetLength(0);i+)for(int j=0;jw.GetLength(1);j+)wi,j+=dwi,j;public void adjustWV(double w,double dw)for(int i=0;iw.Length;i+)wi

12、+=dwi;/数据仿真函数public double sim(double psim)for(int i=0;iinNum;i+)xi= psimi/in_rate;/in_rate为归一化系数for(int j=0;jhideNum;j+)o1j=0.0;for(int i=0;iinNum;i+)o1j=o1j+wi,j*xi;x1j=1.0/(1.0+Math.Exp(-o1j-b1j);for(int k=0;koutNum;k+)o2k=0.0;for(int j=0;jhideNum;j+)o2k=o2k+vj,k*x1j;x2k=1.0/(1.0+Math.Exp(-o2k-b2

13、k);x2k=in_rate*x2k;return x2; /end sim/保存矩阵w,vpublic void saveMatrix(double , w,string filename)StreamWriter sw=File.CreateText(filename);for(int i=0;iw.GetLength(0);i+)for(int j=0;jw.GetLength(1);j+)sw.Write(wi,j+ );sw.WriteLine();sw.Close();/保存矩阵b1,b2public void saveMatrix(double b,string filename

14、)StreamWriter sw=File.CreateText(filename);for(int i=0;ib.Length;i+)sw.Write(bi+ );sw.Close();/读取矩阵W,Vpublic void readMatrixW(double , w,string filename)StreamReader sr;trysr = new StreamReader(filename,Encoding.GetEncoding(gb2312);String line;int i=0;while (line = sr.ReadLine() != null)string s1=li

15、ne.Trim().Split( );for(int j=0;js1.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)StreamReader sr;trysr =

16、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(The file could not be read:);Console.WriteLine(e.Message);/end bpnet /end na

17、mespace/主调用程序namespace BpANNet/ / Class1 的摘要说明。/ class Class1/ / 应用程序的主入口点。/ STAThreadstatic void Main(string args)/0.1399,0.1467,0.1567,0.1595,0.1588,0.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.9

18、0,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,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组，输入层为六个节点,p156double , p1=new double,0.1399,0.1467,0.1567,0.1595,0

19、.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,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=

20、0;dostudy+;bp.train(p1,t1);/ bp.rate=0.95-(0.95-0.3)*study/50000;/ Console.Write(第 + study+次学习： );/ Console.WriteLine( 均方差为 +bp.e);while(bp.e0.001 & study 50000);Console.Write(第 + study+次学习： );Console.WriteLine( 均方差为 +bp.e);bp.saveMatrix(bp.w,w.txt);bp.saveMatrix(bp.v,v.txt);bp.saveMatrix(bp.b1,b1.t

21、xt);bp.saveMatrix(bp.b2,b2.txt);/ double , p2=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 , p2=new double,0.1399,0.1467,0.1567,0.1595,0.1588,0.1622,0.1622,0.1611,0.1615,0.1685,0.1789,0.1790;int aa=bp.inNum;int bb=bp.outNum;int cc=p2.GetLength(0);double p21=new doubleaa;double t2=new doublebb;for(int n=0;ncc;n+)for(int i=0;iaa;i+)p21i=p2n,i;t2=bp.sim(p21);for(int i=0;it2.Length;i+)Console.WriteLine(t2i+ );Console.ReadLine();