对偶单纯形法C语言实现.docVIP

下载本文档

35
0
约3.82千字
约 7页
2018-06-17 发布于河南
举报
版权申诉

对偶单纯形法C语言实现.doc

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

对偶单纯形法C语言实现

程序源代码： #includestdio.h #includemath.h #define MAX_N 100 #define M 1000000.0; int m, n; double A[MAX_N][MAX_N], C[MAX_N], b[MAX_N], seta[MAX_N], CZ[MAX_N]; int num[MAX_N]; int GetOutCow() //求换出基 { int i,k; int flag; double min = 0; for (i = 0; i m; i ++) { if (b[i] = 0) flag = 1; else { flag = 0; break; } } if (flag == 1) return -1; for (i = 0; i m; i ++) { if (b[i] 0 min b[i]) { min = b[i]; k = i; } } return k; } int GetInCow(int p) //求换入基 { int i, j; int flag = 0; double min; for (j = 0; j n; j ++) { if (A[p][j] = 0) flag = 1; else { flag = 0; break; } } if (flag == 1) { printf(\n原线性规划问题无可行解!\n); return -1; } for (j = 0; j n; j ++) { if (A[p][j] 0) seta[j] = CZ[j] / A[p][j]; else seta[j] = M; } min = M; for (j = 0; j n; j ++) { if (min = seta[j]) { min = seta[j]; i = j; } } num[p] = i + 1; return i; } void change(int p, int q) //计算新的单纯形表 { int i, j; double temp1, temp2, temp3; temp1 = A[p][q]; for (i = 0; i m; i ++) { if (i != p) { if (A[i][q] != 0) { temp2 = A[i][q] / temp1; for (j = 0; j n; j ++) A[i][j] = A[i][j] - A[p][j] * temp2; b[i] = b[i] - b[p] * temp2; } } } temp3 = CZ[q] / temp1; for (i = 0; i n; i ++) CZ[i] = CZ[i] - A[p][i] * temp3; for (j = 0; j n; j ++) A[p][j] = A[p][j] / temp1; b[p] = b[p] / temp1; } void print1() { int i; printf(\n--------------------------------------------------------------------------\n); printf(\t); for(i = 1; i = n; i ++) printf(X(%d)\t, i); printf(RHS\n); printf(--------------------------------------------------------------------------\n); printf(\t); for (i = 0; i n; i ++) { printf(%.3lf\t, -C[i]); } printf(\n--------------------------------------------------------------------------\n); } void print2() { int i, j; printf(\n--------------------------------------------------------------------------\n); for (i = 0; i m; i ++) { prin