?A=[0?1?0?0?1?0?0?1?0?0;??????%?不等式約束????????????
????3?0?0?1?0?0?2?0?0?0;
????0?3?0?0?1?0?0?2?0?0;
????0?0?3?0?0?1?0?0?2?0;
????1?1?1?0?0?0?0?0?0?0;
????0?0?0?1?1?1?0?0?0?0];
b=[650;650;850;400;650;555];
AA=[0?0?0?1?1?1?0?0?0?-1????%?等式約束
????1?0?0?1?0?0?1?0?0?0;
????0?0?1?0?0?1?0?0?1?0;
????0?0?0?0?0?0?1?1?1?0];
?bb=[185;450;350;400];
c=?[-1150?-1150?-1150?-760?-760?-760?-350?-350?-350?0?];???%?目標(biāo)函數(shù)系數(shù)

A=input（‘請(qǐng)輸入不等式約束方程組系數(shù)矩陣：‘）;?????????????%?單純形法求解標(biāo)準(zhǔn)形線性規(guī)劃問題:?min?c*x?s.t.?A*x<=b??AA*x=bb??
b=input（‘請(qǐng)輸入不等式約束方程組常數(shù)項(xiàng)向量b（b>=0）：‘）;???%其中bb>=0??b>0，??x>=0
AA=input（‘請(qǐng)輸入等式約束方程組系數(shù)矩陣：‘）;???????????????
bb=input（‘請(qǐng)輸入等式約束方程組常數(shù)項(xiàng)向量bb（bb>=0）：‘）;
c=input（‘請(qǐng)輸入目標(biāo)函數(shù)系數(shù)向量:‘）;
model=input（‘請(qǐng)輸入采用的單純形法，選擇大M法請(qǐng)輸入1，選擇兩階段法請(qǐng)輸入2：‘）;


if?isempty（AA）&&isempty（A）???????????????????%矩陣全為空
????disp（‘矩陣全為空，無法計(jì)算退出程序‘）;
end


if?（~isempty（AA））&&（~isempty（A））????????????%矩陣全不空，
????[e1f1]=size（A）;
????A=[Aeye（e1）];
????[e2f2]=size（AA）;
????AA=[AAzeros（e2e1）];
????A=[A;AA];??????????????????????????????%得到了添加了松弛變量的系數(shù)矩陣?。
????b=[b;bb];????????????????????????%b為約束方程組常數(shù)項(xiàng)。??
????c1=[czeros（1e1）];??????????%到這，已經(jīng)化成了標(biāo)準(zhǔn)型。
??????if?model==2
?????????[xzstage1_sturct_STstage2_sturct_STST]?=?twophase（c1Ab）;
??????????x=x（1:f1）????%?輸出參數(shù):?x最優(yōu)解?z最優(yōu)目標(biāo)函數(shù)值ST存儲(chǔ)單純形表數(shù)據(jù)
??????????z
??????????ST?%stage1_sturct_ST是第一階段的單純形表的集合。%stage2_sturct_ST是第二階段的單純形表的集合。若需顯示，你自行寫入
????????????????%ST是最終的單純形表。
??????elseif?model==1
?????????[xz~~~sturct_STST]?=?MySimplex_method（c1Ab）?;
??????????x=x（1:f1）
??????????z
??????????ST?????????????????%sturct_ST是單純形表的集合（大M法只需要一個(gè)單純形表）%ST是最終的單純形表。
??????else
?????????disp（‘方法選擇錯(cuò)誤，退出程序。‘）
?
??????end
???
end




if?isempty（AA）&&（~isempty（A））????%當(dāng)AA為空，A不為空時(shí)。
?????[e1f1]=size（A）;
??????A=[Aeye（e1）];
??????b=b;
??????c1=[czeros（1e1）];?%到這，已經(jīng)化成了標(biāo)準(zhǔn)型。
??????if?model==2
??????[xzstage1_sturct_STstage2_sturct_STST]?=?twophase（c1Ab）;
??????z
???????x=x（1:f1）
???????ST
??????
?????elseif?model==1
??????[xz~~~sturct_STST]?=?MySimplex_method（c1Ab）?;
???????x=x（1:f1）
???????z
???????ST
?????else
????????disp（‘方法選擇錯(cuò)誤，退出程序。‘）
?
?????end

end



if?（~isempty（AA））&&isempty（A）????%當(dāng)A為空，AA不為空時(shí)。
????A=AA;
????b=bb;
????c1=c;??%本身已經(jīng)是標(biāo)準(zhǔn)型
????if?model==2
?????[xzstage1_sturct_STstage2_sturct_STST]?=?twophase（c1Ab）;
??????x
??????z
??????ST
???%stage1_sturct_ST是第一階段的單純形表的集合。%stage2_sturct_ST是第二階段的單純形表的集合。
???elseif?model==1
??????[xz~~~sturct_STST]?=?MySimplex_method（c1Ab）?;
???????x
???????z
???????ST
?????else
????????disp（‘方法選擇錯(cuò)誤，退出程序。‘）
????end
?

end