防灾数学分享 http://blog.sciencenet.cn/u/fzmath 防灾科技学院数学教研室

博文

无约束最优化问题求解

已有 5864 次阅读 2017-9-8 22:29 |个人分类:教学辅导|系统分类:科研笔记


无约束最优化问题是最简单的一类最优化问题,其一般数学描述为

              $\min\limits_{\textbfsymbol{x}\in \mathbb{R}^n} f(\textbfsymbol{x})$

其中 $\textbfsymbol{x}=[x_1,x_2,\cdots,x_n]^T$ 称为优化变量, $f(\cdot)$ 函数称为目标函数。

1、解析法和图解法

令梯度 $\nabla f=\textbfsymbol{0}$ 得到驻点,即为极值可疑点。  图解法只是适合低维的一元、二元和三元函数。


例1 求一元函数最值

$f(t)=\text{e}^{-3t}\sin(4t+2)+4\text{e}^{-0.5t}\cos(2t)-0.5$


解析解方法

>> syms t; y = exp(-3*t)*sin(4*t+2)+4*exp(-0.5*t)*cos(2*t)-0.5;

>> y1 = diff(y,t);t0 = solve(y1),ezplot(y,[0,4])

警告: Cannot solve symbolically. Returning a numeric approximation instead.

> In solve (line 303)


t0 =


-0.22133627468652136234203676475842


>> y1 = diff(y,t);t0 = solve(y1),ezplot(y,[-4,4])

警告: Cannot solve symbolically. Returning a numeric approximation instead.

> In solve (line 303)


t0 =


-0.22133627468652136234203676475842


>> y1 = diff(y,t);t0 = solve(y1),ezplot(y,[-1,4])

警告: Cannot solve symbolically. Returning a numeric approximation instead.

> In solve (line 303)


t0 =


-0.22133627468652136234203676475842


>> y2 = diff(y1);b = subs(y2,t,t0)


b =


-51.712061718150102502340461429533

这说明二阶导数 b<0,

t0 =-0.22133 为极大值点







https://wap.sciencenet.cn/blog-292361-1074984.html

上一篇:LaTex编辑器WinEDT10.2与TeXLive2016套装配合
下一篇:MATLAB中的两个常用数值优化函数fminsearch()和fminunc()
收藏 IP: 124.238.132.*| 热度|

1 杨正瓴

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-5-9 10:50

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部