惩罚函数程序（Penalty function program）.docVIP

惩罚函数程序（Penalty function program） #包含iostream #包括 H / /记时间必须要加的头文件窗口。 使用名称空间； 双func0（双*） { 返回（- 3 * [ 1 ] + 2 * 0 [ 2 - 0.8 ] /（（2 * A [ 1 ] + A [ 2 ]）+（- 4 * * 2 [ 0 ] * 2 [ * ] [ 1 ] - [ 2 ] /（7 * 0 [ 0 ] + 3 * A [ 1 ] ]； } 双func1（双*，双R） { 返回（- 3 * [ 1 ] + 2 * 0 [ 2 - 0.8 ] /（（2 * A [ 1 ] + A [ 2 ]）+（- 4 * * 2 [ 0 ] * 2 [ * ] [ 1 ] - [ 2 ] /（7 * 0 [ 0 ] + 3 * A [ 1 ] ]— R *（log（- [ 0 ] + [ 1 ] - A - 2 - 1））+ log（-（- [ 0 ] + [ 1 ] - 2 1 +）+日志（-（12 * 34.8 0 + 5 * 1 + 12 * 2 ]）） 日志（-（12 * 2 [ 0 * 12 ] + 7 * 1 [ 2 - 29.1））+日志（-（- 6 * 0 ] + 1 2 + 4.1））+日志（A [ 0 ] ] [ 1 ] ]； } 双func11（双*）/ / 8个约束条件 { 返回一个[ 0 ] + [ 1 ] - A [ 2 ] - 1； } 双func22（双*） { 返回[ 0 ] + A [ 1 ] - A [ 2 ] + 1； } 双func33（双*） { 返回12 * [ 0 ] + 5 * * [ 1 ] + 12 * A [ 2 ] - 34.8； } 双func44（双*） { 返回12 * [ 0 ] + 12 * * [ 1 ] + 7 * A [ 2 ] - 29.1； } 双func55（双*） { 返回- 6 * 4.1 [ 0 ] + [ 1 ] + [ 2 ] +； } 双func66（双*） { 返回[ 0 ]； } 双func77（双*） { 返回[ 1 ]； } 双func88（双*） { 返回[ 2 ]； } 国际xunhuancishu = 0； 双搜索（双AA，BB双，双R，int n）/参数：搜索方向，初始点，影响系数/用于求解最优步长 { 双L = 1e-2； 双A1 = 0，B1 = 0，C1 = 0； C1＝+1； 双*，*，*； =（双*）malloc（sizeof（双）* N）； B =（双*）malloc（sizeof（双）* N）； C =（双*）malloc（sizeof（双）* N）； 对于（int = i 0；i ；n；+ +） *（C++）； / /，（func11（C） 0 | | func22（C） 0 | | func33（C） 0 | | func44（C） 0）/先从负方向进行搜索，保证从负方向进行搜索不会越界 / { / / cout 777 endl； // 2； //； //（int = i 0；i ；n；+ +） //； / } / / L = L /再从正方向进行搜索； //； //（int = i 0；i ；n；+ +） //； 而（func11（C） 0 | | func22（C） 0 | | func33（C） 0 | | func44（C） 0 | | func55（C） 0 | | func66（C） 0 | | func77（C） 0 | | func88（C） 0）/再从正方向进行搜索，防止跑出范围 { l＝2； C1＝+1； 对于（int = i 0；i ；n；+ +） *（C++）； } 对于（int = i 0；i ；n；+ +） *（+我）= BB [我]；//确定L以后，就进行初始点的计算，此时L是一个正值 如果（晶圆厂（func1（A，R）- func1（C，R）） 1e-6） { 如果（func1（A，R） func1（C，R））/负方向进行搜索 { L = L； C1＝+1； 对于（int = i 0；i ；n；+ +） *（C++）； 而（func11（C） 0 | | func22（C） 0 | | func33（C） 0 | | func44（C） 0 | | func55（C） 0 | | func66（C） 0 | | func77（C） 0 | | func88（C） 0）/用于防止搜索点跑出约束范围 { L = 2；//这个地方尽量是除以2，具体原因还不清楚 C1＝+1； 对于（int = i 0；i ；n；+ +） *（C++）； } B1 = C1 + L； 对于（int