开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 数值算法/人工智能 > 数学计算 > 查看源码
InterpolaterDlg.cpp
资源名称：科学与工程数值计算算法配套源码vc++.zip　[点击查看]
上传用户：weigute
上传日期：2007-03-02
资源大小：1287k
文件大小：22k
源码类别：

数学计算

开发平台：

Visual C++

InterpolaterDlg.cpp：源码内容
  1. // InterpolaterDlg.cpp : implementation file
  2. //
  4. #include "stdafx.h"
  5. #include "Interpolater.h"
  6. #include "InterpolaterDlg.h"
  8. #include "Interpolate.h"
  9. #include "Spline3Dlg.h"
  11. #ifdef _DEBUG
  12. #define new DEBUG_NEW
  13. #undef THIS_FILE
  14. static char THIS_FILE[] = __FILE__;
  15. #endif
  17. /////////////////////////////////////////////////////////////////////////////
  18. // CAboutDlg dialog used for App About
  20. class CAboutDlg : public CDialog
  21. {
  22. public:
  23. CAboutDlg();
  25. // Dialog Data
  26. //{{AFX_DATA(CAboutDlg)
  27. enum { IDD = IDD_ABOUTBOX };
  28. //}}AFX_DATA
  30. // ClassWizard generated virtual function overrides
  31. //{{AFX_VIRTUAL(CAboutDlg)
  32. protected:
  33. virtual void DoDataExchange(CDataExchange* pDX);    // DDX/DDV support
  34. //}}AFX_VIRTUAL
  36. // Implementation
  37. protected:
  38. //{{AFX_MSG(CAboutDlg)
  39. //}}AFX_MSG
  40. DECLARE_MESSAGE_MAP()
  41. };
  43. CAboutDlg::CAboutDlg() : CDialog(CAboutDlg::IDD)
  44. {
  45. //{{AFX_DATA_INIT(CAboutDlg)
  46. //}}AFX_DATA_INIT
  47. }
  49. void CAboutDlg::DoDataExchange(CDataExchange* pDX)
  50. {
  51. CDialog::DoDataExchange(pDX);
  52. //{{AFX_DATA_MAP(CAboutDlg)
  53. //}}AFX_DATA_MAP
  54. }
  56. BEGIN_MESSAGE_MAP(CAboutDlg, CDialog)
  57. //{{AFX_MSG_MAP(CAboutDlg)
  58. // No message handlers
  59. //}}AFX_MSG_MAP
  60. END_MESSAGE_MAP()
  62. /////////////////////////////////////////////////////////////////////////////
  63. // CInterpolaterDlg dialog
  65. CInterpolaterDlg::CInterpolaterDlg(CWnd* pParent /*=NULL*/)
  66. : CDialog(CInterpolaterDlg::IDD, pParent)
  67. {
  68. //{{AFX_DATA_INIT(CInterpolaterDlg)
  69. m_strXi = _T("");
  70. m_strYi = _T("");
  71. m_nNumNodes = 0;
  72. m_nNumTargetNodes = 1;
  73. m_strTargetXi = _T("");
  74. m_strStepOrBorderCondition = _T("");
  75. //}}AFX_DATA_INIT
  76. // Note that LoadIcon does not require a subsequent DestroyIcon in Win32
  77. m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
  78. }
  80. void CInterpolaterDlg::DoDataExchange(CDataExchange* pDX)
  81. {
  82. CDialog::DoDataExchange(pDX);
  83. //{{AFX_DATA_MAP(CInterpolaterDlg)
  84. DDX_Text(pDX, IDC_EDIT1, m_strXi);
  85. DDX_Text(pDX, IDC_EDIT2, m_strYi);
  86. DDX_Text(pDX, IDC_EDIT3, m_nNumNodes);
  87. DDX_Text(pDX, IDC_EDIT5, m_nNumTargetNodes);
  88. DDX_Text(pDX, IDC_EDIT6, m_strTargetXi);
  89. DDX_Text(pDX, IDC_EDIT4, m_strStepOrBorderCondition);
  90. //}}AFX_DATA_MAP
  91. }
  93. BEGIN_MESSAGE_MAP(CInterpolaterDlg, CDialog)
  94. //{{AFX_MSG_MAP(CInterpolaterDlg)
  95. ON_WM_SYSCOMMAND()
  96. ON_WM_PAINT()
  97. ON_WM_QUERYDRAGICON()
  98. ON_BN_CLICKED(IDC_BUTTON1, OnButton1)
  99. ON_BN_CLICKED(IDC_BUTTON2, OnButton2)
  100. ON_BN_CLICKED(IDC_BUTTON3, OnButton3)
  101. ON_BN_CLICKED(IDC_BUTTON4, OnButton4)
  102. ON_BN_CLICKED(IDC_BUTTON5, OnButton5)
  103. ON_BN_CLICKED(IDC_BUTTON6, OnButton6)
  104. ON_BN_CLICKED(IDC_BUTTON7, OnButton7)
  105. ON_BN_CLICKED(IDC_BUTTON8, OnButton8)
  106. ON_BN_CLICKED(IDC_BUTTON9, OnButton9)
  107. ON_BN_CLICKED(IDC_BUTTON10, OnButton10)
  108. ON_BN_CLICKED(IDC_BUTTON11, OnButton11)
  109. ON_BN_CLICKED(IDC_BUTTON12, OnButton12)
  110. ON_BN_CLICKED(IDC_BUTTON13, OnButton13)
  111. ON_BN_CLICKED(IDC_BUTTON14, OnButton14)
  112. ON_BN_CLICKED(IDC_BUTTON15, OnButton15)
  113. ON_BN_CLICKED(IDC_BUTTON16, OnButton16)
  114. ON_BN_CLICKED(IDC_BUTTON17, OnButton17)
  115. //}}AFX_MSG_MAP
  116. END_MESSAGE_MAP()
  118. /////////////////////////////////////////////////////////////////////////////
  119. // CInterpolaterDlg message handlers
  121. BOOL CInterpolaterDlg::OnInitDialog()
  122. {
  123. CDialog::OnInitDialog();
  125. // Add "About..." menu item to system menu.
  127. // IDM_ABOUTBOX must be in the system command range.
  128. ASSERT((IDM_ABOUTBOX & 0xFFF0) == IDM_ABOUTBOX);
  129. ASSERT(IDM_ABOUTBOX < 0xF000);
  131. CMenu* pSysMenu = GetSystemMenu(FALSE);
  132. if (pSysMenu != NULL)
  133. {
  134. CString strAboutMenu;
  135. strAboutMenu.LoadString(IDS_ABOUTBOX);
  136. if (!strAboutMenu.IsEmpty())
  137. {
  138. pSysMenu->AppendMenu(MF_SEPARATOR);
  139. pSysMenu->AppendMenu(MF_STRING, IDM_ABOUTBOX, strAboutMenu);
  140. }
  141. }
  143. // Set the icon for this dialog.  The framework does this automatically
  144. //  when the application's main window is not a dialog
  145. SetIcon(m_hIcon, TRUE); // Set big icon
  146. SetIcon(m_hIcon, FALSE); // Set small icon
  148. // TODO: Add extra initialization here
  150. return TRUE;  // return TRUE  unless you set the focus to a control
  151. }
  153. void CInterpolaterDlg::OnSysCommand(UINT nID, LPARAM lParam)
  154. {
  155. if ((nID & 0xFFF0) == IDM_ABOUTBOX)
  156. {
  157. CAboutDlg dlgAbout;
  158. dlgAbout.DoModal();
  159. }
  160. else
  161. {
  162. CDialog::OnSysCommand(nID, lParam);
  163. }
  164. }
  166. // If you add a minimize button to your dialog, you will need the code below
  167. //  to draw the icon.  For MFC applications using the document/view model,
  168. //  this is automatically done for you by the framework.
  170. void CInterpolaterDlg::OnPaint() 
  171. {
  172. if (IsIconic())
  173. {
  174. CPaintDC dc(this); // device context for painting
  176. SendMessage(WM_ICONERASEBKGND, (WPARAM) dc.GetSafeHdc(), 0);
  178. // Center icon in client rectangle
  179. int cxIcon = GetSystemMetrics(SM_CXICON);
  180. int cyIcon = GetSystemMetrics(SM_CYICON);
  181. CRect rect;
  182. GetClientRect(&rect);
  183. int x = (rect.Width() - cxIcon + 1) / 2;
  184. int y = (rect.Height() - cyIcon + 1) / 2;
  186. // Draw the icon
  187. dc.DrawIcon(x, y, m_hIcon);
  188. }
  189. else
  190. {
  191. CDialog::OnPaint();
  192. }
  193. }
  195. // The system calls this to obtain the cursor to display while the user drags
  196. //  the minimized window.
  197. HCURSOR CInterpolaterDlg::OnQueryDragIcon()
  198. {
  199. return (HCURSOR) m_hIcon;
  200. }
  202. /*
  203. 10
  205. 0.10,0.15,0.25,0.40,0.50,0.57,0.70,0.85,0.93,1.00
  207. 0.904837,0.860708,0.778801,0.670320,0.606531,0.565525,0.496585,0.427415,0.394554,0.367879
  209. 1
  211. 0.63
  212. */
  213. // 一元全区间不等距插值
  214. void CInterpolaterDlg::OnButton1() 
  215. {
  216. // 获取数据
  217. UpdateData();
  218. double *x = new double[m_nNumNodes];
  219. double *y = new double[m_nNumNodes];
  220. CInterpolate::GetNodesFromString(m_strXi, m_nNumNodes, x, ",");
  221. CInterpolate::GetNodesFromString(m_strYi, m_nNumNodes, y, ",");
  223. // 插值运算
  224. double *t = new double[m_nNumTargetNodes];
  225. CInterpolate::GetNodesFromString(m_strTargetXi, m_nNumTargetNodes, t, ",");
  227. // 显示结果
  228. m_strTargetYi = "";
  229. for (int i=0; i<m_nNumTargetNodes; ++i)
  230. {
  231. double yt = CInterpolate::GetValueLagrange(m_nNumNodes, x, y, t[i]);
  232. CString s;
  233. s.Format("f(%f) = %frn", t[i], yt);
  234. m_strTargetYi += s;
  235. }
  236. AfxMessageBox(m_strTargetYi, MB_OK|MB_ICONINFORMATION);
  238. // 释放内存
  239. delete[] t;
  240. delete[] x;
  241. delete[] y;
  242. }
  244. /*
  245. 10
  247. 0.1
  249. 0.1;
  251. 0.904837,0.818731,0.740818,0.670320,0.606531,0.548812,0.496585,0.449329,0.406570,0.367879
  253. 3
  255. 0.25,0.63,0.95
  256. */
  257. // 一元全区间等距插值
  258. void CInterpolaterDlg::OnButton2() 
  259. {
  260. // 获取数据
  261. UpdateData();
  262. double x0;
  263. double *y = new double[m_nNumNodes];
  264. CInterpolate::GetNodesFromString(m_strXi, 1, &x0, ",");
  265. CInterpolate::GetNodesFromString(m_strYi, m_nNumNodes, y, ",");
  267. // 插值运算
  268. double *t = new double[m_nNumTargetNodes];
  269. CInterpolate::GetNodesFromString(m_strTargetXi, m_nNumTargetNodes, t, ",");
  271. // 显示结果
  272. double dblStep = atof(m_strStepOrBorderCondition);
  273. m_strTargetYi = "";
  274. for (int i=0; i<m_nNumTargetNodes; ++i)
  275. {
  276. double yt = CInterpolate::GetValueLagrange(m_nNumNodes, x0, dblStep, y, t[i]);
  277. CString s;
  278. s.Format("f(%f) = %frn", t[i], yt);
  279. m_strTargetYi += s;
  280. }
  281. AfxMessageBox(m_strTargetYi, MB_OK|MB_ICONINFORMATION);
  283. // 释放内存
  284. delete[] t;
  285. delete[] y;
  286. }
  288. /*
  289. 5
  291. 1.615,1.634,1.702,1.828,1.921
  293. 2.41450,2.46459,2.65271,3.03035,3.34066
  295. 2
  297. 1.682,1.813
  298. */
  299. // 一元三点不等距插值
  300. void CInterpolaterDlg::OnButton3() 
  301. {
  302. // 获取数据
  303. UpdateData();
  304. double *x = new double[m_nNumNodes];
  305. double *y = new double[m_nNumNodes];
  306. CInterpolate::GetNodesFromString(m_strXi, m_nNumNodes, x, ",");
  307. CInterpolate::GetNodesFromString(m_strYi, m_nNumNodes, y, ",");
  309. // 插值运算
  310. double *t = new double[m_nNumTargetNodes];
  311. CInterpolate::GetNodesFromString(m_strTargetXi, m_nNumTargetNodes, t, ",");
  313. // 显示结果
  314. m_strTargetYi = "";
  315. for (int i=0; i<m_nNumTargetNodes; ++i)
  316. {
  317. double yt = CInterpolate::GetValueLagrange3(m_nNumNodes, x, y, t[i]);
  318. CString s;
  319. s.Format("f(%f) = %frn", t[i], yt);
  320. m_strTargetYi += s;
  321. }
  322. AfxMessageBox(m_strTargetYi, MB_OK|MB_ICONINFORMATION);
  324. // 释放内存
  325. delete[] t;
  326. delete[] x;
  327. delete[] y;
  328. }
  330. /*
  331. 10
  333. 0.1
  335. 0.1
  337. 0.904837,0.818731,0.740818,0.670320,0.606531,0.548812,0.496585,0.449329,0.406570,0.367879
  339. 3
  341. 0.23,0.63,0.95
  342. */
  343. // 一元三点等距插值
  344. void CInterpolaterDlg::OnButton4() 
  345. {
  346. // 获取数据
  347. UpdateData();
  348. double x0;
  349. double *y = new double[m_nNumNodes];
  350. CInterpolate::GetNodesFromString(m_strXi, 1, &x0, ",");
  351. CInterpolate::GetNodesFromString(m_strYi, m_nNumNodes, y, ",");
  353. // 插值运算
  354. double *t = new double[m_nNumTargetNodes];
  355. CInterpolate::GetNodesFromString(m_strTargetXi, m_nNumTargetNodes, t, ",");
  357. // 显示结果
  358. double dblStep = atof(m_strStepOrBorderCondition);
  359. m_strTargetYi = "";
  360. for (int i=0; i<m_nNumTargetNodes; ++i)
  361. {
  362. double yt = CInterpolate::GetValueLagrange3(m_nNumNodes, x0, dblStep, y, t[i]);
  363. CString s;
  364. s.Format("f(%f) = %frn", t[i], yt);
  365. m_strTargetYi += s;
  366. }
  367. AfxMessageBox(m_strTargetYi, MB_OK|MB_ICONINFORMATION);
  369. // 释放内存
  370. delete[] t;
  371. delete[] y;
  372. }
  374. /*
  375. 10
  377. -1.0,-0.8,-0.65,-0.4,-0.3,0.0,0.2,0.45,0.8,1.0
  379. 0.0384615,0.0588236,0.0864865,0.2,0.307692,1.0,0.5,0.164948,0.0588236,0.0384615
  381. 2
  383. 0.85,0.25
  384. */
  385. // 连分式不等距插值
  386. void CInterpolaterDlg::OnButton5() 
  387. {
  388. // 获取数据
  389. UpdateData();
  390. double *x = new double[m_nNumNodes];
  391. double *y = new double[m_nNumNodes];
  392. CInterpolate::GetNodesFromString(m_strXi, m_nNumNodes, x, ",");
  393. CInterpolate::GetNodesFromString(m_strYi, m_nNumNodes, y, ",");
  395. // 插值运算
  396. double *t = new double[m_nNumTargetNodes];
  397. CInterpolate::GetNodesFromString(m_strTargetXi, m_nNumTargetNodes, t, ",");
  399. // 显示结果
  400. m_strTargetYi = "";
  401. for (int i=0; i<m_nNumTargetNodes; ++i)
  402. {
  403. double yt = CInterpolate::GetValuePqs(m_nNumNodes, x, y, t[i]);
  404. CString s;
  405. s.Format("f(%f) = %frn", t[i], yt);
  406. m_strTargetYi += s;
  407. }
  408. AfxMessageBox(m_strTargetYi, MB_OK|MB_ICONINFORMATION);
  410. // 释放内存
  411. delete[] t;
  412. delete[] x;
  413. delete[] y;
  414. }
  416. /*
  417. 11
  419. 0.2
  421. -1.0
  423. 0.0384615,0.0588236,0.1,0.2,0.5,1.0,0.5,0.2,0.1,0.0588236,0.0384615
  425. 2
  427. -0.75, -0.05
  428. */
  429. // 连分式等距插值
  430. void CInterpolaterDlg::OnButton6() 
  431. {
  432. // 获取数据
  433. UpdateData();
  434. double x0;
  435. double *y = new double[m_nNumNodes];
  436. CInterpolate::GetNodesFromString(m_strXi, 1, &x0, ",");
  437. CInterpolate::GetNodesFromString(m_strYi, m_nNumNodes, y, ",");
  439. // 插值运算
  440. double *t = new double[m_nNumTargetNodes];
  441. CInterpolate::GetNodesFromString(m_strTargetXi, m_nNumTargetNodes, t, ",");
  443. // 显示结果
  444. double dblStep = atof(m_strStepOrBorderCondition);
  445. m_strTargetYi = "";
  446. for (int i=0; i<m_nNumTargetNodes; ++i)
  447. {
  448. double yt = CInterpolate::GetValuePqs(m_nNumNodes, x0, dblStep, y, t[i]);
  449. CString s;
  450. s.Format("f(%f) = %frn", t[i], yt);
  451. m_strTargetYi += s;
  452. }
  453. AfxMessageBox(m_strTargetYi, MB_OK|MB_ICONINFORMATION);
  455. // 释放内存
  456. delete[] t;
  457. delete[] y;
  458. }
  460. /*
  461. 10
  463. 0.1,0.15,0.3,0.45,0.55,0.6,0.7,0.85,0.9,1.0
  465. 0.904837,0.860708,0.740818,0.637628,0.576950,0.548812,0.496585,0.427415,0.406570,0.367879
  467. 1
  469. 0.356
  470. */
  471. // 埃尔米特不等距插值
  472. void CInterpolaterDlg::OnButton7() 
  473. {
  474. // 获取数据
  475. UpdateData();
  476. double *x = new double[m_nNumNodes];
  477. double *y = new double[m_nNumNodes];
  478. CInterpolate::GetNodesFromString(m_strXi, m_nNumNodes, x, ",");
  479. CInterpolate::GetNodesFromString(m_strYi, m_nNumNodes, y, ",");
  481. double *dy = new double[m_nNumNodes];
  482. for (int i=0; i<m_nNumNodes; ++i)
  483. dy[i]=-y[i];
  485. // 插值运算
  486. double *t = new double[m_nNumTargetNodes];
  487. CInterpolate::GetNodesFromString(m_strTargetXi, m_nNumTargetNodes, t, ",");
  489. // 显示结果
  490. m_strTargetYi = "";
  491. for (i=0; i<m_nNumTargetNodes; ++i)
  492. {
  493. double yt = CInterpolate::GetValueHermite(m_nNumNodes, x, y, dy, t[i]);
  494. CString s;
  495. s.Format("f(%f) = %frn", t[i], yt);
  496. m_strTargetYi += s;
  497. }
  498. AfxMessageBox(m_strTargetYi, MB_OK|MB_ICONINFORMATION);
  500. // 释放内存
  501. delete[] dy;
  502. delete[] t;
  503. delete[] x;
  504. delete[] y;
  505. }
  507. /*
  508. 10
  510. 0.1
  512. 0.1
  514. 0.904837,0.818731,0.740818,0.670320,0.606531,0.548812,0.496585,0.449329,0.406570,0.367879
  516. 1
  518. 0.752
  519. */
  520. // 埃尔米特等距插值
  521. void CInterpolaterDlg::OnButton8() 
  522. {
  523. // 获取数据
  524. UpdateData();
  525. double x0;
  526. double *y = new double[m_nNumNodes];
  527. CInterpolate::GetNodesFromString(m_strXi, 1, &x0, ",");
  528. CInterpolate::GetNodesFromString(m_strYi, m_nNumNodes, y, ",");
  530. double *dy = new double[m_nNumNodes];
  531. for (int i=0; i<m_nNumNodes; ++i)
  532. dy[i]=-y[i];
  534. // 插值运算
  535. double *t = new double[m_nNumTargetNodes];
  536. CInterpolate::GetNodesFromString(m_strTargetXi, m_nNumTargetNodes, t, ",");
  538. // 显示结果
  539. double dblStep = atof(m_strStepOrBorderCondition);
  540. m_strTargetYi = "";
  541. for (i=0; i<m_nNumTargetNodes; ++i)
  542. {
  543. double yt = CInterpolate::GetValueHermite(m_nNumNodes, x0, dblStep, y, dy, t[i]);
  544. CString s;
  545. s.Format("f(%f) = %frn", t[i], yt);
  546. m_strTargetYi += s;
  547. }
  548. AfxMessageBox(m_strTargetYi, MB_OK|MB_ICONINFORMATION);
  550. // 释放内存
  551. delete[] dy;
  552. delete[] t;
  553. delete[] y;
  554. }
  556. /*
  557. 10
  559. -1.0,-0.8,-0.65,-0.4,-0.3,0.0,0.2,0.4,0.6,0.8,1.0
  561. 0.0384615,0.0588236,0.0864865,0.2,0.307692,1.0,0.5,0.2,0.1,0.0588236,0.0384615
  563. 2
  565. -0.75,0.05
  566. */
  567. // 埃特金不等距逐步插值
  568. void CInterpolaterDlg::OnButton9() 
  569. {
  570. // 获取数据
  571. UpdateData();
  572. double *x = new double[m_nNumNodes];
  573. double *y = new double[m_nNumNodes];
  574. CInterpolate::GetNodesFromString(m_strXi, m_nNumNodes, x, ",");
  575. CInterpolate::GetNodesFromString(m_strYi, m_nNumNodes, y, ",");
  577. // 插值运算
  578. double *t = new double[m_nNumTargetNodes];
  579. CInterpolate::GetNodesFromString(m_strTargetXi, m_nNumTargetNodes, t, ",");
  581. // 显示结果
  582. m_strTargetYi = "";
  583. for (int i=0; i<m_nNumTargetNodes; ++i)
  584. {
  585. double yt = CInterpolate::GetValueAitken(m_nNumNodes, x, y, t[i]);
  586. CString s;
  587. s.Format("f(%f) = %frn", t[i], yt);
  588. m_strTargetYi += s;
  589. }
  590. AfxMessageBox(m_strTargetYi, MB_OK|MB_ICONINFORMATION);
  592. // 释放内存
  593. delete[] t;
  594. delete[] x;
  595. delete[] y;
  596. }
  598. /*
  599. 10
  601. 0.1
  603. 0.1
  605. 0.904837,0.818731,0.740818,0.670320,0.606531,0.548812,0.496585,0.449329,0.406570,0.367879
  607. 2
  609. 0.15,0.55
  610. */
  611. // 埃特金等距逐步插值
  612. void CInterpolaterDlg::OnButton10() 
  613. {
  614. // 获取数据
  615. UpdateData();
  616. double x0;
  617. double *y = new double[m_nNumNodes];
  618. CInterpolate::GetNodesFromString(m_strXi, 1, &x0, ",");
  619. CInterpolate::GetNodesFromString(m_strYi, m_nNumNodes, y, ",");
  621. // 插值运算
  622. double *t = new double[m_nNumTargetNodes];
  623. CInterpolate::GetNodesFromString(m_strTargetXi, m_nNumTargetNodes, t, ",");
  625. // 显示结果
  626. double dblStep = atof(m_strStepOrBorderCondition);
  627. m_strTargetYi = "";
  628. for (int i=0; i<m_nNumTargetNodes; ++i)
  629. {
  630. double yt = CInterpolate::GetValueAitken(m_nNumNodes, x0, dblStep, y, t[i]);
  631. CString s;
  632. s.Format("f(%f) = %frn", t[i], yt);
  633. m_strTargetYi += s;
  634. }
  635. AfxMessageBox(m_strTargetYi, MB_OK|MB_ICONINFORMATION);
  637. // 释放内存
  638. delete[] t;
  639. delete[] y;
  640. }
  642. /*
  643. 11
  645. -1.0,-0.95,-0.75,-0.55,-0.3,0.0,0.2,0.45,0.6,0.8,1.0
  647. 0.0384615,0.0424403,0.06639,0.116788,0.307692,1.0,0.5,0.164948,0.1,0.0588236,0.0384615
  649. 2
  651. -0.85,0.15
  652. */
  653. // 光滑不等距插值
  654. void CInterpolaterDlg::OnButton11() 
  655. {
  656. // 获取数据
  657. UpdateData();
  658. double *x = new double[m_nNumNodes];
  659. double *y = new double[m_nNumNodes];
  660. CInterpolate::GetNodesFromString(m_strXi, m_nNumNodes, x, ",");
  661. CInterpolate::GetNodesFromString(m_strYi, m_nNumNodes, y, ",");
  663. // 插值运算
  664. double *t = new double[m_nNumTargetNodes];
  665. CInterpolate::GetNodesFromString(m_strTargetXi, m_nNumTargetNodes, t, ",");
  667. // 显示结果
  668. m_strTargetYi = "";
  669. for (int i=0; i<m_nNumTargetNodes; ++i)
  670. {
  671. double dblCoef[5];
  672. double yt = CInterpolate::GetValueAkima(m_nNumNodes, x, y, t[i], dblCoef);
  673. CString s;
  674. s.Format("f(%f) = %frns0 = %f  s1 = %f  s2 = %f  s3 = %frnrn", 
  675. t[i], yt, dblCoef[0], dblCoef[1], dblCoef[2], dblCoef[3]);
  676. m_strTargetYi += s;
  677. }
  678. AfxMessageBox(m_strTargetYi, MB_OK|MB_ICONINFORMATION);
  680. // 释放内存
  681. delete[] t;
  682. delete[] x;
  683. delete[] y;
  684. }
  686. /*
  687. 11
  689. 0.2
  691. -1.0
  693. 0.0384615,0.0588236,0.1,0.2,0.5,1.0,0.5,0.2,0.1,0.0588236,0.0384615
  695. 2
  697. -0.65,0.25
  698. */
  699. // 光滑等距插值
  700. void CInterpolaterDlg::OnButton12() 
  701. {
  702. // 获取数据
  703. UpdateData();
  704. double x0;
  705. double *y = new double[m_nNumNodes];
  706. CInterpolate::GetNodesFromString(m_strXi, 1, &x0, ",");
  707. CInterpolate::GetNodesFromString(m_strYi, m_nNumNodes, y, ",");
  709. // 插值运算
  710. double *t = new double[m_nNumTargetNodes];
  711. CInterpolate::GetNodesFromString(m_strTargetXi, m_nNumTargetNodes, t, ",");
  713. // 显示结果
  714. double dblStep = atof(m_strStepOrBorderCondition);
  715. m_strTargetYi = "";
  716. for (int i=0; i<m_nNumTargetNodes; ++i)
  717. {
  718. double dblCoef[5];
  719. double yt = CInterpolate::GetValueAkima(m_nNumNodes, x0, dblStep, y, t[i], dblCoef);
  720. CString s;
  721. s.Format("f(%f) = %frns0 = %f  s1 = %f  s2 = %f  s3 = %frnrn", 
  722. t[i], yt, dblCoef[0], dblCoef[1], dblCoef[2], dblCoef[3]);
  723. m_strTargetYi += s;
  724. }
  725. AfxMessageBox(m_strTargetYi, MB_OK|MB_ICONINFORMATION);
  727. // 释放内存
  728. delete[] t;
  729. delete[] y;
  730. }
  732. /*
  733. 12
  735. 1.86548,-0.046115
  737. 0.52,8.0,17.95,28.65,50.65,104.6,156.6,260.7,364.4,468.0,507.0,520.0
  739. 5.28794,13.84,20.2,24.9,31.1,36.5,36.6,31.0,20.9,7.8,1.5,0.2
  741. 8
  743. 4.0,14.0,30.0,60.0,130.0,230.0,450.0,515.0
  744. */
  745. // 第一种边界条件的三次样条函数插值、微商与积分
  746. void CInterpolaterDlg::OnButton13() 
  747. {
  748. // 获取数据
  749. UpdateData();
  750. double *x = new double[m_nNumNodes];
  751. double *y = new double[m_nNumNodes];
  752. CInterpolate::GetNodesFromString(m_strXi, m_nNumNodes, x, ",");
  753. CInterpolate::GetNodesFromString(m_strYi, m_nNumNodes, y, ",");
  755. double *dy = new double[m_nNumNodes];
  756. double *ddy = new double[m_nNumNodes];
  758. double dblTmp[2];
  759. CInterpolate::GetNodesFromString(m_strStepOrBorderCondition, 2, dblTmp, ",");
  760. dy[0] = dblTmp[0];
  761. dy[m_nNumNodes-1] = dblTmp[1];
  763. // 插值运算
  764. double *t = new double[m_nNumTargetNodes];
  765. double *z = new double[m_nNumTargetNodes];
  766. double *dz = new double[m_nNumTargetNodes];
  767. double *ddz = new double[m_nNumTargetNodes];
  769. CInterpolate::GetNodesFromString(m_strTargetXi, m_nNumTargetNodes, t, ",");
  771. double yt = CInterpolate::GetValueSpline1(m_nNumNodes, x, y, dy, ddy, m_nNumTargetNodes, t, z, dz, ddz);
  773. // 显示结果
  774. m_strTargetYi.Format("积分值 ＝ %frnrn", yt);
  775. m_strTargetYi += "    t[i]t    z[i]t    dz[i]t    ddz[i]rn";
  776. for (int i=0; i<m_nNumTargetNodes; ++i)
  777. {
  778. CString s;
  779. s.Format("%f  %f  %f  %frn", t[i], z[i], dz[i], ddz[i]);
  780. m_strTargetYi += s;
  781. }
  782. AfxMessageBox(m_strTargetYi, MB_OK|MB_ICONINFORMATION);
  784. // 释放内存
  785. delete[] t;
  786. delete[] x;
  787. delete[] y;
  788. delete[] dy;
  789. delete[] ddy;
  790. delete[] z;
  791. delete[] dz;
  792. delete[] ddz;
  793. }
  795. /*
  796. 12
  798. -0.279319,0.011156
  800. 0.52,8.0,17.95,28.65,50.65,104.6,156.6,260.7,364.4,468.0,507.0,520.0
  802. 5.28794,13.84,20.2,24.9,31.1,36.5,36.6,31.0,20.9,7.8,1.5,0.2
  804. 8
  806. 4.0,14.0,30.0,60.0,130.0,230.0,450.0,515.0
  807. */
  808. // 第二种边界条件的三次样条函数插值、微商与积分
  809. void CInterpolaterDlg::OnButton14() 
  810. {
  811. // 获取数据
  812. UpdateData();
  813. double *x = new double[m_nNumNodes];
  814. double *y = new double[m_nNumNodes];
  815. CInterpolate::GetNodesFromString(m_strXi, m_nNumNodes, x, ",");
  816. CInterpolate::GetNodesFromString(m_strYi, m_nNumNodes, y, ",");
  818. double *dy = new double[m_nNumNodes];
  819. double *ddy = new double[m_nNumNodes];
  821. double dblTmp[2];
  822. CInterpolate::GetNodesFromString(m_strStepOrBorderCondition, 2, dblTmp, ",");
  823. ddy[0] = dblTmp[0];
  824. ddy[m_nNumNodes-1] = dblTmp[1];
  826. // 插值运算
  827. double *t = new double[m_nNumTargetNodes];
  828. double *z = new double[m_nNumTargetNodes];
  829. double *dz = new double[m_nNumTargetNodes];
  830. double *ddz = new double[m_nNumTargetNodes];
  832. CInterpolate::GetNodesFromString(m_strTargetXi, m_nNumTargetNodes, t, ",");
  834. double yt = CInterpolate::GetValueSpline2(m_nNumNodes, x, y, dy, ddy, m_nNumTargetNodes, t, z, dz, ddz);
  836. // 显示结果
  837. m_strTargetYi.Format("积分值 ＝ %frnrn", yt);
  838. m_strTargetYi += "    t[i]t    z[i]t    dz[i]t    ddz[i]rn";
  839. for (int i=0; i<m_nNumTargetNodes; ++i)
  840. {
  841. CString s;
  842. s.Format("%f  %f  %f  %frn", t[i], z[i], dz[i], ddz[i]);
  843. m_strTargetYi += s;
  844. }
  845. AfxMessageBox(m_strTargetYi, MB_OK|MB_ICONINFORMATION);
  847. // 释放内存
  848. delete[] t;
  849. delete[] x;
  850. delete[] y;
  851. delete[] dy;
  852. delete[] ddy;
  853. delete[] z;
  854. delete[] dz;
  855. delete[] ddz;
  856. }
  858. /*
  859. 37
  861. 36
  862. */
  863. // 第三种边界条件的三次样条函数插值、微商与积分
  864. void CInterpolaterDlg::OnButton15() 
  865. {
  866. // 获取数据
  867. UpdateData();
  869. double x[37], y[37], dy[37], ddy[37];
  870. double t[36], z[36], dz[36], ddz[36];
  872. for (int i=0; i<37; ++i)
  873. {
  874. x[i]=i*6.2831852/36.0;
  875.         y[i]=sin(x[i]);
  876. }
  878. // 插值运算
  879. for (i=0; i<36; ++i)
  880. {
  881.         t[i]=(0.5+i)*6.2831852/36.0;
  882. }
  884. double yt = CInterpolate::GetValueSpline3(37, x, y, dy, ddy, 36, t, z, dz, ddz);
  886. // 显示结果
  887. CString s, ss;
  888. s.Format("积分值 ＝ %frnrn", yt);
  889.     ss.Format(" x(i) y(i)=sin(x) dy(i)=cos(x) ddy(i)=-sin(x)rnrn");
  890. s += ss;
  891.     ss.Format("%6.1f%10.6f%11.6f%12.6frn",x[0],y[0],dy[0],ddy[0]);
  892. s += ss;
  893.     for (i=0; i<36; i++)
  894.     { 
  895. double u=t[i]*36.0/0.62831852;
  896.         ss.Format("%6.1f%10.6f%11.6f%12.6frn",u,z[i],dz[i],ddz[i]);
  897. s += ss;
  899.         u=x[i+1]*36.0/0.62831852;
  900.         int j=i+1;
  901.         ss.Format("%6.1f%10.6f%11.6f%12.6frn",u,y[j],dy[j],ddy[j]);
  902. s += ss;
  903.     }
  905. CSpline3Dlg dlg;
  907. dlg.m_strResult = s;
  908. dlg.DoModal();
  910. }
  912. // 二元三点插值
  913. void CInterpolaterDlg::OnButton16() 
  914. {
  915. int i,j;
  916.     double u,v,w;
  917.     double x[6],y[5],z[30];
  919.     for (i=0;i<=5;i++) 
  920. x[i]=0.2*i;
  921.     for (j=0;j<=4;j++) 
  922. y[j]=0.25*j;
  923.     for (i=0;i<=5;i++)
  924.     {
  925. for (j=0;j<=4;j++)
  926. z[i*5+j]=exp(-(x[i]-y[j]));
  927. }
  929. // 插值运算
  930. CString s, ss;
  932.     u = 0.9; 
  933. v = 0.8;
  934.     w = CInterpolate::GetValueTqip(6, x, 5, y, z, u, v);
  935.     s.Format("z(%f, %f) = %fn",u,v,w);
  937.     u = 0.3; 
  938. v = 0.9;
  939.     w = CInterpolate::GetValueTqip(6, x, 5, y, z, u, v);
  940.     ss.Format("z(%f, %f) = %fn",u,v,w);
  942. s += ss;
  944. // 显示结果
  945. AfxMessageBox(s, MB_OK|MB_ICONINFORMATION);
  946. }
  948. // 二元全区间插值
  949. void CInterpolaterDlg::OnButton17() 
  950. {
  951. int i,j;
  952.     double u,v,w;
  953.     double x[11],y[11],z[121];
  955.     for (i=0;i<=10;i++)
  956.     { 
  957. x[i]=0.1*i; 
  958. y[i]=x[i];
  959. }
  960.     
  961. for (i=0;i<=10;i++)
  962.     {
  963. for (j=0;j<=10;j++)
  964. z[i*11+j]=exp(-(x[i]-y[j]));
  965. }
  967. // 插值运算
  968. CString s, ss;
  970.     u=0.35; 
  971. v=0.65;
  972.     w = CInterpolate::GetValueLagrange2(11, x, 11, y, z, u, v);
  973.     s.Format("z(%f, %f) = %fn",u,v,w);
  975.     u=0.45; 
  976. v=0.55;
  977.     w = CInterpolate::GetValueLagrange2(11, x, 11, y, z, u, v);
  978.     ss.Format("z(%f, %f) = %fn",u,v,w);
  980. s += ss;
  982. // 显示结果
  983. AfxMessageBox(s, MB_OK|MB_ICONINFORMATION);
  984. }