實验3线性规划的灵敏性分析.docVIP

实验3 线性规划的灵敏性分析 成绩 专业班级 信息123 学号 201212030317 姓名 李帅 报告日期 实验类型：●验证性实验 ○综合性实验 ○设计性实验 实验目的：熟练线性规划的灵敏性分析。 实验内容：线性规划的灵敏性分析4个（题目自选b,c,A灵敏性分析） 实验原理 在线性规划单纯形法求出最优解的情况下，分析b,c,A分别变化对最优解的影响，确定最优解或最有基的变化范围，在变化的情况下能求出最优解。 实验步骤 1 要求上机实验前先编写出程序代码 2 编辑录入程序 3 调试程序并记录调试过程中出现的问题及修改程序的过程 4 经反复调试后，运行程序并验证程序运行是否正确。 5 记录运行时的输入和输出。 实验报告：根据实验情况和结果撰写并递交实验报告。 1.分析b变化对最优解的影响的代码： 参考程序：function yy=lingmiudu(bb,b,b0,AA,k) m=length(b); if k==0 B=AA(:,1:m); if B*bb=0 flg=the base is optimal else flg=the base is changed end else for i=1:m if AA(i,k)0 d1=-b(i)/AA(i,k); end if AA(i,k)0 d2=-b(i)/AA(i,k); end end ld=max(d1)+b0(k); gd=min(d2)+b0(k); flg=[ld,gd]; end yy=flg 2．分析c变化对最优解的影响的代码： function yy=lingc(cc,c,A,AA,sgma,xx,k，i) m=length(c); B=AA(1,1:m); if k==0 if sgma（i）+cc（i）0 flg= the base is optimal else flg= the base is changed end else for j=1:m if (AA(k,j)0)(sgma(j)~=0) c1=sgma(j)/AA(k,j); end if (AA(k,j)0)(sgma(j)~=0) c2=sgma(j)/AA(k,j); end end cl=max(c1)+c(xx(k)); cg=min(c2)+c(xx(k)); flg=[cl,cg]; end yy=flg end 实验原理：1）资源系数变化是指发生变化，及，相应的，只要0,则最优解不变。经检验计算。 运行结果： 资源系数检验 A=[1 0 0 1 2;0 1 0 4 0;0 0 1 0 4] b=[8,16,12] c=[0 0 0 2 3] m1=0 n=5 m=3 [xx,b,fm,sgma,AA,flg]=myprgmh(m1,m,n,A,b,c) xx = 4 3 5 b = 4 4 2 fm = 14 sgma = -1.5000 -0.1250 0 0 0 AA = 0 0.2500 0 1.0000 0 -2.0000 0.5000 1.0000 0 0 0.5000 -0.1250 0 0 1.0000 flg = There is only one solution bb=[]； k=2； b0=[8 16 12]； b=[4 4 2]； lingmiudu(bb,b,b0,AA,k) yy = 8 32 ans= 8 32 再次输入：  bb=[8 34 12] b0=[8 16 12] k=0 b=[4 4 2] lingmiudu(bb,b,b0,AA,k) flg = the base is changed yy = the base is chang