1 單純形法

（1） 單純形法是解線性規劃問題的一個重要方法。
其原理的基本框架為：
第一步：將LP線性規劃變標準型，確定一個初始可行解（頂點）。
第二步：對初始基可行解最優性判別，若最優，停止；否則轉下一步。
第三步：從初始基可行解向相鄰的基可行解（頂點）轉換，且使目標值有所改善—目標函式值增加，重複第二和第三步直到找到最優解。
（2） 用程式進行運算前，要將目標函式及約束方程變成標準形式。

於非標準形式須作如下變換：
a) 目標函式為極小值min z=CX時，轉換為max z=-CX形式；
b) 在約束方程中有 “≤”時，在加上一個鬆弛變數；
c) 在約束方程中有 “≥”時，採用減去一個鬆弛變數，再加上一個人工變數；
d) 在約束方程中有 “=”時，加上一個人工變數；
e) 所有的人工變數，鬆弛變數的目標函式係數置為0。
（3） 對於標準形式的線性規劃問題。用單純形法計算步驟的框圖。

2 程式測試及結果：

線性規劃問題如下:
max z=2*x1-3*x2+3x3;
x1+ x2 -x3<=7;
x1- x2 +x3<=-7;
x1-2*x2 +2*x3<=4;
x1,x2,x3>=0;

3 C++實現程式碼

// Simplex.cpp : 定義控制檯應用程式的入口點。
//
//
/********************************* 
----------------------------------- 
單純形法求解線性規劃問題（C++實現程式碼）
----------------------------------- 
Author:牧之丶  Date：2014年
Email:
 
[email protected] 
**********************************/  
#include "stdafx.h"
#include<iostream> 
#include<math.h> 
using namespace std;
#define M 10000        //全域性變數大M
 float juzhen[11][31];//核心矩陣表 
 int m=0,n=0,t=0;//m:結構向量的個數 //n:約束不等式個數  //t:目標函式型別：－1代表求求最小值，1代表求最大值                          
void input() //輸入介面函式
 

{   
    int i,j; 
    cout<<"－－－－－－－－－－單純形法的參 數 輸 入－－－－－－－－－－－"<<endl; 
    cout<<"請按提示輸入下列引數："<<endl<<endl; 
    cout<<"  結構向量數m:  "<<"   m=  "; 
    cin>>m; 
    cout<<endl<<"  約束不等式數n:"<<"   n=  "; 
    cin>>n; 
    for (i=0;i<=n+1;i++) 
    for (j=0;j<=m+n+n;j++)  
         juzhen [i][j]=0;      //初始化矩陣，所有元素均為0
   //讀入約束條件 
   cout<<endl<<"  約束方程矩陣的係數及不等式方向（1代表<=,-1代表>=）:"<<endl<<endl<<"          "; 
    for (i=1;i<=m;i++) 
       cout<<"   x"<<i; 
       cout<<"  不等式方向 "<<" 常數項"<<endl; 
    for (i=1;i<=n;i++)  
    {  
      cout<<"    不等式"<<i<<"  "; 
      for (j=1;j<=m+2;j++)  
          cin>>juzhen [i][j]; 
    } 
   for (i=1;i<=n;i++) 
     {  
       juzhen [i][0]=juzhen [i][m+2]; 
       juzhen [i][m+2]=0; 
     } 
  //讀入目標條件 
        cout<<endl<<endl<<" 目標函式的係數及型別（求最小值:1；求最大值:-1）:"<<endl<<endl<<"                "; 
   for(i=1;i<=m;i++) 
        cout<<"x"<<i<<"   "; 
        cout<<"型別"<<endl<<"  "; 
        cout<<"  目標函式:   "; 
   for (i=1;i<=m;i++) 
     cin>>juzhen [0][i]; 
     cin>>t; 
  //矩陣調整 
  if(t==-1)  
      for(i=1;i<=m;i++)  
          juzhen [0][i]=(-1)*juzhen [0][i]; 
  for(i=1;i<=n;i++) 
  {  
      juzhen [i][m+i]=juzhen [i][m+1]; 
         if(i!=1) 
             juzhen [i][m+1]=0; 
  } 
} 
//演算法函式 
void comput() 
{ 
    int i,j,flag,temp1,temp2,h,k=0,temp3[10]; 
    float a,b[11],temp,temp4[11],temp5[11],f=0,aa,d,c; 
    //初始化 
  for(i=1;i<=n;i++)  
     temp3[i]=0; 
   for(i=0;i<11;i++) 
   {    temp4[i]=0; 
        temp5[i]=0; 
   } 
    for(i=1;i<=n;i++) 
  { 
     if(juzhen [i][m+i]==-1) 
      { 
       juzhen [i][m+n+i]=1; 
       juzhen [0][m+n+i]=M; 
       temp3[i]=m+n+i; 
      } 
     else  
         temp3[i]=m+i; 
  } 
     for(i=1;i<=n;i++) 
           temp4[i]=juzhen [0][temp3[i]]; 
     //迴圈求解 
do{    
    for(i=1;i<=m+n+n;i++) 
     {   
        a=0; 
        for(j=1;j<=n;j++) 
          a+=juzhen [j][i]*temp4[j]; 
          juzhen [n+1][i]=juzhen [0][i]-a; 
        } 
    for(i=1;i<=m+n+n;i++) 
         { 
            if(juzhen [n+1][i]>=0) flag=1; 
            else  
            {   
                flag=-1; 
                 break; 
             } 
          } 
    if(flag==1) 
    {     for(i=1;i<=n;i++) 
             {  
                  if(temp3[i]<=m+n) temp1=1; 
                           else 
                              {  
                                  temp1=-1; 
                                  break; 
                               } 
              } 
    //輸出結果 
    cout<<endl<<endl; 
    cout<<"－－－－－－－－－－結 果 輸 出－－－－－－－－－－－"<<endl<<endl; 
        if(temp1==1) 
          {  
             cout<<" 此線性規劃的最優解存在!"<<endl<<endl<<"  最優解為："<<endl<<endl<<"     "; 
              for(i=1;i<=n;i++) 
                 temp5[temp3[i]]=juzhen [i][0]; 
              for(i=1;i<=m;i++) 
                 f+=t*juzhen [0][i]*temp5[i];          
              for(i=1;i<=m;i++)  
              { 
                 cout<<"x"<<i<<" = "<<temp5[i]; 
                 if(i!=m) 
                     cout<<"， "; 
              } 
              cout<<" ;"<<endl<<endl<<"     最優目標函式值f= "<<f<<endl<<endl; 
              return ; 
            } 
            else  
            {    
                cout<<" 此線性規劃無解"<<endl<<endl; 
                return ; 
            } 
      } 
   if(flag==-1) 
   {    
       temp=100000; 
       for(i=1;i<=m+n+n;i++) 
          if(juzhen [n+1][i]<temp)     
            { 
              temp=juzhen [n+1][i]; 
              h=i;                     
            }               
        for(i=1;i<=n;i++) 
          {  
           if(juzhen [i][h]<=0) temp2=1; 
                else {  
                            temp2=-1; 
                              break; 
                    }                        
          } 
   } 
    if(temp2==1)   
     {  
        cout<<"此線性規劃無約束"; 
         return  ; 
     } 
   if(temp2==-1)  
   {    
      c=100000; 
      for(i=1;i<=n;i++) 
        { 
           if(juzhen [i][h]!=0)  b[i]=juzhen [i][0]/juzhen [i][h]; 
           if(juzhen [i][h]==0)  b[i]=100000; 
           if(b[i]<0)     b[i]=100000; 
           if(b[i]<c)  
           { 
              c=b[i]; 
               k=i; 
          } 
      } 
         temp3[k]=h; 
         temp4[k]=juzhen[0][h]; 
         d=juzhen [k][h]; 
      for(i=0;i<=m+n+n;i++) 
        juzhen [k][i]=juzhen [k][i]/d; 
      for(i=1;i<=n;i++) 
      { if(i==k)  
        continue; 
         aa=juzhen [i][h]; 
         for(j=0;j<=m+n+n;j++) 
           juzhen [i][j]=juzhen [i][j]-aa*juzhen [k][j];  
      } 
   }     
} 
while(1); 
     return ; 
}

int _tmain(int argc, _TCHAR* argv[])
{   cout<<"-------------------單純形演算法程式----------------------"<<endl<<endl; 
  input(); 
  comput();
  system("pause");
    return 0;
}