灵敏度分析使用MATLAB编写

1、实验二、线性规划的灵敏度分析(一) 实验目的1. 线性规划求解的单纯形法的灵敏度分析的编程实现2掌握使用 matlab 、Lingo 、 Excel 的规划求解功能求解，并利用“敏感性报告”进行 分析。(二)实验内容课本例 1 解的灵敏度分析 (1)：调用单纯形程序： function x,z,flg,sgma=simplexfun(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 numb

2、ers of row in A and n is the con number in A% n1 is the nubers of artificial variables,and artificial variables are default as 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 matr

3、ix 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)&(ll0 thita(i)=A1(i,k)b(i);endend r8,c8=find(thita999);if sum(c8)mmith,rr=min(thita);r=rr(1);aa=A1(r,k);for i=1

4、:mif i=rb(r)=b(r)/aa;for j=1:nA1(r,j)=A1(r,j)/aa ; endendendfor i=1:mif i=rcc=A1(i,k) b(i)=b(i)-b(r)*cc;for j=1:nA1(i,j)=A1(i,j)-A1(r,j)*cc; end endendcb(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);elseflg=3;masg=-1;x=unbound solution;z=inf;ende

5、nd if flg=3if n1=0sgma1=c-(cb*B1)*A rc,ccc=find(sgma1-0.0000000001);if sum(rc)=n-mflg=1;elseflg=2;endx=zeros(n,1);for i=1:mx(xx(i)=b(i);endelsex=zeros (n ,1); for i=1:mx(xx(i)=b(i);endxa=x( (n_n 1+1): n,:); ra=fi nd(xa);if sum(ra)=0sgma仁c-(cb*B1)*A;rc,ccc=fi nd(sgma10 aa(i)=-b(i)/B1(i,r);endif B1(i,r)0 bb(i)=-b