开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > Windows编程 > DNA > 查看源码
OpticalNetworkDlg.cpp
资源名称：jiyinsuanfaVC.rar　[点击查看]
上传用户：gzsenex
上传日期：2019-05-12
资源大小：79k
文件大小：24k
源码类别：

DNA

开发平台：

Visual C++

OpticalNetworkDlg.cpp：源码内容
  1. // OpticalNetworkDlg.cpp : implementation file
  2. //
  4. #include "stdafx.h"
  5. #include "OpticalNetwork.h"
  6. #include "OpticalNetworkDlg.h"
  7. #include "math.h"
  8. #include "stdlib.h"
  9. #ifdef _DEBUG
  10. #define new DEBUG_NEW
  11. #undef THIS_FILE
  12. static char THIS_FILE[] = __FILE__;
  13. #endif
  14. #define  MAX_GROUPSIZE 100       // 最大群体规模
  15. #define  MAX_OBDNUM    10        // 最大OBD数目
  16. double TotalF=0;                 // 总适应度值
  17. int    *GeneSerial[MAX_GROUPSIZE];     // 当前一代个体序列   
  18. int    *NextGeneSerial[MAX_GROUPSIZE]; // 下一代个体序列
  19. double Cost[MAX_GROUPSIZE];      // 总费用
  20. double F[MAX_GROUPSIZE];         // 每代个体适应度值
  21. double Pm[MAX_GROUPSIZE];        // 辅助变量
  22. double *OBD_OLT;                 // OBD_OLT之间路径长度
  23. double *ONU_OBD;                 // ONU_OBD之间路径长度
  24. double temp[MAX_GROUPSIZE][20];
  25. int    nCurGANum    = 0;         // 基因算法当前进行的代数
  26. int    nGACrossNum  = 0;         // 交叉个体数目
  27. int    nGAVariNum   = 0;         // 变异个体数目
  28. int    nOBDLimit[MAX_OBDNUM];    // 每个OBD的最大连接数         
  29. int    nCurOBDLimit[MAX_GROUPSIZE][MAX_OBDNUM];   // 当前一代中每个OBD的连接数
  30. int    nNextOBDLimit[MAX_GROUPSIZE][MAX_OBDNUM];  // 进行交叉、变异等操作后，每个OBD的连接数
  31. int    colornum=0;               // 绘图时采用不同颜色的数目
  32. COLORREF COLOR[16]={0x008080F0,0x0000D080,0x00800080,0x008080E0,  // 用于绘制OLT、ONU、OBD连线的颜色
  33.                     0x0080E000,0x00E08000,0x0000E080,0x003020D0,
  34. 0x0080F080,0x00F08080,0x008080F0,0x00F0F080,
  35. 0x0080F020,0x00208080,0x008020F0,0x00A02080,};
  36.                     
  38. /////////////////////////////////////////////////////////////////////////////
  39. // CAboutDlg dialog used for App About
  41. class CAboutDlg : public CDialog
  42. {
  43. public:
  44. CAboutDlg();
  46. // Dialog Data
  47. //{{AFX_DATA(CAboutDlg)
  48. enum { IDD = IDD_ABOUTBOX };
  49. //}}AFX_DATA
  51. // ClassWizard generated virtual function overrides
  52. //{{AFX_VIRTUAL(CAboutDlg)
  53. protected:
  54. virtual void DoDataExchange(CDataExchange* pDX);    // DDX/DDV support
  55. //}}AFX_VIRTUAL
  57. // Implementation
  58. protected:
  59. //{{AFX_MSG(CAboutDlg)
  60. //}}AFX_MSG
  61. DECLARE_MESSAGE_MAP()
  62. };
  64. CAboutDlg::CAboutDlg() : CDialog(CAboutDlg::IDD)
  65. {
  66. //{{AFX_DATA_INIT(CAboutDlg)
  67. //}}AFX_DATA_INIT
  68. }
  70. void CAboutDlg::DoDataExchange(CDataExchange* pDX)
  71. {
  72. CDialog::DoDataExchange(pDX);
  73. //{{AFX_DATA_MAP(CAboutDlg)
  74. //}}AFX_DATA_MAP
  75. }
  77. BEGIN_MESSAGE_MAP(CAboutDlg, CDialog)
  78. //{{AFX_MSG_MAP(CAboutDlg)
  79. // No message handlers
  80. //}}AFX_MSG_MAP
  81. END_MESSAGE_MAP()
  83. /////////////////////////////////////////////////////////////////////////////
  84. // COpticalNetworkDlg dialog
  86. COpticalNetworkDlg::COpticalNetworkDlg(CWnd* pParent /*=NULL*/)
  87. : CDialog(COpticalNetworkDlg::IDD, pParent)
  88. {
  89. memset(m_OBDX,0,sizeof(double)*100);
  90. memset(m_OBDY,0,sizeof(double)*100);
  91. memset(m_ONUX,0,sizeof(double)*100);
  92. memset(m_ONUY,0,sizeof(double)*100);
  93. m_ONUX[0] =98;   m_ONUY[0] =70;
  94. m_ONUX[1] =16;   m_ONUY[1] =48;
  95. m_ONUX[2] =50;   m_ONUY[2] =15;
  96. m_ONUX[3] =94;   m_ONUY[3] =2;
  97. m_ONUX[4] =8;    m_ONUY[4] =8;
  98. m_ONUX[5] =88;   m_ONUY[5] =12;
  99. m_ONUX[6] =80;   m_ONUY[6] =30;
  100. m_ONUX[7] =18;   m_ONUY[7] =66;
  101. m_ONUX[8] =86;   m_ONUY[8] =45;
  102. m_ONUX[9] =22;   m_ONUY[9] =77;
  103. m_ONUX[10]=62;   m_ONUY[10]=86;
  104. m_ONUX[11]=24;   m_ONUY[11]=2;
  105. m_ONUX[12]=66;   m_ONUY[12]=32;
  106. m_ONUX[13]=87;   m_ONUY[13]=59;
  107. m_ONUX[14]=33;   m_ONUY[14]=88;
  108. m_ONUX[15]=44;   m_ONUY[15]=89;
  109. m_ONUX[16]=97;   m_ONUY[16]=8;
  110. m_ONUX[17]=18;   m_ONUY[17]=12;
  111. m_ONUX[18]=42;   m_ONUY[18]=58;
  112. m_ONUX[19]=73;   m_ONUY[19]=66;
  114. m_OBDX[0] =45;   m_OBDY[0] =55;
  115. m_OBDX[1] =78;   m_OBDY[1] =78;
  116. m_OBDX[2] =25;   m_OBDY[2] =25;
  117. m_OBDX[3] =30;   m_OBDY[3] =50;
  118. m_OBDX[4] =90;   m_OBDY[4] =70;
  119. m_OBDX[5] =30;   m_OBDY[5] =40;
  120. m_OBDX[6] =80;   m_OBDY[6] =10;
  121. m_OBDX[7] =10;   m_OBDY[7] =70;
  122. m_OBDX[8] =85;   m_OBDY[8] =40;
  123. m_OBDX[9] =80;   m_OBDY[9] =20;
  124. m_CurGANum = 0;
  125. m_MiniCost = 0.0;
  126. m_bIsGa = FALSE;
  127. for(int c=0;c<16;c++)
  128.     PenONU[c].CreatePen(PS_SOLID,1,COLOR[c]);
  129. PenOBD.CreatePen(PS_SOLID,2,0x00100000);
  130. //{{AFX_DATA_INIT(COpticalNetworkDlg)
  131. m_OLTX = 50.0;
  132. m_OLTY = 50.0;
  133. m_nONUNum = 20;
  134. m_GAVariProb = 0.05;
  135. m_GACrossProb = 0.6;
  136. m_GANum = 100;
  137. m_nGroupSize = 50;
  138. m_nOBDLimit = 4;
  139. m_nOBDNum = 10;
  140. m_SleepTime = 500;
  141. //}}AFX_DATA_INIT
  142. // Note that LoadIcon does not require a subsequent DestroyIcon in Win32
  143. m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
  144. }
  146. void COpticalNetworkDlg::DoDataExchange(CDataExchange* pDX)
  147. {
  148. CDialog::DoDataExchange(pDX);
  149. //{{AFX_DATA_MAP(COpticalNetworkDlg)
  150. DDX_Text(pDX, IDC_EDIT_OLTX, m_OLTX);
  151. DDV_MinMaxDouble(pDX, m_OLTX, 0., 100.);
  152. DDX_Text(pDX, IDC_EDIT_OLTY, m_OLTY);
  153. DDV_MinMaxDouble(pDX, m_OLTY, 0., 100.);
  154. DDX_Text(pDX, IDC_EDIT_ONUNUM, m_nONUNum);
  155. DDV_MinMaxInt(pDX, m_nONUNum, 1, 100);
  156. DDX_Text(pDX, IDC_EDIT_VARIPROB, m_GAVariProb);
  157. DDV_MinMaxDouble(pDX, m_GAVariProb, 0., 1.);
  158. DDX_Text(pDX, IDC_EDIT_CROSSPROB, m_GACrossProb);
  159. DDV_MinMaxDouble(pDX, m_GACrossProb, 0., 1.);
  160. DDX_Text(pDX, IDC_EDIT_GANUM, m_GANum);
  161. DDX_Text(pDX, IDC_EDIT_GROUPSIZE, m_nGroupSize);
  162. DDX_Text(pDX, IDC_EDIT_OBDLIMIT, m_nOBDLimit);
  163. DDX_Text(pDX, IDC_EDIT_OBDNUM, m_nOBDNum);   
  164. DDX_Text(pDX, IDC_EDIT_SLEEPTIME, m_SleepTime);
  165. DDV_MinMaxInt(pDX, m_SleepTime, 0, 10000);
  166. //}}AFX_DATA_MAP
  167. int i=0;
  168. int Start=0;
  169. for(Start=IDC_EDIT_OBDX1;Start<IDC_EDIT_OBDX1+10;Start++)
  170. {
  171. DDX_Text(pDX, IDC_EDIT_OBDX1+i, m_OBDX[i]);
  172.     i++;
  173. }
  174. i=0;
  175. for(Start=IDC_EDIT_OBDY1;Start<IDC_EDIT_OBDY1+10;Start++)
  176. {
  177. DDX_Text(pDX, IDC_EDIT_OBDY1+i, m_OBDY[i]);
  178.     i++;
  179. }
  180. i=0;
  181. for(Start=IDC_EDIT_ONUX1;Start<IDC_EDIT_ONUX1+20;Start++)
  182. {
  183. DDX_Text(pDX, IDC_EDIT_ONUX1+i, m_ONUX[i]);
  184.     i++;
  185. }
  186.     i=0;
  187. for(Start=IDC_EDIT_ONUY1;Start<IDC_EDIT_ONUY1+20;Start++)
  188. {
  189. DDX_Text(pDX, IDC_EDIT_ONUY1+i, m_ONUY[i]);
  190.     i++;
  191. }
  192. }
  194. BEGIN_MESSAGE_MAP(COpticalNetworkDlg, CDialog)
  195. //{{AFX_MSG_MAP(COpticalNetworkDlg)
  196. ON_WM_SYSCOMMAND()
  197. ON_WM_PAINT()
  198. ON_WM_QUERYDRAGICON()
  199. ON_BN_CLICKED(IDC_BUTTON_GA, OnButtonGa)
  200. ON_BN_CLICKED(IDABOUT, OnAbout)
  201. //}}AFX_MSG_MAP
  202. END_MESSAGE_MAP()
  204. /////////////////////////////////////////////////////////////////////////////
  205. // COpticalNetworkDlg message handlers
  207. BOOL COpticalNetworkDlg::OnInitDialog()
  208. {
  209. CDialog::OnInitDialog();
  211. // Add "About..." menu item to system menu.
  213. // IDM_ABOUTBOX must be in the system command range.
  214. ASSERT((IDM_ABOUTBOX & 0xFFF0) == IDM_ABOUTBOX);
  215. ASSERT(IDM_ABOUTBOX < 0xF000);
  217. CMenu* pSysMenu = GetSystemMenu(FALSE);
  218. if (pSysMenu != NULL)
  219. {
  220. CString strAboutMenu;
  221. strAboutMenu.LoadString(IDS_ABOUTBOX);
  222. if (!strAboutMenu.IsEmpty())
  223. {
  224. pSysMenu->AppendMenu(MF_SEPARATOR);
  225. pSysMenu->AppendMenu(MF_STRING, IDM_ABOUTBOX, strAboutMenu);
  226. }
  227. }
  229. // Set the icon for this dialog.  The framework does this automatically
  230. //  when the application's main window is not a dialog
  231. SetIcon(m_hIcon, TRUE); // Set big icon
  232. SetIcon(m_hIcon, FALSE); // Set small icon
  234. // TODO: Add extra initialization here
  235. CAboutDlg about;
  236. about.DoModal();
  238. return TRUE;  // return TRUE  unless you set the focus to a control
  239. }
  241. void COpticalNetworkDlg::OnSysCommand(UINT nID, LPARAM lParam)
  242. {
  243. if ((nID & 0xFFF0) == IDM_ABOUTBOX)
  244. {
  245. CAboutDlg dlgAbout;
  246. dlgAbout.DoModal();
  247. }
  248. else
  249. {
  250. CDialog::OnSysCommand(nID, lParam);
  251. }
  252. }
  254. // If you add a minimize button to your dialog, you will need the code below
  255. //  to draw the icon.  For MFC applications using the document/view model,
  256. //  this is automatically done for you by the framework.
  258. void COpticalNetworkDlg::OnPaint() 
  259. {
  260. if (IsIconic())
  261. {
  262. CPaintDC dc(this); // device context for painting
  264. SendMessage(WM_ICONERASEBKGND, (WPARAM) dc.GetSafeHdc(), 0);
  266. // Center icon in client rectangle
  267. int cxIcon = GetSystemMetrics(SM_CXICON);
  268. int cyIcon = GetSystemMetrics(SM_CYICON);
  269. CRect rect;
  270. GetClientRect(&rect);
  271. int x = (rect.Width() - cxIcon + 1) / 2;
  272. int y = (rect.Height() - cyIcon + 1) / 2;
  274. // Draw the icon
  275. dc.DrawIcon(x, y, m_hIcon);
  276. }
  277. else
  278. {   
  279. if(!m_bIsGa)
  280. {
  281.    int NetLeft,NetBottom;
  282.    int x_coordinate=0,y_coordinate=0;
  283.        char txt[4];
  284.        CString info;
  285.        CPen Pen;
  286.        CPen* OldPen;  
  287.        CDC* pDC=GetDC();
  288.    CRect RectClient,Workarea;
  289.            GetClientRect(RectClient);
  290.            Workarea.left=RectClient.left+RectClient.right/3-15;
  291.            Workarea.right=RectClient.right;
  292.            Workarea.top=RectClient.top;
  293.            Workarea.bottom=RectClient.bottom;
  294.            pDC->Rectangle(Workarea);
  295.            pDC->SetBkColor(0x00FFFFFF);
  296.        NetLeft=Workarea.left+25;
  297.        NetBottom=Workarea.bottom-20;
  298.        pDC->MoveTo(NetLeft,NetBottom);
  299.            pDC->LineTo(NetLeft+500,NetBottom);
  300.        pDC->MoveTo(NetLeft,NetBottom);
  301.        pDC->LineTo(NetLeft,NetBottom-500);  
  302.        pDC->TextOut(NetLeft+500,NetBottom-12,"X");
  303.        pDC->TextOut(NetLeft+5,NetBottom-505,"Y");
  304.        pDC->SetTextColor(0x000000FF);
  305.        info.Format("已进化代数 %d",0);//+1);//wsprintf(txt,"%d",m_CurGANum+1);
  306.        pDC->TextOut(NetLeft+100,NetBottom-520,info);
  307.        info="最小费用 ";//_MiniCost);//wsprintf(txt,"%4.6f",m_MiniCost);
  308.        pDC->TextOut(NetLeft+250,NetBottom-520,info);
  309.  
  310.        Pen.CreatePen(PS_SOLID,2,0x00100000);
  311.        OldPen=pDC->SelectObject(&PenOBD);
  312.        for(int b1=0;b1<m_nOBDNum;b1++)
  313.    {   
  314.        pDC->MoveTo(NetLeft+m_OBDX[b1]*5,NetBottom-m_OBDY[b1]*5);
  315.        pDC->LineTo(NetLeft+m_OLTX*5,NetBottom-m_OLTY*5);
  316.    }
  317.        Pen.DeleteObject();
  318.        pDC->SelectObject(OldPen);
  319.        pDC->SetTextColor(GetSysColor(COLOR_WINDOWTEXT));
  320.        for(x_coordinate=0;x_coordinate<=10;x_coordinate++)
  321.    {
  322.            char xn[4];
  323.          wsprintf(xn,"%d",x_coordinate*10);
  324.        pDC->TextOut(NetLeft+x_coordinate*50-8,NetBottom+2,xn);
  325.    }
  326.        for(y_coordinate=1;y_coordinate<=10;y_coordinate++)
  327.    {
  328.            char yn[4];
  329.        wsprintf(yn,"%d",y_coordinate*10);
  330.        if(y_coordinate==10)
  331.        pDC->TextOut(NetLeft-24,NetBottom-y_coordinate*50,yn);
  332.        else
  333.        pDC->TextOut(NetLeft-22,NetBottom-y_coordinate*50,yn);
  334.    }
  335.        pDC->SetTextColor(0x000000FF);
  336.        pDC->TextOut(NetLeft+m_OLTX*5-7,NetBottom-m_OLTY*5-7,"★");
  337.        pDC->SetTextColor(0x00FF0000);
  338.        for(int a1=0;a1<m_nOBDNum;a1++)
  339.    {
  340.        wsprintf(txt,"%d",a1+1);
  341.            pDC->TextOut(NetLeft+m_OBDX[a1]*5-3,NetBottom-m_OBDY[a1]*5-9,txt);
  342.    }
  343.        pDC->SetTextColor(0x0000FF00);
  344.        for(int a2=0;a2<m_nONUNum;a2++)
  345.    {
  346.                wsprintf(txt,"%d",a2+1);
  347.                pDC->TextOut(NetLeft+m_ONUX[a2]*5-3,NetBottom-m_ONUY[a2]*5-9,txt);
  348.    }
  349.        pDC->SetTextColor(0x00FFFFFF);
  350.        pDC->TextOut(NetLeft+507,NetBottom-508,"c");
  351.        pDC->SetTextColor(GetSysColor(COLOR_WINDOWTEXT));
  352.        pDC->SetBkColor(GetSysColor(COLOR_WINDOW));
  353. }
  354.     else
  355.  DrawNetwork();
  356.  CDialog::OnPaint();
  357. }
  358. }
  360. // The system calls this to obtain the cursor to display while the user drags
  361. //  the minimized window.
  362. HCURSOR COpticalNetworkDlg::OnQueryDragIcon()
  363. {
  364. return (HCURSOR) m_hIcon;
  365. }
  367. // 搜索最小C值
  368. double COpticalNetworkDlg::FindMiniC(double *data, int len)
  369. { 
  370. double mini=data[0];
  371. int mininum=0;
  372. for(int n=0;n<len;n++)
  373. {
  374. if(data[n]<mini)
  375. {
  376. mini=data[n];
  377.     mininum=n;
  378. }
  379. }
  380. return mini;
  381. }
  383. // 基因算法的具体实现
  384. void COpticalNetworkDlg::OnButtonGa() 
  385. {     
  386.  // 取得基因算法所需参数 
  387.  UpdateData(TRUE);
  388.      // 分配内存
  389.      OBD_OLT=(double *)malloc(sizeof(double)*m_nOBDNum);
  390.      ONU_OBD=(double *)malloc(sizeof(double)*m_nONUNum*m_nOBDNum);
  391.      // 根据群体规模和交叉概率计算出交叉个体个数
  392.  nGACrossNum = m_nGroupSize*m_GACrossProb/2;
  393.  // 根据群体规模和变异概率计算出变异个体个数
  394.      nGAVariNum  = m_nGroupSize*m_GAVariProb;
  395.  // 设置已进行过基因算法标志
  396.  m_bIsGa=TRUE;
  397.  // 为当前一代个体分配内存
  398.      for(int k=0;k<m_nGroupSize;k++)
  399.           GeneSerial[k]=(int *) malloc(sizeof(int)*m_nONUNum);
  400.   // 为当前一代个体交叉和变异后的新产生的个体分配内存
  401.      for(int u=0;u<nGACrossNum*2+nGAVariNum;u++)
  402.           NextGeneSerial[u]=(int *)malloc(sizeof(int)*m_nONUNum);
  403.      // 参数有效性检查
  404.  if(m_nONUNum>=m_nOBDLimit*m_nOBDNum)
  405.  {
  406.   AfxMessageBox("ONU数目过多，或者OBD数目过少，或OBD最大连接数目过小!",MB_OK|MB_ICONINFORMATION,NULL);
  407.           return; 
  408.  }
  409.  if(m_nGroupSize<nGACrossNum*2+nGAVariNum)
  410.  {
  411.   AfxMessageBox("群体规模过小，或者交叉概率过大，或者变异概率过大!",MB_OK|MB_ICONINFORMATION,NULL);
  412.           return; 
  413.  }
  414.      for(int n1=0;n1<m_nOBDNum;n1++)
  415.      {
  416.           OBD_OLT[n1]=sqrt((m_OBDX[n1]-m_OLTX)*(m_OBDX[n1]-m_OLTX)+(m_OBDY[n1]-m_OLTY)*(m_OBDY[n1]-m_OLTY));
  417.           nOBDLimit[n1]=m_nOBDLimit;
  418.      } 
  419.      int LimitTotal=0;
  420.      int CurGroupSize=1;
  421.      for(int p=0;p<m_nOBDNum;p++)
  422.           LimitTotal+=nOBDLimit[p];
  423.      for(int q=0;q<m_nONUNum;q++)
  424.      {
  425.           if(CurGroupSize>MAX_GROUPSIZE*10)
  426.           break;
  427.           else
  428.           {
  429.              CurGroupSize=LimitTotal*CurGroupSize;
  430.              LimitTotal--;
  431.           }
  432.      }
  433.      if(m_nGroupSize>CurGroupSize||m_nGroupSize>MAX_GROUPSIZE*10)
  434.           return;
  435.  // 计算ONU和OBD之间距离
  436.      for(int i=0;i<m_nONUNum;i++)
  437.           for(int j=0;j<m_nOBDNum;j++)                                            
  438.      *(ONU_OBD+i*m_nOBDNum+j)=sqrt((m_ONUX[i]-m_OBDX[j])*(m_ONUX[i]-m_OBDX[j])+(m_ONUY[i]-m_OBDY[j])*(m_ONUY[i]-m_OBDY[j]));
  439.      // 产生初始基因序列   
  440.  for(int i1=0;i1<m_nGroupSize;i1++)
  441.      {          
  442.           for(int k=0;k<m_nOBDNum;k++)
  443.               nCurOBDLimit[i1][k]=0;
  444.           for(int j1=0;j1<m_nONUNum;j1++)
  445.           {
  446.              int RandInt;
  447.       Random:RandInt=RandomInt(0,m_nOBDNum-1);
  448.              if(nCurOBDLimit[i1][RandInt]<=nOBDLimit[RandInt]-1)
  449.              {
  450.                 *(GeneSerial[i1]+j1)=RandInt;
  451.                 nCurOBDLimit[i1][RandInt]++; 
  452.              }
  453.              else
  454.                 goto Random;   // 不满足约束条件2时，重新产生基因序列
  455.           } 
  456.      }
  457.  // 基因算法具体实现
  458.      for( nCurGANum=0;nCurGANum<m_GANum;nCurGANum++)
  459.      {   
  460.           // 计算当前一代群体中每个个体的适应度数值F
  461.   for(int i=0;i<m_nGroupSize;i++)
  462.           {    
  463.              Cost[i]=0;
  464.  for(int j=0;j<m_nONUNum;j++)
  465.  {   
  466.  Cost[i]=Cost[i]+*(ONU_OBD+j*m_nOBDNum+*(GeneSerial[i]+j))+*(OBD_OLT+*(GeneSerial[i]+j));             
  467.          temp[i][j]=*(ONU_OBD+j*m_nOBDNum+*(GeneSerial[i]+j));
  468.  }  
  469.  F[i]=1/Cost[i]; 
  470.  TotalF=F[i]+TotalF;
  471.           }
  472.           // 归一化F值
  473.   for(int a=0;a<m_nGroupSize;a++)
  474.              Pm[a]=F[a]/TotalF*100;
  475.           // 将当前一代群体中的个体按F值从大到小排序
  476.   for(int b=0;b<m_nGroupSize-1;b++)
  477.   {
  478.   FindMiniF(Pm,m_nGroupSize-b);
  479.   }
  480.   m_CurGANum=nCurGANum;
  481.        // 查找当前一代中的最小费用个体
  482.   m_MiniCost=FindMiniC(Cost,m_nGroupSize);
  483.   // 绘制图形
  484.   DrawNetwork();
  485.           // 交叉操作
  486.   for(int c=0;c<nGACrossNum;c++)
  487.           {  
  488.  int nExchangePos,temp=0;
  489.              int nNextOBD1;
  490.              int nNextOBD2;
  491.  int LimitLoop=0;
  492.              int new1=0;     // 发生交叉个体1在当前一代群体中位置
  493.  int new2=0;     // 发生交叉个体2在当前一代群体中位置
  494.              //int temp=0;
  495.       Limit1:new1=RandomInt(0,m_nGroupSize-1);   
  496.  temp=RandomInt(0,m_nGroupSize-1);
  497.              // 由于当前一代群体中个体已经按照F从大到小排序
  498.  // 如果new1>temp,则位置new1处个体F值小于temp处个体F值
  499.  // 所以选定new1=temp处个体进行交叉，体现了F值较大的个体有较多机会进行交叉
  500.  if(new1>temp)  
  501.                 new1=temp;
  502.              new2=RandomInt(0,m_nGroupSize-1);
  503.              temp=RandomInt(0,m_nGroupSize-1);
  504.              if(new2>temp)
  505.                 new2=temp;             
  506.              for(int g=0;g<m_nOBDNum;g++)
  507.              { 
  508.                  nNextOBDLimit[2*c][g]=nCurOBDLimit[new1][g];
  509.                  nNextOBDLimit[2*c+1][g]=nCurOBDLimit[new2][g];
  510.              } 
  511.             
  512.       Limit2:nExchangePos=RandomInt(0,m_nONUNum-1);  // 被选中交叉个体的基因序列中用于交叉的基因位
  513.              nNextOBD1=nNextOBDLimit[2*c][*(GeneSerial[new2]+nExchangePos)]+1;
  514.              nNextOBD2=nNextOBDLimit[2*c+1][*(GeneSerial[new1]+nExchangePos)]+1;
  515.  LimitLoop++;
  516.              if((nNextOBD1>nOBDLimit[*(GeneSerial[new2]+nExchangePos)])||(nNextOBD2>nOBDLimit[*(GeneSerial[new1]+nExchangePos)]))
  517.                  if(LimitLoop>50)   // 跳转到Limit2重新产生交叉位满50次，仍未满足约束条件2，
  518. goto Limit1;    // 则跳转到Limit1，重新产生交叉个体在当前一代群体中位置
  519.  else
  520.     goto Limit2;    // 交叉后，不满足约束条件2，跳转到Limit2重新产生交叉位
  521.              
  522.      // 交叉成功，产生新的个体
  523.  for(int m=0;m<m_nONUNum;m++)
  524.              {   
  525.                  if(m==nExchangePos)
  526.                  {
  527.                     *(NextGeneSerial[2*c]+nExchangePos)=*(GeneSerial[new2]+nExchangePos);
  528.                     *(NextGeneSerial[2*c+1]+nExchangePos)=*(GeneSerial[new1]+nExchangePos);
  529.                     nNextOBDLimit[2*c][*(GeneSerial[new2]+nExchangePos)]++;
  530.                     nNextOBDLimit[2*c][*(GeneSerial[new1]+nExchangePos)]--;
  531.                     nNextOBDLimit[2*c+1][*(GeneSerial[new2]+nExchangePos)]--;
  532.                     nNextOBDLimit[2*c+1][*(GeneSerial[new1]+nExchangePos)]++;                   
  533.                  }
  534.                  else
  535.                  {
  536.                     *(NextGeneSerial[2*c]+m)=*(GeneSerial[new1]+m);
  537.                     *(NextGeneSerial[2*c+1]+m)=*(GeneSerial[new2]+m);
  538.                  }
  539.               }            
  540.            } 
  542.            // 变异操作
  543.    for(int d=0;d<nGAVariNum;d++)
  544.            {
  545.                int nVariIndividual=0;
  546.                int nVariPos=0;
  547.    int VariLoop=0;
  548.          Vari1:VariLoop=0;
  549.    nVariIndividual=RandomInt(0,m_nGroupSize-1); // 发生变异个体在当前一代群体中位置
  550.    nVariPos=RandomInt(0,m_nONUNum-1);           // 发生变异个体的基因序列中变异的基因位
  551.                for(int z=0;z<m_nOBDNum;z++)
  552.                    nNextOBDLimit[2*nGACrossNum+d][z]=nCurOBDLimit[nVariIndividual][z];
  553.                for(int h=0;h<m_nONUNum;h++)
  554.                {
  555.                    if(h==nVariPos)
  556.                    {
  557.                       int nVariValue;
  558.                       int nNextVari;
  559.                Vari2: nVariValue=RandomInt(0,m_nOBDNum-1);  // 产生变异位变异后的值
  560.                       if(nVariValue==*(GeneSerial[nVariIndividual]+nVariPos))
  561.   {   
  562.   VariLoop++;
  563.   if(VariLoop==100) // 跳转到Vari2满50次，随机产生的变异位的值仍等于待变异位的当前值
  564.   goto Vari1;   // 则跳转到Vair1，重新产生变异个体在当前一代群体中位置
  565.   else
  566.       goto Vari2; // 变异位的值等于待变异位的当前值，跳转到Vari2，重新产生变异位变异后的值
  567.   }
  568.                       nNextVari=nNextOBDLimit[2*nGACrossNum+d][nVariValue]+1;
  569.                       if(nNextVari>nOBDLimit[nVariValue])
  570.                       {   
  571.   VariLoop++;
  572.   if(VariLoop==100)// 跳转到Vari2重新产生变异位满50次，仍未满足约束条件2，
  573.   goto Vari1;  // 则跳转到Vair1，重新产生变异个体在当前一代群体中位置
  574.   else
  575.       goto Vari2;  // 变异后，不满足约束条件2，跳转到Vari2重新产生变异位
  576.   }
  577.   // 变异成功
  578.                       *(NextGeneSerial[2*nGACrossNum+d]+h)=nVariValue;
  579.                       nNextOBDLimit[2*nGACrossNum+d][nVariValue]++;
  580.                       nNextOBDLimit[2*nGACrossNum+d][*(GeneSerial[nVariIndividual]+nVariPos)]--;
  581.                    }
  582.                    else
  583.                       *(NextGeneSerial[2*nGACrossNum+d]+h)=*(GeneSerial[nVariIndividual]+h);
  584.                 }
  585.            } 
  586.            
  587.            // 产生下一代群体
  588.    for(int e=0;e<2*nGACrossNum+nGAVariNum;e++)
  589.            {
  590.                for(int f=0;f<m_nONUNum;f++)
  591.                    *(GeneSerial[m_nGroupSize-e-1]+f)=*(NextGeneSerial[e]+f);
  592.                for(int x=0;x<m_nOBDNum;x++)
  593.                    nCurOBDLimit[m_nGroupSize-e-1][x]=nNextOBDLimit[e][x];
  594.            }
  595.    Sleep(m_SleepTime);
  596.       }
  597. }
  599. // 产生1～100之间的随机整数
  600. int COpticalNetworkDlg::RandomInt(int low,int high)
  601. {
  602.     int result;
  603. result=rand();
  604.     return result%(high+1);
  605. }
  607. // 搜索最小F值
  608. void COpticalNetworkDlg::FindMiniF(double *data,int len)
  609. {
  610.    double mini=data[0];
  611.    int    mininum=0;
  612.    int    CurOBDTemp;
  613.    int    *GeneSerTemp;
  614.    double datatemp;   
  615.    for(int n=0;n<len;n++)
  616.       if(mini>data[n])
  617.       {
  618.          mini=data[n];
  619.          mininum=n;
  620.       }
  621.    if(mininum!=len-1)
  622.    {
  623.       datatemp=*(data+len-1);
  624.       *(data+len-1)=*(data+mininum);
  625.       *(data+mininum)=datatemp;
  626.       GeneSerTemp=GeneSerial[len-1];
  627.   GeneSerial[len-1]=GeneSerial[mininum];
  628.   GeneSerial[mininum]=GeneSerTemp;   
  629.   for(int n2=0;n2<m_nOBDNum;n2++)
  630.   {
  631.   CurOBDTemp=nCurOBDLimit[mininum][n2];
  632.       nCurOBDLimit[mininum][n2]=nCurOBDLimit[len-1][n2];
  633.   nCurOBDLimit[len-1][n2]=CurOBDTemp; 
  634.   }
  635.    }
  636. }
  638. void COpticalNetworkDlg::DrawNetwork()
  639. {
  640.      int NetLeft,NetBottom;
  641.  int x_coordinate=0,y_coordinate=0;
  642.  char txt[4];
  643.  CString info;  
  644.  CDC* pDC=GetDC();
  645.      CRect RectClient,Workarea;
  646.      GetClientRect(RectClient);
  647.      Workarea.left=RectClient.left+RectClient.right/3-15;
  648.      Workarea.right=RectClient.right;
  649.      Workarea.top=RectClient.top;
  650.      Workarea.bottom=RectClient.bottom;
  651.      pDC->Rectangle(Workarea);
  652.      pDC->SetBkColor(0x00FFFFFF);
  653.  NetLeft=Workarea.left+25;
  654.  NetBottom=Workarea.bottom-20;
  655.  pDC->MoveTo(NetLeft,NetBottom);
  656.      pDC->LineTo(NetLeft+500,NetBottom);
  657.  pDC->MoveTo(NetLeft,NetBottom);
  658.  pDC->LineTo(NetLeft,NetBottom-500);  
  659.  pDC->TextOut(NetLeft+500,NetBottom-12,"X");
  660.  pDC->TextOut(NetLeft+5,NetBottom-505,"Y");
  661.  pDC->SetTextColor(0x000000FF);
  662.  info.Format("已进化代数 %d",m_CurGANum+1);
  663.  pDC->TextOut(NetLeft+100,NetBottom-520,info);
  664.  info.Format("最小费用 %5.5f",m_MiniCost);
  665.  pDC->TextOut(NetLeft+250,NetBottom-520,info);
  666.  pDC->SetTextColor(GetSysColor(COLOR_WINDOWTEXT));
  667.  for(x_coordinate=0;x_coordinate<=10;x_coordinate++)
  668.  {
  669.     char xn[4];
  670. wsprintf(xn,"%d",x_coordinate*10);
  671. pDC->TextOut(NetLeft+x_coordinate*50-8,NetBottom+2,xn);
  672.  }
  673.  for(y_coordinate=1;y_coordinate<=10;y_coordinate++)
  674.  {
  675.     char yn[4];
  676. wsprintf(yn,"%d",y_coordinate*10);
  677. if(y_coordinate==10)
  678. pDC->TextOut(NetLeft-24,NetBottom-y_coordinate*50,yn);
  679. else
  680. pDC->TextOut(NetLeft-22,NetBottom-y_coordinate*50,yn);
  681.  }
  683.  PtrOldPen=pDC->SelectObject(&PenOBD);
  684.  for(int b1=0;b1<m_nOBDNum;b1++)
  685.  {   
  686.  pDC->MoveTo(NetLeft+m_OBDX[b1]*5,NetBottom-m_OBDY[b1]*5);
  687.  pDC->LineTo(NetLeft+m_OLTX*5,NetBottom-m_OLTY*5);
  688.  }
  689.  for(int b2=0;b2<m_nONUNum;b2++)
  690.  {  
  691.  pDC->SelectObject(&PenONU[colornum]);
  692.  pDC->MoveTo(NetLeft+m_ONUX[b2]*5,NetBottom-m_ONUY[b2]*5);
  693.  pDC->LineTo(NetLeft+m_OBDX[*(GeneSerial[0]+b2)]*5,NetBottom-m_OBDY[*(GeneSerial[0]+b2)]*5);
  694.  colornum++;
  695.  if(colornum==16)
  696. colornum=0;  
  697.  }
  698.  pDC->SelectObject(PtrOldPen);    
  699.  pDC->SetTextColor(0x000000FF);
  700.  pDC->TextOut(NetLeft+m_OLTX*5-7,NetBottom-m_OLTY*5-7,"★");
  701.  pDC->SetTextColor(0x00FF0000);
  702.  for(int a1=0;a1<m_nOBDNum;a1++)
  703.  {
  704.  wsprintf(txt,"%d",a1+1);
  705.      pDC->TextOut(NetLeft+m_OBDX[a1]*5-3,NetBottom-m_OBDY[a1]*5-9,txt);
  706.  }
  707.  pDC->SetTextColor(0x0000FF00);
  708.  for(int a2=0;a2<m_nONUNum;a2++)
  709.  {
  710.          wsprintf(txt,"%d",a2+1);
  711.          pDC->TextOut(NetLeft+m_ONUX[a2]*5-3,NetBottom-m_ONUY[a2]*5-9,txt);
  712.  }
  713.  pDC->SetTextColor(0x00FFFFFF);
  714.  pDC->TextOut(NetLeft+507,NetBottom-508,"c");
  715.  pDC->SetTextColor(GetSysColor(COLOR_WINDOWTEXT));
  716.  pDC->SetBkColor(GetSysColor(COLOR_WINDOW));
  717.  
  718. }
  720. void COpticalNetworkDlg::OnAbout() 
  721. {
  722.      CAboutDlg dlg;
  723.      dlg.DoModal();
  724. }
  726. void COpticalNetworkDlg::OnOK() 
  727. {
  728. CDialog::OnOK();
  729. }