中科大《优化设计》课程大作业之约束优化实验报告.docx

中科大《优化设计》课程大作业之约束优化实验报告
中科大《优化设计》课程大作业之约束优化实验报告
中科大《优化设计》课程大作业之约束优化实验报告
约束优化设计实验报告
力学系
型号：联想ｙ4７0
ＣPU：i5－24５0Ｍ
内存:2ＧＢ
系统:win７-６4位
求解问题：
如上是以下三个约束方法共同需要求解的问题，预估结果:在(x1，ｘ２，x3）≈(２3,１3,１2)点附近存在极值。其中,每个方法对应的初始条件分别为:
（1）随机试验法
设计变量范围:
随机试验点数：N=10０0
精度:eｐｓ=０。001
（2）随机方向法
初始点：ｘ０=（25，15,５）
初始步长：a0＝
精度:ｅｐｓ=０。0０1
(3）线性规划单纯形法
中科大《优化设计》课程大作业之约束优化实验报告
中科大《优化设计》课程大作业之约束优化实验报告
中科大《优化设计》课程大作业之约束优化实验报告
初始复合形：X＝［2０ 23　2５ 3０;10　13 1５ ２０;1０ 9 ５ 0］
顶点个数:n＝4
精度：eｐs=
计算结果：
约束方法
所需时间
迭代次数
极值点
极值
随机试验法

0
（，１2。６3，１2。06）
3４４5。５0
随机方向法
０。6５
４９
(２2。６7，，11。９8）
34４１。9９
线性规划单纯形法
6。２３
5
（2２。59,，1２。１2)
3445。6１
程序说明：主程序为mａin,运行maiｎ后按提示即可得到相应约束方法的求解结果。
程序如下:
主程序
cｌｅar;
glｏbaｌ kk；
ｋk＝0；
disp(＇’）;
ｄisp('2。随机方向法＇）;
disp（'3。线性规划单纯形法’);
wｈiｌe　1
n０=iｎｐｕt('请输入上面所想选择约束优化方法的编号(1、２、3）：’）;
ｉｆ　n0==１||ｎ0==２|｜ｎ0==3
　 　ｂrｅak；
ｅnd
　disp（'此次输入无效．')；
eｎｄ　
ｄisp(’　’)；
ｄiｓp（’～～～～~~~～~～~～～～～～～~~~～~～~~～~~~~～~～~～～～～～～~~～～～～’）;
[ｘｘ，yy]=ｆminｓ(n0)；
fprｉnｔf(’迭代次数为:　%8。0f\n', ｋk）；
ｄisp（'所求极值点的坐标向量为：’)；
中科大《优化设计》课程大作业之约束优化实验报告
中科大《优化设计》课程大作业之约束优化实验报告
中科大《优化设计》课程大作业之约束优化实验报告
ｆｐrintf（’ 　 　 　 ％１6。5ｆ\n＇, ｘx)；
fｐrｉnｔf（’所求函数的极值为： 　%１6．５f\n', ｙy）；
调用函数
ｆunctioｎ [xｘ，yｙ］=fmｉns(n０）
iｆ n0＝＝1
　 　ｔｉｃ;［xx,yｙ]=suiｊｉsｈiyan()；ｔoc；
elsｅｉf ｎ0==2
　 tic；［xｘ,ｙy]=suｉjiｆangxiａng();toc；
eｌｓeｉｆ n0==3
　　 tｉｃ;[xx,yy］＝daｎ