基于塞尔维斯特准则的二次型标量函数的正定性判断分解.doc

窗体顶端 窗体顶端 窗体底端 Leverrier-Faddeeva递推算法的程序设计范例，编写程序实现基于塞尔维斯特准则的二次型标量函数的正定性判断。要求：编写程序设计报告，含：任务需求分析、程序流程、实现方法、测试验收、源代码等， 根据赛尔维斯特(Sylvester)，可以通过判断对称矩阵P的定号性来确定二次型函数的定号性。 则 矩阵P正定的充要条件是所有的顺子主子式都是正的，即； 矩阵P负定的充要条件是 （三）程序设计流程图： （四） 程序设计实现 输入：矩阵A，n阶方阵 输出：二次型标量函数V(x)是正定的、负定的、既不是正定的也不是负定的 结果 程序验证与测试 包含三个cpp文件 Mat.cpp包含类Mat的各成员函数、 Base.cpp实现矩阵的一些基本操作 包括：求行列式的值、 矩阵的余子式、 判断正定性 Mat.cpp #include Mat.h /********类Mat成员函数************/ void Mat::CreateMat(int r,int c,int type,double *A) { int i; this-rows=r; this-cols=c; this-type=type; data = new double[r*c]; //为矩阵申请（r*c）大小的内存空间。 if(A!=NULL)switch(type) { case INT: for(i=0; ir*c; i++)//循环赋值 data[i] = (int)A[i]; break; case DOUBLE: for(i=0; ir*c; i++) data[i] = (double)A[i]; break; default: data = NULL; break; } else { for(i=0; ir*c; i++)//循环赋值 data[i] = 0; } } void Mat::PrintMat() { int j; if(data==NULL)coutendl还没给该矩阵赋值endl; else { for(int i=0;irows;i++) { switch(type) { case INT: case DOUBLE: for(j=0;jcols;j++)cout”\t”data[i*cols+j]; default:break; } printf(\n); } } } double Mat::Read(int i,int j) { return this-data[i*cols+j]; } （2）Main.cpp #include base.h void main() { Mat A; int n; double* a; cout请输入方阵的维数：; cinn; a=new double[n*n]; cout请输入方阵：endl; for(int i=0;in*n;i++) cina[i]; A.CreateMat(n,n,DOUBLE,a); //创建矩阵 A 并赋值 cout原矩阵：endl; A.PrintMat(); //输出原矩阵 coutendl; Judge(A); //判断是否为实对数矩阵、正定性 // system(pause); } （3）Base.cpp #include base.h #include mat.h //============用递归法求行列式的值（并返回）================= double Fun( int n, double *a ) { if(1==n)return a[0]; double* b; double sum = 0; //i，j为行与列，sum为行列式的值 int x=0,num=0; //用x判断（-1）^（i+j）的正负 b=new double[(n-1)*(n-1)]; //数组用来存放代数余子式 for(int c=0;cn;c++) { for(int i=1;in;i++) //实现将余子式存入数组b中 { for(int j=0;jn;j++) { if(j==c)continue; b[num++]