高效求逆矩阵方法

}
}
//消元
for (i = k - 1;i >= 0;i--)
{
//获得倍数
b = (*(m.m_pfarr+i*m.m_icolumn+k)) / (*(m.m_pfarr+k*m.m_icolumn+k));
//行变换
for (j = 0;j < m.m_icolumn;j++)
{
temp = *(m.m_pfarr+i*m.m_icolumn+j) - b*(*(m.m_pfarr+k*m.m_icolumn+j));
{
if (j <= A->m_icolumn - 1)
{
*(m.m_pfarr+i*m.m_icolumn+j) = *(A->m_pfarr+i*A->m_icolumn+j);
}
else
{
if (i == j - A->m_icolumn)
{
*(m.m_pfarr+i*m.m_icolumn+j) = 1;
}
}
}
}
//高斯消元
//变换下三角
for (k = 0;k < m.m_irow - 1;k++)
{
//如果坐标为k,k的数为0,则行变换
//if (m.Read(k,k) == 0)
if (*(m.m_pfarr+k*m.m_icolumn+k)== 0)
{
for (i = k + 1;i < m.m_irow;i++)
*(m.m_pfarr+i*m.m_icolumn+j) = temp;
}
}
}
//将左边方阵化为单位矩阵
for (i = 0;i < m.m_irow;i++)
{
if (*(m.m_pfarr+i*m.m_icolumn+i) != 1 && *(m.m_pfarr+i*m.m_icolumn+i) != 0)
{
if (*(m.m_pfarr+i*m.m_icolumn+k) != 0)
{
break;
}
}
if (i >= m.m_irow)
{
return FALSE;
}
else
{
memcpy(pDtemp, m.m_pfarr+k*m.m_icolumn, sizeof(double)*m.m_icolumn);
memcpy(m.m_pfarr+k*m.m_icolumn, m.m_pfarr+(k+1)*m.m_icolumn, sizeof(double)*m.m_icolumn);
memcpy(m.m_pfarr+(k+1)*m.m_icolumn, pDtemp, sizeof(double)*m.m_icolumn);
}
}
//消元
for (i = k + 1;i < m.m_irow;i++)
{
//获得倍数
b = (*(m.m_pfarr+i*m.m_icolumn+k))/(*(m.m_pfarr+k*m.m_icolumn+k));
//行变换
for (j = 0;j < m.m_icolumn;j++)
{
temp = *(m.m_pfarr+i*m.m_icolumn+j) - b*(*(m.m_pfarr+k*m.m_icolumn+j));
{
return FALSE;
}
/*计算跌是否为0，如果为0，则是奇异矩阵
if (CMatrix_Calc::Det(A) <= 0 )
{
return FALSE;
} */
//增广矩阵m = A | B初始化
m.InitMatrix();
for (i = 0;i < m.m_irow;i++)
{
for (j = 0;j < m.m_icolumn;j++)
}
}
}
//求得逆矩阵
for (i = 0; i< B->m_irow; i++)
{
memcpy(B->m_pfarr+i*B->m_icolumn, m.m_pfarr+i*m.m_icolumn+m.m_irow, sizeof(double)*m.m_irow);
}
delete []pDtemp;
{
if (*(m.m_pfarr+i*m.m_icolumn+k) != 0)
{
break;
}
}
if (i >= m.m_irow)
{
return FALSE;
}
else
{
memcpy(pDtemp, m.m_pfarr+k*m.m_icolumn, sizeof(double)*m.m_icolumn);
/*函数名:Inverse
*作用：求出矩阵的逆矩阵,消元法
*输入参数：要求逆矩阵的矩阵A，目标矩阵B
*返回值：TRUE
*失败返回：FALSE;
* */
/************************************************************************/
bool CMatrix_Calc::Inverse(CMatrix *A, CMatrix *B)
{
int i = 0;
int j = 0;
int k = 0;
CMatrix m(A->m_irow,2 * A-&g 0;
double b = 0;
double *pDtemp = new double [m.m_icolumn];
//判断是否可以运算
if (A->m_irow != A->m_icolumn || B->m_irow != B->m_icolumn || A->m_irow != B->m_irow)
*(m.m_pfarr+i*m.m_icolumn+j) = temp;
}
}
}
//变换上三角
for (k = m.m_irow - 1;k > 0;k--)
{
//如果坐标为k,k的数为0,则行变换
if (*(m.m_pfarr+k*m.m_icolumn+k) == 0)
{
for (i = k - 1;i < m.m_irow;i--)
return TRUE;
}
//这代码非最优，不能判断是否存在奇异矩阵，因为用于高阶矩阵求行列值会消耗时间，所以不必判定。
//以下是类的对象，还有求逆矩阵的方法。只供学习参考，对于生产应用需要进一步优化。
//用高斯消元的方法，经过测试200阶矩阵运行需要几十毫秒。
class CMatrix
{
public:
int m_irow;
{
//获得倍数
b = 1 / (*(m.m_pfarr+i*m.m_icolumn+i));
//行变换
for (j = 0;j < m.m_icolumn;j++)
{
temp = (*(m.m_pfarr+i*m.m_icolumn+j))*b;
*(m.m_pfarr+i*m.m_icolumn+j) = temp;
int m_icolumn;
double *m_pfarr;//初始化需要赋值0
//构造
CMatrix(int mm = 0,int nn = 0);
//析构
~CMatrix(void);
};
/************************************************************************/
memcpy(m.m_pfarr+k*m.m_icolumn, m.m_pfarr+(k-1)*m.m_icolumn, sizeof(double)*m.m_icolumn);
memcpy(m.m_pfarr+(k-1)*m.m_icolumn, pDtemp, sizeof(double)*m.m_icolumn);