广西大学优化设计课程作业.docVIP

优化设计大作业 专 业： 指导教师 学 生： 日 期： 目录 第一部分 优化设计大作业 ……………………………………1 课程练习题之一 ………………………………………1 课程练习题之二 ………………………………………4 课程练习题之三 ………………………………………8 第二部分 优化设计课后作业 …………………………………13 机械优化的基本概念与数学模型……………………13 优化设计的最优性条件………………………………14 无约束规划的解法 …………………………………22 线性规划与二次规划 ………………………………23 第七章 机械结构优化 ………………………………………25 第八章 结构优化的灵敏度分析………………………………27 第一部分 优化设计大作业 课程练习题之一 3-2.求函数 f(X)= x1^4-5*x1^2*x2-2*x1*x2+4*x1^2+6*x2^2+4.5*x1-4*x2+5的极小值，初始点为X(0)=[-2，2]T，误差ε不大于0.001。 解.此问题为无约束非线性优化问题，用牛顿法进行迭代求解。 1.建立数学模型： Find x1 , x2 Min f=x1^4-5*x1^2*x2-2*x1*x2+4*x1^2+6*x2^2+4.5*x1-4*x2+5 初始点X(0)=[-2，2], ε≤0.001 2.运行结果： 迭代值 迭代次数 x1 x2 f 1 -0.50000000000000 -0.58333333333333 8.33333333333333 2 -0.56606217616580 0.37068221070812 3.00457834092293 3 -0.65287432092310 0.39898283419758 2.97857158805445 4 -0.65084893915893 0.40135860990393 2.97849714375511 5 -0.65083910731569 0.40135662524523 2.97849714338108 从中我们可以看出：通过牛顿法迭代5次可得出结果，当x1=-0.65083910731569，x2=0.40135662524523时，目标函数最优值fmin=2.97849714338108，且满足。 3.迭代曲线： 4.检验结果： 用matlab自带优化程序检验 程序为： x0=[-2,2]; [x,fval]=fminsearch(x(1)^4-5*x(1)^2*x(2)-2*x(1)*x(2)+4*x(1)^2+6*x(2)^2+4.5*x(1)-4*x(2)+5,x0) x =-0.65086658687466 0.40137142333985 fval =2.97849714628600 经检验用牛顿法进行迭代优化结果是正确的，优化结果达到精度要求，ε≤0.001。 5.讨论： （1）由以上迭代曲线可知，牛顿法迭代收敛速度很快，本优化经过迭代3次后目标函数值趋于平稳。也可采用黄金分割法，变尺度法等其他方法优化。由于本题比较简单，不必采用变尺度法来优化。 （2）采用matlab编程解决了求导和计算海森阵比较复杂的难题，编程简单方便。 5.题3-2的matlab源程序： function ZY32 format long syms x1 x2 %定义符号变量x1,x2 f=x1^4-2*x1^2*x2-5*x1*x2+6*x1^2+7*x2^2+4.5*x1-4*x2+5 %定义函数f df=[diff(f,x1);diff(f,x2)] %diff(f,x1)用于对函数f中变量x1求偏导 %diff(f,x2)用于对函数f中变量x2求偏导 f1=diff(f,x1,2); %diff(f,x1,2)用于对函数f中变量x1求而二次偏导 f2=diff(diff(f,x1),x2); %diff(diff(f,x1),x2)用于对函数f中变量下x1,x2求偏导 f3=diff(diff(f,x2),x1); f4=diff(f,x2,2); %diff(f,x2,2)用于对函数f中变量x2求而二次偏导 ddf=[f1,f2;f3,f4]; %求函数f的海森阵 x1=-