目标规划非线性单纯形算法.docxVIP

#include stdio.h #include stdlib.h #include iostream.h #define FILENAME in_25.txt float a[10][15]; //记录约束条件函数的变量系数 float b[15]; 〃记录约束条件的函数值 float c[15]; 〃记录目标函数的变量系数及检验数rj int m; 〃记录约束条件的个数 int n; 〃记录变量的个数 int or_n; 〃记录原始的变量个数 int flag; 〃记录标志，0为非规范型，1为规范型 //(规范型flag= 1时通过单纯形法直接求解) 〃(非规范型flag=0时通过二阶段单纯型法求解) int vi=0; 〃引入人工变量的个数 int manual_var_raw; //如果基变量中包含人工变量，则记录人工变量所在的行 int in_lie; 〃确定进基变量所在的列 int out_hang; 〃确定出基变量所在的行 void distill_target() 〃提取目标函数的系数并存放到数组c中 { FILE *fp; char ch,prior,middle,next; char ch_prior,ch_prior_prior; char ch_next,ch_next_next; int i; if((fp=fopen(FILENAME,r))==NULL) printf(Can not open this file!\n); ch=fgetc(fp); prior=middle=next=ch; while(ch!=#) { while(ch!=# ch!=x) { prior=middle; middle=next; ch=fgetc(fp); next=ch; }//while(ch!=# ch!=x) if(next==x) { ch=fgetc(fp); i=int(ch)-48; switch(middle) { case +: c[i]=1; break; case =: c[i]=1; break; case -: c[i]=-1;break; default: { if(prior==+|| prior===) c[i]=float(middle)-48; else if(prior=0 prior=9) { fseek(fp,-5l,1); ch_prior=fgetc(fp); if(ch_prior==+ || ch_prior===) c[i]=(float(prior)-48)*10+float(middle)-48; if(ch_prior==-) c[i]=-((float(prior)-48)*10+float(middle)-48); fseek(fp,4l,1); if(ch_prior=0 ch_prior=9) { fseek(fp,-6l,1); ch_prior_prior=fgetc(fp); if(ch_prior_prior==+ || ch_prior_prior===) c[i]=(float(ch_prior)-48)*100+(float(prior)-48)*10+float(middle)-48; if(ch_prior_prior==-) c[i]=-(float(ch_prior)-48)*100+((float(prior)-48)*10+float(middle)-48); fseek(fp,5l,1); } } else c[i]=-(float(middle)-48); } } if((ch=fgetc(fp))==/) // 若系数是分数 { ch=fgetc(fp); if(ch=0 ch=9) { ch_next=ch; ch=fgetc(fp); if(ch=0 ch=9) { ch_next_next=ch; c[i]=c[i]/((float(ch_next)-48)*10+float(ch_next_next)-48); prior=middle=ch_next; next=ch_next_next; } else { c[i]=c[i]/(float(ch_next)-48); prior=middle=next=ch_next; } } } else fseek(fp,-1l,1); } prior=middle; middle=next; next=ch; }//while(ch!=#) fclose(fp); } 〃函数inout()读取输入文件inout.txt中的约束条件的系数并存放到数组a中 〃确定约束条件的个数m和变量的个数n,并由此构成增广矩阵 void inout() { FIL