sccadlg.cpp

免疫算法求解函数优化问题运行时请选择所要的函数,并且点Button1执行

// SCCADlg.cpp : implementation file
//
/**********************************************************************
二进制编码的简单克隆选择算法源程序
武汉大学  软件工程国家重点实验室 林志毅
2006.10.17
**********************************************************************/
#include "stdafx.h"
#include "SCCA.h"
#include "math.h"
#include "stdlib.h"
#include "stdio.h"
#include "time.h"
#include "SCCADlg.h"

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

/////////////////////////////////////////////////////////////////////////////
// CAboutDlg dialog used for App About

class CAboutDlg : public CDialog
{
public:
	CAboutDlg();

// Dialog Data
	//{{AFX_DATA(CAboutDlg)
	enum { IDD = IDD_ABOUTBOX };
	//}}AFX_DATA

	// ClassWizard generated virtual function overrides
	//{{AFX_VIRTUAL(CAboutDlg)
	protected:
	virtual void DoDataExchange(CDataExchange* pDX);    // DDX/DDV support
	//}}AFX_VIRTUAL

// Implementation
protected:
	//{{AFX_MSG(CAboutDlg)
	//}}AFX_MSG
	DECLARE_MESSAGE_MAP()
};

CAboutDlg::CAboutDlg() : CDialog(CAboutDlg::IDD)
{
	//{{AFX_DATA_INIT(CAboutDlg)
	//}}AFX_DATA_INIT
}

void CAboutDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CAboutDlg)
	//}}AFX_DATA_MAP
}

BEGIN_MESSAGE_MAP(CAboutDlg, CDialog)
	//{{AFX_MSG_MAP(CAboutDlg)
		// No message handlers
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CSCCADlg dialog

CSCCADlg::CSCCADlg(CWnd* pParent /*=NULL*/)
	: CDialog(CSCCADlg::IDD, pParent)
{
	//{{AFX_DATA_INIT(CSCCADlg)
	m_SP2 = _T("");
	m_N = 0;
	m_c = 0;
	m_ge = 0;
	m_zz = 0.0;
	m_x2 = 0.0;
	m_x3 = 0.0;
	m_x4 = 0.0;
	m_x5 = 0.0;
	m_y1 = 0.0;
	m_y2 = 0.0;
	m_y3 = 0.0;
	m_y4 = 0.0;
	m_y5 = 0.0;
	m_Nb = 0.0;
	m_fev = 0;
	m_now = 0;
	m_x1 =0;
	//}}AFX_DATA_INIT
	// Note that LoadIcon does not require a subsequent DestroyIcon in Win32
	m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}

void CSCCADlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CSCCADlg)
	DDX_CBString(pDX, IDC_COMBO1, m_SP2);
	DDX_Text(pDX, IDC_EDIT1, m_N);
	DDX_Text(pDX, IDC_EDIT2, m_c);
	DDX_Text(pDX, IDC_EDIT3, m_ge);
	DDX_Text(pDX, IDC_EDIT4, m_zz);
	DDX_Text(pDX, IDC_EDIT5, m_x2);
	DDX_Text(pDX, IDC_EDIT6, m_x3);
	DDX_Text(pDX, IDC_EDIT7, m_x4);
	DDX_Text(pDX, IDC_EDIT8, m_x5);
	DDX_Text(pDX, IDC_EDIT9, m_y1);
	DDX_Text(pDX, IDC_EDIT10, m_y2);
	DDX_Text(pDX, IDC_EDIT11, m_y3);
	DDX_Text(pDX, IDC_EDIT12, m_y4);
	DDX_Text(pDX, IDC_EDIT13, m_y5);
	DDX_Text(pDX, IDC_EDIT14, m_Nb);
	DDX_Text(pDX, IDC_EDIT15, m_fev);
	DDX_Text(pDX, IDC_EDIT16, m_now);
	DDX_Text(pDX, IDC_EDIT17, m_x1);
	//}}AFX_DATA_MAP
}

BEGIN_MESSAGE_MAP(CSCCADlg, CDialog)
	//{{AFX_MSG_MAP(CSCCADlg)
	ON_WM_SYSCOMMAND()
	ON_WM_PAINT()
	ON_WM_QUERYDRAGICON()
	ON_BN_CLICKED(IDC_BUTTON1, OnButton1)
	ON_CBN_EDITCHANGE(IDC_COMBO1, OnEditchangeCombo1)
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CSCCADlg message handlers

BOOL CSCCADlg::OnInitDialog()
{
	CDialog::OnInitDialog();

	// Add "About..." menu item to system menu.

	// IDM_ABOUTBOX must be in the system command range.
	ASSERT((IDM_ABOUTBOX & 0xFFF0) == IDM_ABOUTBOX);
	ASSERT(IDM_ABOUTBOX < 0xF000);

	CMenu* pSysMenu = GetSystemMenu(FALSE);
	if (pSysMenu != NULL)
	{
		CString strAboutMenu;
		strAboutMenu.LoadString(IDS_ABOUTBOX);
		if (!strAboutMenu.IsEmpty())
		{
			pSysMenu->AppendMenu(MF_SEPARATOR);
			pSysMenu->AppendMenu(MF_STRING, IDM_ABOUTBOX, strAboutMenu);
		}
	}

	// Set the icon for this dialog.  The framework does this automatically
	//  when the application's main window is not a dialog
	SetIcon(m_hIcon, TRUE);			// Set big icon
	SetIcon(m_hIcon, FALSE);		// Set small icon
	
	// TODO: Add extra initialization here
	
	return TRUE;  // return TRUE  unless you set the focus to a control
}

void CSCCADlg::OnSysCommand(UINT nID, LPARAM lParam)
{
	if ((nID & 0xFFF0) == IDM_ABOUTBOX)
	{
		CAboutDlg dlgAbout;
		dlgAbout.DoModal();
	}
	else
	{
		CDialog::OnSysCommand(nID, lParam);
	}
}

// If you add a minimize button to your dialog, you will need the code below
//  to draw the icon.  For MFC applications using the document/view model,
//  this is automatically done for you by the framework.

void CSCCADlg::OnPaint() 
{
	if (IsIconic())
	{
		CPaintDC dc(this); // device context for painting

		SendMessage(WM_ICONERASEBKGND, (WPARAM) dc.GetSafeHdc(), 0);

		// Center icon in client rectangle
		int cxIcon = GetSystemMetrics(SM_CXICON);
		int cyIcon = GetSystemMetrics(SM_CYICON);
		CRect rect;
		GetClientRect(&rect);
		int x = (rect.Width() - cxIcon + 1) / 2;
		int y = (rect.Height() - cyIcon + 1) / 2;

		// Draw the icon
		dc.DrawIcon(x, y, m_hIcon);
	}
	else
	{
		CDialog::OnPaint();
	}
}

// The system calls this to obtain the cursor to display while the user drags
//  the minimized window.
HCURSOR CSCCADlg::OnQueryDragIcon()
{
	return (HCURSOR) m_hIcon;
}

void CSCCADlg::OnButton1() 
{
	// TODO: Add your control notification handler code here
	Window();
	//-----------------------------------------
    UpdateData();
	//-----------------------------------------
	if(m_SP2=="")
	{
		m_SP2="F001";
		UpdateData(false);
	}
	if(m_N==0||m_ge==0||m_Nb==0)
	{
		m_N=50;
		m_ge=100;
		m_Nb=5;
        UpdateData(false);
	}
    //设置算法默认参数
	//--------------------------------------------
	int N=m_N;
	int funn;
	int c;
	double Nb=m_Nb;
	int ge=m_ge;
	double vup,vdown,ting;
	int i,j,k,L;
	//定义算法变量
	//============================================
	//定义测试问题的变量区间
	if(m_SP2=="F001")
	{
		funn=1;
		vup=0.5;
		vdown=-0.5;
		c=1;
		ting=100000000000000000;
	}
	else if(m_SP2=="F003")
	{
		funn=3;
		vup=1;
		vdown=-1;
		c=2;
		ting=100000000000000000;
	}
	else if(m_SP2=="F004")
	{
		funn=4;
		vup=10;
		vdown=-10;
		c=2;
		ting=100000000000000000;
	}
	else if(m_SP2=="F005")
	{
		funn=5;
		vup=1;
		vdown=-1;
		c=2;
		ting=100000000000000000;
	}
	else if(m_SP2=="F006")
	{
		funn=6;
		vup=1;
		vdown=-1;
		c=2;
		ting=100000000000000000;
	}
	else if(m_SP2=="F007")
	{
		funn=7;
		vup=5.12;
		vdown=-5.12;
		c=2;
		ting=100000000000000000;
	}
	else if(m_SP2=="F008")
	{
		funn=8;
		vup=5.12;
		vdown=-5.12;
		c=2;
		ting=100000000000000000;
	}
    else if(m_SP2=="F009")
	{
		funn=9;
		vup=5.12;
		vdown=-5.12;
		c=2;
		ting=100000000000000000;
	}
	else if(m_SP2=="F011")
	{
		funn=11;
		vup=5.12;
		vdown=-5.12;
		c=2;
		ting=100000000000000000;
	}
	else if(m_SP2=="F013")
	{
		funn=13;
		vup=10;
		vdown=-10;
		c=2;
		ting=100000000000000000;
	}
	else if(m_SP2=="F014")
	{
		funn=14;
		vup=512;
		vdown=-512;
		c=2;
		ting=100000000000000000;
	}
	m_c=c;
	UpdateData(false);
    //-----------------------------------------------------------------------------------
	//算法程序（输入变量，N种群规模，funn测试问题，c变量维数,Nb克隆比例，ge最大进化代数
	//          vup变量上限，vdown变量下限，ting终止亲合度值（如果没有请设置成足够大的数）
	//输出变量：it进化代数，guiji[]最优值进化轨迹，fev评价次数
	//-----------------------------------------------------------------------------------
	double Nc;
	Nc=N*Nb;
    int CL=20;  //编码长度
	int(*v)[100];
	v=new int[300][100];
	double(*v1)[5];
	v1 = new double[300][5];

	double fal[5];
	double fit1[50],fit[300],guiji[10000];

	double z=32768,ab,max;
	double Ns,yop;
	int fev=0;
	int it=0;

	double sollution =-1*10000000;
	int sub[50][100];