线性规划求解对矩阵.docVIP

昆 明 理 工 大 学 理 学 院 信息与计算科学专业 操作性实验报告 年级： 09级 姓名： 陈龙飞 学号： 200911101111 指导教师： 范敏 实验课程名称： 运筹学 开课实验室： 理216 实验成绩： 学风(5) 观察能力(15) 操作能力(30) 调试能力(50) 其它 总分 实验内容： 实验/作业题目:用线性规划的方法求解矩阵对策 3 3 2 9 8 2 4 3 2 6 6 5 6 4 4 2 2．实验/作业课时：2课时 3．实验过程(包括实验环境、实验内容的描述、完成实验要求的知识或技能)： 实验环境：多媒体计算机一台 实验内容： 、对矩阵进行优超处理。 、求解对策矩阵。 程序结构(程序中的函数调用关系图) 优超算法： C++代码： #include iostream using namespace std; int main() { void youchao(); //int youchao_a[4][4]={{3,3,2,9},{8,2,4,3},{2,6,6,5},{6,4,4,2}}; int a[4][4]={{1,0,3,4},{-1,4,0,1},{2,2,2,3},{0,4,1,5}}; coutthe array is:endl; for (int i=0;i4;i++) { for (int j=0;j4;j++) { couta[i][j] ; } if (j==4) cout\n; } // int tt=1; for(int k=0;k4;k++) /* for(int h=0;h4;h++) { if(a[k][h]=a[k][h+tt]) a[k][h]=NULL; }*/ { if (a[0][k]=a[1][k]a[0][k]=a[2][k]a[0][k]=a[3][k]) { a[0][k]=NULL; } else if (a[1][k]=a[0][k]a[1][k]=a[2][k]a[1][k]=a[3][k]) { a[1][k]=NULL; } else if (a[2][k]=a[0][k]a[2][k]=a[1][k]a[2][k]=a[3][k]) { a[2][k]=NULL; } else if (a[3][k]=a[0][k]a[3][k]=a[1][k]a[3][k]=a[2][k]) { a[3][k]=NULL; } } for (int e=0;e4;e++) { if (a[e][0]=a[e][1]a[e][0]=a[e][2]a[e][0]=a[e][3]) { a[e][0]=NULL; } else if (a[e][1]=a[e][0]a[e][1]=a[e][2]a[e][1]=a[e][3]) { a[e][1]=NULL; } else if (a[e][2]=a[e][0]a[e][2]=a[e][1]a[e][2]=a[e][3]) { a[e][2]=NULL; } else if (a[e][3]=a[e][0]a[e][3]=a[e][1]a[e][3]=a[e][2]) { a[e][3]=NULL; } } coutendlendl优超后的矩阵为endl; for(int k2=0;k24;k2++) { for (int k1=0;k14;k1++) { couta[k2][k1] ; } if(k1==4) coutendl; } return 0; } 结果： 当然c++里面的那个矩阵可以换成作业的里面这个，但是作业里面的那个不可以优超，所以效果不明显，我有重新找了一个矩阵，原来都是一样的 Matlab代码： %请输入你要优超的矩阵 L=[3 3 2 9; 8 2 4 3; 2 6 6 5; 6 4 4 2]; %调用main()函数 mian(L) function E0=mian(a) global a; global Y; global k2; global T; global T1;