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

实验3 线性规划的灵敏性分析 成绩 专业班级 信息121 班 学号 201212030120 姓名 刘帅 报告日期 实验类型：●验证性实验 ○综合性实验 ○设计性实验 实验目的：熟练线性规划图解法的灵敏性分析。 实验内容：线性规划的灵敏性分析4个（题目自选b,c灵敏性分析） 实验原理 在线性规划图解法求出最优解的情况下，分析b,c分别变化对最优解的影响，确定最优解的变化范围，在变化的情况下能求出最优解。 实验步骤 1 要求上机实验前先编写出程序代码 2 编辑录入程序 3 调试程序并记录调试过程中出现的问题及修改程序的过程 4 经反复调试后，运行程序并验证程序运行是否正确。 5 记录运行时的输入和输出。 预习编写程序代码： 实验报告：根据实验情况和结果撰写并递交实验报告。 唯一最优解： max z=x1+x2 建立simplex.m文件 function [x,z,flg,sgma]=simplex(A,A1,b,c,m,n,n1,cb,xx) % A,b are the matric in A*x=b % c is the matrix in max z=c*x % A1 is the matric in simplex table % m is the numbers of row in A and n is the column number in A % n1 is the nubers of artificial variables,and artificial variables are default at the last % n1 variables in x. % cb is the worth coefficient matrix for basic variables % xx is the index matrix for basic variables % B1 is the invers matrix for the basic matrix in simplex table.The initial % matrix is default as the last m con in the matrix A. x=zeros(n,1); z=0; B1=A1(:,n-m+1:n) sgma1=c-(cb*B1)*A; [masg,kk]=max(sgma1); k=kk(1); flg=0; ll=0; while (masg0)(ll20) ll=ll+1; thita=1000+zeros(m,1); for i=1:m if A1(i,k)0 thita(i)=A1(i,k)\b(i); end end [r8,c8]=find(thita999); if sum(c8)m [mith,rr]=min(thita); r=rr(1); aa=A1(r,k); for i=1:m if i==r b(r)=b(r)/aa; for j=1:n A1(r,j)=A1(r,j)/aa; end end end for i=1:m if i~=r cc=A1(i,k) b(i)=b(i)-b(r)*cc; for j=1:n A1(i,j)=A1(i,j)-A1(r,j)*cc; end end end cb(r)=c(k); xx(r)=k; B1=A1(:,n-m+1:n); sgma1=c-(cb*B1)*A; [masg,kk]=max(sgma1); k=kk(1); thita=100+zeros(m,1); else flg=3; masg=-1; x=unbound solution; z=inf; end end if flg~=3 if n1==0 sgma1=c-(cb*B1)*A [rc,ccc]=find(sgma1-0.0000000001); if sum(rc)==n-m flg=1; else flg=2; end x=zeros(n,1); for i=1:m x(xx(i))=b(i); end z=c*x; else x=zeros(n,1); for i=1:m x(xx(i))=b(i); end xa=x((n-n1+1):n,:); ra=find(xa); if sum(ra)==0 sgma1=c-(cb*B1)*A; [rc,ccc]=find(sgma1-0.0000000