贝叶斯决策dlg.cpp

贝叶斯决策解决城市最短路径问题,人工只能实验项目

// 贝叶斯决策Dlg.cpp : implementation file
//

#include "stdafx.h"
#include "贝叶斯决策.h"
#include "贝叶斯决策Dlg.h"
#include "math.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()

/////////////////////////////////////////////////////////////////////////////
// CMyDlg dialog

CMyDlg::CMyDlg(CWnd* pParent /*=NULL*/)
	: CDialog(CMyDlg::IDD, pParent)
{
	f=i=k=0;
	move=true;
	select=true;
	//{{AFX_DATA_INIT(CMyDlg)
	m_Lei = 0;
	m_TeZheng = 0;
	m_X = _T("0");
	m_f = 2.0;
	//}}AFX_DATA_INIT
	// Note that LoadIcon does not require a subsequent DestroyIcon in Win32
	m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}

void CMyDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CMyDlg)
	DDX_Control(pDX, IDC_LIST4, m_List4);
	DDX_Control(pDX, IDC_LIST3, m_List3);
	DDX_Control(pDX, IDC_LIST2, m_List2);
	DDX_Control(pDX, IDC_LIST1, m_List1);
	DDX_Text(pDX, IDC_EDIT1, m_Lei);
	DDX_Text(pDX, IDC_EDIT2, m_TeZheng);
	DDX_Text(pDX, IDC_EDIT3, m_X);
	DDX_Text(pDX, IDC_EDIT4, m_f);
	DDV_MinMaxDouble(pDX, m_f, 1.e-011, 99999999999.);
	//}}AFX_DATA_MAP
}

BEGIN_MESSAGE_MAP(CMyDlg, CDialog)
	//{{AFX_MSG_MAP(CMyDlg)
	ON_WM_SYSCOMMAND()
	ON_WM_PAINT()
	ON_WM_QUERYDRAGICON()
	ON_BN_CLICKED(IDC_BUTTON1, OnButton1)
	ON_LBN_SELCHANGE(IDC_LIST1, OnSelchangeList1)
	ON_LBN_SELCHANGE(IDC_LIST2, OnSelchangeList2)
	ON_LBN_SELCHANGE(IDC_LIST3, OnSelchangeList3)
	ON_BN_CLICKED(IDC_BUTTON2, OnButton2)
	ON_BN_CLICKED(IDC_RADIO1, OnRadio1)
	ON_BN_CLICKED(IDC_RADIO2, OnRadio2)
	ON_BN_CLICKED(IDC_BUTTON3, OnButton3)
	ON_BN_CLICKED(IDC_BUTTON4, OnButton4)
	ON_WM_MOUSEMOVE()
	ON_WM_SETCURSOR()
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CMyDlg message handlers

BOOL CMyDlg::OnInitDialog()
{
	CDialog::OnInitDialog();
	((CButton*)GetDlgItem(IDC_RADIO1))->SetCheck(1);
	// 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 CMyDlg::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 CMyDlg::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 CMyDlg::OnQueryDragIcon()
{
	return (HCURSOR) m_hIcon;
}

void CMyDlg::OnOK() 
{
	// TODO: Add extra validation here

	UpdateData(true);
	double *E,*e,*pj,*E1,*g;
	double *m1;
	double *m2;
	double *m3;
	double m4=0;
	double p3;
	double p1=0;
	int p2;
	CString str;

	int result;
	int m,j,a,b,c,num=0;
	if (f<k||f==0) 
	{
		MessageBox("样本不合格！");
		return;
	}
	E=new double [k*k];
	g=new double [i];
	E1=new double [k*k];
	e=new double [k*k*i];
	pj=new double [k*i];
	m1=new double [k];
	m2=new double [k];
	m3=new double [k];
	x=new int [k];

	
	for(m=0;m<i;m++)
		for(a=0;a<k;a++)
		{
			pj[m*k+a]=0;
			num=0;
			for(c=0;c<m;c++)
			{
				num+=n[c]*k;
			}
			num+=n[m]*a;
			for(j=0;j<n[m];j++)
				pj[m*k+a]+=master[num+j];
			pj[m*k+a]/=n[m];
	
		}
		
		for(j=0;j<i;j++)
		{
			num=0;
			for(c=0;c<j;c++)
			{
				num+=n[c]*k;
			}
			for(m=0;m<k;m++)
			{
				for(a=0;a<k;a++)
				{
					e[j*k*k+m*k+a]=0;
					for(c=0;c<n[j];c++)
					{
						e[j*k*k+m*k+a]+=(pj[j*k+m]-master[num+n[j]*m+c])*(pj[j*k+a]-master[num+n[j]*a+c]);
					}
					e[j*k*k+m*k+a]/=n[m]-1;
				}
			}
			
		}

///////////////////////////////////
		for(a=0;a<k;a++)
		{
			m_List4.GetText(a,str);
			sscanf(str,"%d",&p2);
			x[a]=p2;
		}
			
///////////////////////////////////////
		num=0;
		for(b=0;b<i;b++)
			num+=n[b];
///////协方差相等//////////////////////////
	if(select)
	{		
		for(j=0;j<k;j++)
		for(m=0;m<k;m++)
		{
			E[j*k+m]=0;
			for(a=0;a<i;a++)
			E[j*k+m]+=e[a*k*k+j*k+m];
		}
	
		if (GetValue(E,k)==0) 
		{
			MessageBox("协方差矩阵不可逆!");
			delete m1;delete m2;delete m3;delete e;
			delete E;delete E1;delete g;delete pj;
			return;
		}
		
		ivy(E,k);
		for(m=0;m<i;m++)
		{
			for(a=0;a<k;a++)
			m1[a]=x[a]-pj[m*k+a];
			
			transpose(m1,m2,1,k);
			mult(m1,E,m3,1,k,k);
			mult(m3,m2,&m4,1,k,1);
			
	        p1=1.0*n[m]/num;
			
			g[m]=-0.5*(m4)+log(p1);

		}
	}
/////////////两类协方差不等/////////
			else
			{
				for(m=0;m<i;m++)
					{
						for(a=0;a<k*k;a++)
						{
							E[a]=e[m*k*k+a];
							E1[a]=e[m*k*k+a];
						}
						if (GetValue(E,k)==0) 
						{
							MessageBox("协方差矩阵不可逆!");
							delete m1;delete m2;delete m3;delete e;
							delete E;delete E1;delete g;delete pj;
							return;
						}							
						ivy(E1,k);

						for(a=0;a<k;a++)
						m1[a]=x[a]-pj[m*k+a];
						
						transpose(m1,m2,1,k);
						mult(m1,E1,m3,1,k,k);
						mult(m3,m2,&m4,1,k,1);
							
						p1=1.0*n[m]/num;
						p3=GetValue(E,k);
						g[m]=-0.5*(m4)+log(p1)-0.5*log(p3);
					}
			}
			result=Max(g,i);
			
			CFont font;
			font.CreateFont(30,12,0,0,700,0,0,0,0,0,0,0,0,"宋体");


			((CStatic*)GetDlgItem(IDC_STATIC_RESULT))->SetFont(&font);	
			if (result!=-1) 
			{
				str.Format("该样本属于第 %d 类",result);
				CWnd::SetDlgItemText(IDC_STATIC_RESULT,str);
			}
			else
			{
			CWnd::SetDlgItemText(IDC_STATIC_RESULT,"无法判别该样本!");
			}
//////////////////////////////////////////////////////////////////////////

			CStatic *out=((CStatic*)GetDlgItem(IDC_STATIC_MAP));
			CDC *pDC=out->GetDC();
			
			
			CRect rect;
			out->GetClientRect(&rect);
			x0=rect.Width()/2;
			y0=rect.Height()/2;
			
			if(k==2&&m_f!=0)
			{	
				move=false;
				pDC->MoveTo(x0-rect.Width()/2,y0);
				pDC->LineTo(x0+rect.Width()/2,y0);
				
				pDC->MoveTo(x0,y0-rect.Height()/2);
				pDC->LineTo(x0,y0+rect.Height()/2);

				CPen redPen,greenPen,bluePen,*oldPen,blackpen,*pen;
				pen=new CPen[i];
				CStatic *out=((CStatic*)GetDlgItem(IDC_STATIC_MAP));
				CDC *pDC=out->GetDC();
				redPen.CreatePen(PS_SOLID,1,RGB(255,0,0));
				greenPen.CreatePen(PS_SOLID,1,RGB(0,255,0));
				bluePen.CreatePen(PS_SOLID,1,RGB(0,0,255));
				blackpen.CreatePen(PS_SOLID,1,RGB(0,0,0));
				for(a=0;a<i;a++)
				{
					if (a%3<1) 
					{
						pen[a].CreatePen(PS_SOLID,1,RGB(10*a%256,150*a%256,255));
					}
					else if (a%3<2) 
					{
						pen[a].CreatePen(PS_SOLID,1,RGB(255,150*a%256,10*a%256));
					}
					else if (a%3<3) 
					{
						pen[a].CreatePen(PS_SOLID,1,RGB(10*a%256,255,150*a%256));
					}
				}
			
				oldPen=pDC->SelectObject(&redPen);
				
					for(b=10;b<rect.Width()-10;b++)
					for(a=10;a<rect.Height()-10;a++)
					{
						c=Select((1.0*a-rect.Width()/2.0)/m_f,(-1.0*b+rect.Height()/2.0)/m_f);

					if (c!=-1)
					{
					
							pDC->SelectObject(&pen[c-1]);
							pDC->MoveTo(a,b);
							pDC->LineTo(a+1,b);
					}
					else
					{
						pDC->SelectObject(&blackpen);
						pDC->MoveTo(a,b);
						pDC->LineTo(a+1,b);
					}

					}
	
					pDC->SelectObject(&blackpen);
					pDC->MoveTo(x0-rect.Width()/2,y0);
					pDC->LineTo(x0+rect.Width()/2,y0);
					
					pDC->MoveTo(x0,y0-rect.Height()/2);
					pDC->LineTo(x0,y0+rect.Height()/2);
					
					pDC->SetBkMode(TRANSPARENT);
					pDC->TextOut(x0+105,y0+2,"X");
					pDC->TextOut(x0+2,y0-118,"Y");


			}
			else if(m_f!=0&&k>=2)
			{
				pDC->SetBkMode(TRANSPARENT);
				pDC->SetTextColor((0,0,255));
				pDC->TextOut(x0-50,y0-50,"注  意！！！");
				pDC->SetTextColor((0,0,255));
				pDC->TextOut(x0-110,y0+10,"程序只提供2特征数的分类界线图");
			}
			out->ReleaseDC(pDC);

//////////////////////////////////////////////////////////////////////////
			delete m1;delete m2;delete m3;delete e;
			delete E;delete E1;delete g;delete pj;delete x;

//	CDialog::OnOK();
}

void CMyDlg::OnButton1() 
{
	// TODO: Add your control notification handler code here
	CFileDialog filedialog(true);
	FILE *fp;
	CString str;
	
	CString pFileName;
	int m,j;