最优控制大作业.doc

《最优控制大作业.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、单纯行法第一题和第二题采用了单纯形法进行解决，单纯形法的理论依据是：线形规划问题的可行域是n维向量空间Rn中的多面凸集，其最优值如果存在必在该凸集的某顶点处达到。顶点所对应的可行解称为基本可行解。单纯形法的基本思想是：先找出一个基本可行解，对它进行鉴别，看是否是最优解；若不是，则按照一定法则转换到另一改进的基本可行解,再鉴别；若仍不是，则再转换,按此重复进行。因基本可行解的个数有限，故经有限次转换必能得出问题的最优解。如果问题无最优解也可用此法判别。单纯形法的一般解题步骤可归纳如下：1.把线性规划问题的约束方程组表达成典范型方程组，找出基本可行解作为初始基本可行解。2.

2、若基本可行解不存在，即约束条件有矛盾，则问题无解。3.若基本可行解存在，从初始基本可行解作为起点，根据最优性条件和可行性条件，引入非基变量取代某一基变量，找出目标函数值更优的另一基本可行解。4.按步骤3进行迭代,直到对应检验数满足最优性条件（这时目标函数值不能再改善），即得到问题的最优解。5.若迭代过程中发现问题的目标函数值无界，则终止迭代。用单纯形法求解线性规划问题所需的迭代次数主要取决于约束条件的个数。1某工厂生产A和B两种产品。已知制造A产品，每公斤要用煤9吨、电力4千万、劳力3个；制造产品B，每公斤要用煤4吨、电力5千瓦-劳力10个。又知制成产品A每公斤的产值是

3、7万元；B每公斤的产值是12万元。现该厂只有煤360吨、电力200千瓦、劳力300个。问在这种条件下，应该生产A、B产品各多少才能使产值为最高。试写出其数学模型，即约束方程和目标函数，并利用单纯形法求解该线性规划问题。解：设生产A、B产品各吨使引入附加变量，使不等式约束变为等式约束程序清单如下：#include#includeusingnamespacestd;intvarIn(doubledelta[5]);//计算进基变量intvarOut(doublesita[3],doublea[3][5],doubleb[3],doubl

4、edelta[5]);//计算出基变量voidcal(doublesita[3],doublea[3][5],doubleb[3],doubledelta[5]);//计算方程组系数和判别数double*A;//产品A的产量double*B;//产品B的产量intmain(){doublemax=0;//最高产值doubledelta[5]={-7,-12,0,0,0};//判别数doubleb[3]={360,200,300};//基可行解doublea[3][5]={9,4,1,0,0,//系数矩阵A4,5,0,1,0,3,10,0,0,1};doublesita[

5、3];//出基变量判别数intin;//进入基变量位置intout;//出入基变量位置intflag=0;//判别数是否全部非负的标志while(!flag){in=varIn(delta);out=varOut(sita,a,b,delta);cal(sita,a,b,delta);flag=1;for(inti=0;i<5;i++){if(delta[i]<0)//是否所有的判别数都大于0flag=0;}}cout<<"生产A和B产品的产量和最大产值："<

6、(*A)+12*(*B);cout<<"产值最大时A产品的产量:x="<<*A<<","<

7、[i];k=i;}}returnk;}intvarOut(doublesita[3],doublea[3][5],doubleb[3],doubledelta[5]){intl=0;doubleminsita=1000;for(inti=1;i<3;i++){if(a[i][varIn(delta)]>0)sita[i]=b[i]/a[i][varIn(delta)];if(sita[i]