优化设计习题.docx
- 文档编号:26431383
- 上传时间:2023-06-19
- 格式:DOCX
- 页数:14
- 大小:182.29KB
优化设计习题.docx
《优化设计习题.docx》由会员分享,可在线阅读,更多相关《优化设计习题.docx(14页珍藏版)》请在冰豆网上搜索。
优化设计习题
1.设函数f(x)=2x;+4x“2+4x2+4无+8x2+50,将/(x)写成矩阵形式,并求其梯度矢量和Hesse矩阵,并证明该函数为凸函数。
矩阵形式:
/(X)=-XrAX+B*X+C
2
X.
afl
44
「b「
'4'
式中X二
1
A=
ii
12
=
B=
i
=
C=50o
_X2-
_a21
322_
48
8
在点X(k)处的梯度矢量为
V/-(X X. atl 44 「b「 4 式中x= 1 A= ii 12 = B= i = 0 _X2- _321 322_ 48 8 因为几=4,f: iX2 =4, 总= : 4,f: iX2=8 所以原函数的Hessian矩阵为 z"44 HX= 48 因为原函数的定义域是实数集,属于非空凸集,在其定义域内,对于任意自变量都有Aj=4>0,亠="j=16>0,即原函数Hessian矩阵的各阶顺序 48| 主子式均大于零,可以说明Hessian矩阵是正定矩阵,所以原函数为凸函数<> 2.约束优化问题的数学模型为 minf(x)=(X)-2)2+(x2一2)2 s.t.g©)=彳+£-1SO ^2(x)=-x2<0 ^3(x)=-Xj<0 用作图法求该问题的极小点“,并验证该点满足血加-Tucker条 件。 (1)该数学模型的可行域如图1阴影(包括边界)所示 以点(2,2)为圆心以R为半径做圆,”的值即为/(劝的值,可见其最小 —U可行点; (R厅\ F二,、二的目标函数梯度是所有紧约束梯度的线性组合,所以该点满足KuhnI22丿 -Tucker条件° 3.对于极小化/(X),而受限于约束g“(x)no("=i,2,..-M)的优化问题,写出其内点和外点罚函数表达式,并说明内点罚函数法和外点罚函数法的特点。 (I)内点罚函数表达式为p (2))=/(x)+a£计刃或 p(x,川))=/(X)-川)》ln(g“(x)),罚因子严)为递减正数列。 (|>1 (2)外点罚函数表达式为"(X,严)卜/(X)+严£{max[O,-g“(X)]f,罚因子 W=1 7•⑹为递增正数列。 (3)内点罚函数的特点: [ 初始点可行;迭代在可行域内部;最终点可行;只含不等式约束。 (4)外点罚函数的特点: 初始点任意选取;迭代在可行域外部;最终点不可行;可含不等式与等式约束。 4.在长为350cm、宽为260cm的长方形不锈钢板的四角,各剪去一个小正方形,做成一个无盖的储水箱,试确定正方形的边长,使储水箱的容积最大。 建立优化设计模型,并用法求解(迭代3次)。 1)设剪去的正方形边长为X,如图2所示 350-2x 260—2x 图2第四题图 储水箱体积/(x)=x•(350-2x)-(260-2a)优化模型为 nin/(x)=-x•(350-2a)•(260—2x) s.t.0 2)首先用进退法确定单峰区间: a)取初始步长h=30,置初始值u3=0,f3=f(u3)=0,置k=0; b)置u=u3+h=30,f=f(u)=-1740000,k=k+l=l; c)因为f d). o>置u=u3+h=90,f=f(u)=-1224000,k=k+l=3; f)因为f>f3,kH1,所以置ul=u2=0, f1=f2=0,u2=u3=-l740000,u3=u=90,f3=f=-l224000。 g)令a=min{ul,u3)=0,b=max{ul,u3)=90,所以该函数的一个单峰区间为[0,90]o 3)法求解: a)置初始搜索区间[0,90],则a二0,b=90,左右试探点al=a+(b-a)=, a2=a+(b-a)=,相应的函数值1=(al)=+006,2=(a2)=+006; b)因为1>2,所以a=al=,al=a2=,a2=a+(b-a)=。 迭代第一次完成。 c)1=(al)=+006,2=(a2)=+006,因为1<2,置b=a2=, a2=al=,2=1=+006,并计算al=a+(b~a)=。 第二次迭代完成。 d)1=(al)=+006,2=(a2)=+006,因为1>2,所以a=al=, al=a2=,a2=a+(b~a)=。 第三次次迭代完成。 e)/ 三次迭代后,al=,a2=,1=(al)=+006,2=(a2)=+006, 因为1>2,所以置a*=a2=,a*即为问题的解。 所以用法迭代三次的最优解为x*二,优化模型最小值为f(x*)=+006,即正方形的边长为52.5222cm时,储水箱的容积最大,最大为1993600cm2<> 5.用最速下降法求下列无约束优化问题极小点: Minf(A)二打+4卅+60 设初始点取为/0)={2,2}T,迭代3次。 3) •n求解一维问题 17 min0(a)=/(X⑹+^/(0>)=(2-4a)2+4(2-16a)2+60,当6/=—时取得 /、设计空间 /x.<5行域 1- XI +X,<4 1- 0" X,<5 •]■ / /0" 1''/« >1254 6 设计空间非劣解 图3设计空间 F.(x) 图4目标空间 7.简述遗传算法、拓扑优化算法的思想。 遗传算法是模拟生物在自然环境中的遗传和进化过程而形成的一种自适应全局优化概率搜索算法,最早由美国密歇根大学的Holland教授提出,起源于20世纪60年代对于自然和人工自适应系统的研究。 遗传算法是从代表问题可能潜在的解集的一个种群开始的,而一个种群则由经过基因编码的一定数目的个体组成。 每个个体实际上是染色体带有特征的实体。 染色体作为遗传物质的主要载体,即多个基因的集合,其内部表现(即基因型)是某种基因组合,它决定了个体的形状的外部表现。 因此,在一开始需要实现从表现型到基因型的映射即编码工作。 由于仿照基因编码的工作很复杂,我们往往进行简化,如二进制编码,初代种群产生之后,按照适者生存和优胜劣汰的原理,逐代演化产生出越来越好的近似解,在每一代,根据问题域中个体的适应度大小选择个体,并借助于自然遗传学的遗传算子进行组合交叉和变异,产生出代表新的解集的种群。 这个过程将导致种群像自然进化一样的后生代种群比前代更加适应于环境,末代种群中的最优个体经过解码,可以作为问题近似最优解。 拓扑优化以材料分布为优化对象,通过拓扑优化,可以在均匀分布材料的设计空间中找到最佳的分布方案。 具体来说,是在一个给定的空间区域里,依据已知的负载或支撑等约束条件,解决材料的分布问题,从而使结构的刚度达到最大或使输出位移、应力等目标达到规定要求的一种结构设计方法,是有限元分析和优化方法有机结合的方法。 拓扑优化的研究领域主要分为连续体拓扑优化和离散结构拓扑优化。 不论哪个领域,都要依赖于有限元方法。 连续体拓扑优化是把优化空间的材料离散成有限个单元(壳单元或者体单元),离散结构拓扑优化是在设计空间内建立一个由有限个梁单元组成的基结构,然后根据算法确定设计空间内单元的去留,保留下来的单元即构成最终的拓扑方案,从而实现拓扑优化。 拓扑优化的建模方法有变密度法,水平集法,独立映射法,计划法等。 求解方法即算法主要有: QC法(优化准则法),MMA法(移动渐进法),SLP(序列线性规划法),SQP(序列二次规划法)。 8.依据自己工作,提出优化设计问题,并写出优化设计数学模型,编写优化程序,求出优化结果并分析。 ftincludeOftincludeO#include<>#include<>#include<>^defineFILEPATH//***〃数据文件名 ^defineNCOLONY500^defineCITY500 int xColony=80; int xCity; double edgeSpeed=5000; double •> probabl=; double probab2=;; long N0CHANGE-20000; long maxGen=80000; intcolony[NCOLONY][CITY];doublecityXY[CITY][2]; double citydisECITY][CITY]; double dis』[NCOLONY]; 讶 double sumbest,sumTemp; double speed; int temp[CITY],ibest; clockt tiineStartttimeNow,timeTemp; long GenNum,Ni; voidmain() registerintCl.j,k,posC,posCl; intkltk2,11,12tposflag; registerdoubledisChange; statici=0; timeStart=timeNow=timeTemp=c1ock();init(); for(;;) { for(j=0;j disChange=O;posflag=0; posC=rand()%xCity; for(: ;) { if((rand()/ { doposCl=rand()%xCity; while(posCl==posC): Cl=colony[i][posCl]; } else { doj=rand()%xCo1ony;whi1e(j==i); k=posilion(colony[j],temp[posC]); Cl=colony[j][(k+l)%xCity]: posCl=position(temp,Cl); } if(speed>edgeSpeed&&posCKposC+2)break; /////////////////////// if((posC+l)%xCity==posCl(posC-l+xCity)%xCity==posCl)break;kl=temp[posC]: k2=temp[(posC+l)%xCity]; 1l=temp[posCl];12=temp[(posCl+l)%xCity]; disChange+=citydis[kl][ll]+citydis[k2][12]-citydis[kl][k2]-citydistil][ 12]; invert(posC.posCl); posflag++; if(posflag>xCity-1)break;//////////// posC++; if(posC>=xCity)posC=0; if(speed { disp[i]+=disChange; disChange=0; tempTest(i); } ) if(speed>=edgeSpeed&&disChange<0) { disp[i]+=disChange; k disChange=O; tempTest(i); ) i++; if(i>=xColony) { Ni++;GenNum++;i=0; probabl=probab1*(1-GenNum*maxGen): » ■ if(speed { mapped(); } if(GenNurn>-maxGen||Ni>二NOCHANGE) { LastCPO; printBest(GenNuin,Ni);exit (1); else for(j=kl;j {tl=j%xCity;t2=(j+l)%xCity; tempdis+=citydis[tmp[tl]][tmp[t2]];) returntempdis; 课后习题: 1: 用黄金分割法编程求解函数/以)=10(召+冬-5)2+(召-兀2)2由点 X(0)=[3.86933.3797]出发,沿方向=[-0・6,0・5]卩的极小值点。 解: X(0)=[3.8693,1.3797]出发,沿=[-0・6,0・5]卩方向,可得 x2=-0.8333^+4.6041,带入原函数得: /(X)=3.6389x)2-9.1007x,+22.7651, 容易确定其搜索区间为[0.3.8693]o黄金分割法C程序代码如下。 Doublea,b,e,k,u,al,a2,f1,f2,xljieguo,x2jieguo; e=;k=(sprt(5)-l)/2; a=0;b=;u二b-a; while(u>e) { al=a+(l-k)*(b-a);a2=a+k**(b~a); fl=*al**al+; f2=*a2**a2+; if(fl b=a2; else a=a2; u=b~a; } xlJieguo=(b+a)/2;x2Jieguo=*x1jieguo+; 结果显示X⑴=[1.25,3.562475]。 2: 用最速下降法求些列无约束优化问题极小点。 Mi心X)=#+4£_160),设初始点取为X(0)=[2,2f,迭代三次。 解: 依题意以及最速下降法可知: ,其中/=-W(X⑹)。 由题意可知V/-(Xa))=[2xf,8xJ]r,代入原函数可得/(Xa+,))=f(X(k)-akVf(Xik}))=/(X⑷一/[2彳,迟丁),是函数达到极小值则 满足: 川X;])=2(屮)-)•(—2屮>)+8(垮〉_8a・『)・(-8护)=0,即满足 宀空')二16({广。 根据上式,从x⑹出发求得,a0=—=0.13077,所 2(x;*))2+128(x^))2130 以X")=[2,2]7+0.13007[4,16]7=[2.52308,4.08112]「;刃=0.12723,所以 X (2)=[2.5230&4.08112]7+0.12723(2x2.5230&8x4.08112]f=[2.16510.8.23505]7«2=0.11620,所以,更具前面计算可得: 乂⑶二[2.16510,8.23505]7+0.11620[2x2.1651058x8.2350习7=[2.66827,15.89035”3: 证明DFP法具有二次终止性。 考虑函数/(X)=+x"GX+刃X+C,H为n阶正定矩阵,由DPF方法可以得出 2 y汀=W(xz)-巧(X。 依题意可知只需证明对于心0丄…冲-1如下软式子 成立即可: 且当迭代nT次时,即m二nT时,有Hn=G】 H3j=Sj0 s'Hsj=00 证明: 当时,结论显然成立。 设结论对j成立,现证明i+1时结论亦成立。 注意到由精确一维搜索及归纳法假设,对于J\有 gj+Pjhgj^j+工(MWI&)Sj A-j+1 =g;+円+£y[sj=0+乞s: GSj=0 &.J+lj+1 由也=一4+血=%厲+禺+i及Gsj=兀,得 smGsj=-0+心: 汕+1兀=%品円=0 这就证明了定理中的第一组式子。 下证第二组式子,即 Hdj=j=0,l^..j+lo 由DFP公式立即可得而对于j $二兀=s: Gs]=0 yM+y=>^sj=smGsj=° 右rrMiSmj乩十i)jy[]Hmy 有比小=H^y)+-r—->z—-=H加)=s. S+iX+i刃+1门井i〉j+i 定理证毕。 由上证明可知,DFP拟牛顿法是共馳方向法。 若取%=1,则初始方向为负梯度方向,此时方法变成共龜梯度法。 DFP算法是一个在理论分析和实际应用中都起重要作用的算法。 4: 用坐标轮换法编程求解函数的极小点X•: /(X)=10(x,+x2-5)2+(x1-x2)2,由点x(°)=[o,or出发,迭代一次(或精度为io”,编程时使用)。 解: 根据坐标轮换法可得,X;“=Xf+ah,xy=x『+a;e2,依题意知,X;°)=[0,0]7,弓=[1,0]7,e2=[0jfo则,min/(Xj)=min(10(a: —5)2+a;2),求解得£=罟,X|⑴=[罟,0卩,则,min/(XI<2))=min(10(^-«*)2+(^_^)2)求解得罟,则皆=[晋挣,一次迭代完毕。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 优化 设计 习题
![提示](https://static.bdocx.com/images/bang_tan.gif)