单纯形法matlab程序

function Z=dcxf(c,A,N) %定义函数名称为dcxf。

l=length(N);

CB=c(N(1):N(l))

[m,n]=size(A);

b=A(:,n);

A=A(:,1:n-1);

%参数包括目标函数系数(C)，约束条件的系数矩阵(A)，

%其中A的最后一列为约束条件的右端值b，初始基向量的位置(N)。

sigma=c-CB*A;%计算检验数sigma。

display('初始单纯形表为:');%输出初始的单纯形表table=[nan,nan,nan,c;CB',N',b,A;nan,nan,nan,si gma]

opt=1;step=0;

while opt

step=step+1;%定义循环，直到第“step”步找到最优解(opt=0)。

if sum(sigma>0)==0 %利用检验数判断是否得到最优解，并给出提示。

display('没有得到最优解,继续迭代．'); opt=0;

else

inb=find(sigma==max(sigma)); %利用单纯形方法找到入基变量的位置

num=length(inb);

Inb=inb(num)

flag=0;

for i=1:m %利用单纯形方法找出出基变量的位置

if A(i,inb)>0

theta(i)=b(i)/A(i,inb);

else

theta(i)=inf;

end

end

outb=find(theta==min(theta));

num=length(outb); %判断足否出现退化现象，如出现退化，给il{语言提示，并取最后出现的符合出基条件的变量为出基变量。

if num~=1

display('出现退化情况．');

end

outb=outb(num);

for i=1:m %将单纯形表进行“转轴”运算，得到新的单纯形表。

for j=1:n-1

if i==outb

Anew(i,j)=A(outb,j)/A(outb,inb);

bnew(i)=b(outb)/A(outb,inb);

else

Anew(i,j)=A(i,j)-A(outb,j)/A(outb,inb)*A(i,inb );

bnew(i)=b(i)-b(outb)/A(outb,inb)*A(i,inb);

end

end

end

display('主元素为:'),a=[A(outb,inb),outb,inb] %输出主元素，计算新单纯形表的检验数。

A=Anew;b=bnew;

N(outb)=inb;

for i=1:l

CB(i)=c(N(i));

end

sigma=c-CB*A;

end

display('迭代得到的单纯形表为:'); %输出得到的新单纯形表。并给出提示语句。

table=[nan,nan,nan,c;CB',N',b',A;nan,nan,nan,s igma]

end

display('得到最优解:');

Z=CB*b'; %计算检验数