资源描述：
《化实验报告(单纯形法的matlab程序lingo程序)》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、南昌航空大学数学与信息科学学院实验报告实验一：线性规划单纯形算法一、实验目的通过实验熟悉单纯形法的原理，掌握Matlab循环语句的应用,提高编程的能力和技巧。二、实验用仪器设备、器材或软件环境WindowsXp操作系统,Matlab6.5，计算机三、算法对于一般的标准形式线性规划问题(求极小问题),首先给定一个初始基本可行解。设初始基为B,然后执行如下步骤：(1).解,求得，(2).计算单纯形乘子,,得到,对于非基变量，计算判别数,令,R为非基变量集合若判别数,则得到一个最优基本可行解，运算结束；否则，转到下一步(3).解,得到；若
2、,即的每个分量均非正数，则停止计算，问题不存在有限最优解，否则，进行步骤(4).(4).确定下标r,使为离基变量。为进基变量，用替换，得到新的基矩阵，返回步骤（1）。对于极大化问题，可以给出完全类似的步骤，只是确定进基变量的准则不同。对于极大化问题，应令四、计算框图第7页南昌航空大学数学与信息科学学院实验报告开始初始可行解令计算单纯形乘子，计算判别数（非基变量）令是得到最优解解方程，得到。否是不存在有限最优解确定下标，是否为进基变量，用替换，得到新的基矩阵五、计算程序第7页南昌航空大学数学与信息科学学院实验报告function[x,
3、f]=zuiyouhua(A,b,c)size(A)=[m,n];i=n+1:n+m;%基变量集合,后面m个松弛变量为初始基变量;N=1:n;%初始非基变量;B=eye(m,m);xb=b';xn=zeros(m,1);f1=0;w=zeros(1,m);z=-c;%初始判别数;flag=1;while(1)[a,k]=max(z);%x(k)为进基变量;ifa<=0flag=0;breakelsey=inv(B)*A(:,k)ify<=0flag=0;fprintf('不存在最优解')breakend第7页南昌航空大学数学与信息科
4、学学院实验报告t=find(y>0);[a,r1]=min(b1(t)./y(t))r=t(r1);%基变量中第r个变量为退基变量;i(:,r)=kB(:,r)=A(:,k);%换基，即将原基中第ｒ个变量换成第ｋ个变量;cb=c(:,i);%新的价值系数;xb=inv(B)*b;b0=xb;x=zeros(1,n+m)x(:,i)=xb'f=cb*xbz=cb*inv(B)*A-c;%可用z=cb*(BA)-c,判别数.endend六、数值实验及结果分析求解线性规划问题:在工作区输入:A=[3,3,1,0;-4,-4,0,1;2,
5、-1,0,0];b=[30,16,12]';第7页南昌航空大学数学与信息科学学院实验报告c=[-3,1,0,0];[x,f]=zuiyouhua(A,b,c)x=7.33332.666700056.00000f=-19.3333检验结果正确七、心得体会通过这次试验，使我对单纯形法的计算有了更进一步的了解。但是在编程过程中由于对matlab不是很熟悉还是遇到了很多麻烦，所以我觉得老师在让我们编程的时候不能只是简单的介绍一下算法，更要着重说明一下软件的使用方法。这样我们在编程的时候就能更加的得心应手。本次完全仿照老师给的程序，没有能够形
6、成自己的东西。自己编程的能力还是很差的，对于这种已经给出算法的程序也不能正确的编写出来。所以在今后要加强这方面的学习。实验二：Lingo求解动态规划问题第7页南昌航空大学数学与信息科学学院实验报告一、实验目的通过本实验熟悉动态规划的原理，了解动态规划的应用，并能利用数学软件（Lingo）求解动态规划模型。二、问题重述某公司打算向他的营业区增设4个销售点，各区赚取的利润与增设的销售点个数有关，其数据为：销售店增加数第一区利润（万元）第二区利润（万元）第三区利润（万元）第四区利润（万元）0160190200250131022529830
7、8254144539948736005176016554705632721674试求各区应分配几个增设的销售书店，才能使利润最大？其值是多少？三、数学模型设为第区增设销售点的个数，为增设第个点所得到的盈利。故问题模型为：四、计算编程model:sets:quyu/1..4/;zl/0..4/;lirun(quyu,zl):g,c;endsets第7页南昌航空大学数学与信息科学学院实验报告data:g=160310541600705,190225445517632,200298399601721,250308487655674;end
8、datamax=@sum(lirun(i,j):g(i,j)*c(i,j));@for(quyu(i):@sum(lirun(i,j):c(i,j))<=1);@for(lirun:@bin(c));@sum(lirun(i,j):(
