newtondlg.cpp

最优化方法newton法求极值的vc++实现 环境：vs2008

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

#include "stdafx.h"
#include "Newton.h"
#include "NewtonDlg.h"

#include "qd.h"
#include "Calculation.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()


// CNewtonDlg 对话框




CNewtonDlg::CNewtonDlg(CWnd* pParent /*=NULL*/)
	: CDialog(CNewtonDlg::IDD, pParent)
	, m_Exprssion(_T("x*x*x*x-4*x*x*x-6*x*x-16*x+4"))
	, m_Error(0.0001)
	, m_FirstValue(0)
	, m_Times(0)
	, m_OutValue(0)
	, m_OutValue2(0)
{
	m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}

void CNewtonDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	DDX_Text(pDX, IDC_EDIT1, m_Exprssion);
	DDX_Text(pDX, IDC_EDIT2, m_Error);
	DDX_Text(pDX, IDC_EDIT3, m_FirstValue);
	DDX_Text(pDX, IDC_EDIT4, m_Times);
	DDX_Text(pDX, IDC_EDIT5, m_OutValue);
	DDX_Text(pDX, IDC_EDIT6, m_OutValue2);
}

BEGIN_MESSAGE_MAP(CNewtonDlg, CDialog)
	ON_WM_SYSCOMMAND()
	ON_WM_PAINT()
	ON_WM_QUERYDRAGICON()
	//}}AFX_MSG_MAP
	ON_BN_CLICKED(IDOK, &CNewtonDlg::OnBnClickedOk)
END_MESSAGE_MAP()


// CNewtonDlg 消息处理程序

BOOL CNewtonDlg::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 CNewtonDlg::OnSysCommand(UINT nID, LPARAM lParam)
{
	if ((nID & 0xFFF0) == IDM_ABOUTBOX)
	{
		CAboutDlg dlgAbout;
		dlgAbout.DoModal();
	}
	else
	{
		CDialog::OnSysCommand(nID, lParam);
	}
}

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

void CNewtonDlg::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 CNewtonDlg::OnQueryDragIcon()
{
	return static_cast<HCURSOR>(m_hIcon);
}

CCalculation cacu;
void CNewtonDlg::OnBnClickedOk()
{
	// TODO: 在此添加控件通知处理程序代码
	this->UpdateData();
	m_Times=0;
	double f1,f2,firstvalue;
	CString ef1,ef2,t;
	ef1=GetExpretion(m_Exprssion);
	ef2=GetExpretion(ef1);
	firstvalue=m_FirstValue;
	f1=cacu.StoN(cacu.MainPro(ef1,firstvalue));
	f2=cacu.StoN(cacu.MainPro(ef2,firstvalue));
	m_OutValue=firstvalue-f1/f2;
	if(Abs(m_OutValue-firstvalue)<=m_Error)
	{
		m_Times=0;
		m_OutValue2=cacu.StoN(cacu.MainPro(m_Exprssion,m_OutValue));
	}
	else
	{
		while(Abs(m_OutValue-firstvalue)>m_Error)
		{
			firstvalue=m_OutValue;
			t.Format("%f",firstvalue);
			firstvalue=cacu.StoN(t);
			f1=cacu.StoN(cacu.MainPro(ef1,firstvalue));
			f2=cacu.StoN(cacu.MainPro(ef2,firstvalue));
			m_OutValue=firstvalue-f1/f2;
			//firstvalue=m_OutValue;
			m_Times++;
			
		}
	}
		
		t.Format("%f",m_OutValue);
		m_OutValue=cacu.StoN(t);
		m_OutValue2=cacu.StoN(cacu.MainPro(m_Exprssion,m_OutValue));
	
		UpdateData(FALSE);	
	//OnOK();
}