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

实验3 线性规划的灵敏性分析 成绩 专业班级 数学122 学号 201212010218 姓名 王怡 报告日期 . 实验类型：●验证性实验 ○综合性实验 ○设计性实验 实验目的：熟练线性规划的灵敏性分析。 实验内容：线性规划的灵敏性分析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; end end 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)=m1 vv=vv+2; end end if vv=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 f