3fga33.rar
3fga33.rar
代码片段和文件信息
属性 大小 日期 时间 名称
----------- --------- ---------- ----- ----
文件 2506752 2009-05-13 11:40 最优化作业凹印优化设计.doc
文件 19968 2009-06-27 10:15 最优化作业数据公式表.xls
文件 4859 2009-06-27 10:14 最优化作业最优化程序Test est.cpp
文件 4260 2009-06-04 11:23 最优化作业最优化程序TestTest.dsp
文件 516 2009-06-04 11:23 最优化作业最优化程序TestTest.dsw
文件 50176 2009-06-30 10:04 最优化作业最优化程序TestTest.ncb
文件 74752 2009-06-30 10:04 最优化作业最优化程序TestTest.opt
文件 242 2009-06-30 10:01 最优化作业最优化程序TestTest.plg
文件 4096 2009-06-04 11:26 最优化作业最优化程序TestUntitled Project.IAB
文件 368 2009-06-04 11:26 最优化作业最优化程序TestUntitled Project.IAD
文件 4096 2009-06-04 11:26 最优化作业最优化程序TestUntitled Project.IMB
文件 368 2009-06-04 11:26 最优化作业最优化程序TestUntitled Project.IMD
文件 4 2009-06-04 11:54 最优化作业最优化程序TestUntitled Project.PFI
文件 776 2009-06-04 11:54 最优化作业最优化程序TestUntitled Project.PO
文件 824 2009-06-04 11:27 最优化作业最优化程序TestUntitled Project.PR
文件 2536 2009-06-04 11:28 最优化作业最优化程序TestUntitled Project.PRI
文件 7048 2009-06-04 12:08 最优化作业最优化程序TestUntitled Project.PS
文件 19779 2009-06-04 12:08 最优化作业最优化程序TestUntitled Project.WK3
目录 0 2009-06-30 10:05 最优化作业最优化程序Test
目录 0 2009-06-30 10:04 最优化作业最优化程序
目录 0 2009-06-30 10:04 最优化作业
----------- --------- ---------- ----- ----
2701420 21
#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<<“
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 est.cpp
文件 4260 2009-06-04 11:23 最优化作业最优化程序TestTest.dsp
文件 516 2009-06-04 11:23 最优化作业最优化程序TestTest.dsw
文件 50176 2009-06-30 10:04 最优化作业最优化程序TestTest.ncb
文件 74752 2009-06-30 10:04 最优化作业最优化程序TestTest.opt
文件 242 2009-06-30 10:01 最优化作业最优化程序TestTest.plg
文件 4096 2009-06-04 11:26 最优化作业最优化程序TestUnti
文件 368 2009-06-04 11:26 最优化作业最优化程序TestUnti
文件 4096 2009-06-04 11:26 最优化作业最优化程序TestUnti
文件 368 2009-06-04 11:26 最优化作业最优化程序TestUnti
文件 4 2009-06-04 11:54 最优化作业最优化程序TestUnti
文件 776 2009-06-04 11:54 最优化作业最优化程序TestUnti
文件 824 2009-06-04 11:27 最优化作业最优化程序TestUnti
文件 2536 2009-06-04 11:28 最优化作业最优化程序TestUnti
文件 7048 2009-06-04 12:08 最优化作业最优化程序TestUnti
文件 19779 2009-06-04 12:08 最优化作业最优化程序TestUnti
目录 0 2009-06-30 10:05 最优化作业最优化程序Test
目录 0 2009-06-30 10:04 最优化作业最优化程序
目录 0 2009-06-30 10:04 最优化作业
----------- --------- ---------- ----- ----
2701420 21
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌抄袭侵权/违法违规的内容, 请发送邮件举报,一经查实,本站将立刻删除。
评论列表(条)