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

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

数学计算

开发平台：

Visual C++

LECalculatorDlg.cpp：源码内容
  1. // LECalculatorDlg.cpp : implementation file
  2. //
  4. #include "stdafx.h"
  5. #include "LECalculator.h"
  6. #include "LECalculatorDlg.h"
  7. #include "LEquations.h"
  8. #include "ImagDataDlg.h"
  9. #include "Complex.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. // CLECalculatorDlg dialog
  65. CLECalculatorDlg::CLECalculatorDlg(CWnd* pParent /*=NULL*/)
  66. : CDialog(CLECalculatorDlg::IDD, pParent)
  67. {
  68. //{{AFX_DATA_INIT(CLECalculatorDlg)
  69. m_strCoefficient = _T("");
  70. m_strConstant = _T("");
  71. //}}AFX_DATA_INIT
  72. // Note that LoadIcon does not require a subsequent DestroyIcon in Win32
  73. m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
  74. }
  76. void CLECalculatorDlg::DoDataExchange(CDataExchange* pDX)
  77. {
  78. CDialog::DoDataExchange(pDX);
  79. //{{AFX_DATA_MAP(CLECalculatorDlg)
  80. DDX_Text(pDX, IDC_EDIT1, m_strCoefficient);
  81. DDX_Text(pDX, IDC_EDIT2, m_strConstant);
  82. //}}AFX_DATA_MAP
  83. }
  85. BEGIN_MESSAGE_MAP(CLECalculatorDlg, CDialog)
  86. //{{AFX_MSG_MAP(CLECalculatorDlg)
  87. ON_WM_SYSCOMMAND()
  88. ON_WM_PAINT()
  89. ON_WM_QUERYDRAGICON()
  90. ON_BN_CLICKED(IDC_BUTTON1, OnButton1)
  91. ON_BN_CLICKED(IDC_BUTTON2, OnButton2)
  92. ON_BN_CLICKED(IDC_BUTTON3, OnButton3)
  93. ON_BN_CLICKED(IDC_BUTTON4, OnButton4)
  94. ON_BN_CLICKED(IDC_BUTTON5, OnButton5)
  95. ON_BN_CLICKED(IDC_BUTTON6, OnButton6)
  96. ON_BN_CLICKED(IDC_BUTTON7, OnButton7)
  97. ON_BN_CLICKED(IDC_BUTTON8, OnButton8)
  98. ON_BN_CLICKED(IDC_BUTTON9, OnButton9)
  99. ON_BN_CLICKED(IDC_BUTTON10, OnButton10)
  100. ON_BN_CLICKED(IDC_BUTTON11, OnButton11)
  101. ON_BN_CLICKED(IDC_BUTTON12, OnButton12)
  102. ON_BN_CLICKED(IDC_BUTTON13, OnButton13)
  103. ON_BN_CLICKED(IDC_BUTTON14, OnButton14)
  104. ON_BN_CLICKED(IDC_BUTTON15, OnButton15)
  105. //}}AFX_MSG_MAP
  106. END_MESSAGE_MAP()
  108. /////////////////////////////////////////////////////////////////////////////
  109. // CLECalculatorDlg message handlers
  111. BOOL CLECalculatorDlg::OnInitDialog()
  112. {
  113. CDialog::OnInitDialog();
  115. // Add "About..." menu item to system menu.
  117. // IDM_ABOUTBOX must be in the system command range.
  118. ASSERT((IDM_ABOUTBOX & 0xFFF0) == IDM_ABOUTBOX);
  119. ASSERT(IDM_ABOUTBOX < 0xF000);
  121. CMenu* pSysMenu = GetSystemMenu(FALSE);
  122. if (pSysMenu != NULL)
  123. {
  124. CString strAboutMenu;
  125. strAboutMenu.LoadString(IDS_ABOUTBOX);
  126. if (!strAboutMenu.IsEmpty())
  127. {
  128. pSysMenu->AppendMenu(MF_SEPARATOR);
  129. pSysMenu->AppendMenu(MF_STRING, IDM_ABOUTBOX, strAboutMenu);
  130. }
  131. }
  133. // Set the icon for this dialog.  The framework does this automatically
  134. //  when the application's main window is not a dialog
  135. SetIcon(m_hIcon, TRUE); // Set big icon
  136. SetIcon(m_hIcon, FALSE); // Set small icon
  138. // TODO: Add extra initialization here
  140. return TRUE;  // return TRUE  unless you set the focus to a control
  141. }
  143. void CLECalculatorDlg::OnSysCommand(UINT nID, LPARAM lParam)
  144. {
  145. if ((nID & 0xFFF0) == IDM_ABOUTBOX)
  146. {
  147. CAboutDlg dlgAbout;
  148. dlgAbout.DoModal();
  149. }
  150. else
  151. {
  152. CDialog::OnSysCommand(nID, lParam);
  153. }
  154. }
  156. // If you add a minimize button to your dialog, you will need the code below
  157. //  to draw the icon.  For MFC applications using the document/view model,
  158. //  this is automatically done for you by the framework.
  160. void CLECalculatorDlg::OnPaint() 
  161. {
  162. if (IsIconic())
  163. {
  164. CPaintDC dc(this); // device context for painting
  166. SendMessage(WM_ICONERASEBKGND, (WPARAM) dc.GetSafeHdc(), 0);
  168. // Center icon in client rectangle
  169. int cxIcon = GetSystemMetrics(SM_CXICON);
  170. int cyIcon = GetSystemMetrics(SM_CYICON);
  171. CRect rect;
  172. GetClientRect(&rect);
  173. int x = (rect.Width() - cxIcon + 1) / 2;
  174. int y = (rect.Height() - cyIcon + 1) / 2;
  176. // Draw the icon
  177. dc.DrawIcon(x, y, m_hIcon);
  178. }
  179. else
  180. {
  181. CDialog::OnPaint();
  182. }
  183. }
  185. // The system calls this to obtain the cursor to display while the user drags
  186. //  the minimized window.
  187. HCURSOR CLECalculatorDlg::OnQueryDragIcon()
  188. {
  189. return (HCURSOR) m_hIcon;
  190. }
  192. /*
  193. 0.2368,0.2471,0.2568,1.2671
  194. 0.1968,0.2071,1.2168,0.2271
  195. 0.1581,1.1675,0.1768,0.1871
  196. 1.1161,0.1254,0.1397,0.1490
  198. 1.8471
  199. 1.7471
  200. 1.6471
  201. 1.5471
  202. */
  203. // 全选主元高斯消去法
  204. void CLECalculatorDlg::OnButton1() 
  205. {
  206. // 获取数据
  207. UpdateData();
  209. // 构造矩阵
  210. CMatrix mtxA, mtxB, mtxC;
  212. // 初始化系数矩阵
  213. mtxA.FromString(m_strCoefficient, ",");
  214. // 初始化常数数矩阵
  215. mtxB.FromString(m_strConstant, ",");
  217. // 构造线性方程组
  218. CLEquations leqs(mtxA, mtxB);
  220. // 全选主元高斯消去法
  221. if (leqs.GetRootsetGauss(mtxC))
  222. m_strResult = mtxC.ToString(" ");
  223. else
  224. m_strResult = "求解失败";
  226. // 显示数据
  227. AfxMessageBox(m_strResult, MB_OK|MB_ICONINFORMATION);
  228. }
  230. /*
  231. 1.0,3.0,2.0,13.0
  232. 7.0,2.0,1.0,-2.0
  233. 9.0,15.0,3.0,-2.0
  234. -2.0,-2.0,11.0,5.0
  236. 9.0,0.0
  237. 6.0,4.0
  238. 11.0,7.0
  239. -2.0,-1.0
  240. */
  241. // 全选主元高斯－约当消去法
  242. void CLECalculatorDlg::OnButton2() 
  243. {
  244. // 获取数据
  245. UpdateData();
  247. // 构造矩阵
  248. CMatrix mtxA, mtxB, mtxC;
  250. // 初始化系数矩阵
  251. mtxA.FromString(m_strCoefficient, ",");
  252. // 初始化常数数矩阵
  253. mtxB.FromString(m_strConstant, ",");
  255. // 构造线性方程组
  256. CLEquations leqs(mtxA, mtxB);
  258. // 全选主元高斯－约当消去法
  259. if (leqs.GetRootsetGaussJordan(mtxC))
  260. m_strResult = mtxC.ToString(" ");
  261. else
  262. m_strResult = "求解失败";
  264. // 显示数据
  265. AfxMessageBox(m_strResult, MB_OK|MB_ICONINFORMATION);
  266. }
  268. /*
  269. 1.0,3.0,2.0,13.0
  270. 7.0,2.0,1.0,-2.0
  271. 9.0,15.0,3.0,-2.0
  272. -2.0,-2.0,11.0,5.0
  274. 3.0,-2.0,1.0,6.0
  275. -2.0,7.0,5.0,8.0
  276. 9.0,-3.0,15.0,1.0
  277. -2.0,-2.0,7.0,6.0
  279. 2.0
  280. 7.0
  281. 3.0
  282. 9.0
  284. 1.0
  285. 2.0
  286. -2.0
  287. 3.0
  288. */
  289. // 复系数方程组的全选主元高斯消去法
  290. void CLECalculatorDlg::OnButton3() 
  291. {
  292. // 获取系数和常数矩阵的虚部矩阵数据
  293. CImagDataDlg dlg;
  294. if (dlg.DoModal() != IDOK)
  295. return;
  297. // 获取数据
  298. UpdateData();
  300. // 构造矩阵
  301. CMatrix mtxA, mtxB, mtxAI, mtxBI, mtxC, mtxCI;
  303. // 初始化系数矩阵实部矩阵
  304. mtxA.FromString(m_strCoefficient, ",");
  305. // 初始化系数矩阵实部矩阵
  306. mtxB.FromString(m_strConstant, ",");
  308. // 初始化系数矩阵虚部矩阵
  309. mtxAI.FromString(dlg.m_strCoefficient, ",");
  310. // 初始化系数矩阵虚部矩阵
  311. mtxBI.FromString(dlg.m_strConstant, ",");
  313. // 构造线性方程组
  314. CLEquations leqs(mtxA, mtxB);
  316. // 复系数方程组的全选主元高斯消去法
  317. if (leqs.GetRootsetGauss(mtxAI, mtxBI, mtxC, mtxCI))
  318. {
  319. for (int i=0; i<mtxA.GetNumColumns(); ++i)
  320. {
  321. CComplex cpxS(mtxC.GetElement(i, 0), mtxCI.GetElement(i, 0));
  322. m_strResult += cpxS.ToString() + "rn";
  323. }
  324. }
  325. else
  326. m_strResult = "求解失败";
  328. // 显示数据
  329. AfxMessageBox(m_strResult, MB_OK|MB_ICONINFORMATION);
  330. }
  332. /*
  333. 1.0,3.0,2.0,13.0
  334. 7.0,2.0,1.0,-2.0
  335. 9.0,15.0,3.0,-2.0
  336. -2.0,-2.0,11.0,5.0
  338. 3.0,-2.0,1.0,6.0
  339. -2.0,7.0,5.0,8.0
  340. 9.0,-3.0,15.0,1.0
  341. -2.0,-2.0,7.0,6.0
  343. 2.0,-2.0
  344. 7.0,3.0
  345. 3.0,2.0
  346. 9.0,1.0
  348. 1.0,3.0
  349. 2.0,7.0
  350. -2.0,9.0
  351. 3.0,2.0
  352. */
  353. // 复系数方程组的全选主元高斯－约当消去法
  354. void CLECalculatorDlg::OnButton4() 
  355. {
  356. // 获取系数和常数矩阵的虚部矩阵数据
  357. CImagDataDlg dlg;
  358. if (dlg.DoModal() != IDOK)
  359. return;
  361. // 获取数据
  362. UpdateData();
  364. // 构造矩阵
  365. CMatrix mtxA, mtxB, mtxAI, mtxBI, mtxC, mtxCI;
  367. // 初始化系数矩阵实部矩阵
  368. mtxA.FromString(m_strCoefficient, ",");
  369. // 初始化系数矩阵实部矩阵
  370. mtxB.FromString(m_strConstant, ",");
  372. // 初始化系数矩阵虚部矩阵
  373. mtxAI.FromString(dlg.m_strCoefficient, ",");
  374. // 初始化系数矩阵虚部矩阵
  375. mtxBI.FromString(dlg.m_strConstant, ",");
  377. // 构造线性方程组
  378. CLEquations leqs(mtxA, mtxB);
  380. // 复系数方程组的全选主元高斯－约当消去法
  381. if (leqs.GetRootsetGaussJordan(mtxAI, mtxBI, mtxC, mtxCI))
  382. {
  383. for (int i=0; i<mtxA.GetNumColumns(); ++i)
  384. {
  385. CComplex cpxS1(mtxC.GetElement(i, 0), mtxCI.GetElement(i, 0));
  386. CComplex cpxS2(mtxC.GetElement(i, 1), mtxCI.GetElement(i, 1));
  387. m_strResult += cpxS1.ToString() + " , " + cpxS2.ToString() + "rn";
  388. }
  389. }
  390. else
  391. m_strResult = "求解失败";
  393. // 显示数据
  394. AfxMessageBox(m_strResult, MB_OK|MB_ICONINFORMATION);
  395. }
  397. /*
  398. 13,12,0,0,0
  399. 11,10,9,0,0
  400.  0, 8,7,6,0
  401.  0, 0,5,4,3
  402.  0, 0,0,2,1
  404. 3.0
  405. 0.0
  406. -2.0
  407. 6.0
  408. 8.0
  409. */
  410. // 求解三对角线方程组的追赶法
  411. void CLECalculatorDlg::OnButton5() 
  412. {
  413. // 获取数据
  414. UpdateData();
  416. // 构造矩阵
  417. CMatrix mtxA, mtxB, mtxC;
  419. // 初始化系数矩阵
  420. mtxA.FromString(m_strCoefficient, ",");
  421. // 初始化常数数矩阵
  422. mtxB.FromString(m_strConstant, ",");
  424. // 构造线性方程组
  425. CLEquations leqs(mtxA, mtxB);
  427. // 求解三对角线方程组的追赶法
  428. if (leqs.GetRootsetTriDiagonal(mtxC))
  429. m_strResult = mtxC.ToString(" ");
  430. else
  431. m_strResult = "求解失败";
  433. // 显示数据
  434. AfxMessageBox(m_strResult, MB_OK|MB_ICONINFORMATION);
  435. }
  437. /*
  438. 3.0,-4.0,1.0,0.0,0.0,0,0,0
  439. -2.0,-5.0,6.0,1.0,0.0,0,0,0
  440. 1.0,3.0,-1.0,2.0,-3.0,0,0,0
  441. 0,2.0,5.0,-5.0,6.0,-1.0,0,0
  442. 0,0,-3.0,1.0,-1.0,2.0,-5.0,0
  443. 0,0,0,6.0,1.0,-3.0,2.0,-9.0
  444. 0,0,0,0,-4.0,1.0,-1.0,2.0
  445. 0,0,0,0,0,5.0,1.0,-7.0
  447. 13.0,29.0,-13.0
  448. -6.0,17.0,-21.0
  449. -31.0,-6.0,4.0
  450. 64.0,3.0,16.0
  451. -20.0,1.0,-5.0
  452. -22.0,-41.0,56.0
  453. -29.0,10.0,-21.0
  454. 7.0,-24.0,20.0
  455. */
  456. // 一般带型方程组的求解
  457. void CLECalculatorDlg::OnButton6() 
  458. {
  459. // 获取数据
  460. UpdateData();
  462. // 构造矩阵
  463. CMatrix mtxA, mtxB, mtxC;
  465. // 初始化系数矩阵
  466. mtxA.FromString(m_strCoefficient, ",");
  467. // 初始化常数数矩阵
  468. mtxB.FromString(m_strConstant, ",");
  470. // 构造线性方程组
  471. CLEquations leqs(mtxA, mtxB);
  473. // 一般带型方程组的求解
  474. if (leqs.GetRootsetBand(5, mtxC))
  475. m_strResult = mtxC.ToString(" ");
  476. else
  477. m_strResult = "求解失败";
  479. // 显示数据
  480. AfxMessageBox(m_strResult, MB_OK|MB_ICONINFORMATION);
  481. }
  483. /*
  484. 5.0,7.0,6.0,5.0,1.0
  485. 7.0,10.0,8.0,7.0,2.0
  486. 6.0,8.0,10.0,9.0,3.0
  487. 5.0,7.0,9.0,10.0,4.0
  488. 1.0,2.0,3.0,4.0,5.0
  490. 24.0,96.0
  491. 34.0,136.0
  492. 36.0,144.0
  493. 35.0,140.0
  494. 15.0,60.0
  495. */
  496. // 求解对称方程组的分解法
  497. void CLECalculatorDlg::OnButton7() 
  498. {
  499. // 获取数据
  500. UpdateData();
  502. // 构造矩阵
  503. CMatrix mtxA, mtxB, mtxC;
  505. // 初始化系数矩阵
  506. mtxA.FromString(m_strCoefficient, ",");
  507. // 初始化常数数矩阵
  508. mtxB.FromString(m_strConstant, ",");
  510. // 构造线性方程组
  511. CLEquations leqs(mtxA, mtxB);
  513. // 求解对称方程组的分解法
  514. if (leqs.GetRootsetDjn(mtxC))
  515. m_strResult = mtxC.ToString(" ");
  516. else
  517. m_strResult = "求解失败";
  519. // 显示数据
  520. AfxMessageBox(m_strResult, MB_OK|MB_ICONINFORMATION);
  521. }
  523. /*
  524. 5.0,7.0,6.0,5.0
  525. 7.0,10.0,8.0,7.0
  526. 6.0,8.0,10.0,9.0
  527. 5.0,7.0,9.0,10.0
  529. 23.0,92.0
  530. 32.0,128.0
  531. 33.0,132.0
  532. 31.0,124.0
  533. */
  534. // 求解对称正定方程组的平方根法
  535. void CLECalculatorDlg::OnButton8() 
  536. {
  537. // 获取数据
  538. UpdateData();
  540. // 构造矩阵
  541. CMatrix mtxA, mtxB, mtxC;
  543. // 初始化系数矩阵
  544. mtxA.FromString(m_strCoefficient, ",");
  545. // 初始化常数数矩阵
  546. mtxB.FromString(m_strConstant, ",");
  548. // 构造线性方程组
  549. CLEquations leqs(mtxA, mtxB);
  551. // 求解对称正定方程组的平方根法
  552. if (leqs.GetRootsetCholesky(mtxC))
  553. m_strResult = mtxC.ToString(" ");
  554. else
  555. m_strResult = "求解失败";
  557. // 显示数据
  558. AfxMessageBox(m_strResult, MB_OK|MB_ICONINFORMATION);
  559. }
  561. /*
  562. 0.0,0.0,-1.0,0.0,0.0,0.0,2.0,0.0
  563. 0.0,6.0,0.0,0.0,0.0,-6.0,0.0,0.0
  564. 0.0,0.0,0.0,2.0,0.0,0.0,0.0,-4.0
  565. 3.0,0.0,0.0,0.0,-2.0,0.0,1.0,0.0
  566. 0.0,0.0,6.0,0.0,0.0,0.0,5.0,0.0
  567. 1.0,0.0,0.0,0.0,-3.0,0.0,0.0,2.0
  568. 0.0,4.0,0.0,-1.0,0.0,0.0,0.0,0.0
  569. 0.0,0.0,1.0,0.0,-1.0,0.0,0.0,-2.0
  571. 4.0
  572. 6.0
  573. -8.0
  574. -2.0
  575. 27.0
  576. -9.0
  577. 2.0
  578. -4.0
  579. */
  580. // 求解大型稀疏方程组的全选主元高斯－约去消去法
  581. void CLECalculatorDlg::OnButton9() 
  582. {
  583. // 获取数据
  584. UpdateData();
  586. // 构造矩阵
  587. CMatrix mtxA, mtxB, mtxC;
  589. // 初始化系数矩阵
  590. mtxA.FromString(m_strCoefficient, ",");
  591. // 初始化常数数矩阵
  592. mtxB.FromString(m_strConstant, ",");
  594. // 构造线性方程组
  595. CLEquations leqs(mtxA, mtxB);
  597. // 求解大型稀疏方程组的全选主元高斯－约去消去法
  598. if (leqs.GetRootsetGgje(mtxC))
  599. m_strResult = mtxC.ToString(" ");
  600. else
  601. m_strResult = "求解失败";
  603. // 显示数据
  604. AfxMessageBox(m_strResult, MB_OK|MB_ICONINFORMATION);
  605. }
  607. /*
  608. 6.0,5.0,4.0,3.0,2.0,1.0
  609. 5.0,6.0,5.0,4.0,3.0,2.0
  610. 4.0,5.0,6.0,5.0,4.0,3.0
  611. 3.0,4.0,5.0,6.0,5.0,4.0
  612. 2.0,3.0,4.0,5.0,6.0,5.0
  613. 1.0,2.0,3.0,4.0,5.0,6.0
  615. 11.0
  616. 9.0
  617. 9.0
  618. 9.0
  619. 13.0
  620. 17.0
  621. */
  622. // 求解托伯利兹方程组的列文逊方法
  623. void CLECalculatorDlg::OnButton10() 
  624. {
  625. // 获取数据
  626. UpdateData();
  628. // 构造矩阵
  629. CMatrix mtxA, mtxB, mtxC;
  631. // 初始化系数矩阵
  632. mtxA.FromString(m_strCoefficient, ",");
  633. // 初始化常数数矩阵
  634. mtxB.FromString(m_strConstant, ",");
  636. // 构造线性方程组
  637. CLEquations leqs(mtxA, mtxB);
  639. // 求解托伯利兹方程组的列文逊方法
  640. if (leqs.GetRootsetTlvs(mtxC))
  641. m_strResult = mtxC.ToString(" ");
  642. else
  643. m_strResult = "求解失败";
  645. // 显示数据
  646. AfxMessageBox(m_strResult, MB_OK|MB_ICONINFORMATION);
  647. }
  649. /*
  650. 7.0,2.0,1.0,-2.0
  651. 9.0,15.0,3.0,-2.0
  652. -2.0,-2.0,11.0,5.0
  653. 1.0,3.0,2.0,13.0
  655. 4.0
  656. 7.0
  657. -1.0
  658. 0.0
  659. */
  660. // 高斯－赛德尔迭代法
  661. void CLECalculatorDlg::OnButton11() 
  662. {
  663. // 获取数据
  664. UpdateData();
  666. // 构造矩阵
  667. CMatrix mtxA, mtxB, mtxC;
  669. // 初始化系数矩阵
  670. mtxA.FromString(m_strCoefficient, ",");
  671. // 初始化常数数矩阵
  672. mtxB.FromString(m_strConstant, ",");
  674. // 构造线性方程组
  675. CLEquations leqs(mtxA, mtxB);
  677. // 高斯－赛德尔迭代法
  678. if (leqs.GetRootsetGaussSeidel(mtxC))
  679. m_strResult = mtxC.ToString(" ");
  680. else
  681. m_strResult = "求解失败";
  683. // 显示数据
  684. AfxMessageBox(m_strResult, MB_OK|MB_ICONINFORMATION);
  685. }
  687. /*
  688. 5.0,7.0,6.0,5.0
  689. 7.0,10.0,8.0,7.0
  690. 6.0,8.0,10.0,9.0
  691. 5.0,7.0,9.0,10.0
  693. 23.0
  694. 32.0
  695. 33.0
  696. 31.0
  697. */
  698. // 求解对称正定方程组的共轭梯度法
  699. void CLECalculatorDlg::OnButton12() 
  700. {
  701. // 获取数据
  702. UpdateData();
  704. // 构造矩阵
  705. CMatrix mtxA, mtxB, mtxC;
  707. // 初始化系数矩阵
  708. mtxA.FromString(m_strCoefficient, ",");
  709. // 初始化常数数矩阵
  710. mtxB.FromString(m_strConstant, ",");
  712. // 构造线性方程组
  713. CLEquations leqs(mtxA, mtxB);
  715. // 求解对称正定方程组的共轭梯度法
  716. leqs.GetRootsetGrad(mtxC);
  717. m_strResult = mtxC.ToString(" ");
  719. // 显示数据
  720. AfxMessageBox(m_strResult, MB_OK|MB_ICONINFORMATION);
  721. }
  723. /*
  724. 1.0,1.0,-1.0
  725. 2.0,1.0,0.0
  726. 1.0,-1.0,0.0
  727. -1.0,2.0,1.0
  729. 2.0
  730. -3.0
  731. 1.0
  732. 4.0
  733. */
  734. // 求解线性最小二乘问题的豪斯荷尔德变换法
  735. void CLECalculatorDlg::OnButton13() 
  736. {
  737. // 获取数据
  738. UpdateData();
  740. // 构造矩阵
  741. CMatrix mtxA, mtxB, mtxC, mtxQ, mtxR;
  743. // 初始化系数矩阵
  744. mtxA.FromString(m_strCoefficient, ",");
  745. // 初始化常数数矩阵
  746. mtxB.FromString(m_strConstant, ",");
  748. // 构造线性方程组
  749. CLEquations leqs(mtxA, mtxB);
  751. // 求解线性最小二乘问题的豪斯荷尔德变换法
  752. if (leqs.GetRootsetMqr(mtxC, mtxQ, mtxR))
  753. {
  754. m_strResult = "解矩阵为：rn";
  755. m_strResult += mtxC.ToString(" ");
  756. m_strResult += "rnrnQ矩阵的为：rn";
  757. m_strResult += mtxQ.ToString(" ");
  758. m_strResult += "rnrnR矩阵的为：rn";
  759. m_strResult += mtxR.ToString(" ");
  760. }
  761. else
  762. m_strResult = "求解失败";
  764. // 显示数据
  765. AfxMessageBox(m_strResult, MB_OK|MB_ICONINFORMATION);
  766. }
  768. /*
  769. 1.0,1.0,-1.0
  770. 2.0,1.0,0.0
  771. 1.0,-1.0,0.0
  772. -1.0,2.0,1.0
  774. 2.0
  775. -3.0
  776. 1.0
  777. 4.0
  778. */
  779. // 求解线性最小二乘问题的广义逆法
  780. void CLECalculatorDlg::OnButton14() 
  781. {
  782. // 获取数据
  783. UpdateData();
  785. // 构造矩阵
  786. CMatrix mtxA, mtxB, mtxC, mtxAP, mtxU, mtxV;
  788. // 初始化系数矩阵
  789. mtxA.FromString(m_strCoefficient, ",");
  790. // 初始化常数数矩阵
  791. mtxB.FromString(m_strConstant, ",");
  793. // 构造线性方程组
  794. CLEquations leqs(mtxA, mtxB);
  796. // 求解线性最小二乘问题的广义逆法
  797. if (leqs.GetRootsetGinv(mtxC, mtxAP, mtxU, mtxV))
  798. {
  799. m_strResult = "解矩阵为：rn";
  800. m_strResult += mtxC.ToString(" ");
  801. m_strResult += "rnrn广义逆矩阵的为：rn";
  802. m_strResult += mtxAP.ToString(" ");
  803. }
  804. else
  805. m_strResult = "求解失败";
  807. // 显示数据
  808. AfxMessageBox(m_strResult, MB_OK|MB_ICONINFORMATION);
  809. }
  811. /*
  812. 3.4336,-0.5238,0.67105,-0.15272
  813. -0.5238,3.28326,-0.73051,-0.2689
  814. 0.67105,-0.73051,4.02612,0.01835
  815. -0.15272,-0.2689,0.01835,2.75702
  817. -1.0
  818. 1.5
  819. 2.5
  820. -2.0
  821. */
  822. // 病态方程组的求解
  823. void CLECalculatorDlg::OnButton15() 
  824. {
  825. // 获取数据
  826. UpdateData();
  828. // 构造矩阵
  829. CMatrix mtxA, mtxB, mtxC;
  831. // 初始化系数矩阵
  832. mtxA.FromString(m_strCoefficient, ",");
  833. // 初始化常数数矩阵
  834. mtxB.FromString(m_strConstant, ",");
  836. // 构造线性方程组
  837. CLEquations leqs(mtxA, mtxB);
  839. // 病态方程组的求解
  840. if (leqs.GetRootsetMorbid(mtxC))
  841. m_strResult = mtxC.ToString(" ");
  842. else
  843. m_strResult = "求解失败";
  845. // 显示数据
  846. AfxMessageBox(m_strResult, MB_OK|MB_ICONINFORMATION);
  847. }