2011-2012-2实验2线性规划的灵敏性分析.docVIP

实验2 线性规划的灵敏性分析成绩专业班级 数学101班 学号 201012010103 姓名李旭红 报告日期 2012-4-25 实验类型：验证性实验 综合性实验 设计性实验实验目的：熟练线性规划的灵敏性分析。实验内容：线性规划的灵敏性分析4个（题目自选b,c,A灵敏性分析）实验原理 在线性规划单纯形法求出最优解的情况下，分析b,c,A分别变化对最优解的影响，确定最优解或最有基德变化范围，在变化的情况下能求出最优解。实验步骤1 要求上机实验前先编写出程序代码 2 编辑录入程序3 调试程序并记录调试过程中出现的问题及修改程序的过程4 经反复调试后，运行程序并验证程序运行是否正确。5 记录运行时的输入和输出。 预习编写程序代码：一、单纯性法的程序：function xx,b,fm,sgma,AA,flg=myprgmh(m1,m,n,A,b,c)B0=A(:,1:m); % B0 初始可行基矩阵（单位矩阵）;cb=c(:,1:m);xx=1:m; % xx 变量的下标;sgma=c-(cb*B0)*A; % sgma 检验数; h=-1;sta=ones(m,1);for i=1:n if sgma(i)0 h=1; endend vv=0; while h0 msg,mk=max(sgma); for i=1:m if A(i,mk)0 sta(i)=b(i)/A(i,mk); else sta(i)=10000; end end mst,mr=min(sta); if mst=10000 flg=unbounded solution; fm=inf; xx=; b=; h=-1; vv=1; AA=; else zy=A(mr,mk) for i=1:m if i=mr for j=1:n A(i,j)=A(i,j)/zy; end b(i)=b(i)/zy; else end end for i=1:m if i=mr amk=A(i,mk); b(i)=b(i)-amk*b(mr); for j=1:n A(i,j)=A(i,j)-amk*A(mr,j); end else end A; B1=A(:,1:m); % B1 新基的逆矩阵; cb(mr)=c(mk); xx(mr)=mk; sgma=c-cb*A; h=-1; for i=1:n if sgma(i)0 h=1; end end end cb b fm=sum(cb*b); if (h=-1)&(vv=1) vv=0; for i=1:m if xx(i)=2; flg=nofeasibel; xx=; fm=; b=; vv=1; AA=; end if vv=1 AA=A; ss=size(find(sgma) ww=ss(2) if ww=n-m flg=There is only one solution; else flg=There are many solutions; end end end end end end二、分析b对最优解的影响的程序代码：1、function Y=bfanwei(b,b0,AA)s,t=size(b);A=zeros(s,1);Y=zeros(s,2);B=AA(:,1:s);Z=zeros(s,2);for i=1:sA=B(:,i);for j=1:s if (A(j)0) Z(j,1)=-b(j)/A(j); Z(j,2)=10000; elseif (A(j)=0 flg=The base is optimal still else flg=The base is not optimal endelse ld=max()+b0(k) gd=min()+b0(k) flg=ld,gdendyy=flg实验原理：资源数量变化是指系数发生变化，即。并假设规划问题的其他系数不变，实验报告：根据实验情况和结果撰写并递交实验报告。（一） A=0 0 1 1 -2 1 0;1 0 0 -4 1 2 -1;0 1 0 -2 0 1 0; b=11;3;1; c=-100 -100 0 3 -1 -1 0; xx,b,fm,sgma,AA,flg=myprgmh(2,3,7,A,b,c)zy = 1cb = -100 -100 -1b = 10 1 1zy = 3cb = 3 -100 -1b = 3.3333 1.0000 7.6667zy = 1cb = 3 -1 -1b = 4 1 9ss = 1 4ww = 4xx = 4 5 6b = 4 1 9fm = 2sgma = -99.6667 -99.3333 -0.3333 0 0 0 -0.3333AA = 0.6667 -1.6667 0.3333 1.0000 0 0 -0.6667 1.0000 -2.0000 0 0 1.0000 0 -1.0000 1.3333 -2.3333 0.6667 0 0 1.0000 -1.3333flg =There is only one solution b0=11;3;1; b=4;1;9; Y=bfanwei(b,b0,AA)Y = 1.0e+004 * 0.0010 1.0011 -0.9997 0.0004 -0.0011 1.0001（二） A=0 0 1 0 5 10;1 0 0 -1 1 1;0 1 0 0 0 1; b=50;1;4; c=-1000 -1000 0 0 1 3; xx,b,fm,sgma,AA,flg=myprgmh(2,3,6,A,b,c)xx = 5 6 4b = 2 4 5fm = 14sgma = 1.0e+003 * -1.0000 -1.0010 -0.0002 0 0 0AA = 0 -2.0000 0.2000 0 1.0000 0 0 1.0000 0 0 0 1.0000 -1.0000 -1.0000 0.2000 1.0000 0 0flg =There is only one solution b0=50;1;4; b=2;4;5; Y=bfanwei(b,b0,AA)Y = -9950 55 -3 2 -6 10004 detlb=1;1;1; b0=50;1;4; b=2;4;5; k=0; m=3; yy=lingmindf(detlb,b0,b,AA,k,m)B = 0 -2.0000 0.2000 0 1.0000 0 -1.0000 -1.0000 0.2000nb = 51 2 5flg =The base is not optimalyy =The base is not optimalyy =T