#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）;//若產生的起始點不在約束范圍內，剛繼續(xù)重新生成，直到在約束范圍內為止
	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;//最優(yōu)解，跳出while循環(huán)
		}

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

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????2506752??2009-05-13?11:40??最優(yōu)化作業(yè)\凹印優(yōu)化設計.doc

?????文件??????19968??2009-06-27?10:15??最優(yōu)化作業(yè)\數據公式表.xls

?????文件???????4859??2009-06-27?10:14??最優(yōu)化作業(yè)\最優(yōu)化程序\Test\test.cpp

?????文件???????4260??2009-06-04?11:23??最優(yōu)化作業(yè)\最優(yōu)化程序\Test\Test.dsp

?????文件????????516??2009-06-04?11:23??最優(yōu)化作業(yè)\最優(yōu)化程序\Test\Test.dsw

?????文件??????50176??2009-06-30?10:04??最優(yōu)化作業(yè)\最優(yōu)化程序\Test\Test.ncb

?????文件??????74752??2009-06-30?10:04??最優(yōu)化作業(yè)\最優(yōu)化程序\Test\Test.opt

?????文件????????242??2009-06-30?10:01??最優(yōu)化作業(yè)\最優(yōu)化程序\Test\Test.plg

?????文件???????4096??2009-06-04?11:26??最優(yōu)化作業(yè)\最優(yōu)化程序\Test\Untitled?Project.IAB

?????文件????????368??2009-06-04?11:26??最優(yōu)化作業(yè)\最優(yōu)化程序\Test\Untitled?Project.IAD

?????文件???????4096??2009-06-04?11:26??最優(yōu)化作業(yè)\最優(yōu)化程序\Test\Untitled?Project.IMB

?????文件????????368??2009-06-04?11:26??最優(yōu)化作業(yè)\最優(yōu)化程序\Test\Untitled?Project.IMD

?????文件??????????4??2009-06-04?11:54??最優(yōu)化作業(yè)\最優(yōu)化程序\Test\Untitled?Project.PFI

?????文件????????776??2009-06-04?11:54??最優(yōu)化作業(yè)\最優(yōu)化程序\Test\Untitled?Project.PO

?????文件????????824??2009-06-04?11:27??最優(yōu)化作業(yè)\最優(yōu)化程序\Test\Untitled?Project.PR

?????文件???????2536??2009-06-04?11:28??最優(yōu)化作業(yè)\最優(yōu)化程序\Test\Untitled?Project.PRI

?????文件???????7048??2009-06-04?12:08??最優(yōu)化作業(yè)\最優(yōu)化程序\Test\Untitled?Project.PS

?????文件??????19779??2009-06-04?12:08??最優(yōu)化作業(yè)\最優(yōu)化程序\Test\Untitled?Project.WK3

?????目錄??????????0??2009-06-30?10:05??最優(yōu)化作業(yè)\最優(yōu)化程序\Test

?????目錄??????????0??2009-06-30?10:04??最優(yōu)化作業(yè)\最優(yōu)化程序

?????目錄??????????0??2009-06-30?10:04??最優(yōu)化作業(yè)

-----------?---------??----------?-----??----

??????????????2701420????????????????????21