numerical_analysisdlg.cpp

在程序设计部分主要讨论了在读写函数表达式时

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

#include "stdafx.h"
#include "numerical_analysis.h"
#include "numerical_analysisDlg.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()


// Cnumerical_analysisDlg 对话框



Cnumerical_analysisDlg::Cnumerical_analysisDlg(CWnd* pParent /*=NULL*/)
	: CDialog(Cnumerical_analysisDlg::IDD, pParent)
	, m_integration_start(2)
	, m_iter_max_num(16)
	, m_integration_end(5)
	, m_precision(0.0001)
	, m_expression(_T("sin(x*cos(x))*exp(x)"))
	, m_output_filename(_T(""))
	, m_var_name(_T("x"))
	, mb_out_value(TRUE)
	, mb_out_preci(FALSE)
{
	m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}

void Cnumerical_analysisDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	DDX_Text(pDX, IDC_EDIT_START, m_integration_start);
	DDX_Text(pDX, IDC_EDIT_ITER_MAX, m_iter_max_num);
	DDX_Text(pDX, IDC_EDIT_END, m_integration_end);
	DDX_Text(pDX, IDC_EDIT_PRECISION, m_precision);
	DDX_Text(pDX, IDC_EDIT_EXPRESSION, m_expression);
	DDX_Text(pDX, IDC_EDIT_OUT_PUT_FILE, m_output_filename);
	DDX_Control(pDX, IDC_LIST_RESULT, m_result);
	DDX_Text(pDX, IDC_EDIT_VAR_NAME, m_var_name);
	DDX_Control(pDX, IDC_RADIO_OUT_VALUE, CB_out_value);
	DDX_Control(pDX, IDC_COMBO_METHOD, m_method_list);
}

BEGIN_MESSAGE_MAP(Cnumerical_analysisDlg, CDialog)
	ON_WM_SYSCOMMAND()
	ON_WM_PAINT()
	ON_WM_QUERYDRAGICON()
	//}}AFX_MSG_MAP
	ON_BN_CLICKED(IDC_BUTTON_COMPUTE, OnBnClickedButtonCompute)
	ON_BN_CLICKED(IDC_RADIO_OUT_VALUE, OnBnClickedRadioOutValue)
	ON_BN_CLICKED(IDC_RADIO_OUT_PRECI, OnBnClickedRadioOutPreci)
	ON_BN_CLICKED(IDC_BUTTON_OUTPUT_FILE_SELECT, OnBnClickedButtonOutputFileSelect)
	ON_BN_CLICKED(IDC_BUTTON_OUTPUT, OnBnClickedButtonOutput)
	ON_WM_CLOSE()
	ON_BN_CLICKED(IDOK, OnBnClickedOk)
	ON_BN_CLICKED(IDC_BUTTON_GET_VARIABLE, OnBnClickedButtonGetVariable)
	ON_EN_KILLFOCUS(IDC_EDIT_PRECISION, OnEnKillfocusEditPrecision)
	ON_EN_KILLFOCUS(IDC_EDIT_ITER_MAX, OnEnKillfocusEditIterMax)
	ON_EN_KILLFOCUS(IDC_EDIT_EXPRESSION, OnEnKillfocusEditExpression)
END_MESSAGE_MAP()


// Cnumerical_analysisDlg 消息处理程序

BOOL Cnumerical_analysisDlg::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：在此添加额外的初始化代码
	memset(&num_inte_result, 0, sizeof(Iter_result));
	num_inte_result.status = ROMBERG;
	num_inte_result.final_precision = 1e-6;

	memset(&name_list, 0, sizeof(Var_name_list));

	name_list.var_name = m_var_name.GetBuffer();

	m_method_list.InsertString(METHOD_TRAPEZOIDAL, "梯形精度终止迭代");
	m_method_list.InsertString(METHOD_SIMPSON_VALUE, "Simpson精度终止迭代");
	m_method_list.InsertString(METHOD_COTES_VALUE, "Cotes精度终止迭代");
	m_method_list.InsertString(METHOD_ROMBERG_VALUE, "Romberg精度终止迭代");

	m_method_list.SetCurSel(METHOD_ROMBERG_VALUE);
	
	CB_out_value.SetCheck(1);

	
	return TRUE;  // 除非设置了控件的焦点，否则返回 TRUE
}

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

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

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

void Cnumerical_analysisDlg::OnBnClickedButtonCompute()
{
	// TODO: 在此添加控件通知处理程序代码
	UpdateData(true);

	if (analysis_variable() != 0)
	{
		return;
	}

	if (num_inte_result.t_value != NULL)
	{
		free_iter_result(&num_inte_result);
	}

	num_inte_result.integration_start = m_integration_start;
	num_inte_result.integration_end = m_integration_end;
	num_inte_result.final_precision = m_precision;
	num_inte_result.iter_max = m_iter_max_num;

	switch (m_method_list.GetCurSel())
	{
	case METHOD_TRAPEZOIDAL :
		num_inte_result.status = TRAPEZOIDAL;
		break;

	case METHOD_SIMPSON_VALUE :
		num_inte_result.status = SIMPSON;
		break;

	case METHOD_COTES_VALUE :
		num_inte_result.status = COTES;
		break;

	case METHOD_ROMBERG_VALUE :
		num_inte_result.status = ROMBERG;
		break;

	default :
		break;
	}

	int return_value;

	return_value = numerical_integration(m_expression.GetBuffer(), &name_list, &num_inte_result);
	if (return_value != OK)
	{
		MessageBox("计算积分错误！", "计算错误", MB_ICONWARNING | MB_OK);
		return;
	}

	set_user_result_data(num_inte_result);
	
}


void Cnumerical_analysisDlg::set_user_result_data(Iter_result &result)
{
	if (result.iter_t == 0)
	{
		return;
	}

	LVITEM cur_lviterm;

	cur_lviterm.mask = LVIF_TEXT|LVIF_STATE;
	cur_lviterm.iImage = NULL;
	cur_lviterm.state = INDEXTOSTATEIMAGEMASK(1);
	cur_lviterm.stateMask = LVIS_STATEIMAGEMASK;
	cur_lviterm.lParam = 0;
	//cur_lviterm.pszText

	
	m_result.DeleteAllItems();

	//delete all column
	for (; m_result.DeleteColumn(0););

	unsigned long int i;
	CString line_infor;
	
	if (mb_out_value)
	{
		m_result.InsertColumn(LINE_INFOR, "迭代次数", LVCFMT_CENTER, 68);
		m_result.InsertColumn(TRAPEZOIDAL_VALUE, "梯度值", LVCFMT_CENTER, 88);
		m_result.InsertColumn(SIMPSON_VALUE, "simpson值", LVCFMT_CENTER, 88);
		m_result.InsertColumn(COTES_VALUE, "cotes值", LVCFMT_CENTER, 88);
		m_result.InsertColumn(ROMBERG_VALUE, "romberg值", LVCFMT_CENTER, 88);
		for (i = 0; i<result.iter_rom; i++)
		{
			cur_lviterm.iItem = i;

			line_infor.Format(TEXT("%d次迭代"), i+1);
			m_result.InsertItem(i, line_infor);

			line_infor.Format(TEXT("%.4g"), result.t_value[i]);
			cur_lviterm.iSubItem = TRAPEZOIDAL_VALUE;
			cur_lviterm.pszText = line_infor.GetBuffer();
			m_result.SetItem(&cur_lviterm);

			line_infor.Format(TEXT("%.4g"), result.sim_value[i]);
			cur_lviterm.iSubItem = SIMPSON_VALUE;
			cur_lviterm.pszText = line_infor.GetBuffer();
			m_result.SetItem(&cur_lviterm);

			line_infor.Format(TEXT("%.4g"), result.cote_value[i]);
			cur_lviterm.iSubItem = COTES_VALUE;
			cur_lviterm.pszText = line_infor.GetBuffer();
			m_result.SetItem(&cur_lviterm);

			line_infor.Format(TEXT("%.4g"), result.rom_value[i]);
			cur_lviterm.iSubItem = ROMBERG_VALUE;
			cur_lviterm.pszText = line_infor.GetBuffer();
			m_result.SetItem(&cur_lviterm);
		}

		if (result.iter_cote > 0)
		{
			cur_lviterm.iItem = result.iter_cote - 1;
			
			line_infor.Format(TEXT("%d次迭代"), result.iter_cote);
			m_result.InsertItem(result.iter_cote - 1, line_infor);

			line_infor.Format(TEXT("%.4g"), result.t_value[result.iter_cote - 1]);
			cur_lviterm.iSubItem = TRAPEZOIDAL_VALUE;
			cur_lviterm.pszText = line_infor.GetBuffer();
			m_result.SetItem(&cur_lviterm);

			line_infor.Format(TEXT("%.4g"), result.sim_value[result.iter_cote - 1]);
			cur_lviterm.iSubItem = SIMPSON_VALUE;
			cur_lviterm.pszText = line_infor.GetBuffer();
			m_result.SetItem(&cur_lviterm);

			line_infor.Format(TEXT("%.4g"), result.cote_value[result.iter_cote - 1]);
			cur_lviterm.iSubItem = COTES_VALUE;
			cur_lviterm.pszText = line_infor.GetBuffer();
			m_result.SetItem(&cur_lviterm);

		}
		
		if (result.iter_sim > 0)
		{
			cur_lviterm.iItem = result.iter_sim - 1;
			
			line_infor.Format(TEXT("%d次迭代"), result.iter_sim);
			m_result.InsertItem(result.iter_sim - 1, line_infor);

			line_infor.Format(TEXT("%.4g"), result.t_value[result.iter_sim - 1]);
			cur_lviterm.iSubItem = TRAPEZOIDAL_VALUE;
			cur_lviterm.pszText = line_infor.GetBuffer();
			m_result.SetItem(&cur_lviterm);

			line_infor.Format(TEXT("%.4g"), result.sim_value[result.iter_sim - 1]);
			cur_lviterm.iSubItem = SIMPSON_VALUE;
			cur_lviterm.pszText = line_infor.GetBuffer();
			m_result.SetItem(&cur_lviterm);

		
		}

		if (result.iter_t > 0)
		{
			cur_lviterm.iItem = result.iter_t - 1;

			line_infor.Format(TEXT("%d次迭代"), result.iter_t);
			m_result.InsertItem(result.iter_t - 1, line_infor);
			
			line_infor.Format(TEXT("%.4g"), result.t_value[result.iter_t - 1]);
			cur_lviterm.iSubItem = TRAPEZOIDAL_VALUE;
			cur_lviterm.pszText = line_infor.GetBuffer();
			m_result.SetItem(&cur_lviterm);
		}
	
	}
	else 
	{
		m_result.InsertColumn(LINE_INFOR, "迭代次数", LVCFMT_CENTER, 68);
		m_result.InsertColumn(TRAPEZOIDAL_PRECI, "梯度精度", LVCFMT_CENTER, 88);
		m_result.InsertColumn(SIMPSON_PRECI, "simpson精度", LVCFMT_CENTER, 88);
		m_result.InsertColumn(COTES_PRECI, "cotes精度", LVCFMT_CENTER, 88);
		m_result.InsertColumn(ROMBERG_PRECI, "romberg精度", LVCFMT_CENTER, 88);

		for (i = 0; i<result.iter_rom; i++)
		{
			cur_lviterm.iItem = i;

			line_infor.Format(TEXT("%d次迭代"), i+1);
			m_result.InsertItem(i, line_infor);

			line_infor.Format(TEXT("%.4g"), result.t_precision[i]);
			cur_lviterm.iSubItem = TRAPEZOIDAL_PRECI;
			cur_lviterm.pszText = line_infor.GetBuffer();
			m_result.SetItem(&cur_lviterm);

			line_infor.Format(TEXT("%.4g"), result.sim_precision[i]);
			cur_lviterm.iSubItem = SIMPSON_PRECI;
			cur_lviterm.pszText = line_infor.GetBuffer();
			m_result.SetItem(&cur_lviterm);
			line_infor.Format(TEXT("%.4g"), result.cote_precision[i]);
			cur_lviterm.iSubItem = COTES_PRECI;
			cur_lviterm.pszText = line_infor.GetBuffer();
			m_result.SetItem(&cur_lviterm);

			line_infor.Format(TEXT("%.4g"), result.rom_precision[i]);
			cur_lviterm.iSubItem = ROMBERG_PRECI;
			cur_lviterm.pszText = line_infor.GetBuffer();
			m_result.SetItem(&cur_lviterm);

		}

		if (result.iter_cote > 0)
		{