CMatrix

传说中的CMatrix,支持实矩阵和复数矩阵，支持矩阵的加减乘除求逆运算，矩阵的字符串输入输出，以及：实矩阵求逆的全选主元高斯－约当法复矩阵求逆的全选主元高斯－约当法对称正定矩阵的求逆托伯利兹矩阵求逆的埃兰特方法求行列式值的全选主元高斯消去法求矩阵秩的全选主元高斯消去法对称正定矩阵的乔里斯基分解与行列式的求值矩阵的三角分解一般实矩阵的QR分解一般实矩阵的奇异值分解求广义逆的奇异值分解法约化对称矩阵为对称三对角阵的豪斯荷尔德变换法实对称三对角阵的全部特征值与特征向量的计算约化一般实矩阵为赫申伯格矩阵的初等相似变换法求赫申伯格矩阵全部特征值的QR方法求实对称矩阵特征值与特征向量的雅可比法求实对称矩阵特征值与特征向量的雅可比过关法

// Matrix.cpp: implementation of the CMatrix class.
//
//////////////////////////////////////////////////////////////////////

#include “stdafx.h“
//#include “ldi.h“
#include “Matrix.h“

#ifdef _DEBUG
#undef THIS_FILE
static char THIS_FILE[]=__FILE__;
#define new DEBUG_NEW
#endif

//////////////////////////////////////////////////////////////////////
// Construction/Destruction
//////////////////////////////////////////////////////////////////////

//////////////////////////////////////////////////////////////////////
// 基本构造函数
//////////////////////////////////////////////////////////////////////
CMatrix::CMatrix（）
{
	m_nNumColumns = 1;
	m_nNumRows = 1;
	m_pData = NULL;
	BOOL bSuccess = Init（m_nNumRows m_nNumColumns）;
	ASSERT（bSuccess）;
}

//////////////////////////////////////////////////////////////////////
// 指定行列构造函数
//
// 参数：
// 1. int nRows - 指定的矩阵行数
// 2. int nCols - 指定的矩阵列数
//////////////////////////////////////////////////////////////////////
CMatrix::CMatrix（int nRows int nCols）
{
	m_nNumRows = nRows;
	m_nNumColumns = nCols;
	m_pData = NULL;
	BOOL bSuccess = Init（m_nNumRows m_nNumColumns）;
	ASSERT（bSuccess）;
}

//////////////////////////////////////////////////////////////////////
// 指定值构造函数
//
// 参数：
// 1. int nRows - 指定的矩阵行数
// 2. int nCols - 指定的矩阵列数
// 3. double value[] - 一维数组，长度为nRows*nCols，存储矩阵各元素的值
//////////////////////////////////////////////////////////////////////
CMatrix::CMatrix（int nRows int nCols double value[]）
{
	m_nNumRows = nRows;
	m_nNumColumns = nCols;
	m_pData = NULL;
	BOOL bSuccess = Init（m_nNumRows m_nNumColumns）;
	ASSERT（bSuccess）;

	SetData（value）;
}

//////////////////////////////////////////////////////////////////////
// 方阵构造函数
//
// 参数：
// 1. int nSize - 方阵行列数
//////////////////////////////////////////////////////////////////////
CMatrix::CMatrix（int nSize）
{
	m_nNumRows = nSize;
	m_nNumColumns = nSize;
	m_pData = NULL;
	BOOL bSuccess = Init（nSize nSize）;
	ASSERT （bSuccess）;
}

//////////////////////////////////////////////////////////////////////
// 方阵构造函数
//
// 参数：
// 1. int nSize - 方阵行列数
// 2. double value[] - 一维数组，长度为nRows*nRows，存储方阵各元素的值
//////////////////////////////////////////////////////////////////////
CMatrix::CMatrix（int nSize double value[]）
{
	m_nNumRows = nSize;
	m_nNumColumns = nSize;
	m_pData = NULL;
	BOOL bSuccess = Init（nSize nSize）;
	ASSERT （bSuccess）;

	SetData（value）;
}

//////////////////////////////////////////////////////////////////////
// 拷贝构造函数
//
// 参数：
// 1. const CMatrix& other - 源矩阵
//////////////////////////////////////////////////////////////////////
CMatrix::CMatrix（const CMatrix& other）
{
	m_nNumColumns = other.GetNumColumns（）;
	m_nNumRows = other.GetNumRows（）;
	m_pData = NULL;
	BOOL bSuccess = Init（m_nNumRows m_nNumColumns）;
	ASSERT（bSuccess）;

	// copy the pointer
	memcpy（m_pData other.m_pData sizeof（double）*m_nNumColumns*m_nNumRows）;
}

//////

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      75899  2005-09-12 22:04  Matrix.cpp

     文件       3212  2005-09-12 16:59  Matrix.h

----------- ---------  ---------- -----  ----

                79111                    2