小衣服博客:矩阵求逆的算法怎么写
来源:百度文库 编辑:偶看新闻 时间:2024/05/03 15:18:04
//***************************
//求任何一个矩阵的逆矩阵
//***************************
#include
#include
void main( void )
{
float *buffer,*p; //定义数组首地址指针变量
short int row,num; //定义矩阵行数row及矩阵元素个数
short int i,j;
float determ; //定义矩阵的行列式
float comput_D(float *p,short int n); //求矩阵的行列式
float Creat_M(float *p, short int m,short int n,short int k); //求代数余子式
void Print( float *p,short int n); //打印n×n的矩阵
printf("\nPlease input the number of rows: ");
scanf("%d",&row);
num=2 * row * row;
buffer = (float *)calloc(num, sizeof(float)); //分配内存单元
p=buffer;
if(p != NULL)
{
for(i=0;i {
printf("Input the number of %d row ",i+1);
for(j=0;j {
scanf("%f",p++);
}
}
}
else
printf( "Can't allocate memory\n" );
printf("\nThe original matrix is:\n");
Print(buffer,row); //打印该矩阵
determ=comput_D(buffer,row); //求整个矩阵的行列式
p=buffer + row * row;
if (determ != 0)
{
for (i=0;i for (j=0; j *(p+j*row+i)= Creat_M(buffer,i,j,row)/determ;
printf("The determinant is %G\n",determ);
p=buffer + row * row;
printf("\nThe inverse matrix is:\n");
Print(p,row); //打印该矩阵
}
else
printf("The determnant is 0, and there is no inverse matrix !\n");
free( buffer );
}
//--------------------------------------------------------
//功能:求矩阵 n X n 的行列式
//入口参数:矩阵首地址 p;矩阵行数 n
//返回值:矩阵的行列式值
//--------------------------------------------------------
float comput_D(float *p,short int n)
{
short int i,j,m; //i--row; j--column
short int lop=0;
float result=0;
float mid=1;
if (n!=1)
{
lop=(n==2)?1:n; //控制求和循环次数,若为2阶,则循环1次,否则为n次
for(m=0;m {
mid=1; //顺序求和
for(i=0,j=m;i mid = mid * ( *(p+i*n+j%n) );
result+=mid;
}
for(m=0;m {
mid=1; //逆序相减
for(i=0,j=n-1-m+n; i mid=mid * ( *(p+i*n+j%n));
result-=mid;
}
}
else result=*p;
return(result);
}
//----------------------------------------------------
//功能:求k×k矩阵中元素A(mn)的代数余子式
//入口参数:k×k矩阵首地址;元素A的下标m,n; 矩阵行数 k
//返回值: k×k矩阵中元素A(mn)的代数余子式
//----------------------------------------------------
float Creat_M(float *p, short int m,short int n,short int k)
{
short int len;
short int i,j;
float mid_result=0;
short int quo=1;
float *p_creat,*p_mid;
len=(k-1)*(k-1);
p_creat = (float *)calloc(len, sizeof(float)); //分配内存单元
p_mid=p_creat;
for(i=0;i for(j=0;j {
if (i!=m && j!=n)
*p_mid++ =* (p+i*k+j);
}
// Print(p_creat,k-1);
quo = (m + n) %2==0 ? 1:-1;
mid_result = (float ) quo * comput_D(p_creat,k-1);
free(p_creat);
return(mid_result);
}
//-------------------------------------------
//功能:打印n×n的矩阵
//入口参数:n×n矩阵的首地址;该矩阵的行数 n
//返回值: 无
//-------------------------------------------
void Print( float *p,short int n)
{
int i,j;
for (i=0;i {
for (j=0; j printf("%10G ",*p++);
printf("\n");
}
printf("--------------\n");
}
//求任何一个矩阵的逆矩阵
//***************************
#include
#include
void main( void )
{
float *buffer,*p; //定义数组首地址指针变量
short int row,num; //定义矩阵行数row及矩阵元素个数
short int i,j;
float determ; //定义矩阵的行列式
float comput_D(float *p,short int n); //求矩阵的行列式
float Creat_M(float *p, short int m,short int n,short int k); //求代数余子式
void Print( float *p,short int n); //打印n×n的矩阵
printf("\nPlease input the number of rows: ");
scanf("%d",&row);
num=2 * row * row;
buffer = (float *)calloc(num, sizeof(float)); //分配内存单元
p=buffer;
if(p != NULL)
{
for(i=0;i
printf("Input the number of %d row ",i+1);
for(j=0;j
scanf("%f",p++);
}
}
}
else
printf( "Can't allocate memory\n" );
printf("\nThe original matrix is:\n");
Print(buffer,row); //打印该矩阵
determ=comput_D(buffer,row); //求整个矩阵的行列式
p=buffer + row * row;
if (determ != 0)
{
for (i=0;i
printf("The determinant is %G\n",determ);
p=buffer + row * row;
printf("\nThe inverse matrix is:\n");
Print(p,row); //打印该矩阵
}
else
printf("The determnant is 0, and there is no inverse matrix !\n");
free( buffer );
}
//--------------------------------------------------------
//功能:求矩阵 n X n 的行列式
//入口参数:矩阵首地址 p;矩阵行数 n
//返回值:矩阵的行列式值
//--------------------------------------------------------
float comput_D(float *p,short int n)
{
short int i,j,m; //i--row; j--column
short int lop=0;
float result=0;
float mid=1;
if (n!=1)
{
lop=(n==2)?1:n; //控制求和循环次数,若为2阶,则循环1次,否则为n次
for(m=0;m
mid=1; //顺序求和
for(i=0,j=m;i
result+=mid;
}
for(m=0;m
mid=1; //逆序相减
for(i=0,j=n-1-m+n; i
result-=mid;
}
}
else result=*p;
return(result);
}
//----------------------------------------------------
//功能:求k×k矩阵中元素A(mn)的代数余子式
//入口参数:k×k矩阵首地址;元素A的下标m,n; 矩阵行数 k
//返回值: k×k矩阵中元素A(mn)的代数余子式
//----------------------------------------------------
float Creat_M(float *p, short int m,short int n,short int k)
{
short int len;
short int i,j;
float mid_result=0;
short int quo=1;
float *p_creat,*p_mid;
len=(k-1)*(k-1);
p_creat = (float *)calloc(len, sizeof(float)); //分配内存单元
p_mid=p_creat;
for(i=0;i
if (i!=m && j!=n)
*p_mid++ =* (p+i*k+j);
}
// Print(p_creat,k-1);
quo = (m + n) %2==0 ? 1:-1;
mid_result = (float ) quo * comput_D(p_creat,k-1);
free(p_creat);
return(mid_result);
}
//-------------------------------------------
//功能:打印n×n的矩阵
//入口参数:n×n矩阵的首地址;该矩阵的行数 n
//返回值: 无
//-------------------------------------------
void Print( float *p,short int n)
{
int i,j;
for (i=0;i
for (j=0; j
printf("\n");
}
printf("--------------\n");
}
求旋转矩阵算法
已知逆矩阵,怎么求原矩阵?
怎样求一个矩阵的逆矩阵?
三角矩阵的TC算法
求逆矩阵的公式
谁知道用C程序写个 矩阵连乘的算法,
转置矩阵的C++算法
怎么求矩阵的秩啊
写一个矩阵的类时怎么重载[ ]运算符?
高手进!求生成数独矩阵的算法
含有分块toeplitz矩阵系统的matlab算法
用Matlab求矩阵中的未知数--求位姿矩阵逆解的Matlab方法
如何求解矩阵的逆矩阵
矩阵和逆矩阵的概念
求解逆矩阵的c++程序怎么写的?哪位高手帮下忙!不胜感激!!我的邮箱smartgyg@tom.com
逆矩阵的求法
二阶矩阵的伴随矩阵如何求的
写者问题的算法
求矩阵的特征值和特征向量
求矩阵特征值和特征向量的软件.
求个矩阵的秩!!谢谢
如何求矩阵的本征值啊?
写编程总结的时候的算法报告和功能说明要怎么写
求逆矩阵用什么软件最好??