?? matrix.cpp
字號:
/////////////////////////////////////////////////////////////////////////////
// Matrix.cpp : Implementation of the class Matrix
//
/////////////////////////////////////////////////////////////////////////////
#include "stdafx.h"
#include "Matrix.h"
#include <math.h>
#include <stdlib.h>
#include <conio.h>
#include <stdio.h>
#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif
//////////////////////////////////////////////////////////////////////
// Construction/Destruction
//////////////////////////////////////////////////////////////////////
CMatrix::CMatrix()
{
m_nRow = 0;
m_nCol = 0;
m_pTMatrix.resize (m_nRow);
for(unsigned int i=0; i < m_nRow; i++)
{
for(unsigned int j=0; j < m_nCol; j++)
{
m_pTMatrix[i].resize (m_nCol);
m_pTMatrix[i][j] = (double) 0;
}
}
}
CMatrix::~CMatrix()
{
}
CMatrix::CMatrix(unsigned int nRow,unsigned int nCol)
{
// 動態(tài)分配二維數(shù)組
TMatrix tMatrix;
tMatrix.resize (nRow);
for(unsigned int i=0; i < nRow; i++)
{
for(unsigned int j=0; j < nCol; j++)
{
tMatrix[i].resize(nCol);
tMatrix[i][j] = (double) 0;
}
}
// 對對象變量賦值
m_nRow = nRow;
m_nCol = nCol;
m_pTMatrix = tMatrix;
}
CMatrix::CMatrix(CMatrix& cMatrixB)
{
// Initialize the variable
m_nRow = cMatrixB.m_nRow ;
m_nCol = cMatrixB.m_nCol ;
m_pTMatrix = cMatrixB.m_pTMatrix ;
// Copy Data
for(unsigned int i=0; i< cMatrixB.m_nRow; i++)
{
for(unsigned int j=0; j < cMatrixB.m_nCol; j++)
{
m_pTMatrix [i][j] = cMatrixB.m_pTMatrix [i][j];
}
}
}
/////////////////////////////////////////////////////////////////////////////
// CMatrix member functions
//
CMatrix CMatrix::operator +(CMatrix& cMatrixB)
{
// 要滿足矩陣相加的條件: 行列數(shù)目相等!
if(m_nRow != cMatrixB.m_nRow || m_nCol != cMatrixB.m_nCol )
{
::AfxMessageBox (TEXT("執(zhí)行相加的兩個矩陣維數(shù)不相等!"),MB_OK | MB_ICONERROR);
}
CMatrix cMatrix = *this;
for(unsigned int i=0; i < m_nRow; i++)
{
for(unsigned int j=0; j < m_nCol; j++)
{
cMatrix.m_pTMatrix [i][j] = m_pTMatrix [i][j] + cMatrixB.m_pTMatrix [i][j];
}
}
return cMatrix;
}
CMatrix CMatrix::operator -(CMatrix& cMatrixB)
{
// 要滿足矩陣相加的條件: 行列數(shù)數(shù)目相等!
if(m_nRow != cMatrixB.m_nRow || m_nCol != cMatrixB.m_nCol )
{
::AfxMessageBox (TEXT("執(zhí)行相減的兩個矩陣維數(shù)不相等!"),MB_OK | MB_ICONERROR);
}
CMatrix cMatrix = *this;
for(unsigned int i=0; i < m_nRow; i++)
{
for(unsigned int j=0; j < m_nCol; j++)
{
cMatrix.m_pTMatrix [i][j] = m_pTMatrix [i][j] - cMatrixB.m_pTMatrix [i][j];
}
}
return cMatrix;
}
CMatrix CMatrix::operator *(CMatrix& cMatrixB)
{
if( m_nCol != cMatrixB.m_nRow )
{
::AfxMessageBox (TEXT("執(zhí)行相乘的兩個矩陣維數(shù)不滿足相乘的條件!"),MB_OK | MB_ICONERROR);
}
CMatrix cResultMatrix(m_nRow,cMatrixB.m_nCol);
for(unsigned int i=0; i < m_nRow; i++)
{
for(unsigned int j=0; j < cMatrixB.m_nCol; j++)
{
for(unsigned int m=0; m < m_nCol; m++)
{
cResultMatrix.m_pTMatrix [i][j] += m_pTMatrix [i][m] * cMatrixB.m_pTMatrix [m][j];
}
}
}
return cResultMatrix;
}
CMatrix CMatrix::operator * (double nValue)
{
CMatrix cMatrix = *this;
for(unsigned int i=0; i < m_nRow; i++)
{
for(unsigned int j=0; j < m_nCol; j++)
{
cMatrix.m_pTMatrix [i][j] =m_pTMatrix [i][j] * nValue;
}
}
return cMatrix;
}
CMatrix& CMatrix::operator =(CMatrix& cMatrixB)
{
if( (m_nRow != cMatrixB.m_nRow) || (m_nCol != cMatrixB.m_nCol) )
{
::AfxMessageBox(TEXT("等號左右兩邊的矩陣的維數(shù)不相等!"),MB_OK | MB_ICONERROR);
return *this; // return invalid value
}
// 給變量賦值
m_nRow = cMatrixB.m_nRow ;
m_nCol = cMatrixB.m_nCol ;
m_pTMatrix = cMatrixB.m_pTMatrix ;
// 賦值操作
for(unsigned int i=0; i < cMatrixB.m_nRow; i++)
{
for(unsigned int j=0; j< cMatrixB.m_nCol; j++)
{
m_pTMatrix [i][j] = cMatrixB.m_pTMatrix [i][j];
}
}
return *this;
}
CMatrix& CMatrix::operator += (CMatrix& cMatrixB)
{
if(m_nRow != cMatrixB.m_nRow || m_nCol != cMatrixB.m_nCol )
{
//printf("錯誤!執(zhí)行相加的兩個矩陣維數(shù)不相等!\n");
::AfxMessageBox (TEXT("運算符的兩邊矩陣的維數(shù)不相等!"),MB_OK | MB_ICONERROR);
return *this; // return invalid value
}
// 賦值操作
for(unsigned int i=0; i < cMatrixB.m_nRow; i++)
{
for(unsigned int j=0; j< cMatrixB.m_nCol; j++)
{
m_pTMatrix [i][j] += cMatrixB.m_pTMatrix [i][j];
}
}
return *this;
}
CMatrix CMatrix::Transpose()
{
CMatrix cMatrix(m_nCol,m_nRow);
for(unsigned int i=0; i < m_nRow; i++)
{
for(unsigned int j=0; j < m_nCol; j++)
{
cMatrix.m_pTMatrix [j][i] = m_pTMatrix [i][j];
}
}
return cMatrix;
}
/////////////////////////////////////////////////////////////////////////////
// 將矩陣的所有的元素按列合成一列
// 例如：
// matrix = [
// 1 2 3
// 4 5 6
// 7 8 9
// ]
// CMatrix cMatrix = matrix.MergeColumnsToColumnVector();
// cMatrix =
// [ 1
// 4
// 7
// 2
// 5
// 8
// 3
// 6
// 9 ]
/////////////////////////////////////////////////////////////////////////////
CMatrix CMatrix::MergeColumnsToColumnVector()
{
CMatrix cMatrix(m_nRow * m_nCol,(unsigned int)1);
// 對矩陣賦值
for(unsigned int j=0; j < m_nCol; j++)
{
for(unsigned int i=0; i < m_nRow; i++)
{
cMatrix.m_pTMatrix [i + j * m_nRow][(unsigned int)0] = m_pTMatrix [i][j];
}
}
return cMatrix;
}
/////////////////////////////////////////////////////////////////////////////
// Get the total value of the matrix
/////////////////////////////////////////////////////////////////////////////
double CMatrix::GetTotalElementValue()
{
double nTotalValue = 0;
for(unsigned int i=0; i < m_nRow; i++)
{
for( unsigned int j=0; j < m_nCol; j++)
{
nTotalValue += m_pTMatrix [i][j];
}
}
return nTotalValue;
}
/////////////////////////////////////////////////////////////////////////////
// Get System Error
/////////////////////////////////////////////////////////////////////////////
double CMatrix::GetSystemError() const
{
double nSystemError = 0;
for(unsigned int i=0; i < m_nRow; i++)
{
for( unsigned int j=0; j < m_nCol; j++)
{
nSystemError += m_pTMatrix [i][j] * m_pTMatrix [i][j];
}
}
return nSystemError;
}
/////////////////////////////////////////////////////////////////////////////
// Make all the matrix elements to be changed into absolute value
/////////////////////////////////////////////////////////////////////////////
CMatrix CMatrix::AbsoluteValue ()
{
CMatrix cMatrix = *this;
for(unsigned int i=0; i < m_nRow; i++)
{
for(unsigned int j=0; j < m_nCol; j++)
{
cMatrix.m_pTMatrix [i][j] = fabs( m_pTMatrix [i][j]);
}
}
return cMatrix;
}
CMatrix CMatrix::Inverse()
{
/////////////////////////////////////////////////////////////////////////
// Using Gauss - Jordan Method
// 參考書目: 計算機數(shù)值方法 --->施吉林 陳桂枝
/////////////////////////////////////////////////////////////////////////
/////////////////////////////////////////////////////////////////////////
// 判斷是否是可逆陣:
// 可逆陣一定是方陣!!!
if ( m_nRow != m_nCol)
{
//printf("錯誤!矩陣的行列數(shù)不相等,是非可逆陣!\n");
::AfxMessageBox (TEXT("矩陣的行列數(shù)不相等,是非可逆陣!"),MB_OK | MB_ICONERROR);
}
CMatrix cMatrix = *this;
//***********************************************************************
// 思路:(非常規(guī)思維!)
// 動態(tài)分配整型數(shù)組(2*m_nCol)來存儲每次交換的行坐標的值
// 不論有沒有行交換都記錄在數(shù)組中,
// 1.沒進行行交換的兩個數(shù)據(jù)相等,在SwapMatrixRow()函數(shù)中
// 檢測到兩個值相等立即返回,在SwapMatrixCol()函數(shù)中也一樣,
// 檢測到兩個值相等立即返回,不占用系統(tǒng)資源;
// 2.不相等的就交換
//***********************************************************************
// 分配內(nèi)存
int *pIntArray = new int [2*m_nCol];
// nSetp --- 約化步數(shù),按列展開
for(unsigned int k=0; k < cMatrix.m_nCol; k++)
{
/////////////////////////////////////////////////////////////////////
// 進行行交換 ---> 游戲規(guī)則:
// 為保證計算過程的數(shù)值穩(wěn)定性,在第k步約化時,先在{a(ik)}|i=k->n中選按
// 模最大者作為約化主元素,并交換矩陣相應(yīng)的行
// 標記主元素
double nMaxElement = cMatrix.m_pTMatrix [k][k];
// 標記主元素所在的行數(shù)
unsigned int nMainRow = k;
for(unsigned int nCount = k+1; nCount < cMatrix.m_nCol; nCount++)
{
if( fabs(nMaxElement) < fabs(cMatrix.m_pTMatrix [nCount][k]) )
{
nMaxElement = cMatrix.m_pTMatrix [nCount][k];
nMainRow = nCount;
}
}
// 將欲交換的行數(shù)存在數(shù)組中
pIntArray [2*k] = k;
pIntArray [2*k+1] = nMainRow;
// 交換行
cMatrix.SwapMatrixRow(k,nMainRow);
//Display();
// 判斷是否是可逆陣
if(cMatrix.m_pTMatrix [k][k] == 0)
{
//printf("錯誤!此矩陣為非可逆陣!\n");
::AfxMessageBox (TEXT("此矩陣為非可逆陣,沒有逆矩陣!"),MB_OK | MB_ICONERROR);
}
cMatrix.m_pTMatrix [k][k] = 1/(cMatrix.m_pTMatrix [k][k]);
// 算主列
for(unsigned int i=0; i < cMatrix.m_nRow; i++)
{
if( i != k)
cMatrix.m_pTMatrix [i][k] = -(cMatrix.m_pTMatrix [k][k]) * (cMatrix.m_pTMatrix [i][k]);
//int nTempValue = m_pTMatrix [i][k];
}
//printf("\n");
// 約化非主行
for(unsigned int m=0; m < cMatrix.m_nRow; m++)
{
if ( m == k)
continue;
for(unsigned int n=0; n < cMatrix.m_nCol; n++)
{
if ( n == k)
continue;
cMatrix.m_pTMatrix [m][n] += cMatrix.m_pTMatrix [m][k] * cMatrix.m_pTMatrix [k][n];
//printf("%10f ",m_pTMatrix [m][n]);
}
//printf("\n");
}
// 算主行
for(unsigned int j=0; j < cMatrix.m_nCol; j++)
{
if( j != k)
cMatrix.m_pTMatrix [k][j] = (cMatrix.m_pTMatrix [k][k]) * (cMatrix.m_pTMatrix [k][j]);
}
}
/////////////////////////////////////////////////////////////////////////
// 進行列交換 ---> 對交換行后的矩陣進行列交換 ---> 還原矩陣
// 游戲規(guī)則:
// 將開始矩陣中進行的行交換 ---> 現(xiàn)用相對應(yīng)的列交換進行還原,即可得到所求的
// 逆矩陣
for(int i=2*m_nCol-1; i > 0; i--)
{
cMatrix.SwapMatrixCol(pIntArray[i],pIntArray[i-1]);
i--;
}
delete []pIntArray;
return cMatrix;
}
void CMatrix::SwapMatrixRow(unsigned int nRow1,unsigned int nRow2)
{
if( nRow1 == nRow2)
return;
double *pArray = new double;
for(unsigned int i=0; i < m_nCol; i++)
{
// Swap the datum of the two rows
pArray[0] = m_pTMatrix [nRow1][i];
m_pTMatrix [nRow1][i] = m_pTMatrix [nRow2][i];
m_pTMatrix [nRow2][i] = pArray[0];
}
delete pArray;
}
void CMatrix::SwapMatrixCol(unsigned int nCol1,unsigned int nCol2)
{
if( nCol1 == nCol2)
return;
double *pArray = new double;
for(unsigned int i=0; i < m_nRow; i++)
{
// Swap the datum of the two columns
pArray[0] = m_pTMatrix [i][nCol1];
m_pTMatrix [i][nCol1] = m_pTMatrix [i][nCol2];
m_pTMatrix [i][nCol2] = pArray[0];
}
delete pArray;
}
?? 快捷鍵說明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -