文件名称：区间约束Powell方法寻优代码
文件大小：144KB
文件格式：7Z
更新时间：2013-03-25 14:43:47
Powell优化 约束 无约束 BOBYQA算法
powell方法是计算最值的有效方法，一般情况下是无约束的，目前Powell已发展了一种称为BOBYQA的带区间约束算法。 压缩文件里面optimization.h是无约束powell寻优的一个类，以前在baidu上找到的。powell.h和powell.cpp是带区间约束的powell算法BOBYQA的C++包装,由于BOBYQA算法只能计算二维以上，一维算法是俺自己写的。调用时，一维算法采用ExecuteBrent，二维以上采用BoundedPowell。 BOBYQA目录下是powell方法祖师爷M. J. D. Powell写的FORTRAN代码，powell.lib是俺把Fortran代码编译成的静态库。由于这个lib是Fortran Power Station 4.0做的，msfrt40.dll就成了一个甩不掉的累赘。 如果使用Fortran编程的话，就不需要俺的包装了，直接用Fortran代码吧。 使用方法请参考Fortran代码中main.f，有问题请联系baita00@yahoo.com.cn。 另外链接时好像还需要fps4.0中的msfrt.lib，原来没有加进去，现在加不进去了。需要的话请邮件联系。 鉴于很多同志不清楚用法，下面略加说明。 BoundedPowell函数的参数说明。 第一个参数，待优化函数的指针，这个函数必须定义成_stdcall类型，其本身带有三个参数，第一个参数是优化变量个数，第二个参数是优化变量数组，第三个参数就是函数的值，由于使用Fortran的关系，这几个参数都必须传递地址。 举例说，假定待优化函数为f(x)=x^2，那么只有一个参数，这个函数应该定义为 void _stdcall objfun(int *n, double *para, double *f) { *f=para[0]*para[0]; } 第二个参数int n，待优化计算变量的个数 第三个参数double *x, 待优化变量，这是一个数组，长度为n 第四个参数double *xlb,变量的下界，数组 第五个参数double *xub,变量的上界，数组 第五个参数double rhobeg,第六个参数double rhoend，这两个参数是Powell同志定义的两个半径，具体我也说不清楚，你想搞清楚的话可以搜Powerll同志的文献。一般你把rhobeg设置为1，rhoend设置为计算的精度，比如你希望精度为万分之一，rhoend就为1e-4。 第七个参数int maxfun，最大迭代次数 第八个参数，BOBYQA计算时的返回代码，具体如下 // //BOBYQA iflag返回值的含义： //IFLAG=1, Return from BOBYQA because NPT is not in the required interval //IFLAG=2, Return from BOBYQA because one of the differences XU(I)-XL(I)/6X is less than 2*RHOBEG. //IFLAG=3, Return from BOBYQA because FCN has been called MAXFUN times. //IFLAG=4, Return from BOBYQA because of much cancellation in a denominator. //IFLAG=5, Return from BOBYQA because a trust region step has failed to reduce Q. //
【文件预览】：
MSFRT40.DLL
Optimization.h
bobyqa
----Modified()
--------calfun.f(268B)
--------powell.h(213B)
--------Powell.mdp(40KB)
--------Powell.mak(6KB)
--------trsbox.f(13KB)
--------prelim.f(5KB)
--------main.f(2KB)
--------altmov.f(9KB)
--------bobyqa.f(5KB)
--------rescue.f(13KB)
--------update.f(2KB)
--------README.txt(3KB)
--------bobyqb.f(22KB)
--------Makefile(206B)
----Original()
--------calfun.f(268B)
--------trsbox.f(13KB)
--------prelim.f(5KB)
--------main.f(2KB)
--------altmov.f(9KB)
--------bobyqa.f(5KB)
--------rescue.f(13KB)
--------update.f(2KB)
--------README.txt(3KB)
--------SampleOutput.txt(9KB)
--------bobyqb.f(22KB)
--------Makefile(206B)
Powell.cpp
Powell.lib
Powell.h