许久没来过了,求论坛的大神帮忙解决下矩阵的问题?希望杨大哥能看到。
// 问题是这样的:// 设 int m,n; (3 <= m < n <= 16 ),
// int x[m][n]; 数组的每个数为区间[-17,17]的随机整数,srand 、rand配合产生,
// 用数组 x[m][n] 构造一个m行n列矩阵,依次排列(先排列行),
// 请问?
// 怎样才能快速构造出一个线性无关的m行n列矩阵,
// 如已构造的矩阵线性相关则重新产生随机数数组。
#include <stdio.h> #define M 2 //矩阵的行 #define N 3 //矩阵的列 int Fun(int n, int a[M][M] ); /*函数声明*/ int main() { int x[M][N];//此处假设m=2,n=3 int a[M][M];//取自x[2][3]的方阵 int k=0;//统计行列式不为0的阶数 int i,j; while(k!=M)//如果最大不为0的行列式不为M,则不是线性无关 { for (i=0;i<=M-1;i++) { for(j=0;j<=N-1;j++) { x[i][j]=rand()%35;//[-17,17]之间共有35个数字,rand()%35产生0-34的随机数字 x[i][j]=x[i][j]-17; } } for (i=0;i<=M-1;i++) { for(j=0;j<=M-1;j++) { a[i][j]=x[i][j]; } } for (i=0;i<=M-1;i++) { if(Fun(i,a)!=0) k +=1; } } for (i=0;i<=M-1;i++) { for(j=0;j<=N-1;j++) { printf("%d,",x[i][j]); } printf("\n"); } return 0; } int Fun( int n, int a[M][M] ) { int b[M][M] = {{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; /*将值返回*/ }