矩陣求逆（c++）

標(biāo)簽（空格分隔）： 技術(shù)博客

---

簡(jiǎn)要過(guò)程介紹

方法的名稱是“Gauss-Jordan (or reduced row) elimination method”。

設(shè)單位對(duì)角矩陣為I，則MM?1=I

主要過(guò)程為，擺一個(gè)相同大小的對(duì)角矩陣在旁邊，將原矩陣變成對(duì)角矩陣的過(guò)程中，對(duì)對(duì)角矩陣施以相同的變化。原理為，對(duì)矩陣施以特定變化等同于對(duì)矩陣進(jìn)行線性計(jì)算。

實(shí)現(xiàn)過(guò)程

第一步：

準(zhǔn)備階段：進(jìn)行 行與行的變換，使矩陣對(duì)角位的數(shù)值非0。

過(guò)程如下：
+ 按順序我們先從第一行開始。
+ 查看后面所有行中位于第一個(gè)位置的元素的絕對(duì)值，找到絕對(duì)值最大的那一行，將其與第一行位置交換。
+ 如果絕對(duì)值最大為0，此矩陣不可逆，退出。
+ 緊接著做第二行，依舊查看后續(xù)行中位于第二列的元素中絕對(duì)值，將絕對(duì)值最大的行與第二行交換。

代碼為一個(gè)4*4的矩陣求逆（4*4矩陣在圖形學(xué)中用途最廣）

int i, j, k;Matrix44<T> s;Matrix44<T> t(*this);for (i = 0; i < 3; i++) {//找到下三角每列的絕對(duì)最大值int pivot = i;T pivotsize = t[i][i];if (pivotsize < 0)pivotsize = -pivotsize;for (j = i + 1; j < 4; j++) {T tmp = t[j][i];if (tmp < 0)tmp = -tmp;if (tmp > pivotsize) {pivot = j;pivotsize = tmp;}}if (pivotsize == 0) {//can not inversereturn Matrix44();}if (pivot != i)//交換兩行，使對(duì)角位的值為該列最大{for (j = 0; j < 4; j++) {T tmp;tmp = t[i][j];t[i][j] = t[pivot][j];t[pivot][j] = tmp;tmp = s[i][j];s[i][j] = s[pivot][j];s[pivot][j] = tmp;}}

第二步：

將下三角所有數(shù)值置為0。

對(duì)于交換后的每一行，從它的下一行開始進(jìn)行操作。
+ 對(duì)于第 i 行，那么從 i+1行開始，對(duì)于每一行，設(shè)定一個(gè)因子。
+ 該行-（第i行*因子），使該行的第 i 列的值為 0 。
結(jié)果為，做完第 i 行，后續(xù)所有行的第 i 列都為 0 。

//將下三角值設(shè)定為0for (j = i + 1; j < 4; j++) {T f = t[j][i] / t[i][i];for (k = 0; k < 4; k++) {t[j][k] -= t[i][k] * f;s[j][k] -= s[i][k] * f;}} }

第三步：

• 
  
• 先判斷現(xiàn)在是否有對(duì)角位為0 的情況，如果有，則證明矩陣不可逆。因?yàn)槿绻藭r(shí)對(duì)角位為0，則該行一定可以被其他行表示。
• 再將對(duì)角位置為1：
  

每一行都乘以一個(gè)相同因子使對(duì)角位都為1。
（現(xiàn)在就可以做這一步，因?yàn)楹罄m(xù)步驟并不會(huì)改變對(duì)角位了。）
后續(xù)需要把上三角都置為0，過(guò)程與第二步類似。

for (i = 3; i >= 0; i--) {T f;f = t[i][i];if (f == 0)return Matrix44<T>();for (j = 0; j < 4; j++) {//將對(duì)角位置為1t[i][j] /= f;s[i][j] /= f;}for (j = 0; j < i; j++) {f = t[j][i];for (k = 0; k < 4; k++) {t[j][k] -= f*t[i][k];s[j][k] -= f*s[i][k];}}}

完整代碼如下：

int i, j, k;Matrix44<T> s;Matrix44<T> t(*this);for (i = 0; i < 3; i++) {//找到下三角每列的絕對(duì)最大值int pivot = i;T pivotsize = t[i][i];if (pivotsize < 0)pivotsize = -pivotsize;for (j = i + 1; j < 4; j++) {T tmp = t[j][i];if (tmp < 0)tmp = -tmp;if (tmp > pivotsize) {pivot = j;pivotsize = tmp;}}if (pivotsize == 0) {//can not inversereturn Matrix44();}if (pivot != i)//交換兩行，使對(duì)角位的值為該列最大{for (j = 0; j < 4; j++) {T tmp;tmp = t[i][j];t[i][j] = t[pivot][j];t[pivot][j] = tmp;tmp = s[i][j];s[i][j] = s[pivot][j];s[pivot][j] = tmp;}}//將下三角值設(shè)定為0for (j = i + 1; j < 4; j++) {T f = t[j][i] / t[i][i];for (k = 0; k < 4; k++) {t[j][k] -= t[i][k] * f;s[j][k] -= s[i][k] * f;}} }for (i = 3; i >= 0; i--) {T f;f = t[i][i];if (f == 0)return Matrix44<T>();for (j = 0; j < 4; j++) {//將對(duì)角位置為1t[i][j] /= f;s[i][j] /= f;}for (j = 0; j < i; j++) {f = t[j][i];for (k = 0; k < 4; k++) {t[j][k] -= f*t[i][k];s[j][k] -= f*s[i][k];}}}return s; }

代碼來(lái)源：

https://www.scratchapixel.com

注：這是一個(gè)圖形學(xué)學(xué)習(xí)的網(wǎng)站

總結(jié)

以上是生活随笔為你收集整理的矩阵求逆（c++）的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。