資源簡介
大M法matlab代碼 用于幫助學生學習在不能直接找到初始解的情況下如何使用大M法
代碼片段和文件信息
clc;
A=[111100100;-120010010;-10200-1001];
b=[6;4;2];
c=[-2;1;-1;0;0;0;99999;99999;99999];
[mn]=size(A);
x0=[0;0;0;0;0;0;6;4;2];
v=[789123456];
x0=x0(v:);
while?1
????m_B=A(:v(1:m));
????m_N=A(:v(m+1:n));
????m_I=eye(n-mn-m);
????d=[-inv(m_B)*m_N;m_I];
????cp=c(v:);
????r=cp‘*d;
????if?all(r>=0)
????????display(sprintf(‘x?is?optimal‘));
????????cp=c(v:);
????????z=cp‘*x0;
????????return;
????end
????[rqq]=min(r)
評論
共有 條評論