运筹学解对偶单纯形法.docx
- 文档编号:24755032
- 上传时间:2023-06-01
- 格式:DOCX
- 页数:13
- 大小:276.79KB
运筹学解对偶单纯形法.docx
《运筹学解对偶单纯形法.docx》由会员分享,可在线阅读,更多相关《运筹学解对偶单纯形法.docx(13页珍藏版)》请在冰豆网上搜索。
运筹学解对偶单纯形法
题目:
对偶单纯形法解线性规划问题
小组成员:
摘要:
运筹学是辅助人们进行科学管理的一种数学方法.而对偶单纯形法是线性规划中重要的数学方法,在简化运算,解决实际问题中具有重要的应用。
它是解决研究线性约束条件下线性目标函数的极值问题的数学理论和方法,广泛应用于军事作战、经济分析、经营管理和工程技术等方面。
为合理地利用有限的人力、物力、财力等资源作出的最优决策,提供科学的依据。
在经济管理、交通运输、工农业生产等经济活动中,提高经济效果是人们不可缺少的要求。
关键词:
对偶单纯形法线性规划最优解
正文:
单纯形法和对偶单纯形法的基本思想:
给出一个线性规划问题:
Maxz=CX
AX≤b
X≥0
其对偶问题是:
Minw=Yb
YA≥C
Y≥0
单纯形法解决线性规划问题的思想是:
从问题
(1)的一个基解X0开始迭代到另一个基解,在迭代过程中保持基解的可行性,同时它对应的对偶问题
(2)的基解Y0=CBB-1的不可行性逐步消失,直到Y0是问题
(2)的可行解时,X0就是问题
(1)的最优解了。
对偶单纯形法正是基于对称的想法,从一个基解X0开始,X0不是基可行解,但它的检验数全部非正,对应的对偶问题的基解Y0=CBB-1是基可行解;从X0迭代到另一个基解X1,在迭代过程中保持它对应的对偶问题的基解是基可行解,逐步消除原问题基解的不可行性,最终达到两者同时为可行解,也就同时是最优解了。
这就是对偶单纯形法的基本思想。
算法:
用对偶单纯形法解决生产资料分配问题的步骤:
Step1 找出一组以定基元素x0i和人工变量为基变量的正则解X0,若X0是可行的,则X0是最优解,
停止,否则转向STEP2;
Step2 确定换出变量x0l,其中x0l=min{x0r;x0r<0};
Step3 如果对所有非基变量x0j,βlj≥0,则该问题无可行解,运算停止,否则转向STEP4;
Step4 确定换入变量x0k,其中σkβlk=minσtβlt;βlt<0;1≤t≤n+m;
Step5 取x0l为换出变量,x0k为换入变量进行迭代,然后重复上过程直到得到最优解。
程序:
#include
#include
intm,n;
floatM=1000000.0;
floatA[100][100];
floatC[100];
floatb[100];
floatseta[100];
intnum[100];
floatz=0;
voidinput();
voidprint();
intduioudanchunxing1();
intduioudanchunxing2(inta);
voidduioudanchunxing3(inta,intb);
voidinput()
{printf("请输入方程组的系数矩阵维数,m行n列:
\n");
scanf("%d%d",&m,&n);
inti,j;
printf("请输入方程组的系数矩阵A(%d行%d列):
\n",m,n);
for(i=0;i for(j=0;j scanf("%f",&A[i][j]); printf("\n请输入初始基变量的数字代码num矩阵: \n"); for(i=0;i scanf("%d",&num[i]); printf("\n请输入方程组右边的值矩阵b: \n"); for(i=0;i scanf("%f",&b[i]); printf("\n请输入目标函数各个变量的系数所构成的系数阵C: \n"); for(i=0;i scanf("%f",&C[i]); } intduioudanchunxing1() { inti,k; intflag; floatmin=0; for(i=0;i if(b[i]>=0) flag=1; else{flag=0;break;} if(flag==1) return-1; for(i=0;i { if(min>b[i]) {min=b[i];k=i;} } returnk; } intduioudanchunxing2(inta) { inti,j; intflag=0; floatmin; for(j=0;j if(A[a][j]>=0) flag=1; else{flag=0;break;} if(flag==1) {printf("\n该线性规划无最优解! \n");return-1;} for(j=0;j { if(A[a][j]<0) seta[j]=-C[j]/A[a][j]; elseseta[j]=M; } min=M; for(j=0;j { if(min>=seta[j]) {min=seta[j];i=j;} } num[a]=i+1; returni; } voidduioudanchunxing3(intp,intq) { inti,j,c,l; floattemp1,temp2,temp3; c=q; l=p; temp1=A[c][l]; b[c]=b[c]/temp1; for(j=0;j A[c][j]=A[c][j]/temp1; for(i=0;i { if(i! =c) if(A[c][l]! =0) { temp2=A[i][l]; b[i]=b[i]-b[c]*temp2; for(j=0;j A[i][j]=A[i][j]-A[c][j]*temp2; } } temp3=C[l]; for(i=0;i C[i]=C[i]-A[c][i]*temp3; z=z+b[c]*temp3; } voidprint() { inti,j; printf("\n--------------------------------------------------------------------------\n"); printf("\t"); for(i=0;i { printf("%.3f\t",-C[i]); } printf("%.3f",z); printf("\n--------------------------------------------------------------------------\n"); for(i=0;i { printf("x(%d)\t",num[i]); for(j=0;j printf("%.3f\t",A[i][j]); printf("%.3f\n",b[i]); } printf("\n--------------------------------------------------------------------------\n"); } main() { inti,j=0; intp,q; input(); for(i=0;i { if(A[i][num[i]-1]<=0) { b[i]=-b[i]; for(j=0;j A[i][j]=-A[i][j]; } } printf("\n--------------------------------------------------------------------------\n"); printf("\t"); for(i=0;i printf("X(%d)\t",i+1);printf("RHS\n"); while (1) { q=duioudanchunxing1(); if(q==-1) { printf("\n所得解已经是最优解! \n"); print(); for(i=0;i {printf("x(%d)=%.3f\t",num[i],b[i]);} printf("z=%.3f",z); break; } print(); p=duioudanchunxing2(q); if(q==-1)break; duioudanchunxing3(p,q); } } 流程图: duioudanchunxing1(); duioudanchunxing2(inta); duioudanchunxing3(inta,intb); 参考文献: [1]胡运权,甘应爱.运筹学教程[M].北京: 清华大学出版社,2009.[2]王周宏.运筹学基础[M].北京: 清华大学出版社,北京交通大学出版社,2010.[3]何钦铭,颜晖.C语言程序设计[M].浙江: 浙江科学技术出版社,2003.[4]吕凤煮.C++语言程序设计教程[M].北京: 人民邮电出版社.2009.
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 运筹学 对偶 单纯