机械设计s优化设计程序.docxVIP

优化设计与应用Matlab编程作业1、试用牛顿法求函数的极小点。解：牛顿法一维搜索函数文件minNewton如下：function [x,minf]=minNewton(f,x0,eps)%目标函数：f;%初始点：x0;%精度：eps;%目标函数取最小值时的自变量：x;%目标函数的最小值：minf;format long if nargin==2 eps=1.0e-6; enddf=diff(f);%求目标函数的一阶导数d2f=diff(df);%求目标函数的二阶导数k=0;tol=1;while toleps dfx=subs(df,findsym(df),x0);%目标函数的一阶导数值 d2fx=subs(d2f,findsym(d2f),x0);%目标函数的二阶导数值 x1=x0-dfx/d2fx; k=k+1; tol=abs(dfx); x0=x1;endx=x1;minf=subs(f,findsym(f),x);format short;命令行输入：syms t;f=t^4-4*t^3-6*t^2-16*t+4;[x,fx]=minNewton(f,3)M函数文件的运行结果如下：x =4fx = -1562、利用梯度法求目标函数的极小值，设初始点为，收敛精度。解：梯度法函数文件minFD如下：function [x,minf] = minFD(f,x0,var,eps