单纯形法(C程序).docVIP

下载本文档

32
0
约 5页
2018-10-12 发布于重庆
举报
版权申诉

单纯形法(C程序).doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

单纯形法(C程序)

#include stdio.h int main() { int i,j,r,k,l,jj[4],m=4,n=7,maxjj,mini,count=0; float a[4][7]={{0,2,-1,1,0,0,0},{60,3,1,1,1,0,0},{10,1,-1,2,0,1,0},{20,1,1,-1,0,0,1}}; float b[4][4],e[4][4],t[4][4],y[4],maxcj,tb[4],tp[4],cb[4],cj[7],theta,lk,z; jj[1]=4;jj[2]=5;jj[3]=6; for(i=0;im;i++) for(j=0;jm;j++) b[i][j]=0.0; for(i=1;im;i++) b[i][i]=1.0; printf(*********************************\n); for(i=0;im;i++) { for(j=0;jn;j++) printf(%6.1f,a[i][j]);printf(\n);} for(i=1;im;i++) cb[i]=a[0][jj[i]]; for(j=1;jm;j++) { y[j]=0.0; for(i=1;im;i++) y[j]+=cb[i]*b[i][j]; } for(j=1;jn;j++) { for(r=1;rm;r++) if(jj[r]==j){cj[j]=0.0;break; } cj[j]=a[0][j]; for(i=1;im;i++) cj[j]=cj[j]-y[i]*a[i][j]; } maxcj=cj[1]; maxjj=1; for(j=2;jn;j++) if(maxcjcj[j]) { maxcj=cj[j]; maxjj=j;} k=maxjj; for(j=1;jn;j++)printf(%6.1f,cj[j]);printf( k=%d\n,k); while(maxcj0) { count++; for(i=1;im;i++) { tb[i]=0.0; for(r=1;rm;r++)tb[i]+=b[i][r]*a[r][0]; } for(i=1;im;i++) { tp[i]=0.0; for(r=1;rm;r++)tp[i]+=b[i][r]*a[r][k]; } theta=1000.0; for(j=1;jm;j++) { if(tp[j]0) if (thetatb[j]/tp[j]) { theta=tb[j]/tp[j]; mini=j;} } l=mini; lk=tp[l]; jj[l]=k; printf(*********************************\n); for(i=0;im;i++) for(j=0;jm;j++) e[i][j]=0.0; for(i=1;im;i++) e[i][i]=1.0; for(i=1;im;i++)e[i][l]=(-1)*tp[i]/lk; e[l][l]=1.0/lk; for(i=1;im;i++) { for(j=1;jm;j++)printf(%6.1f,e[i][j]);printf(\n); } for(i=1;im;i++) for(j=1;jm;j++) { t[i][j]=0.0; for(r=1;rm;r++)t[i][j]+=e[i][r]*b[r][j];} for(i=1;im;i++) { for(j=1;jm;j++) { b[i][j]=t[i][j]; printf(%6.1f,b[i][j]);} printf( x%d\n,jj[i]); } for(i=1;im;i++) cb[i]=a[0][jj[i]]; for(j=1;jm;j++) { y[j]=0.0; for(i=1;im;i++) y[j]+=cb[i]*b[i][j]; } for(j=1;jn;j++) { for(r=1;rm;r++) if(jj[r]==j){cj[j]=0.0;break; } cj[j]=a[0][j]; for(i=1;im;i++) cj[j]=cj[j]-y[i]*a[i][j]; } maxcj=cj[1]; maxj