单纯形法.txt

单纯形法 主要进行最优求解

a=[-1 1 1 0 0 2;1 2 0 1 0 10;3 1 0 0 1 15];
b=[-2 -3 0 0 0 0];
c=[a;b];
b1=b(1,1:5);
[k,j]=min(b);

for i=1:3
o(i)=c(i,6)/c(i,j);
end
[d,e]=min(o);
k1=1/c(e,j)*c(e, : );
k2=-c(2,j)/c(e,j)*c(e, : )+c(2, : );
k3=-c(3,j)/c(e,j)*c(e, : )+c(3, : );
k4=-c(4,j)/c(e,j)*c(e, : )+c(4, : );
k=[k1;k2;k3;k4];

b=k4(1,1:5);
[k1,j]=min(b);
for i=1:3 
t(i)=k(i,6)/k(i,j);
if t(i)<0
    t(i)=100;
 else
    t(i)=t(i);
  end
end
[d,e]=min(t);
switch e
    case 1
        q1=1/k(1,j)*k(1, : );
        q2=-k(1,j)/k(1,j)*k(1, : )+k(2, : );
        q3=-k(3,j)/k(1,j)*k(1, : )+k(3, : );
        q4=-k(4,j)/k(1,j)*k(1, : )+k(4, : );
        q=[q1;q2;q3;q4];
                
    case 2
        q2=1/k(2,j)*k(2, : );
        q1=-k(1,j)/k(2,j)*k(2, : )+k(1, : );
        q3=-k(3,j)/k(2,j)*k(2, : )+k(3, : );
        q4=-k(4,j)/k(2,j)*k(2, : )+k(4, : );
        q=[q1;q2;q3;q4];
        
        
        
    case 3
          q3=1/k(3,j)*k(3, : );
          q1=-k(1,j)/k(3,j)*k(3, : )+k(1, : );
          q2=-k(2,j)/k(2,j)*k(3, : )+k(2, : );
          q4=-k(4,j)/k(4,j)*k(3, : )+k(4, : );
          q=[q1;q2;q3;q4];
    otherwise
        disp('不可能!!!')
end


b=q4(1,1:5);
[k2,j]=min(b);
for i=1:3 
t(i)=q(i,6)/q(i,j);
if t(i)<0
    t(i)=100;
 else
    t(i)=t(i);
  end
end
[d,e]=min(t);
switch e
    case 1
        y1=1/q(1,j)*q(1, : );
        y2=-q(1,j)/q(1,j)*q(1, : )+q(2, : );
        y3=-q(3,j)/q(1,j)*q(1, : )+q(3, : );
        y4=-q(4,j)/q(1,j)*q(1, : )+q(4, : );
        y=[y1;y2;y3;y4];
                
    case 2
        y2=1/q(2,j)*q(2, : );
        y1=-q(1,j)/q(2,j)*q(2, : )+q(1, : );
        y3=-q(3,j)/q(2,j)*q(2, : )+q(3, : );
        y4=-q(4,j)/q(2,j)*q(2, : )+q(4, : );
        y=[y1;y2;y3;y4];
        
        
        
    case 3
          y3=1/q(3,j)*q(3, : );
          y1=-q(1,j)/q(3,j)*q(3, : )+q(1, : );
          y2=-q(2,j)/q(3,j)*q(3, : )+q(2, : );
          y4=-q(4,j)/q(3,j)*q(3, : )+q(4, : );
          y=[y1;y2;y3;y4];
    otherwise
        disp('不可能!!!')
end