混合惩罚函数程序.docVIP

混合惩罚函数程序 #include iostream #include windows.h//记时间必须要加的头文件 using namespace std; double func0(double *a)//原目标函数 { return 2*pow(a[0],2)+3*pow(a[1],2)+5*pow(a[2],2)+a[0]+2*a[1]-3*a[2]; } double func1(double *a,double r)//其中r是一个递减的惩罚因子 { return 2*pow(a[0],2)+3*pow(a[1],2)+5*pow(a[2],2)+a[0]+2*a[1]-3*a[2] -r*log(a[0]*a[1]*a[2])+1/sqrt(r)*(pow(a[0]+a[1]-5,2)+pow(a[1]+a[2]-10,2)); } double func11(double *a)//不等式1约束 { return -a[0]; } double func22(double *a)//不等式2约束 { return -a[1]; } double func33(double *a)//不等式3约束 { return -a[2]; } double func44(double *a)//等式1约束 { return a[0]+a[1]-5; } double func55(double *a)//等式2约束 { return a[1]+a[2]-10; } double sousuo(double*aa,double *bb,double r,int n)//参数：搜索方向,初始点,影响系数/用于求解最优步长 { double l=1e-2; double a1,b1,c1; a1=0; c1=a1+l; double *a,*b,*c; a=(double *)malloc(sizeof(double)*n); b=(double *)malloc(sizeof(double)*n); c=(double *)malloc(sizeof(double)*n); for(int i=0;in;i++) *(c+i)=bb[i]+c1*aa[i]; while(func11(c)0||func22(c)0||func33(c)0)//再从正方向进行搜索,防止跑出范围 { l/=2; c1=a1+l; for(int i=0;in;i++) *(c+i)=bb[i]+c1*aa[i]; } for(int i=0;in;i++) { *(a+i)=bb[i]+a1*aa[i]; *(c+i)=bb[i]+c1*aa[i]; }//确定l以后，就进行初始点的计算,此时l是一个正值 if(fabs(func1(a,r)-func1(c,r))1e-8) { if(func1(a,r) func1(c,r))//负方向进行搜索 { l=-l; c1=a1+l; b1=c1+l; for(int i=0;in;i++) { *(c+i)=bb[i]+c1*aa[i]; *(b+i)=bb[i]+b1*aa[i]; } while(func11(b)0||func22(b)0||func33(b)0)//用于防止搜索点跑出约束范围 { l/=2;//这个地方尽量是除以2，具体原因还不清楚 b1=b1-l; for(int i=0;in;i++) *(b+i)=bb[i]+b1*aa[i]; } while(func1(b,r) func1(c,r))//确定求值区间 { l*=2; a1=c1; c1=b1; b1=b1+l; for(int i=0;in;i++) { *(c+i)=bb[i]+c1*aa[i]; *(b+i)=bb[i]+b1*aa[i]; } while(func11(b)0||func22(b)0||func33(b)0)//用于防止搜索点跑出约束范围 { l/=2; b1=b1-l; for(int i=0;in;i++) *(b+i)=bb[i]+b1*aa[i]; } }//。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。这个地方需要修改 double temp;//交换区间边界两个值 temp=a1; a