优化设计习习题

1、1. 设函数，将写成矩阵形式，并求其梯度矢量和Hesse矩阵，并证明该函数为凸函数。矩阵形式：式中，。在点处的梯度矢量为 ，式中，。因为 ，所以 原函数的Hessian矩阵为因为 原函数的定义域是实数集，属于非空凸集，在其定义域内，对于任意自变量都有，即原函数Hessian矩阵的各阶顺序主子式均大于零，可以说明Hessian矩阵是正定矩阵，所以原函数为凸函数。2. 约束优化问题的数学模型为 用作图法求该问题的极小点，并验证该点满足Kuhn Tucker 条件。(1) 该数学模型的可行域如Error! Reference source not found.阴影（包括边界）所示以点（2,2）为圆心

2、以R为半径做圆， 的值即为的值，可见其最小值产生于与可行域边缘相切的圆，其切点即为最小值点。图 1 第二题图 (2) 因为，所以极值点是可行点；又因为是该模型的紧约束，而，可以看出，即极值点的目标函数梯度是所有紧约束梯度的线性组合，所以该点满足Kuhn Tucker 条件。3. 对于极小化，而受限于约束的优化问题，写出其内点和外点罚函数表达式，并说明内点罚函数法和外点罚函数法的特点。(1) 内点罚函数表达式为或，罚因子为递减正数列。(2) 外点罚函数表达式为，罚因子为递增正数列。(3) 内点罚函数的特点：初始点可行；迭代在可行域内部；最终点可行；只含不等式约束。(4) 外点罚函数的特点：初始点

3、任意选取；迭代在可行域外部；最终点不可行；可含不等式与等式约束。4. 在长为350cm、宽为260cm的长方形不锈钢板的四角，各剪去一个小正方形，做成一个无盖的储水箱，试确定正方形的边长，使储水箱的容积最大。建立优化设计模型，并用法求解（迭代3次）。1) 设剪去的正方形边长为x，如Error! Reference source not found.所示图 2 第四题图储水箱体积优化模型为2) 首先用进退法确定单峰区间：a) 取初始步长h=30，置初始值u3=0，f3=f(u3)=0，置k=0；b) 置u=u3+h=30，f=f（u）= -1740000，k=k+1=1；c) 因为f<f3

4、，所以置u2=u3=0，f2=f3=0，u3=u=30，h=2h=60，k=k+1=2；d) 置u=u3+h=90，f=f（u）= -1224000，k=k+1=3；e) 因为f>f3, k1，所以置u1=u2=0，f1=f2=0,u2=u3=-1740000,u3=u=90,f3=f=-1224000。f) 令a=minu1，u3=0，b= maxu1，u3=90，所以该函数的一个单峰区间为0,90。3) 法求解：a) 置初始搜索区间0,90，则a=0，b=90，左右试探点al=a+(b-a)= ，a2=a+(b-a)= ，相应的函数值jl= j(al)= +006， j2= j(a2

5、)= +006；b) 因为jl>j2，所以a=al=，al=a2=，a2= a+(b-a)= 。迭代第一次完成。c) jl= j(al)= +006， j2= j(a2)= +006，因为jl<j2，置b=a2=，a2=a1=， j2j1=+006，并计算al=a+(b-a)= 。第二次迭代完成。d) jl=j(al) = +006，j2= j(a2)= +006，因为jl>j2，所以a=al=，al=a2=，a2= a+(b-a)= 。第三次次迭代完成。e) 三次迭代后，a1=，a2=，jl=j(al) = +006，j2= j(a2)= +006，因为jl>j2，所

6、以置a*=a2=，a*即为问题的解。所以用法迭代三次的最优解为x*=，优化模型最小值为f（x*）=+006，即正方形的边长为52.5222cm时，储水箱的容积最大，最大为1993600cm2。5. 用最速下降法求下列无约束优化问题极小点： Min f(X)=x12+4x22+60设初始点取为X(0)=2, 2T，迭代3次。1) 初始点X(0)=2, 2T。2) ，令。3) 求解一维问题，当时取得最小值，所以令，第一次迭代完成。4) 令。5) 求解一维问题当时取得最小值，所以令，第二次迭代完成。 6) 令。7) 求解一维问题当时取得最小值，所以令，第三次迭代完成。6、对下列双目标约束优化设计问题

7、分别在设计空间和目标空间中用图形表示可行域D和非劣解集F。图3 设计空间图4 目标空间7. 简述遗传算法、拓扑优化算法的思想。遗传算法是模拟生物在自然环境中的遗传和进化过程而形成的一种自适应全局优化概率搜索算法，最早由美国密歇根大学的Holland教授提出，起源于20世纪60年代对于自然和人工自适应系统的研究。遗传算法是从代表问题可能潜在的解集的一个种群开始的，而一个种群则由经过基因编码的一定数目的个体组成。每个个体实际上是染色体带有特征的实体。染色体作为遗传物质的主要载体，即多个基因的集合，其内部表现（即基因型）是某种基因组合，它决定了个体的形状的外部表现。因此，在一开始需要实现从表现型到基

8、因型的映射即编码工作。由于仿照基因编码的工作很复杂，我们往往进行简化，如二进制编码，初代种群产生之后，按照适者生存和优胜劣汰的原理，逐代演化产生出越来越好的近似解，在每一代，根据问题域中个体的适应度大小选择个体，并借助于自然遗传学的遗传算子进行组合交叉和变异，产生出代表新的解集的种群。这个过程将导致种群像自然进化一样的后生代种群比前代更加适应于环境，末代种群中的最优个体经过解码，可以作为问题近似最优解。拓扑优化以材料分布为优化对象，通过拓扑优化，可以在均匀分布材料的设计空间中找到最佳的分布方案。具体来说，是在一个给定的空间区域里，依据已知的负载或支撑等约束条件，解决材料的分布问题，从而使结构的

9、刚度达到最大或使输出位移、应力等目标达到规定要求的一种结构设计方法，是有限元分析和优化方法有机结合的方法。拓扑优化的研究领域主要分为连续体拓扑优化和离散结构拓扑优化。不论哪个领域，都要依赖于有限元方法。连续体拓扑优化是把优化空间的材料离散成有限个单元（壳单元或者体单元），离散结构拓扑优化是在设计空间内建立一个由有限个梁单元组成的基结构，然后根据算法确定设计空间内单元的去留，保留下来的单元即构成最终的拓扑方案，从而实现拓扑优化。拓扑优化的建模方法有变密度法，水平集法，独立映射法，计划法等。求解方法即算法主要有：QC法（优化准则法），MMA法（移动渐进法），SLP（序列线性规划法），SQP（序列二

10、次规划法）。8. 依据自己工作，提出优化设计问题，并写出优化设计数学模型，编写优化程序，求出优化结果并分析。#include<>#include<>#include<>#include<>#include<>#define FILE_PATH "","r" /*/数据文件名#define N_COLONY 500 #define CITY 500 int xColony=80; int xCity; double edgeSpeed=5000; double probab1=; double p

11、robab2=; ;long NOCHANGE=20000; long maxGen=80000; int colonyN_COLONYCITY;double cityXYCITY2;double city_disCITYCITY;double dis_pN_COLONY;double sumbest,sumTemp;double speed;int tempCITY,ibest;clock_t timeStart,timeNow,timeTemp;long GenNum,Ni;void main() register int C1,j,k,pos_C,pos_C1; int k1,k2,l1

12、,l2,pos_flag; register double disChange; static i=0; timeStart=timeNow=timeTemp=clock(); init(); for(;) for(j=0;j<xCity;j+)tempj=colonyij; disChange=0;pos_flag=0; pos_C=rand()%xCity; for(;) if(rand()/<probab1) /内变异算子 do pos_C1=rand()%xCity; while (pos_C1=pos_C); C1=colonyipos_C1; else do j=ran

13、d()%xColony;while(j=i); k=position(colonyj,temppos_C); C1=colonyj(k+1)%xCity; pos_C1=position(temp,C1); if(speed>edgeSpeed && pos_C1<pos_C+2)break; / if(pos_C+1)%xCity=pos_C1 | (pos_C-1+xCity)%xCity=pos_C1 )break; k1=temppos_C; k2=temp(pos_C+1)%xCity; l1=temppos_C1; l2=temp(pos_C1+1)%x

14、City; disChange+=city_disk1l1+city_disk2l2-city_disk1k2-city_disl1l2; invert(pos_C,pos_C1); pos_flag+; if(pos_flag>xCity-1)break; / pos_C+; if(pos_C>=xCity)pos_C=0; if( speed<edgeSpeed && disChange<0) dis_pi+=disChange; disChange=0; tempTest(i); if( speed>=edgeSpeed &&

15、 disChange<0 ) dis_pi+=disChange; disChange=0; tempTest(i); i+; if(i>=xColony) Ni+; GenNum+;i=0; probab1=probab1*(1-GenNum*maxGen); if( speed<edgeSpeed && (rand()/<probab2) ) mapped(); if(GenNum>=maxGen | Ni>=NOCHANGE ) LastCP(); printBest(GenNum,Ni); exit(1); else for(j=k1

16、;j<k2+xCity;j+) t1=j%xCity; t2=(j+1)%xCity; temp_dis+=city_distmpt1tmpt2;return temp_dis;_课后习题：1：用黄金分割法编程求解函数由点出发，沿方向的极小值点。解：出发，沿方向，可得，带入原函数得：，容易确定其搜索区间为。黄金分割法C程序代码如下。Double a,b,e,k,u,a1,a2,f1,f2,x1jieguo,x2jieguo;e=;k=(sprt(5)-1)/2;a=0;b=;u=b-a;while(u>e) a1=a+(1-k)*(b-a);a2=a+k*(b-a); f1=*a1*a1+; f2=*a2*a2+; if(f1<f2) b=a2;else a=a2;u=b-a;x1Jieguo=(b+a)/2;x2Jieguo=*x1jieguo+;结果显示。2：用最速下降法求些列无约束优化问题极小点。，设初始点取为，迭代三次。解：依题意以及最速下降法可知：，其中。由题意可知，代入原函数可得，是函数达到极小值则满足：，即满足。根据上式，从出发求得，所以；，所以，所以，更具前面计算可得：3：证明DFP法具有二次终止性。考虑函数，H为n阶正定矩阵，由DPF方法可以得出，其中，依题意可知只需证明对于如下软式子成立即可：，且当迭代n-1次时，即m=n-1时，有Hn