对称矩形C语言,c语言编程求任意对称正定矩阵的逆。
满意答案
玖林国际999
推荐于 2016.03.21
采纳率:58% 等级:13
已帮助:6009人
#ifndef __MATRIX_DOT_H__
#define __MATRIX_DOT_H__
template
void Swap(T& a, T& b)
{
T temp;
temp = a;
a = b;
b = temp;
}
class CMatrix
{
static double ZERO;//极小值
public:
CMatrix(int row = 3, int col = 3); //
CMatrix(const CMatrix& right); //拷贝构造函数
~CMatrix();
void Show(const char* pre = NULL)const; //输出矩阵
void Free();
int Resize(int row, int col); //重新定义矩阵大小
int GetRow()const{ return m_iRow; } //返回矩阵行数
int GetCol()const{ return m_iCol; } //返回矩阵列数
int RowSwap(int x, int y); //行交换,成功返回1,否则0
int ColSwap(int x, int y); //列交换
static void SetZero(double z); //设定ZERO的值,所有CMatrix实例精度都将改变
const CMatrix Transpose()const; //返回转置矩阵
const CMatrix Adjoint()const; //伴随矩阵
const CMatrix Residue(int row, int col)const;//求对应元素的余子式
const CMatrix Contrary()const;//逆矩阵
const CMatrix Gauss_Jordan(double* pDet = NULL)const;//逆矩阵(高斯-约旦法),pDet为行列式,
//此法精度较低,但效率较高
double Residue_a(int row, int col)const;//求对应元素的代数余子式
double Determinant()const; //返回方阵的行列式
double Det_Recursion()const; //返回方阵的行列式(递归)
int IsZero()const; //判断元素是否全为0(零矩阵)
int IsPhalanx()const; //判断是否为方阵
int IsReverse()const; //判断矩阵是否可逆
int IsNonfunnyPhalanx()const; //判断是否非奇异方阵
double* operator[](int i)const; //操作单个元素
CMatrix& operator=(const CMatrix& right);
CMatrix& operator=(const double* pRight);
CMatrix& operator=(const double** ppRight);
const CMatrix& operator+()const; //一元操作符
const CMatrix operator-()const; //一元操作符
const CMatrix operator+(const CMatrix& right)const;
const CMatrix operator-(const CMatrix& right)const;
const CMatrix operator*(const CMatrix& right)const;
const CMatrix operator*(const double& right)const;
const CMatrix operator/(const double& right)const;
CMatrix& operator+=(const CMatrix& right);
CMatrix& operator-=(const CMatrix& right);
CMatrix& operator*=(const CMatrix& right);
CMatrix& operator*=(const double& right);
CMatrix& operator/=(const double& right);
int operator==(const CMatrix& right)const;
int operator!=(const CMatrix& right)const;
private:
int m_iRow; //行数
int m_iCol; //列数
double** m_ppData; //数据
};
#endif //__MATRIX_DOT_H__
//matrix.cpp
//
#include
#include
#include
#include
#include
#include "matrix.h"
double CMatrix::ZERO = 1e-10;
CMatrix::CMatrix(int row/*=3*/, int col/*=3*/)
{
m_ppData = NULL;
Resize(row, col);
}
//拷贝构造函数
CMatrix::CMatrix(const CMatrix& right)
{
m_ppData = NULL;//一定要加这句初始化(一个对象不会同时调用构造函数和拷贝构造函数)
Resize(right.GetRow(), right.GetCol());
for(int i = 0; i < right.GetRow(); i++)
{
for(int j = 0; j < right.GetCol(); j++)
m_ppData[i][j] = right[i][j];
}
}
CMatrix::~CMatrix()
{
Free();
}
void CMatrix::Free()
{
if(m_ppData != NULL){
for(int i = 0; i < m_iRow; i++)
{
if(m_ppData[i] != NULL)
delete[] m_ppData[i];
m_ppData[i] = NULL;
}
m_ppData = NULL;
}
}
int CMatrix::Resize(int row, int col)
{
assert(row > 0 && col > 0);
//释放空间
Free();
//申请空间
m_iRow = row;
m_iCol = col;
m_ppData = new double*[m_iRow];
assert(m_ppData != NULL);
for(int i = 0; i < m_iRow; i++)
{
m_ppData[i] = new double[m_iCol];
assert(m_ppData[i] != NULL);
}
//初始化
for(i = 0; i < m_iRow; i++)
{
for(int j = 0; j < m_iCol; j++)
m_ppData[i][j] = 0;
}
return 1;
}
//zero
void CMatrix::SetZero(double z)
{
double zero = fabs(z);
if(zero > 1.0f) return;
ZERO = zero;
}
//show
void CMatrix::Show(const char* pre/*=NULL*/)const
{
int i, j;
#ifdef _WINDOWS
if(m_iRow > 10 || m_iCol > 10)
MessageBox(NULL, "矩阵数据量太大,不能输出", "警告", MB_OK);
char buf[4096];
char temp[256];
strcpy(buf, "");
if(pre != NULL)
{
strcpy(buf, pre);
strcat(buf, "\n");
}
for(i = 0; i < m_iRow; i++)
{
for(j = 0; j < m_iCol; j++)
{
sprintf(temp, "%.3f\t", m_ppData[i][j]);
strcat(buf, temp);
}
strcat(buf, "\n");
}
MessageBox(NULL, buf, "提示信息", MB_OK);
#else
if(pre != NULL)
puts(pre);
for(i = 0; i < m_iRow; i++)
{
for(j = 0; j < m_iCol; j++)
printf("%f\t", m_ppData[i][j]);
printf("\n");
}
#endif //_WINDOWS
}
///计算//
//行交换
int CMatrix::RowSwap(int x, int y)
{
if(x < 0 || x >= m_iRow || y < 0 || y >= m_iRow)
return 0;
if(x == y)
return 1;
for(int i = 0; i < m_iCol; i++)
{
Swap(m_ppData[x][i], m_ppData[y][i]);
}
return 1;
}
//列交换
int CMatrix::ColSwap(int x, int y)
{
if(x < 0 || x >= m_iCol || y < 0 || y >= m_iCol)
return 0;
if(x == y)
return 1;
for(int i = 0; i < m_iRow; i++)
{
Swap(m_ppData[i][x], m_ppData[i][y]);
}
return 1;
}
//转置
const CMatrix CMatrix::Transpose()const
{
CMatrix tr(m_iCol, m_iRow);
for(int i = 0; i < m_iRow; i++)
{
for(int j = 0; j < m_iCol; j++)
tr[j][i] = m_ppData[i][j];
}
return tr;
}
//计算方阵的行列式(精度不高)
double CMatrix::Determinant()const
{
assert(m_iRow == m_iCol);
CMatrix temp = *this;
int i,j,m,n,s,t,k=1;
double f=1,c,x,sn;
for (i=0,j=0; i
{
if (temp[i][j]==0)
{
for (m = i;m < m_iRow; m++)
{
if(fabs(temp[m][j]) > ZERO)//0
break;
}
if (m == m_iRow)
{
return 0;
}
else
for (n = j; n < m_iRow; n++)
{
c = temp[i][n];
temp[i][n] = temp[m][n];
temp[m][n] = c;
}
k *= (-1);
}
for (s = m_iRow-1; s>i; s--)
{
x = temp[s][j];
for (t = j; t < m_iRow; t++)
temp[s][t] -= temp[i][t] * (x/temp[i][j]);
}
}
for (i = 0; i < m_iRow; i++)
f *= temp[i][i];
sn = k * f;
return sn;
}
//行列式(递归,精度较高)
double CMatrix::Det_Recursion()const
{
assert(m_iRow == m_iCol);
CMatrix temp;
double ans = 0;
if(m_iRow == 1)
{
return m_ppData[0][0];
}
else if(m_iRow == 2)
{
return m_ppData[0][0]*m_ppData[1][1] - m_ppData[1][0]*m_ppData[0][1];
}
else
{
for(int i = 0; i < m_iRow; i++)
{
temp = Residue(i, 0);//this->Residue(i, 0)
ans += temp.Det_Recursion()*m_ppData[i][0]*pow(-1, i);
}
}
return ans;
}
//计算方阵的余子式
const CMatrix CMatrix::Residue(int row, int col)const
{
CMatrix re;
int index = 0;
assert(m_iRow == m_iCol);
assert(m_iRow >= 2);
assert(row < m_iRow && col < m_iCol);
assert(row >= 0 && col >= 0);
double* pData = NULL;
pData = new double[(m_iRow-1)*(m_iCol-1)];
assert(pData != NULL);
re.Resize(m_iRow-1, m_iCol-1);
for(int i = 0; i < m_iRow; i++)
{
for(int j = 0; j < m_iCol; j++)
if(i != row && j != col)
pData[index++] = m_ppData[i][j];
}
re = pData;
delete[] pData;
pData = NULL;
return re;
}
//计算方阵的代数余子式
double CMatrix::Residue_a(int row, int col)const
{
return (Residue(row, col)).Det_Recursion()*pow(-1, row+col);
}
//伴随矩阵
const CMatrix CMatrix::Adjoint()const
{
CMatrix ad(m_iRow, m_iCol);
for(int i = 0; i < m_iRow; i++)
{
for(int j = 0; j < m_iCol; j++)
ad[j][i] = Residue_a(i, j);
}
return ad;
}
//逆矩阵
const CMatrix CMatrix::Contrary()const
{
assert(IsReverse());
CMatrix co(m_iRow, m_iCol);
co = Adjoint();//this
co /= Det_Recursion();//this
return co;
}
//高斯-约旦法求逆矩阵(全选主元), pDet为原方阵的行列式
const CMatrix CMatrix::Gauss_Jordan(double* pDet/*=NULL*/)const
{
assert(IsReverse());
double fDet = 1.0f;
int flag = 1;
int k = 0, i = 0, j = 0;
CMatrix out(m_iRow, m_iCol);//逆
CMatrix m = *this;//原
CMatrix rhs(2, m_iRow);//保存主元素位置,0 i, 1 j;
for(k = 0; k < m_iRow; k++)
{
//第一步:全选主元
double fMax = 0.0f;
for(i = 0; i < m_iRow; i++)
{
for(j = 0; j < m_iCol; j++)
{
const double f = fabs(m[i][j]);
if(f > fMax)
{
fMax = f;
rhs[0][k] = i;
rhs[1][k] = j;
}
}
}
//if(fMax < 0.00001)//元素全为0
//{
// fDet = 0.0f;
// return out;
//}
if((int)rhs[0][k] != k)
{
flag = -flag;
m.RowSwap((int)rhs[0][k], k);
}
if((int)rhs[1][k] != k)
{
flag = -flag;
m.ColSwap((int)rhs[1][k], k);
}
//计算行列值
fDet *= m[k][k];
//计算逆矩阵
//第二步
m[k][k] = 1.0f/m[k][k];
//第三步
for(j = 0; j < m_iCol; j++)
{
if(j != k)
m[k][j] *= m[k][k];
}
//第四步
for(i = 0; i < m_iRow; i++)
{
if(i != k)
{
for(j = 0; j < m_iCol; j++)
{
if(j != k)
m[i][j] = m[i][j] - m[i][k]*m[k][j];
}
}
}
//第五步
for(i = 0; i < m_iRow; i++)
{
if(i != k)
{
m[i][k] *= -m[k][k];
}
}
}//end for(k);
for(k = m_iRow-1; k >= 0; k--)
{
if((int)rhs[1][k] != k)
{
m.RowSwap((int)rhs[1][k], k);
}
if((int)rhs[0][k] != k)
{
m.ColSwap((int)rhs[0][k], k);
}
}
fDet *= flag;
if(pDet != NULL)
*pDet = fDet;
return m;
}
///判断//
//零矩阵
int CMatrix::IsZero()const
{
for(int i = 0; i < m_iRow; i++)
{
for(int j = 0; j < m_iCol; j++)
if(fabs(m_ppData[i][j]) > ZERO)
return 0;
}
return 1;
}
//方阵
int CMatrix::IsPhalanx()const
{
return (m_iRow == m_iCol);
}
//非奇异方阵
int CMatrix::IsNonfunnyPhalanx()const
{
return (IsPhalanx() && fabs(Det_Recursion()) > ZERO);
}
//可逆矩阵
int CMatrix::IsReverse()const
{
return IsNonfunnyPhalanx();
}
///操作符重载//
// []
double* CMatrix::operator [](int i)const
{
assert(i >= 0 && i < m_iRow);
return m_ppData[i];
}
// =
CMatrix& CMatrix::operator =(const CMatrix& right)
{
if(this == &right) return *this;
if((m_iRow != right.GetRow())
|| (m_iCol != right.GetCol()))// 添加于 2005-11-09
{
Resize(right.GetRow(), right.GetCol());
}
for(int i = 0; i < m_iRow; i++)
{
for(int j = 0; j < m_iCol; j++)
m_ppData[i][j] = right[i][j];
}
return *this;
}
// =
CMatrix& CMatrix::operator =(const double* pRight)
{
assert(pRight != NULL);
for(int i = 0; i < m_iRow; i++)
{
for(int j = 0; j < m_iCol; j++)
m_ppData[i][j] = pRight[m_iCol*i + j];
}
return *this;
}
// =
CMatrix& CMatrix::operator =(const double** ppRight)
{
assert(ppRight != NULL);
for(int i = 0; i < m_iRow; i++)
{
for(int j = 0; j < m_iCol; j++)
m_ppData[i][j] = ppRight[i][j];
}
return *this;
}
// 一元操作符+
const CMatrix& CMatrix::operator +()const
{
return *this;
}
// 一元操作符-
const CMatrix CMatrix::operator -()const
{
CMatrix temp(m_iRow, m_iCol);
for(int i = 0; i < m_iRow; i++)
{
for(int j = 0; j < m_iCol; j++)
temp[i][j] = -m_ppData[i][j];
}
return temp;
}
// +
const CMatrix CMatrix::operator +(const CMatrix& right)const
{
CMatrix temp(m_iRow, m_iCol);
if(m_iRow == right.GetRow()
&& m_iCol == right.GetCol())
{
for(int i = 0; i < m_iRow; i++)
{
for(int j = 0; j < m_iCol; j++)
temp[i][j] = m_ppData[i][j] + right[i][j];
}
}
return temp;
}
// -
const CMatrix CMatrix::operator -(const CMatrix& right)const
{
CMatrix m_temp(m_iRow, m_iCol);
if(m_iRow == right.GetRow()
&& m_iCol == right.GetCol())
{
for(int i = 0; i < m_iRow; i++)
{
for(int j = 0; j < m_iCol; j++)
m_temp[i][j] = m_ppData[i][j] - right[i][j];
}
}
return m_temp;
}
// *
const CMatrix CMatrix::operator *(const CMatrix& right)const
{
double temp = 0;
CMatrix m_temp(m_iRow, right.GetCol());
if(m_iCol != right.GetRow())
return m_temp;
for(int i = 0; i < m_temp.GetRow(); i++)
{
for(int j = 0; j < m_temp.GetCol(); j++)
{
temp = 0;
for(int k = 0; k < right.GetRow(); k++)
temp += m_ppData[i][k] * right[k][j];
m_temp[i][j] = temp;
}
}
return m_temp;
}
// *
const CMatrix CMatrix::operator *(const double& right)const
{
CMatrix m_temp(m_iRow, m_iCol);
for(int i = 0; i < m_iRow; i++)
{
for(int j = 0; j < m_iCol; j++)
m_temp[i][j] = m_ppData[i][j] * right;
}
return m_temp;
}
// /
const CMatrix CMatrix::operator /(const double& right)const
{
CMatrix m_temp(m_iRow, m_iCol);
for(int i = 0; i < m_iRow; i++)
{
for(int j = 0; j < m_iCol; j++)
m_temp[i][j] = m_ppData[i][j] / right;
}
return m_temp;
}
// +=
CMatrix& CMatrix::operator +=(const CMatrix& right)
{
*this = (*this) + right;
return *this;
}
// -=
CMatrix& CMatrix::operator -=(const CMatrix& right)
{
*this = (*this) - right;
return *this;
}
// *=
CMatrix& CMatrix::operator *=(const CMatrix& right)
{
*this = (*this) * right;
return *this;
}
// *=
CMatrix& CMatrix::operator *=(const double& right)
{
*this = (*this) * right;
return *this;
}
// /=
CMatrix& CMatrix::operator /=(const double& right)
{
*this = (*this) / right;
return *this;
}
// ==
int CMatrix::operator ==(const CMatrix& right)const
{
if(this == &right) return 1;
if((m_iRow != right.GetRow())
|| (m_iCol != right.GetCol()))
{
return 0;
}
for(int i = 0; i < m_iRow; i++)
{
for(int j = 0; j < m_iCol; j++)
if(fabs(m_ppData[i][j] - right[i][j]) > ZERO)//0
return 0;
}
return 1;
}
// !=
int CMatrix::operator !=(const CMatrix& right)const
{
return !(*this == right);
}
00分享举报
对称矩形C语言,c语言编程求任意对称正定矩阵的逆。相关推荐
- C语言程序设计之编程求鸡和兔的只数,用穷举法解决
CSDN论坛上网友的问题:有一群鸡和一群兔,两种动物只数相同.两种动物的脚的总数都是三位数,且这两个三位数的六个数字分别是0,1,2,3,4,5.编程求鸡和兔的只数是多少?它们的脚数各是多少? C语言 ...
- C语言实验报告求三角形面积,C语言作业,编程求三角形面积。
满意答案 oikre5602 2014.05.08 采纳率:47% 等级:9 已帮助:363人 123456789#include<stdio.h>#include<math. ...
- c语言编程所得票数,C语言编程求1X2X3····Xn所得的数末尾有多少个零
C语言编程求1X2X3····Xn所得的数末尾有多少个零 发布时间:2020-08-10 02:23:57 来源:51CTO 阅读:312 作者:sonissa 参见大数的阶乘 https://blo ...
- 输入两个数、用python求他们的和_从键盘上输入俩个实型数,编程求它们的和差积商,要求输出时,保留两位小数C语言编程题:从键盘上输入两个...
从键盘上输入俩个实型数,编程求它们的和差积商,要求输出时,保留两位小数 C语言编程题:从键盘上输入两个 www.zhiqu.org 时间: 2020-11-23 参考代码: #includev ...
- c语言求100以内被7整除的最大自然数,编程,求100以内被7整除的最大自然数
java 编程,求1000以内能被17整除数的和 publicclassTest{publicstaticvoidmain(String[]args){intsum=0;for(inti=1;i c# ...
- python输入一个英文句子、统计单词个数_C语言编程求一个英文句子中的单词数和最长单词的位置、长度及输出这个单词。c++编程 从键盘输入一个英文...
C语言编程求一个英文句子中的单词数和最长单词的位置.长度及输出这个单词. c++编程 从键盘输入一个英文 www.zhiqu.org 时间: 2020-11-23 我刚做了一关于英文句子里面每 ...
- c语言输入一批正整数 求其中的偶数和,编程,输入一批整数,先求出其中的偶数和及奇数和,然后输出偶数和与奇数和的差...
设计一个程序,要求输入三个整数,能求出其中最大的数并输出 #includeintmain{inti;inta1,a2,a3;printf("输入三个数:\n");scanf(&qu ...
- a除以b的编程怎么写C语言,输入a,b两个整数,编程求出a除以b得到的商和余数.c语言...
C语言编程问题 编程任意输入两个实数a和b,计算下面函数的值,并输出y值. a-b a>b y= a/b a #includevoidmain(){\x09floata,b,y;\x09prin ...
- c语言坐标海伦公式,C语言:用海伦公式求三角形面积 , C语言编程问题,利用海伦公式求三角形面积...
导航:网站首页 > C语言:用海伦公式求三角形面积 , C语言编程问题,利用海伦公式求三角形面积 C语言:用海伦公式求三角形面积 , C语言编程问题,利用海伦公式求三角形面积 匿名网友: 程序已 ...
- c语言编程求二元一次方程组方程,二元一次方程组练习题 已知二元一次方程的三个系数,用C语言编程求方程的......
导航:网站首页 > 二元一次方程组练习题 已知二元一次方程的三个系数,用C语言编程求方程的... 二元一次方程组练习题 已知二元一次方程的三个系数,用C语言编程求方程的... 相关问题: 匿名网 ...
最新文章
- 混凝土静力受压弹性模量试验计算公式_混凝土静力受压弹性模量试验须进行最少()次预压。A.1B.2C.3D.4...
- 安装ESXi5.5遇到Relocating modules and starting up the kernel的处理
- Hadoop使用常见问题以及解决方法
- 国立大学 计算机,新加坡国立大学 计算机
- UnitTest in .NET(Part 2)
- 错误An entity with the same identity already exists in this EntitySet RIA
- linux导出mysql下ssl证书_Linux系统下生成证书 https证书
- kaggle案例实战
- 开启python之旅
- mac系统安装python(2.7)requests库
- DotNET企业架构应用实践-数据库表记录的唯一性设计的设计兼议主键设定原则
- 超级详细的Vmware下载与安装过程
- Vijos训练计划 1-1 计数问题
- 商家可以通过刷脸支付与顾客建立粘性效应
- 如何在Ubuntu 16.04上使用ProxySQL缓存优化MySQL查询
- 使用 Sqoop 将 30W+ MySQL 数据导入 Hive
- 网站图标开发指南(精)
- Axure第7讲:设置元件文字行距、边距
- 现代信息系统的企业驱动力
- 原来何恺明提出的MAE还是一种数据增强
热门文章
- “那个人样子好怪。”“我也看到了,他好像一条狗。”
- 摄氏度与华氏度之间的转换
- Cadence 17.4 PSpice仿真555定时器输出方波
- (一)外显子组数据分析之软件安装大全
- pbootcms模板,pbootcms模板程序,pboot响应式模板
- 如何用OpenCV给图片加上文字?
- 网站服务器要有固定IP地址吗,你知道动态IP和静态IP服务器有何区别吗?
- 英尺C语言,C语言中关于英尺、英寸、厘米的换算
- python图片标记_用python找出那些被“标记”的照片
- CatBoost快速入门