资源

矩阵平方根示例(omode01)

/***************************************************************//*                                                             *//*S A S S A M P L E L I B R A R Y*//*                                                             *//*姓名:omode01*//*标题:矩阵平方根示例(omode01)*//*产品:或*//*系统:所有*//*关键字:或*//*程序:OPTMODEL*//*数据:*//*                                                             *//*支持:更新:*//*参考编号:*//*MISC:的OPTMODEL章节中的示例1*//*数学编程*//*                                                             *//***************************************************************/过程操作模型;数量n=5;/*矩阵的大小*//*随机原始阵列*/数字A{1..n,1..n}=10-20*rand(“统一”);/*计算的上三角*对称矩阵A*转置(A)*//*除非A是单数,否则应为正def*/数P{i在1..n中,j在i..n}中;对于{i在1..n中,j在i..n}中P[i,j]=和{k在1..n}A[i,k]*A[j,k];/*平方根阵列系数*(对称矩阵的上三角)*/变量q{i在1..n,i..n}中;/*默认初始值q[i,j]=0为*目标的局部最小值,*所以你必须把它移离那个点*/q[1,1]=1;/*最小化q与P的平方差*/最小r=总和{i在1..n中,j在i..n}中(1中的和{k+i+1中的和{k…j}q[i,k]*q[k,j]+j+1中的和{k…n}q[i,k]*q[j,k]-P[i,j])**2;解决;打印q;退出;