mullinedlg.cpp

最小二乘法的圆拟合算法

// mullineDlg.cpp : 实现文件
//

#include "stdafx.h"
#include "mulline.h"
#include "mullineDlg.h"
#include ".\mullinedlg.h"
#include "rp.h"
#include "math.h"

#ifdef _DEBUG
#define new DEBUG_NEW
#endif
 
// 用于应用程序“关于”菜单项的 CAboutDlg 对话框

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

// 对话框数据
	enum { IDD = IDD_ABOUTBOX };

	protected:
	virtual void DoDataExchange(CDataExchange* pDX);    // DDX/DDV 支持

// 实现
protected:
	DECLARE_MESSAGE_MAP()
};

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

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

BEGIN_MESSAGE_MAP(CAboutDlg, CDialog)
END_MESSAGE_MAP()


// CmullineDlg 对话框



CmullineDlg::CmullineDlg(CWnd* pParent /*=NULL*/)
	: CDialog(CmullineDlg::IDD, pParent)
	, m_editnum(0)
	, m_editx(0)
	, m_edity(0)
	, m_editposx(0)
	, m_editposy(0)
	, m_editposz(0)
	, m_editr(0)
	, m_editd(0)
{
	m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}

void CmullineDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	DDX_Text(pDX, IDC_EDIT1, m_editnum);
	DDX_Text(pDX, IDC_EDIT7, m_editx);
	DDX_Text(pDX, IDC_EDIT8, m_edity);
	DDX_Text(pDX, IDC_EDIT2, m_editposx);
	DDX_Text(pDX, IDC_EDIT3, m_editposy);
	DDX_Text(pDX, IDC_EDIT4, m_editposz);
	DDX_Text(pDX, IDC_EDIT5, m_editr);
	DDX_Text(pDX, IDC_EDIT6, m_editd);
}

BEGIN_MESSAGE_MAP(CmullineDlg, CDialog)
	ON_WM_SYSCOMMAND()
	ON_WM_PAINT()
	ON_WM_QUERYDRAGICON()
	//}}AFX_MSG_MAP
	ON_EN_CHANGE(IDC_EDIT5, OnEnChangeEdit5)
	ON_BN_CLICKED(IDC_BUTTON1, OnBnClickedButton1)
	ON_BN_CLICKED(IDC_BUTTON2, OnBnClickedButton2)
	ON_BN_CLICKED(IDOK, OnBnClickedOk)
END_MESSAGE_MAP()


// CmullineDlg 消息处理程序

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

	// 将\“关于...\”菜单项添加到系统菜单中。

	// IDM_ABOUTBOX 必须在系统命令范围内。
	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);
		}
	}

	// 设置此对话框的图标。当应用程序主窗口不是对话框时，框架将自动
	//  执行此操作
	SetIcon(m_hIcon, TRUE);			// 设置大图标
	SetIcon(m_hIcon, FALSE);		// 设置小图标

	// TODO: 在此添加额外的初始化代码
	
	return TRUE;  // 除非设置了控件的焦点，否则返回 TRUE
}

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

// 如果向对话框添加最小化按钮，则需要下面的代码
//  来绘制该图标。对于使用文档/视图模型的 MFC 应用程序，
//  这将由框架自动完成。

void CmullineDlg::OnPaint() 
{
	if (IsIconic())
	{
		CPaintDC dc(this); // 用于绘制的设备上下文

		SendMessage(WM_ICONERASEBKGND, reinterpret_cast<WPARAM>(dc.GetSafeHdc()), 0);

		// 使图标在工作矩形中居中
		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;

		// 绘制图标
		dc.DrawIcon(x, y, m_hIcon);
	}
	else
	{
		CDialog::OnPaint();
	}
}

//当用户拖动最小化窗口时系统调用此函数取得光标显示。
HCURSOR CmullineDlg::OnQueryDragIcon()
{
	return static_cast<HCURSOR>(m_hIcon);
}

void CmullineDlg::OnEnChangeEdit5()
{
	// TODO:  如果该控件是 RICHEDIT 控件，则它将不会
	// 发送该通知，除非重写 CDialog::OnInitDialog()
	// 函数并调用 CRichEditCtrl().SetEventMask()，
	// 同时将 ENM_CHANGE 标志“或”运算到掩码中。

	// TODO:  在此添加控件通知处理程序代码
}

void CmullineDlg::OnBnClickedButton1()
{
	UpdateData(1);
	rp::num=m_editnum;
	rp::n=m_editnum;
	// TODO: 在此添加控件通知处理程序代码
}

void CmullineDlg::OnBnClickedButton2()
{
	op1=0;
	op2=0;
	op3=0;
	op4=0;
	op5=0;
	op6=0;
	op7=0;
	op8=0;
	op9=0;
	d1=0;
	d2=0;
	d3=0;
	d4=0;
	d5=0;
	d6=0;
	d7=0;
	//numi=0;
	r=0;
	UpdateData(1);
	if(rp::n<3)
	{
		AfxMessageBox("can't");
		
	}
	x[rp::k]=m_editx;
	y[rp::k]=m_edity;
	rp::k++;
	rp::num--;
	if(rp::num==0)
	{
		for(i=0;i<rp::n;i++)
		{
			op1+=x[i];
			op2+=x[i]*x[i];
			op3+=x[i]*x[i]*x[i];
			op4+=y[i];
		    op5+=y[i]*y[i];
			op6+=y[i]*y[i]*y[i];
			op7+=x[i]*y[i];
			op8+=x[i]*x[i]*y[i];
			op9+=x[i]*y[i]*y[i];
		}
		numi=(1.0 / rp::n);
		d1=op2-numi*op1*op1;
		d2=op7-numi*op1*op4;
		d3=op5-numi*op4*op4;
		d4=op3-numi*op1*op2;
		d5=op9-numi*op1*op5;
		d6=op8-numi*op2*op4;
		d7=op6-numi*op4*op5;
		if(((d4+d5)*d3-(d6+d7)*d2)<0.0000001)
			Xc=0;
		else
            Xc=((d4+d5)*d3-(d6+d7)*d2)/(2*(d1*d3-d2*d2));
		if(((d6+d7)*(op2-(1/rp::n)*op1*op1)-(d4+d5)*d2)<0.0000001)
			Yc=0;
		else
		    Yc=((d6+d7)*(op2-(1/rp::n)*op1*op1)-(d4+d5)*d2)/(2*(d1*d3-d2*d2));
		for(i=0;i<rp::n;i++)
		{
			r+=(x[i]-Xc)*(x[i]-Xc)+(y[i]-Yc)*(y[i]-Yc);
		}
		Rc=sqrt(r/rp::n);
		if((d1*d3-d2*d2)==0)
		{
			AfxMessageBox("can't");

		}
		
     
	  m_editposx=Xc;
	  m_editposy=Yc;
	  m_editr=Rc;
	  m_editd=2*Rc;
	  UpdateData(0);
	}
	

	// TODO: 在此添加控件通知处理程序代码
}

void CmullineDlg::OnBnClickedOk()
{
	// TODO: 在此添加控件通知处理程序代码
	OnOK();
}