单纯形法MATLAB程序.docx

《单纯形法MATLAB程序.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、。单纯形法（Matlab程序）%%单纯形法（Matlab程序）a=input('inputthemajormatrixA');b=input('inputthematrixb');n=input('inputthejudgement');%%g为计数器(确定循环次数)g=0;whileg<40%%确定非负alength=max(size(n));blength=max(size(b));m=0;fori=1:alengthifn(i)>=0m=m+1;endend;ifm==alengthx=b;breakend;%%找Ks=min

2、(n);fori=1:alengthifn(i)==sk=i;breakend;end;%%a[i,k]的非负性m=0;fori=1:blengthifa(i,k)<0m=m+1;end;end;ifm==blengthdisp('xdoesnotexit');judge=1;breakend;。1。%%找L确定主元cc=100000;fori=1:blengthifa(i,k)>0if(b(i)/a(i,k))

3、a(i,k))==ccl=i;breakendendend;%%计算,a标准化zu=a(l,k);aa=a;fori=1:l-1forj=1:alengthaa(i,j)=a(i,j)-a(l,j)*a(i,k)/a(l,k);endend;fori=l+1:blengthforj=1:alengthaa(i,j)=a(i,j)-a(l,j)*a(i,k)/a(l,k);endend;forj=1:alengthaa(l,j)=a(l,j)/zu;end;%%b的判别bb=b;bb(l)=b(l)/zu;fori=1:l-1bb(i

4、)=b(i)-b(l)*a(i,k)/a(l,k);end;fori=l+1:blengthbb(i)=b(i)-b(l)*a(i,k)/a(l,k);end;b=bb;%%确定判别数。2。tt=n;forj=1:alengthtt(j)=n(j)-a(l,j)*n(k)/a(l,k);end;n=tt;a=aa;%%显示单纯形表sasa=[b'aa;0n];disp('单纯表示例');disp(g+1);disp(sa);g=g+1;judge=2;end;ifjudge==2q=0;result=zeros(alength,2)

5、;forj=1+q:alengthifn(j)==0t=a(:,j);zu=find(t);result(j,1)=j;result(j,2)=x(zu);q=q+1;endifn(j)>0result(j,1)=q+1;q=q+1;endend;disp('最优解');disp(result);disp('循环次数');gend。3。欢迎您的下载，资料仅供参考！致力为企业和个人提供合同协议，策划案计划书，学习资料等等打造全网一站式需求。4