當前位置:首頁 » 編程語言 » 求行列式c語言

求行列式c語言

發布時間: 2022-06-21 02:01:02

① 行列式求值--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;
}

② C語言 計算行列式

利用一定的固定方法把他化成上三角或者下三角,然後用一個新數組儲存對角線上,連乘,得結果

③ 如何用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; /*將值返回*/
}

④ 求行列式的值,用C語言怎麼寫啊

行列式的值的源代碼如下:

#include <stdio.h>

#include <stdlib.h>

void main()

{

int i,j,m,n,s,t,k=1;

double a[n][n],f=1,c,x,sn;

for (i=0;i<n;i++)

for (j=0;j<n;j++)

scanf ("%lf",&a[i][j]);

for (i=0,j=0;i<n&&j<n;i++,j++)

{

if (a[i][j]==0)

{

for (m=i;a[m][j]==0;m++);

if (m==n)

{

sn=0;

printf("deta=%lf ",sn);

exit(0);

}

else

for (n=j;n<n;n++)

{

c=a[i][n];

a[i][n]=a[m][n];

a[m][n]=c;

}

k*=(-1);

}

for (s=n-1;s>i;s--)

{

x=a[s][j];

for (t=j;t<n;t++)

a[s][t]-=a[i][t]*(x/a[i][j]);

}

for (i=0;i<n;i++)

f*=a[i][i];

sn=k*f;

printf ("deta=%lf ",sn);

}

(4)求行列式c語言擴展閱讀

1、每個C語言程序中main函數是有且只有一個。

2、C語言程序寫完後,都是先編譯,後連接,最後運行。(.c-->.obj-->.exe)這個過程中注意.c和.obj文件時無法運行的,只有.exe文件才可以運行

⑤ 用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;
}

⑥ 淺析用c語言編程求解線性代數中行列式值

include <iostream>
#include <stdio.h>
int jc(int m){
if(m!=1) return m*jc(m-1);
else return 1;
}
int c(int m,int n){
if(m>=n) return jc(m)/(jc(n)*jc(m-n));
}
int main(void){
int m,n;
scanf("%d%d",&m,&n);
printf("%d\n",c(m,n));
return 0;
}

⑦ 用c語言編程 求n階行列式

#include<stdio.h>
void main()
{
int n,a,b,y[10],x[10],i,j,temp;
float m[10][10],k=1;
printf("Input \"n\"");
scanf("%d",&n);
printf("input the matrix:\n");
for(b=1;b<=n;b++)
{
printf("next row:\n");
for(a=1;a<=n;a++)scanf("%f",&m[a][b]);
}
for(b=1;b<=n;b++)
{
for(a=1;a<=n;a++)printf("%10f",m[a][b]);
printf("\n");
}
printf("\n");
for(i=1;i<=n;i++){x[i]=1;y[i]=i;}
for(j=1;j<=n;j++)while(m[x[j]][j]==0)x[j]++;
for(b=n;b>=1;b--)
for(a=1;a<=b-1;a++)
if(x[a]>x[a+1])
{
temp=x[a];x[a]=x[a+1];x[a+1]=temp;
temp=y[a];y[a]=y[a+1];y[a+1]=temp;
k*=-1;
}
for(b=1;b<=n;b++)
{
for(a=1;a<=n;a++)printf("%10f",m[a][y[b]]);
printf("\n");
}
printf("\n");
for(a=1;a<=n;a++)
for(b=a+1;b<=n;b++)
if(m[a][y[b]]!=0)
for(i=a+1;i<=n;i++)
m[i][y[b]]-=m[a][y[b]]*m[i][y[b-1]]/m[a][y[b-1]];
for(b=1;b<=n;b++)
{
for(a=1;a<=n;a++)printf("%10f",m[a][y[b]]);
printf("\n");
}
for(i=1;i<=n;i++)k*=m[i][y[i]];
printf("the matrix is %f \n",k);
scanf("%d",i);
}
我在剛學c那會寫的一個程序,演算法不是很優秀,而且有漏洞,希望能給你提供大體思路。
我也幫你找了個別人寫的程序,你也去看看吧。也是網路上的。

⑧ 求矩陣行列式的C語言原代碼

完全c語言的
#include <stdio.h>
#include <math.h>

void getarray(int n);
void showarray(int n);
int getresult(int n);

int array[10][10];

void main()
{
int result,size;
printf("\nPlease input the Array size n(n<10):");
scanf("%d",&size);
getarray(size);
showarray(size);
result=getresult(size);
printf("\nResult=%d",result);

printf("\nPress Any Key to Continue...");
getchar();
getchar();
}

void getarray(int n)
{
int row,col;
for(row=0;row<n;row++)
{
printf("\nPlease input line %d:",row);
for(col=0;col<n;col++)
scanf("%d",&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%d",array[row][col]);
printf("\n");
}
}

int getresult(int n)
{
float temp;
int result=1,switchtime=0;
int row,nextrow,col,stemp;
for(row=0;row<n-1;row++)
{
nextrow=row+1;
if(array[row][row]==0)
{
while(array[nextrow][nextrow]==0)
nextrow++;
if((nextrow==n-1)&&(array[nextrow][nextrow]==0))
return 0;
else
{
switchtime++;
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=(float)array[nextrow][row]/array[row][row];
for(col=0;col<n;col++)
array[nextrow][col]+=-(int)temp*array[row][col];
}
}
showarray(n);
for(row=0;row<n;row++)
result*=array[row][row];
if(switchtime%2)
return -result;
else
return result;
}

熱點內容
如何配置2檸檬酸 發布:2025-02-06 00:57:26 瀏覽:434
積木源碼 發布:2025-02-06 00:55:26 瀏覽:545
變分的運演算法則 發布:2025-02-06 00:55:21 瀏覽:775
x2哪個配置性價比高 發布:2025-02-06 00:40:12 瀏覽:109
豬哥亮訪問張菲 發布:2025-02-06 00:37:52 瀏覽:570
期貨賬戶怎麼改密碼 發布:2025-02-06 00:32:35 瀏覽:279
qq自動上傳群文件 發布:2025-02-06 00:26:25 瀏覽:111
安卓照片放在什麼地方 發布:2025-02-06 00:26:24 瀏覽:988
linux系統鏡像iso 發布:2025-02-06 00:15:39 瀏覽:188
存儲上料模塊的意義 發布:2025-02-06 00:14:14 瀏覽:125