A. c語言求行列式的值
你這樣全部叫別人寫是不可能的,除非把你的代碼發出來,大家給你檢查下錯誤,不要不好意思.
B. 求c語言求矩陣的行列式和逆天天向哪兒
引用型參數應該在能被定義為const的情況下，盡量定義為const 。
#define N 5 /*[注]：修改6為你所要的矩陣階數*/
#include "stdio.h"
#include "conio.h"
/*js()函數用於計算行列式,通過遞歸演算法實現*/
int js(s,n)
int s[][N],n;
{int z,j,k,r,total=0;
int b[N][N];/*b[N][N]用於存放，在矩陣s[N][N]中元素s[0]的餘子式*/
if(n>2) {for(z=0;z<n;z++)
{for(j=0;j<n-1;j++)
for(k=0;k<n-1;k++)
if(k>=z) b[j][k]=s[j+1][k+1];
else b[j][k]=s[j+1][k];
if(z%2==0) r=s[0][z]*js(b,n-1); /*遞歸調用*/
else r=(-1)*s[0][z]*js(b,n-1);
total=total+r;
}
}
else if(n==2) total=s[0][0]*s[1][1]-s[0][1]*s[1][0];
return total;
}
/*n_1()函數用於求原矩陣各元素對應的餘子式，存放在數組b[N][N]中，定義為float型*/
void n_1(s,b,n)
int s[][N],n;
float b[][N];
{int z,j,k,l,m,g,a[N][N];
for(z=0;z<n;z++)
{l=z;
for(j=0;j<n;j++)
{ m=j;
for (k=0;k<n-1;k++)
for(g=0;g<n-1;g++)
{ if(g>=m&&k<l) a[k][g]=s[k][g+1];
else if(k>=l&&g<m) a[k][g]=s[k+1][g];
else if(k>=l&&g>=m) a[k][g]=s[k+1][g+1];
else a[k][g]=s[k][g];
}
b[z][j]=js(a,n-1);
}
}
}
main()
{int a[N][N];
float b[N][N];
int r,z,j;
float temp;
//clrscr();
printf("Input original data:\n");
for(z=0;z<N;z++) /*輸入所需要的數據，為整型數據*/
for(j=0;j<N;j++)
scanf("%d",&a[z][j]);
printf("\nPress Enter continue......");
getchar();
//gotoxy(1,1);
printf("The original matrix is:\n");
for(z=0;z<N;z++)/*列印原矩陣*/
{for(j=0;j<N;j++)
printf("%5d",a[z][j]);
printf("\n");
}
r=js(a,N); /*調用js()函數計算原矩陣的行列式值*/
printf("\nThe original matrix hanglieshi is:|A|==%d\n",r);
if (r==0) printf("Because |A|==0,the original matrix have no nijuzhen!"); /*判斷條件:若|A|==0，則原矩陣無逆矩陣，反之則存在逆矩陣*/
else
{n_1(a,b,N); /*調用n_1()函數，得到原矩陣各元素對應的"餘子式",存放在數組b[N][N]中*/
for(z=0;z<N;z++) /*求代數餘子式，此時b[N][N]中存放的為原矩陣各元素對應的"代數餘子式"*/
for(j=0;j<N;j++)
if((z+j)%2!=0 && b[z][j]!=0) b[z][j]=-b[z][j];
for(z=0;z<N;z++) /*對b[N][N]轉置,此時b[N][N]中存放的為原矩陣的伴隨矩陣*/
for(j=z+2;j<N;j++)
{temp=b[z][j];
b[z][j]=b[j][z];
b[j][z]=temp;
}
printf("Because |A|!=0,the original matrix have nijuzhen!\n");
printf("The bansuijuzhen A* is:\n");
for(z=0;z<N;z++)/* 列印伴隨矩陣A* */
{for(j=0;j<N;j++)
printf("%4.0f\t",b[z][j]);
printf("\n");
}
for(z=0;z<N;z++) /*求逆矩陣，此時b[N][N]中存放的是原矩陣的逆矩陣*/
for(j=0;j<N;j++)
b[z][j]=b[z][j]/r;
printf("\nThe nijuzhen is:(A*)/|A|(|A|=%d)\n",r); /*列印逆矩陣*/
for(z=0;z<N;z++)
{for(j=0;j<N;j++)
printf("%8.3f",b[z][j]);
printf("\n");
}
}
}
C. 用C語言解決：求任意階（n階）矩陣的行列式
很遺憾，上面匿名的程序不正確。
比如n=2 輸入:
3---7
2---1
得出錯誤結果。
而當輸入n=3
0---1---3
3---0---2
5---2---0
時也會得出錯誤結果。
錯誤的原因有2：
1 是數據類型不對，匿名的程序是設定輸入都是整數int，顯然按照行列式的定義結果肯定是整數，但是他程序中使用了整型數的除法，結果是取整整數，雖然他使用了類型強制轉換，但結果顯然不同，有誤差而且有時候這個誤差很大形成錯誤。
2 是演算法有點問題。小可很欣賞匿名的演算法思路，簡潔明快。不過有相當缺陷，這在程序中註明。
下面的程序是在匿名的程序思路上改寫的。考慮到數據類型和精確度問題，程序中行列式數據使用double型。由於tc和win-tc是16位編輯器，對float型和double型數據支持不好，所以程序是在32位編輯器Dev-c++下調試並通過的。
本題的一個完整的c程序如下,程序在Dev-c++下都調試通過，結果正確。
/* 用C語言解決：求任意階（n階）矩陣的行列式值 */
#include <stdio.h>
#include <math.h>
void getarray(int n);
void showarray(int n);
double getresult(int n);
double array[10][10];/*設矩陣不超過10階，可更改*/
int main()
{
int n;
double result;
printf("\nPlease input the Array size n:");
scanf("%d",&n);
getarray(n);
showarray(n);
result=getresult(n);
printf("\nResult=%f\n",result);
system("pause");
return 0;
}
void getarray(int n)
{
int row,col;
for(row=0;row<n;row++)
{
printf("\nPlease input line %d:",row+1);
for(col=0;col<n;col++)
scanf("%lf",&array[row][col]);
}
}
void showarray(int n)
{
int row,col;
printf("\nA=");
for(row=0;row<n;row++)
{
for(col=0;col<n;col++)
printf("\t%f",array[row][col]);
printf("\n");
}
}
double getresult(int n)
{
double temp,result=1.0;
int switchtime=0,flag=0;
int row,nextrow,col,stemp;
for(row=0;row<n-1;row++)
{
nextrow=row+1;
if(array[row][row]==0)/* 開始處理第一列，如果行列式第一行第一個數為零，要交換行 */
{ while(array[nextrow][row]==0)
{
nextrow++; /* 如果行列式第二行第一個數為零，行增加繼續尋找非零數值的行 */
if(nextrow==n)/* 如果遍歷完行列式行列式第一列元素都為零，退出while循環 */
{ flag=1;
break;
}
}
if(flag==1) /* 退出while循環後回到for(row=0;row<n-1;row++)行加1?*/
continue; /* 從array[row][row]==0知列也相應加1，開始處理第二列 */
switchtime++; /* 每交換一次行，行列式符號變化1次，統計變化次數 */
for(col=0;col<n;col++) /* 交換非零行到行列式頂部 */
{
stemp=array[row][col];
array[row][col]=array[nextrow][col];
array[nextrow][col]=stemp;
}
}
for(nextrow=row+1;nextrow<n;nextrow++)
{ /* 類似高斯消去法，消第一行下各行第一列數值到零*/
temp=array[nextrow][row]/array[row][row];
for(col=0;col<n;col++)
array[nextrow][col]+=-temp*array[row][col];/* 化行列式為上三角行列式形式 */
}
}
showarray(n);
for(row=0;row<n;row++)
result*=array[row][row];
if(switchtime%2)
return -result;
else
return result;
}
D. 求矩陣行列式的C語言原代碼
要看是不是用戶輸入行例數，如果是這樣就有些麻煩。。。別問人了，也不要想著別人給你寫，直接去網上找，找得到這種題的代碼的，因為題不很難。。。HOHOHO
就用「求矩陣行列式的C語言原代碼」在網路里找就是了。。。
E. 如何用c語言來求一個四階行列式的值
#include <stdio.h>
#include <stdlib.h>
#define MAX 9 /*定義最大為9階行列式*/
int Fun(int n, int a[MAX][MAX] ); /*函數聲明*/
int main()
{
int n = 0; /*初始化階數n*/
int i = 0, j = 0; /*i,j分別表示行與列*/
int a[MAX][MAX] = {{0}}; /*定義行列式*/
scanf("%d",&n); /*讀入階數*/
while( n != 0) /*輸入為0時退出程序*/
{
for( i = 0; i < n; i++) /*此處2次循環將行列式存入數組中*/
{
for( j = 0; j < n; j++)
scanf("%d", &a[i][j]);
}
printf("%d\n", Fun( n, a ) );
scanf("%d", &n);
}
return 0;
}
/*以下為計算行列式值的遞歸函數*/
int Fun( int n, int a[MAX][MAX] )
{
int b[MAX][MAX] = {{0}}; /*定義數組b並初始化*/
int i = 0, j = 0, sum = 0; /*i，j為行與列，sum為行列式的值*/
int x = 0,c = 0,p=0; /*用x判斷加與減，c,p為中間變數*/
if(n == 1)
return a[0][0];
for(i = 0;i < n; i++) /*此處大循環實現將餘子式存入數組b中*/
{
for(c = 0;c < n-1; c++)
{
for(j = 0;j < n-1;j++)
{
if (c < i){ /*藉助c判斷每行的移動方法*/
p = 0; /*當p=0時，行列式只向左移，即消去對應的第一列的數*/
}
else{ /*否則行列式左移後再上移*/
p = 1;
}
b[c][j] = a[c+p][j+1];
}
}
if(i % 2 == 0){ /*i+j（此時j=0，故只考慮i）為偶數，加法預算*/
x = 1;
}
else{ /*i+j為奇數，減法運算*/
x = (-1);
}
sum += a[i][0] * Fun(n - 1, b ) * x; /*計算行列式的值*/
}
return sum; /*將值返回*/
}
F. c語言計算矩陣行列式問題
簡化代碼，宏定義後，後面的程序直接拿來用就可以了
G. C語言，這是我編的求矩陣行列式的函數，但是出了問題，求解答
你這個程序至少有兩個問題（可能還有別的錯誤，我再看看）
第一個
if(m=1)//應該是m==1
{
return**p;
}//行數為1直接返回
第二個問題，循環嵌套不對，下面是正確的代碼，注意大括弧位置。
for(i=0;i<m;i++)
{
p2[i]=(float**)malloc((m-1)*sizeof(float));
for(j=0;j<m-1;j++)
{
p2[i][j]=(float*)malloc((m-1)*sizeof(float));
}//第一維是p第幾列的餘子式，第二維第三維是餘子式
}
H. C語言求行列式的值
C語言程序如下：
#include<stdio.h>
#include<math.h>
#defineMAX100
doubleD(doublearr[][MAX],intn);
voidDisplay(doublearr[][MAX],intn);
voidmain()
{
doublearr[MAX][MAX];
doublebak[MAX][MAX];
doublevalue;
intn;
inti,j,k,c;
printf("Pleaseentermatrixsizen(1<=n<20):");
scanf("%d",&n);
printf("Pleaseinputmatrixlinebyline: ");
for(i=0;i<n;i++)
{
for(j=0;j<n;j++)
{
scanf("%lf",&arr[i][j]);
}
}
printf("matrixa: ");
Display(arr,n);
/*保存D(0,i)的值到bak中*/
for(i=0;i<n;i++) /*處理第1行的每一列元素*/
{
for(j=1;j<n;j++) /*逐行處理*/
{
for(c=0,k=0;k<n;k++) /*逐列處理*/
{
if(k==i)
{
continue;
}
bak[j-1][c++]=arr[j][k];
}
}
value=D(bak,n-1);
printf("Submatrix: ");
Display(bak,n-1);
printf("DValueoftheSubmatrixis%f ",D(bak,n-1));
}
value=D(arr,n);
printf("result=%f ",value);
}
doubleD(doublearr[][MAX],intn)
{
doublebak[MAX][MAX];
inti,j,k,c,d;
doublesum=0;
if(n==1)
{
returnarr[0][0];
}
for(i=0;i<n;i++) /*處理第1行的每一列元素*/
{
for(j=1;j<n;j++) /*逐行處理*/
{
for(c=0,k=0;k<n;k++) /*逐列處理*/
{
if(k==i)
{
continue;
}
bak[j-1][c++]=arr[j][k];
}
}
/*計算arr(row,col)*D(row,col)的值*/
sum+=(i%2==0?1:-1)*arr[0][i]*D(bak,n-1);
}
returnsum;
}
voidDisplay(doublearr[][MAX],intn)
{
inti,j;
printf("");
for(i=0;i<n;i++)
{
for(j=0;j<n;j++)
{
printf("%10.1lf",arr[i][j]);
}
printf(" ");
}
}
運行測試：
I. 行列式求值--C語言
#include<stdio.h>
#include<stdlib.h>
void main()
{
float f(float *a,int n);
int i,j,n;
scanf("%d",&n);
float *a=(float *)malloc(n*sizeof(float));
for(i=0;i<n;i++)
for(j=0;j<n;j++)
scanf("%f",a+i*n+j);
printf("行列式如下:\n");
for(i=0;i<n;i++)
{
for(j=0;j<n;j++)
printf("%-5.2f ",*(a+i*n+j));
printf("\n");
}
printf("%-8.2f\n",f(a,n));
}
float f(float *a,int n)
{
int i,j,k,p=1;
float *b=(float *)malloc(n*n*sizeof(float));
for(i=0;i<n;i++)
for(j=0;j<n;j++)
*(b+i*n+j)=0;
float s=0;
if(n==1)
return *a;
if(n==2)
return *a**(a+n+1)-*(a+1)**(a+n);
else
for(i=0;i<n;i++)
{
for(j=0;j<n-1;j++)
for(k=0;k<n-1;k++)
if(k>=i)
*(b+n*j+k)=*(a+n*j+k+1);
else
*(b+n*j+k)=*(a+n*j+k);
s+=p**(a+n*(n-1)+i)*f(b,n-1);
p*=-1;
}
return s;
}
J. 求一個 計算矩陣行列式 的C或C++程序！
以下是我原來寫過的一個關於矩陣的類,裡面實現了求方陣的行列式函數:intdet(matrix);
//matrix.h
//*******************************
#ifndefmatrix_h
#definematrix_h
classmatrix
{
private:
intm;
intn;
doublearr[8][8];
public:
matrix();
matrix(double);
voidset();//對矩陣的賦值;
voidshow();//在屏幕上顯示矩陣;
friendinthomotype(matrix,matrix);//是否同型;
friendintmultipliable(matrix,matrix);
//判斷矩陣是否可相乘
friendmatrixoperator+(matrix,matrix);
//矩陣相加,對加號重載;
friendmatrixoperator-(matrix,matrix);
//矩陣機減,對減號重載;
friendmatrixadd(matrix,matrix);
//矩陣相加;
friendmatrixminus(matrix,matrix);
//矩陣相減
friendmatrixmultiply(matrix,matrix);
//矩陣相乘
friendmatrixoperator*(matrix,matrix);
//矩陣相乘,對乘號重載;
frienddoubledet(matrix);
//求方陣的行列式;
intissquare()
{
returnm==n;
}
//判斷是否為方陣;
matrixleft(intx,inty);
//求元素arr[x][y](x,y從0開始)的餘子式;
};
#endif
//*******************************
//matrix.cpp
//下面是對類中成員及友元函數的實現
//*******************************
#include<iostream.h>
#include"matrix.h"
matrix::matrix()
{
m=8;
n=8;
for(inti=0;i<m;i++)
{
for(intj=0;j<n;j++)
arr[i][j]=0;
}
}
matrix::matrix(doublex)
{
m=1;
n=1;
arr[0][0]=x;
}
voidmatrix::set()
{
cout<<"Setthetypeofthematrix(?*?)"<<endl;
cin>>m>>n;
cout<<":"<<endl;
for(inti=0;i<m;i++)
{
for(intj=0;j<n;j++)
cin>>arr[i][j];
}
}
voidmatrix::show()
{
cout<<"Thisisthematrix:"<<endl;
for(inti=0;i<m;i++)
{
for(intj=0;j<n;j++)
{
cout<<arr[i][j]<<ends;
}
cout<<endl;
}
}
inthomotype(matrixx,matrixy)
{
if(x.m==y.m&&x.n==y.n)
return1;
else
return0;
}
intmultipliable(matrixx,matrixy)
{
if(x.n==y.m)
return1;
else
return0;
}
matrixadd(matrixx,matrixy)
{
matrixz;
if(homotype(x,y))
{
z.m=x.m;
z.n=x.n;
for(inti=0;i<z.m;i++)
for(intj=0;j<z.n;j++)
z.arr[i][j]=x.arr[i][j]+y.arr[i][j];
returnz;
}
else
{
cout<<"cannotbeadded!"<<endl;
returnz;
}
}
matrixoperator+(matrixx,matrixy)
{
returnadd(x,y);
}
matrixminus(matrixx,matrixy)
{
matrixz;
if(homotype(x,y))
{
z.m=x.m;
z.n=x.n;
for(inti=0;i<z.m;i++)
for(intj=0;j<z.n;j++)
z.arr[i][j]=x.arr[i][j]-y.arr[i][j];
returnz;
}
else
{
cout<<"cannotbeadded!"<<endl;
returnz;
}
}
matrixoperator-(matrixx,matrixy)
{
returnminus(x,y);
}
matrixmultiply(matrixx,matrixy)
{
matrixz;
if(x.n!=y.m)
{
cout<<"."<<endl;
returnz;
}
for(inti=0;i<x.m;i++)
{
for(intj=0;j<y.n;j++)
{
for(intk=0;k<x.n;k++)
z.arr[i][j]+=x.arr[i][k]*y.arr[k][j];
}
}
z.m=x.m;
z.n=y.n;
returnz;
}
matrixoperator*(matrixx,matrixy)
{
matrixz;
z=multiply(x,y);
returnz;
}
matrixmatrix::left(intx,inty)
{
matrixleftmatrix;
if((x>=m)||(y>=n))
{
cout<<"errer"<<endl;
returnleftmatrix;
}
/*if(!issquare())
{
cout<<"不是方陣!"<<endl;
return*this;
}//*/
leftmatrix.m=m-1;
leftmatrix.n=n-1;
inttestx=0;
inttesty=0;
for(inti=0;i<leftmatrix.m;i++)
{
testy=0;
for(intj=0;j<leftmatrix.n;j++)
{
if(i==x)
testx=1;
if(j==y)
testy=1;
if((!testx)&&(!testy))
leftmatrix.arr[i][j]=this->arr[i][j];
elseif(testx&&(!testy))
leftmatrix.arr[i][j]=this->arr[i+1][j];
elseif((!testx)&&testy)
leftmatrix.arr[i][j]=this->arr[i][j+1];
else
leftmatrix.arr[i][j]=this->arr[i+1][j+1];
}
}
returnleftmatrix;
}
doubledet(matrixx)//遞歸演算法
{
if(!x.issquare())
{
cout<<"不是方陣!"<<endl;
return9.99;
}
if(x.m==0)
return0;
elseif(x.m==1)
returnx.arr[0][0];
elseif(x.m==2)
return(x.arr[0][0]*x.arr[1][1]-x.arr[0][1]*x.arr[1][0]);
else
{
doublenum=0;
inta=1;
for(inti=0;i<x.m;i++)
{
num=num+a*x.arr[0][i]*det(x.left(0,i));
//按第0行展開
a=-a;
}
returnnum;
}
}
//**************************************