线性规划c实现.docxVIP

/* *@aim:线性规划,单纯形算法实现 *@note:这里的实现比常规的实现相比，需要更大的存储空间,但灵活性更强,所需要的代码更少 *@date:2015-6-11 *@Author:小花熊 *@Q:1198253149 */ #ifndef __SIMPLEX_H__ #define __SIMPLEX_H__ #includeCArray2D.h #includeCArray.h//单纯型算法所需要的数据结构//@note:使用的规则/* *@1:松弛约束表达式中第k个基本变量的下标是k,非基本变量的下标为m+i,m为基本变量的最大索引 *@2:系数矩阵的行列(m+n,m+n),m为基本变量数目,n非基本变量的数目 */structCSimplex {//记录基本变量的索引CArrayint *basicVariableIndexes;//记录非基本变量得索引CArrayint *nonbasicVariableIndexes;//记录每个松弛约束中的常数,其长度和nonbasicVariableIndexes+basicVariableIndexes相等CArrayfloat *relaxConstants;//松弛约束表达式的矩形阵列 CArray2Dfloat *relaxExpressRestrict;//目标函数,其内容的长度为nonbasicVariableIndexes-size+bacsicVariableIndexes-size//计算目标函数的最大值CArrayfloat *objectFunc;//目标函数中的常数项floatobjectConstant;//CSimplex();CSimplex(intbasicVariableSize,intnonVariableSize); ~CSimplex(); };enumCSimplexType {CSimplexTypeInvalide,//单纯型输入不可行CSimplexTypeNoBound,//单纯型算法所得到的结果是无界的CSimplexTypeOK,//单纯形算法可能得到一个正常的结果 }; //单纯形算法实现 //@request:调用前simple已经是松弛型的 //@return:返回上面的枚举类型中的一种,如果计算成功，则将结果写入到result中CSimplexTypesimplex_algorithm(CSimplex *simplex,float *result); #endif/* *@aim:单纯形算法实现 *@date:2015-6-11 18:59 */ #includesimplex.h #includeassert.h #includestdio.h //标准型转换为松弛型 //注意,函数可能会修改src的内容staticboolsimplex_normal_switch_relax(CSimplex *src,CSimplex *dst); //注意,函数可能会修改simplex的内容,如果初始现行约束不可行返回CSimplexTypeInvalidestaticCSimplexTypeinit_simplex(CSimplex *simple); //构造函数和析构函数CSimplex::CSimplex() {this-basicVariableIndexes=NULL;this-nonbasicVariableIndexes=NULL;this-relaxConstants=NULL;this-relaxExpressRestrict=NULL;this-objectFunc=NULL;this-objectConstant=0.0f; }CSimplex::CSimplex(intbasicVariableSize,intnonVariableSize) {this-basicVariableIndexes=new CArrayint(basicVariableSize);this-nonbasicVariableIndexes=new CArrayint(nonVariableSize);this-objectFunc=new CArrayfloat(basicVariableSize+nonVariableSize);this-relaxExpressRestrict=new CArray2Dfloat(basicVariableSize+n