?? 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]);
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -