【C++ 程序】 解线性方程组(Cramer法则)
生活随笔
收集整理的這篇文章主要介紹了
【C++ 程序】 解线性方程组(Cramer法则)
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
求解結果為小數。結果為分數和小數兼容的程序詳見我的博客 【C++ 程序】 解線性方程組(Cramer法則)(分數形式結果)。
程序
Source.cpp
#include <iostream> #include "Det.h" using namespace std;int main() {int num;unsigned m = 0;deque<deque<double>> mat;deque<deque<deque<double>>> D;// instructioncout << "Please input like this:\n"<< "a_11 * x_1 + a_12 * x_2 + ... + a_1n * x_n = b_1\n"<< "a_21 * x_1 + a_22 * x_2 + ... + a_2n * x_n = b_2\n"<< "...\n"<< "a_n1 * x_1 + a_n2 * x_2 + ... + a_nn * x_n = b_n"<< "\n" << endl;// inputwhile (1){deque<double> row;while (cin >> num){row.push_back(num);if (getchar() == '\n') break;}mat.push_back(row);m = (row.size() > m) ? row.size() : m;if (mat.size() >= m - 1) break;}// calculate & printfor (unsigned i = 0; i < m; i++){deque<deque<double>> temp = mat; // D[m - 1][m]for (unsigned j = 0; j != m - 1; j++){if (i) temp[j][i - 1] = mat[j][m - 1];temp[j].pop_back();}D.push_back(temp);}double D0 = det_is(D[0]);if (D0 != 0){cout << endl;for (unsigned i = 1; i < m; i++){cout << "x_" << i << " = " << det_is(D[i]) / D0 << endl;}}else cout << "\nNo finite solution!" << endl;return 0; }Det.h
見我的博客 【C++ 程序】 行列式。
輸出示例
分析
- 一開始老是錯:
查了半天發現錯誤在頭文件中:
一開始一直沒加sign = 1的語句,造成了全局變量sign一直處于一個錯誤的位置。
啟示:每一次調用頭文件函數時,并不是重新運行的!運行一次造成的影響會帶給下一次運行。
- 這里面運用到行列式的求法。
- (其他分析以后繼續)
ALL RIGHTS RESERVED ? 2020 Teddy van Jerry
歡迎轉載,轉載請注明出處。
See also
Teddy van Jerry 的導航頁
總結
以上是生活随笔為你收集整理的【C++ 程序】 解线性方程组(Cramer法则)的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: Cramer`s Rule 克莱姆法则(
- 下一篇: java 分批同时处理_java中Lis