interpolaterdlg.cpp

插值类

// InterpolaterDlg.cpp : implementation file
//

#include "stdafx.h"
#include "Interpolater.h"
#include "InterpolaterDlg.h"

#include "Interpolate.h"
#include "Spline3Dlg.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()

/////////////////////////////////////////////////////////////////////////////
// CInterpolaterDlg dialog

CInterpolaterDlg::CInterpolaterDlg(CWnd* pParent /*=NULL*/)
	: CDialog(CInterpolaterDlg::IDD, pParent)
{
	//{{AFX_DATA_INIT(CInterpolaterDlg)
	m_strXi = _T("");
	m_strYi = _T("");
	m_nNumNodes = 0;
	m_nNumTargetNodes = 1;
	m_strTargetXi = _T("");
	m_strStepOrBorderCondition = _T("");
	//}}AFX_DATA_INIT
	// Note that LoadIcon does not require a subsequent DestroyIcon in Win32
	m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}

void CInterpolaterDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CInterpolaterDlg)
	DDX_Text(pDX, IDC_EDIT1, m_strXi);
	DDX_Text(pDX, IDC_EDIT2, m_strYi);
	DDX_Text(pDX, IDC_EDIT3, m_nNumNodes);
	DDX_Text(pDX, IDC_EDIT5, m_nNumTargetNodes);
	DDX_Text(pDX, IDC_EDIT6, m_strTargetXi);
	DDX_Text(pDX, IDC_EDIT4, m_strStepOrBorderCondition);
	//}}AFX_DATA_MAP
}

BEGIN_MESSAGE_MAP(CInterpolaterDlg, CDialog)
	//{{AFX_MSG_MAP(CInterpolaterDlg)
	ON_WM_SYSCOMMAND()
	ON_WM_PAINT()
	ON_WM_QUERYDRAGICON()
	ON_BN_CLICKED(IDC_BUTTON1, OnButton1)
	ON_BN_CLICKED(IDC_BUTTON2, OnButton2)
	ON_BN_CLICKED(IDC_BUTTON3, OnButton3)
	ON_BN_CLICKED(IDC_BUTTON4, OnButton4)
	ON_BN_CLICKED(IDC_BUTTON5, OnButton5)
	ON_BN_CLICKED(IDC_BUTTON6, OnButton6)
	ON_BN_CLICKED(IDC_BUTTON7, OnButton7)
	ON_BN_CLICKED(IDC_BUTTON8, OnButton8)
	ON_BN_CLICKED(IDC_BUTTON9, OnButton9)
	ON_BN_CLICKED(IDC_BUTTON10, OnButton10)
	ON_BN_CLICKED(IDC_BUTTON11, OnButton11)
	ON_BN_CLICKED(IDC_BUTTON12, OnButton12)
	ON_BN_CLICKED(IDC_BUTTON13, OnButton13)
	ON_BN_CLICKED(IDC_BUTTON14, OnButton14)
	ON_BN_CLICKED(IDC_BUTTON15, OnButton15)
	ON_BN_CLICKED(IDC_BUTTON16, OnButton16)
	ON_BN_CLICKED(IDC_BUTTON17, OnButton17)
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CInterpolaterDlg message handlers

BOOL CInterpolaterDlg::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 CInterpolaterDlg::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 CInterpolaterDlg::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 CInterpolaterDlg::OnQueryDragIcon()
{
	return (HCURSOR) m_hIcon;
}

/*
10

0.10,0.15,0.25,0.40,0.50,0.57,0.70,0.85,0.93,1.00

0.904837,0.860708,0.778801,0.670320,0.606531,0.565525,0.496585,0.427415,0.394554,0.367879

1

0.63
*/
// 一元全区间不等距插值
void CInterpolaterDlg::OnButton1() 
{
	// 获取数据
	UpdateData();
	double *x = new double[m_nNumNodes];
	double *y = new double[m_nNumNodes];
	CInterpolate::GetNodesFromString(m_strXi, m_nNumNodes, x, ",");
	CInterpolate::GetNodesFromString(m_strYi, m_nNumNodes, y, ",");

	// 插值运算
	double *t = new double[m_nNumTargetNodes];
	CInterpolate::GetNodesFromString(m_strTargetXi, m_nNumTargetNodes, t, ",");

	// 显示结果
	m_strTargetYi = "";
	for (int i=0; i<m_nNumTargetNodes; ++i)
	{
		double yt = CInterpolate::GetValueLagrange(m_nNumNodes, x, y, t[i]);
		CString s;
		s.Format("f(%f) = %f\r\n", t[i], yt);
		m_strTargetYi += s;
	}
	AfxMessageBox(m_strTargetYi, MB_OK|MB_ICONINFORMATION);

	// 释放内存
	delete[] t;
	delete[] x;
	delete[] y;
}

/*
10

0.1

0.1;

0.904837,0.818731,0.740818,0.670320,0.606531,0.548812,0.496585,0.449329,0.406570,0.367879

3

0.25,0.63,0.95
*/
// 一元全区间等距插值
void CInterpolaterDlg::OnButton2() 
{
	// 获取数据
	UpdateData();
	double x0;
	double *y = new double[m_nNumNodes];
	CInterpolate::GetNodesFromString(m_strXi, 1, &x0, ",");
	CInterpolate::GetNodesFromString(m_strYi, m_nNumNodes, y, ",");

	// 插值运算
	double *t = new double[m_nNumTargetNodes];
	CInterpolate::GetNodesFromString(m_strTargetXi, m_nNumTargetNodes, t, ",");

	// 显示结果
	double dblStep = atof(m_strStepOrBorderCondition);
	m_strTargetYi = "";
	for (int i=0; i<m_nNumTargetNodes; ++i)
	{
		double yt = CInterpolate::GetValueLagrange(m_nNumNodes, x0, dblStep, y, t[i]);
		CString s;
		s.Format("f(%f) = %f\r\n", t[i], yt);
		m_strTargetYi += s;
	}
	AfxMessageBox(m_strTargetYi, MB_OK|MB_ICONINFORMATION);

	// 释放内存
	delete[] t;
	delete[] y;
}

/*
5

1.615,1.634,1.702,1.828,1.921

2.41450,2.46459,2.65271,3.03035,3.34066

2

1.682,1.813
*/
// 一元三点不等距插值
void CInterpolaterDlg::OnButton3() 
{
	// 获取数据
	UpdateData();
	double *x = new double[m_nNumNodes];
	double *y = new double[m_nNumNodes];
	CInterpolate::GetNodesFromString(m_strXi, m_nNumNodes, x, ",");
	CInterpolate::GetNodesFromString(m_strYi, m_nNumNodes, y, ",");

	// 插值运算
	double *t = new double[m_nNumTargetNodes];
	CInterpolate::GetNodesFromString(m_strTargetXi, m_nNumTargetNodes, t, ",");

	// 显示结果
	m_strTargetYi = "";
	for (int i=0; i<m_nNumTargetNodes; ++i)
	{
		double yt = CInterpolate::GetValueLagrange3(m_nNumNodes, x, y, t[i]);
		CString s;
		s.Format("f(%f) = %f\r\n", t[i], yt);
		m_strTargetYi += s;
	}
	AfxMessageBox(m_strTargetYi, MB_OK|MB_ICONINFORMATION);

	// 释放内存
	delete[] t;
	delete[] x;
	delete[] y;
}

/*
10

0.1

0.1

0.904837,0.818731,0.740818,0.670320,0.606531,0.548812,0.496585,0.449329,0.406570,0.367879

3

0.23,0.63,0.95
*/
// 一元三点等距插值
void CInterpolaterDlg::OnButton4() 
{
	// 获取数据
	UpdateData();
	double x0;
	double *y = new double[m_nNumNodes];
	CInterpolate::GetNodesFromString(m_strXi, 1, &x0, ",");
	CInterpolate::GetNodesFromString(m_strYi, m_nNumNodes, y, ",");

	// 插值运算
	double *t = new double[m_nNumTargetNodes];
	CInterpolate::GetNodesFromString(m_strTargetXi, m_nNumTargetNodes, t, ",");

	// 显示结果
	double dblStep = atof(m_strStepOrBorderCondition);
	m_strTargetYi = "";
	for (int i=0; i<m_nNumTargetNodes; ++i)
	{
		double yt = CInterpolate::GetValueLagrange3(m_nNumNodes, x0, dblStep, y, t[i]);
		CString s;
		s.Format("f(%f) = %f\r\n", t[i], yt);
		m_strTargetYi += s;
	}
	AfxMessageBox(m_strTargetYi, MB_OK|MB_ICONINFORMATION);

	// 释放内存
	delete[] t;
	delete[] y;
}

/*
10

-1.0,-0.8,-0.65,-0.4,-0.3,0.0,0.2,0.45,0.8,1.0

0.0384615,0.0588236,0.0864865,0.2,0.307692,1.0,0.5,0.164948,0.0588236,0.0384615

2

0.85,0.25
*/
// 连分式不等距插值
void CInterpolaterDlg::OnButton5() 
{
	// 获取数据
	UpdateData();
	double *x = new double[m_nNumNodes];
	double *y = new double[m_nNumNodes];
	CInterpolate::GetNodesFromString(m_strXi, m_nNumNodes, x, ",");
	CInterpolate::GetNodesFromString(m_strYi, m_nNumNodes, y, ",");

	// 插值运算
	double *t = new double[m_nNumTargetNodes];
	CInterpolate::GetNodesFromString(m_strTargetXi, m_nNumTargetNodes, t, ",");

	// 显示结果
	m_strTargetYi = "";
	for (int i=0; i<m_nNumTargetNodes; ++i)
	{
		double yt = CInterpolate::GetValuePqs(m_nNumNodes, x, y, t[i]);
		CString s;
		s.Format("f(%f) = %f\r\n", t[i], yt);
		m_strTargetYi += s;
	}
	AfxMessageBox(m_strTargetYi, MB_OK|MB_ICONINFORMATION);

	// 释放内存
	delete[] t;
	delete[] x;
	delete[] y;
}

/*
11

0.2

-1.0

0.0384615,0.0588236,0.1,0.2,0.5,1.0,0.5,0.2,0.1,0.0588236,0.0384615

2

-0.75, -0.05
*/
// 连分式等距插值
void CInterpolaterDlg::OnButton6() 
{
	// 获取数据
	UpdateData();
	double x0;
	double *y = new double[m_nNumNodes];
	CInterpolate::GetNodesFromString(m_strXi, 1, &x0, ",");
	CInterpolate::GetNodesFromString(m_strYi, m_nNumNodes, y, ",");

	// 插值运算
	double *t = new double[m_nNumTargetNodes];
	CInterpolate::GetNodesFromString(m_strTargetXi, m_nNumTargetNodes, t, ",");

	// 显示结果
	double dblStep = atof(m_strStepOrBorderCondition);
	m_strTargetYi = "";
	for (int i=0; i<m_nNumTargetNodes; ++i)
	{
		double yt = CInterpolate::GetValuePqs(m_nNumNodes, x0, dblStep, y, t[i]);
		CString s;
		s.Format("f(%f) = %f\r\n", t[i], yt);
		m_strTargetYi += s;
	}
	AfxMessageBox(m_strTargetYi, MB_OK|MB_ICONINFORMATION);

	// 释放内存
	delete[] t;
	delete[] y;
}

/*
10

0.1,0.15,0.3,0.45,0.55,0.6,0.7,0.85,0.9,1.0

0.904837,0.860708,0.740818,0.637628,0.576950,0.548812,0.496585,0.427415,0.406570,0.367879

1

0.356
*/
// 埃尔米特不等距插值
void CInterpolaterDlg::OnButton7() 
{
	// 获取数据
	UpdateData();
	double *x = new double[m_nNumNodes];
	double *y = new double[m_nNumNodes];
	CInterpolate::GetNodesFromString(m_strXi, m_nNumNodes, x, ",");
	CInterpolate::GetNodesFromString(m_strYi, m_nNumNodes, y, ",");

	double *dy = new double[m_nNumNodes];
	for (int i=0; i<m_nNumNodes; ++i)
		dy[i]=-y[i];

	// 插值运算
	double *t = new double[m_nNumTargetNodes];
	CInterpolate::GetNodesFromString(m_strTargetXi, m_nNumTargetNodes, t, ",");

	// 显示结果
	m_strTargetYi = "";
	for (i=0; i<m_nNumTargetNodes; ++i)
	{
		double yt = CInterpolate::GetValueHermite(m_nNumNodes, x, y, dy, t[i]);