对偶单纯形法.docVIP

2023-2023〔1〕专业课程实践论文

对偶单纯形法

范俊，0818180124，R数学08-1班

陈泊宇，0818180116，R数学08-1班

侯建成，0818180128，R数学08-1班

算法理论

对偶单纯形法是解决线性规划问题的一种方法。

当一个线性规划问题具有中心部位具有单位子块，底行相应于单位子块位置元素为0，底行其他元素非负等三个特点而同时又不具有右行列元素非负特点，即满足特点1、3、4时可使用此方法求出最优解。

此方法运用所允许的运算，始终不破坏第1、3、4三个特点，而逐步调出第2个特点的方法，以解出最优解。

其具体的步骤为：

1.从右列负元素中任选一个。

2.从所选行位于中心部位的负元素中确定一个，由于要保持第3、4的特点不被破坏，理应取相应底行的元素与该负元素之比中最大者。

3.进行旋转运算。

如此往复直至第2个特点也被满足，并得到最优解。

该matlab程序可以解决在标准形下满足第1、3、4特点的线性规划问题。

二、算法框图

三、算法程序

建立函数

functionx=lindual(c,A,b)

[n1,n2]=size(A);

A=[-A,eye(n1)];c=[-c,zeros(1,n1)];

x1=[zeros(1,n2),b];lk=[n2+1:n1+n2];

b=-b;

while(1)

x=x1(1:n2);

s1=[lk,b,A];

c;

x1;

cc=[];ci=[];

fori=1:n1

ifb(i)0

cc=[cc,b(i)];

ci=[ci,i];

end

end

nc=length(cc);

ifnc==0

fprintf(到达最优解);

break

end

cliu=cc(1);

cl=ci(1);

forj=1:nc

ifabs(cc(j))abs(cliu)

cliu=cc(j);

cl=j;

end

end

cc1=[];ci1=[];

fori=1:n1+n2

ifA(cl,i)0

cc1=[cc1,A(cl,i)];

ci1=[ci1,i];

end

end

nc1=length(cc1);

ifnc1==0

fprintf(无可行解);

break

end

cliu=c(ci1(1))/cc1(1);

cl1=ci1(1);

forj=1:nc1

ifc(ci1(j))/cc1(j)cliu

cliu=c(ci1(j))/cc1(j);

cl1=ci1(j);

end

end

b(cl)=b(cl)/A(cl,cl1);

A(cl,:)=A(cl,:)/A(cl,cl1);

fork=1:n1

ifk~=cl

b(k)=b(k)-b(cl)*A(k,cl1);

A(k,:)=A(k,:)-A(cl,:).*A(k,cl1);

end

end

c=c-c(cl1).*A(cl,:);

x1(lk(cl))=0;

lk(cl)=cl1;

forkk=1:n1

x1(lk(kk))=b(kk);

end

x=x1(1:n2);

end

四、算法实现

例1.用对偶单纯形法解以下线性规划

min

s.t.6

5

x

解：利用matlab数学软件的解答方法为：

〔1〕输入题中条件：

c=[15245];A=[061;521];b=[2;1];

〔2〕调用函数求出最优解：

x=lindual(c,A,b)

〔3〕求最优值：

z=c*x’

运行结果如下：

例2.用对偶单纯形法解以下线性规划

min

s.t.2

x

x

解：利用matlab数学软件的解答方法为：

〔1〕输入题中条件：

c=[11];A=[21;17];b=[4;7];

〔2〕调用函数求出最优解：

x=lindual(c,A,b)

〔3〕求最优值：

z=c*x’

运行结果如下：