result.cpp

神经网络中的神经元代码

// RESULT.cpp : implementation file
//

#include "stdafx.h"
#include "neuro.h"
#include "RESULT.h"

#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif

/////////////////////////////////////////////////////////////////////////////
// CRESULT dialog


CRESULT::CRESULT(CWnd* pParent /*=NULL*/)
	: CDialog(CRESULT::IDD, pParent)
{
	//{{AFX_DATA_INIT(CRESULT)
		// NOTE: the ClassWizard will add member initialization here
	//}}AFX_DATA_INIT
}


void CRESULT::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CRESULT)
		// NOTE: the ClassWizard will add DDX and DDV calls here
	//}}AFX_DATA_MAP
}


BEGIN_MESSAGE_MAP(CRESULT, CDialog)
	//{{AFX_MSG_MAP(CRESULT)
	ON_BN_CLICKED(ID_MAP, OnMap)
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CRESULT message handlers

void CRESULT::OnMap() 
{
double a[10000];
double  T[10000];
int i;
for(i=0;i<10000;i++)
T[i]=0;
for(i=0;i<10000;i++)
{if(i<=2000)
a[i]=1;
else if(i<=4000)
a[i]=1.5;
else if(i<=6000)
a[i]=2;
else if(i<=8000)
a[i]=2.5;
else
a[i]=3;
}
for(i=0;i<10000;i++)
T[i+1]=T[i]+0.001;
	
  	
    CPaintDC dc(this); // device context for painting
	CDC *pDC=GetDC();

	       //初始化坐标值
	double data_x[10000]; 
	double data_y[10000];
	for(int j=0;j<10000;j++)
	{data_x[j]=T[j];
	data_y[j]=a[j];
	}
	
		//double initX[10]={1.0,2.0,3.0,4.0,5.0,6.0,7.0,8.0,9.0,10.0};
       //double initY[10]={0.5,2.0,2.8,4.0,5.6,6.0,7.0,9.4,13.8,23.4};
       //分别取得X和Y最大值,最小值
       double maxX=data_x[0];
       double minX=data_x[0];
       double maxY=data_y[0];
       double minY=data_y[0];

       for(i=0;i<10000;i++)
       {
              if(data_x[i]>maxX)
                     maxX=data_x[i];
              if(data_x[i]<minX)
                     minX=data_x[i];
              if(data_y[i]>maxY)
                     maxY=data_y[i];
              if(data_y[i]<minY)
                     minY=data_y[i];
       }
       //如果原点必须在X轴上，加上下面2行，否则注释掉
       if(minX>0)
              minX=0;
       //如果原点必须在Y轴上，加上下面2行，否则注释掉
       if(minY>0)
              minY=0;
       //确定图象显示大小
	   	   
       //确定坐标图四周预留的空白大小
	   
     
       int width=500;
       int height=300;
       //确定坐标图四周预留的空白大小
       const int mytop=100;
    const int mybottom=10;
    const int myleft=70;
    const int myright=10;
       //确定X，Y轴每单位显示宽度
       double intervalX=(width-myleft-myright)/(maxX-minX);
       double intervalY=(height-mybottom-mytop)/(maxY-minY);
       //绘制曲线。由于绘图坐标的Y轴是向下延升，所以这里每个点的Y值是用
       //图象高度减去y值大小。
       pDC->MoveTo(int(myleft+(data_x[0]-minX)*intervalX),int(height-(mybottom+(data_y[0]-minY)*intervalY)));
       for(i=0;i<10000;i++)
       {
              pDC->LineTo(int(myleft+(data_x[i]-minX)*intervalX), 
                     int(height-(mybottom+(data_y[i]-minY)*intervalY)));
       }
       //绘制X，Y轴
       //X轴从图形区域最左端到最右端
       double bottomY=0;
       double leftX=0;
       //bottomY表示X轴的y值，leftX表示Y轴的x值
       if(minY>0)
             bottomY=minY;
       if(minX>0)
              leftX=minX;
       pDC->MoveTo(int(myleft),int(height-(mybottom+(bottomY-minY)*intervalY)));
       pDC->LineTo(int(width-myright),int(height-(mybottom+(bottomY-minY)*intervalY)));
       //Y轴从图形区域最底端到最顶端
       pDC->MoveTo(int(myleft+(leftX-minX)*intervalX),int(height-mybottom));
       pDC->LineTo(int(myleft+(leftX-minX)*intervalX),int(mytop));
       //确定显示刻度个数
       const int count=5;const int count1=3;
       //确定每个显示刻度之间的宽度
       double spaceX=(width-myleft-myright)/count;
       double spaceY=(height-mybottom-mytop)/count1;
       //绘制刻度和刻度值
       CString str;
       //X轴
       for(i=0;i<=count;i++)
       {
              str.Format("%.1f",minX+i*(maxX-minX+1)/count);
              pDC->MoveTo(int(myleft+spaceX*i),int(height-(mybottom+(bottomY-minY)*intervalY)));
              pDC->LineTo(int(myleft+spaceX*i),int(height-(mybottom+(bottomY-minY)*intervalY+5)));
              pDC->TextOut(int(myleft+spaceX*i-10),
              int(height-(mybottom+(bottomY-minY)*intervalY-5)),str);
       }
       //Y轴
       for(i=0;i<=count1;i++)
       {
              str.Format("%.1f",minY+i*(maxY-minY)/count1);
              pDC->MoveTo(int(myleft+(leftX-minX)*intervalX),int(height-(mybottom+spaceY*i)));
              pDC->LineTo(int(myleft+(leftX-minX)*intervalX+5),int(height-(mybottom+spaceY*i)));
              pDC->TextOut(int(myleft+(leftX-minX)*intervalX-50),
              int(height-(mybottom+spaceY*i+8)),str);
       }
    //绘制X，Y轴的变量名
       pDC->TextOut((width-myleft)/2,height-5,"t/sec");
       pDC->TextOut(0,height/2+30,"depth/mm");	
}