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

#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