3fga33.rar

3fga33.rar

#include “math.h“
#include “iostream.h“
#include “stdlib.h“


#define E1 0.0001//迭代收敛精度
#define ep 0.0001//复合形法中映射系数
#define n  3        //设计变量的维数
#define k  5     //复合形法中的顶点数


float af=10;//初始映射系数
int ij;
double X0[n]XX[n]X[k][n]FF[k];
void produce （double A[n]double B[n]）;
int p=0;

double a[n]b[n];//
double F（double C[n]）;//目标函数
int cons（double D[n]）;//约束域
void bou（）;//

double *Complex（）;

double F（double C[n]）//目标函数
{
	double F;
	//F=25/（C[0]*C[1]*C[1]*C[1]）;
	//F=4.629+3.297*C[0]*C[1]+（-0.249）*C[0]*C[2];
	F=6.7506*C[0]*C[0]-8.8049*C[0]+（-0.0458）*C[1]*C[1]*C[1]*C[1]+1.1309*C[1]*C[1]*C[1]-9.9148*C[1]*C[1]+35.2889*C[1]+0.0002*C[2]*C[2]+0.0032*C[2]-21.4933;
	return F;
}

int cons（double D[n]）//约束域
{
	//if （（（30/（D[0]*D[1]*D[1]）-50）<=0）&&（（0.0004*D[0]*D[1]-0.001）<=0）&&（（D[0]>=2）&&（D[0]<=4））&&（（D[1]>=0.5）&&（D[1]<=1）））
	if（（（D[0]>=0.5）&&（D[0]<=1.3））&&（（D[1]>=3）&&（D[1]<=9））&&（（D[2]>=60）&&（D[2]<=140）））
		//D[2]为热风温度，D[1]热风速度，D[0]初剂含量
	return 1;//若在约束域内则返回1
	else 
		return 0;
}

void bou（）//计算变量的边界条件
{
	a[0]=0.5;a[1]=3;a[2]=60;
	b[0]=1.3;b[1]=9;b[2]=140;
	
}

void produce（double A[n]double B[n]）//产生复合的K个顶点
{
	double rr;int jj;
	do
	{
		rr=rand（）;  rr=rr/32767;//调用rand随机生成
		for（i=0;i		{
			X[0][i]=A[i]+rr*（B[i]-A[i]）;
			XX[i]=X[0][i];//产生起始点的维坐标
		}
	}
	while （cons（XX）==0）;//若产生的起始点不在约束范围内，刚继续重新生成，直到在约束范围内为止
	for（j=1;j	{	rr=rand（）;rr=rr/32767;
		for（i=0;i		{
			X[j][i]=A[i]+rr*（B[i]-A[i]）;//产生其余的K-1个点
		}
	}
	
	for（j=0;j	{
		for（i=0;i		{
			X0[i]=0;
			for（jj=0;jj			{
				X0[i]+=X[jj][i];//
			}
			X0[i]=（1/（j+1.0））*（X0[i]）;//求出这K个点的形心
			XX[i]=X[j][i];//第J个点
		}
		while（cons（XX）==0）//若其余的K-1个点有某个点不在可行域内
		{
			for（i=0;i			{
				X[j][i]=X0[i]+0.5*（X[j][i]-X0[i]）;//产生新的点，直到这个点在可行域内
				XX[i]=X[j][i];//
			}
		}//至此所有的点均在可行域内
	}
}



double *Complex（）
{
	double EEXc[n]Xh[n]Xg[n]Xl[n]FhFgFlXr[n]*prXl;
	int hgl;

	bou（）;
	produce（ab）;
	EE=E1+1;
	while（EE>E1）//初始迭代条件大于终于迭代条件
	{
		for（i=0;i		{
			Xc[i]=0;
			for（j=0;j			{
				Xc[i]+=X[j][i];
			}
			Xc[i]=1/（k+0.0）*Xc[i];//生成形心
		}
		for（j=0;j		{
			for（i=0;i			{ XX[i]=X[j][i];}
			FF[j]=F（XX）;
		}
		Fh=FF[0];Fg=FF[0];Fl=FF[0];h=0;g=0;l=0;
		for（j=0;j		{
			if（FF[j]>Fh）
			{ Fh=FF[j];h=j;}//最坏点FH
			if（FF[j]			{ Fl=FF[j];l=j;}//最好点
		}
		for（j=0;j			if（（FF[j]>Fg）&&（j!=h））//次坏点FG
			{ Fg=FF[j];g=j; }
		for（i=0;i		{
			Xh[i]=X[h][i];Xg[i]=X[g][i];Xl[i]=X[l][i];//分别在某维时的最坏，次坏，最好点
		}
		EE=0;
		for（j=0;j		{  EE+=pow（（FF[j]-F（Xc））2）;}//顶点与目标函数值之间之差的二次幂的和
		EE=pow（（1/（k+0.0）*EE）0.5）;//和的均方根
		
		
		p++;


		if（EE<=E1）//若满足终止迭代条件
		{
			cout<<“\n F（Xmin）=“<			prXl=&X[l][0];break;//最优解，跳出while循环
		}

		for（i=0;i		{
			Xc[i]=0;
			for（j=0;j			{
				if（j!=h）
					Xc[i]+=X

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件    2506752  2009-05-13 11:40  最优化作业\凹印优化设计.doc

     文件      19968  2009-06-27 10:15  最优化作业\数据公式表.xls

     文件       4859  2009-06-27 10:14  最优化作业\最优化程序\Test\test.cpp

     文件       4260  2009-06-04 11:23  最优化作业\最优化程序\Test\Test.dsp

     文件        516  2009-06-04 11:23  最优化作业\最优化程序\Test\Test.dsw

     文件      50176  2009-06-30 10:04  最优化作业\最优化程序\Test\Test.ncb

     文件      74752  2009-06-30 10:04  最优化作业\最优化程序\Test\Test.opt

     文件        242  2009-06-30 10:01  最优化作业\最优化程序\Test\Test.plg

     文件       4096  2009-06-04 11:26  最优化作业\最优化程序\Test\Untitled Project.IAB

     文件        368  2009-06-04 11:26  最优化作业\最优化程序\Test\Untitled Project.IAD

     文件       4096  2009-06-04 11:26  最优化作业\最优化程序\Test\Untitled Project.IMB

     文件        368  2009-06-04 11:26  最优化作业\最优化程序\Test\Untitled Project.IMD

     文件          4  2009-06-04 11:54  最优化作业\最优化程序\Test\Untitled Project.PFI

     文件        776  2009-06-04 11:54  最优化作业\最优化程序\Test\Untitled Project.PO

     文件        824  2009-06-04 11:27  最优化作业\最优化程序\Test\Untitled Project.PR

     文件       2536  2009-06-04 11:28  最优化作业\最优化程序\Test\Untitled Project.PRI

     文件       7048  2009-06-04 12:08  最优化作业\最优化程序\Test\Untitled Project.PS

     文件      19779  2009-06-04 12:08  最优化作业\最优化程序\Test\Untitled Project.WK3

     目录          0  2009-06-30 10:05  最优化作业\最优化程序\Test

     目录          0  2009-06-30 10:04  最优化作业\最优化程序

     目录          0  2009-06-30 10:04  最优化作业

----------- ---------  ---------- -----  ----

              2701420                    21