判定方阵的正定与负定在C语言中写法.docxVIP

判定方阵的正定与负定在C语言中写法#include#includeint input(int n);float count(int n,float a9999);main() int i,j,m,h,k,l,n=0;float x,a9999=0;printf(这是一个判定方阵(100以内包括100)的正定与负定的程序:n); while(1) /*实现多次运算*/n=input(n); /*输入方阵阶数*/ for(i=0;in;i+) /*逐一输入方阵元素*/ for(j=0;jn;j+) printf(请输入：a%d%d=,i+1,j+1); scanf(%f,&x); aij=x; printf(输入的方阵如下：n); for(i=0;in;i+) /*以矩阵形式输出方阵*/ printf( ); for(j=0;jn;j+) printf(%.2f ,aij); /*输出时保留两位小数*/ printf(n); i=0;j=0;h=0;k=0;m=0; /*i、j(h、k)确定偶(奇)数阶的正负*/ for(l=0;l=0) i+; else j+; else if(x=0) h+; else k+;if(x=0) /*m确定各阶顺序主子式为零的次数*/m+; printf(方阵的%d阶顺序主子式值为:%.2fn,l+1,x);if(n%2=0) /*矩阵为偶数阶时*/if(i=n/2&h=n/2&m=0)printf(该方阵为正定矩阵);else if(j=0&k=0&m!=0)printf(该方阵为半正定矩阵);else if(i=n/2&k=n/2&m=0)printf(该方阵为负定矩阵);elseprintf(该方阵为不定矩阵);else /*矩阵为奇数阶时*/if(i=(n-1)/2&h=(n+1)/2&m=0)printf(该方阵为正定矩阵);else if(j=0&k=0&m!=0)printf(该方阵为半正定矩阵);else if(i=(n-1)/2&k=(n+1)/2&m=0)printf(该方阵为负定矩阵);elseprintf(该方阵为不定矩阵);printf(nnn);return 0;int input(int n)int i;printf(请输入方阵阶数：);scanf(%d,&i);if(i100) /*确保方阵阶数输入正确*/printf(输入错误，请重新输入!nn);i=input(i);return i;float count(int n,float a9999) float b9999,sum =0; /*定义数组b并初始化*/ int i=0,j=0; /*i，j为行与列，sum为行列式的值*/ int x=0,c=0,p=0; /*用x判断加与减，c,p为中间变量*/ if(n=1) return a00; for(i=0;in;i+) /*此处大循环实现将余子式存入数组b中*/ for(c=0;cn-1;c+) for(j=0;jn-1;j+) if(ci) /*借助c判断每行的移动方法*/ p=0; /*当p=0时，行列式只向左移，即消去对应的第一列的数*/ else /*否则行列式左移后再上移*/ p=1; bcj=ac+pj+1; if(i%2=0) /*i+j（此时j=0，故只考虑i）为偶数，加法预算*/ x=1; else /*i+j为奇数，减法运算*