【最优化方法】无约束非线性函数

2023-11-02 11:58

本文主要是介绍【最优化方法】无约束非线性函数,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

向量投影证明

b 1 = d ∗ a 1 ∣ a 1 ∣ d = ∣ a 2 ∣ ∗ c o s θ c o s = ( a 2 , a 1 ) ∣ a 2 ∣ ∣ a 1 ∣ b_1 = d * \frac{a_1}{|a_1|} \\ d = |a_2| * cos\theta \\ cos = \frac{(a_2,a_1)}{|a_2||a_1|} b1=da1a1d=a2cosθcos=a2a1(a2,a1)
由上面3个式子最终得到
b 1 = ( a 2 , a 1 ) ( a 1 , a 1 ) a 1 b1 = \frac{(a_2,a_1)}{(a_1,a_1) } a_1 b1=(a1,a1)(a2,a1)a1

施密特正交化


我们想要构造正交的基坐标系,我们希望各个坐标最好是能够互相正交的,我们让
b 1 = a 1 b_1 = a_1 b1=a1
另一个坐标应该是 b 2 b_2 b2,这个坐标可以根据向量的运算得到
b 2 = a 2 − b 1 b 1 = a 2 , b 1 ( b 1 , b 1 ) b 1 b_2 = a_2 - b_1 \\ b_1 = \frac{a_2,b_1}{(b_1,b_1)} b1 b2=a2b1b1=(b1,b1)a2,b1b1
拓展到高维度也是同样道理,具体可以参考施密特正交化高维度的方法。

最速下降法

二次型

二次型是一个从向量到标量的函数:
f ( x ) = 1 2 x T A x − b x + c f(x) = \frac{1}{2} x^{T} A x - bx + c f(x)=21xTAxbx+c
它的导数是:
A x − b = 0 A x = b Ax-b = 0 \\ Ax = b Axb=0Ax=b
我们可以将计算 A x = b Ax= b Ax=b转为计算二次型的最值问题。
在这里插入图片描述
上图为不同的 A 对 f(x) 的图像的影响。(a)正定矩阵的二次型;(b)负定矩阵的二次型;©奇异矩 阵和非正定矩阵的二次型;(d)不定矩阵的二次型:此时解是一个鞍点,梯度法和 CG 均无法处理该问题。

下降方向

最速下降法也就是梯度法,其中梯度是函数值变化最大的方向。最终通过求导的方式来确定步长。
d = − ∇ f ( x ) x i + 1 = x i + α d d = - \nabla f(\bold x) \\ x_{i+1} = x_i + \alpha d \\ d=f(x)xi+1=xi+αd

证明垂直和最佳步长

现在我们已经确定了下降方向,接下来我们要确定下降的步长 α \alpha α, 步长我们我们使用精准的公式法进行计算。我们令函数为关于
ψ ( α ) = f ( x i + 1 ) ∇ ψ ( α ) = ∂ f ( x i + 1 ) ∂ α = ∂ f ( x i + 1 ) ∂ x i + 1 ∂ x i + 1 ∂ α = ( A x i + 1 − b ) d i = d i + 1 d i = 0 \psi(\alpha) = f(x_{i+1}) \\ \\ \nabla \psi(\alpha)=\frac{\partial f(x_{i+1})}{\partial \alpha} = \frac{\partial f(x_{i+1})}{\partial x_{i+1}} \frac{\partial x_{i+1}}{\partial \alpha} = (A\bold x_{i+1} -b) \bold d_i = d_{i+1} d_i = 0 ψ(α)=f(xi+1)ψ(α)=αf(xi+1)=xi+1f(xi+1)αxi+1=(Axi+1b)di=di+1di=0
所以证明了两次的搜索方向是互相垂直的。下面将计算出 α \alpha α的具体数值:
d i + 1 d i = ( A x i + 1 − b ) T d = ( A ( x i + α d ) − b ) T d i = ( A x i + α A d i − b ) T d i = ( A x i − b + α A d i ) T d i = ( − d i + α A d i ) T d i = − d i T d i + α d i T A d i = 0 \bold d_{i+1} \bold d_i \\ = (A\bold x_{i+1} - b)^T \bold d \\ = (A(x_i + \alpha d) - b)^T d_i \\ = (Ax_i + \alpha A d_i -b)^T d_i \\ = (Ax_i -b + \alpha A d_i)^Td_i \\ = ( -d_i + \alpha A d_i)^T d_i \\= -d_i ^T d_i + \alpha d^T_i A d_{i} =0 di+1di=(Axi+1b)Td=(A(xi+αd)b)Tdi=(Axi+αAdib)Tdi=(Axib+αAdi)Tdi=(di+αAdi)Tdi=diTdi+αdiTAdi=0
化简得到

α = d i T d i d i T A d i \alpha = \frac{d_i^Td_i}{d_i^TAd_i} α=diTAdidiTdi

共轭梯度下降

###

这篇关于【最优化方法】无约束非线性函数的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/330431

相关文章

python获取指定名字的程序的文件路径的两种方法

《python获取指定名字的程序的文件路径的两种方法》本文主要介绍了python获取指定名字的程序的文件路径的两种方法,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要... 最近在做项目,需要用到给定一个程序名字就可以自动获取到这个程序在Windows系统下的绝对路径,以下

JavaScript中的高级调试方法全攻略指南

《JavaScript中的高级调试方法全攻略指南》什么是高级JavaScript调试技巧,它比console.log有何优势,如何使用断点调试定位问题,通过本文,我们将深入解答这些问题,带您从理论到实... 目录观点与案例结合观点1观点2观点3观点4观点5高级调试技巧详解实战案例断点调试:定位变量错误性能分

Python中 try / except / else / finally 异常处理方法详解

《Python中try/except/else/finally异常处理方法详解》:本文主要介绍Python中try/except/else/finally异常处理方法的相关资料,涵... 目录1. 基本结构2. 各部分的作用tryexceptelsefinally3. 执行流程总结4. 常见用法(1)多个e

C++统计函数执行时间的最佳实践

《C++统计函数执行时间的最佳实践》在软件开发过程中,性能分析是优化程序的重要环节,了解函数的执行时间分布对于识别性能瓶颈至关重要,本文将分享一个C++函数执行时间统计工具,希望对大家有所帮助... 目录前言工具特性核心设计1. 数据结构设计2. 单例模式管理器3. RAII自动计时使用方法基本用法高级用法

JavaScript中比较两个数组是否有相同元素(交集)的三种常用方法

《JavaScript中比较两个数组是否有相同元素(交集)的三种常用方法》:本文主要介绍JavaScript中比较两个数组是否有相同元素(交集)的三种常用方法,每种方法结合实例代码给大家介绍的非常... 目录引言:为什么"相等"判断如此重要?方法1:使用some()+includes()(适合小数组)方法2

504 Gateway Timeout网关超时的根源及完美解决方法

《504GatewayTimeout网关超时的根源及完美解决方法》在日常开发和运维过程中,504GatewayTimeout错误是常见的网络问题之一,尤其是在使用反向代理(如Nginx)或... 目录引言为什么会出现 504 错误?1. 探索 504 Gateway Timeout 错误的根源 1.1 后端

GO语言中函数命名返回值的使用

《GO语言中函数命名返回值的使用》在Go语言中,函数可以为其返回值指定名称,这被称为命名返回值或命名返回参数,这种特性可以使代码更清晰,特别是在返回多个值时,感兴趣的可以了解一下... 目录基本语法函数命名返回特点代码示例命名特点基本语法func functionName(parameters) (nam

Python Counter 函数使用案例

《PythonCounter函数使用案例》Counter是collections模块中的一个类,专门用于对可迭代对象中的元素进行计数,接下来通过本文给大家介绍PythonCounter函数使用案例... 目录一、Counter函数概述二、基本使用案例(一)列表元素计数(二)字符串字符计数(三)元组计数三、C

MySQL 表空却 ibd 文件过大的问题及解决方法

《MySQL表空却ibd文件过大的问题及解决方法》本文给大家介绍MySQL表空却ibd文件过大的问题及解决方法,本文给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友参考... 目录一、问题背景:表空却 “吃满” 磁盘的怪事二、问题复现:一步步编程还原异常场景1. 准备测试源表与数据

python 线程池顺序执行的方法实现

《python线程池顺序执行的方法实现》在Python中,线程池默认是并发执行任务的,但若需要实现任务的顺序执行,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋... 目录方案一:强制单线程(伪顺序执行)方案二:按提交顺序获取结果方案三:任务间依赖控制方案四:队列顺序消