最优化方法课程论文

前言PAGE10PAGE11四川理工学院《最优化方法》课程论文题目：线性规划的单纯形算法姓名：专业：统计学班级：2011级1班学号：完成日期：2014年6月27日四川理工学院理学院二O一四年六月摘要线性规划是运筹学中研究较早、发展较快、应用广泛、方法较成熟的一个重要分支，它是辅助人们进行科学管理的一种数学方法。是研究线性约束条件下线性目标函数的极值问题的数学理论和方法。为了得到线性目标函数的极值，我们有多重方法。本文采用单纯性算法求解线性规划问题，并通过Matlab软件编写程序进行求解。关键词：线性规划单纯性算法Matlab编程四川理工学院《最优化方法》课程论文四川理工学院毕业论文目录TOC\o"1-5"\h\z\u24242一、单纯性方法简介 1173431.1单纯性方法提出 1311011.2单纯性方法的基本思想和步骤 1161271.2.1基本思想 1148381.2.2计算步骤 110125二、问题的提出与分析 1165122.1问题提出 163862.2问题分析 23182三、程序设计 2182923.1算法设计 2292213.2算法框图 3189993.3程序编制 414149四、结果分析 6223314.1设计结果 6115584.2进一步讨论和验证 82159五、结束语 892575.1设计的优缺点 8148955.2收获与总结 99903参考文献 1022250附录 11四川理工学院《最优化方法》课程论文四川理工学院《数学建模课程设计》论文单纯性方法简介1.1单纯性方法提出单纯形法，求解线性规划问题的通用方法。单纯形是美国数学家G.B.丹齐克于1947年首先提出来的,这是20世纪数学界最重大的成果之一。由于这一方法的有效性，几十年来一直在几乎所有的领域得到广泛应用。它的理论根据是：线性规划问题的可行域是n维向量空间Rn中的多面凸集，其最优值如果存在必在该凸集的某顶点处达到。顶点所对应的可行解称为基本可行解。1.2单纯性方法的基本思想和步骤1.2.1基本思想单纯形法的基本思想是：先找出一个基本可行解，对它进行鉴别，看是否是最优解；若不是，则按照一定法则转换到另一改进的基本可行解，再鉴别；若仍不是，则再转换，按此重复进行。因基本可行解的个数有限，故经有限次转换必能得出问题的最优解。如果问题无最优解也可用此法判别。1.2.2计算步骤1、对于一般的的线性规划，将其化为标准型；2、求出初始基本可行解；3、先检验其最优性；4、如果不是最优的，则从取负值的非基变量中选取一个最负确定为入基变量；5、选好入基变量后，再在基变量中选取一个出基变量；6、选好入基变量和出基变量后，进行高斯消去，得到新的可行解；7、重复以上过程，直至找到最优解。、二、问题的提出与分析2.1问题提出本文运用单纯性算法求解下列问题：Maxs.t并编写MATLAB程序求解。2.2问题分析在用单纯性算法解决现行规划问题时，我们通常考察标准形现行规划问题，其标准形如下：现在将本文所讨论的线性规划化为标准线性规划的形式：MinS.t.其中A=[230100024010325001]，三、程序设计3.1算法设计解，求得，令，计算目标函数值,以记的第i个分量；计算单纯性乘子w，,得到,对于非基变量，计算判别系数，令，R为非基变量集合，若判别系数，则得到一个最基本可行解，运算结束；否则，转到下一步解，得到；若，即的每一个分量均非正数，则停止计算，问题不存在有限最优解，否则，进行步骤4;确定下标r，使，为出基变量，为入基变量，用替换，得到新的基矩阵B，返回步骤1。3.2算法框图开始开始初始可行解初始可行解令令计算单纯形乘子计算单纯形乘子，计算判别数（非基变量）令 是得到最优解得到最优解解方程，得到解方程，得到。 是不存在有限最优解不存在有限最优解确定下标，是确定下标，是为进基变量，用为进基变量，用替换，得到新的基矩阵3.3程序编制A=input('A=');b=input('b=');c=input('c=');formatrat[m,n]=size(A);E=1:m;E=E';F=n-m+1:n;F=F';D=[E,F];X=zeros(1,n);if(n<m)fprintf('不符合要求需引入松弛变量')flag=0;elseflag=1;B=A(:,n-m+1:n);cB=c(n-m+1:n);whileflagw=cB/B;panbieshu=w*A-c[z,k]=max(panbieshu);fprintf('b''./(B\\A(:,%d))为',k);b'./(B\A(:,k))if(z<0.000000001)flag=0;fprintf('已找到最优解!\n');xB=(B\b')';f=cB*xB';fori=1:nmark=0;forj=1:mif(D(j,2)==i)mark=1;X(i)=xB(D(j,1));endendifmark==0X(i)=0;endendfprintf('基向量为:');Xfprintf('目标函数值为:');felseif(B\A(:,k)<=0)flag=0;fprintf('\n此问题不存在最优解!\n');elseb1=B\b';temp=inf;fori=1:mif((A(i,k)>0)&&(b1(i)/(A(i,k)+eps))<temp)temp=b1(i)/A(i,k);r=i;endendfprintf('x(%d)进基，x(%d)退基\n',k,D(r,2));B(:,r)=A(:,k);cB(r)=c(k);D(r,2)=k;endendendend四、结果分析4.1设计结果在命令窗口中输入：A=[2,3,0,1,0,0;0,2,4,0,1,0;3,2,5,0,0,1]b=[1200,800,2000]得到如下结果:我们可以看到，程序经过4次换基迭代，得到目标函数的最优值为-2600，即目标函数的最小值为-2600。从而，原问题的最大值为2600。4.2进一步讨论和验证对于MATLAB程序的正确性与软件运行的可行性。由于计算量并不是很大，我们通过单纯性表进行手工计算。经过几次换基迭代，我们选取的入基变量和出基变量与以上软件运行过程得到的结果完全相同。由此，我们可以认定目标函数的最小值为-2600，即原问题的最大值为2600。五、结束语5.1设计的优缺点设计优点：设计的程序是根据课本的步骤编写的；程序的编制能得到正确结果；编制的程序得到的结果中具体体现每一步的出基变量与入基变量，清晰明了；设计缺点：不能直观的反应迭代步数，如若迭代次数过多，则想要了解迭代步数则比较麻烦；不能给出完整的单纯性表。5.2收获与总结通过本次课程论文设计，让我对单纯性法有了进一步的了解，明确了它的具体思想理论，算法步骤。此外，通过此次课程设计，初次接触了MATLAB软件，让我对MATLAB软件有了初步的了解，此次论文的完成，主要是通过根据算法设计，编制MATLAB程序，通过MATLAB软件对模型求解。因此，此次设计的最大问题在于怎样设计算法程序，但这对于我们来说难度还是比较大，所以，此次的单纯性算法程序直接利用网上给出的算法程序进行设计。但网上的很多程序也存在很多问题，需要在一次一次的错误中不断的更正问题，直到最后得到模型正确的结果。由于对MATLAB软件的不了解，对于程序设计的优缺点不是很明白。而对于以后，还是希望能多学习一下软件的知识，能够深入了解一下软件的程序设计以及问题分析。参考文献参考文献[1][精通MATLAB最优化计算（第二版）]龚纯等[2]/view/83074402a6c30c2259019e1b.html