拉格朗日函数优化

article/2025/9/20 0:00:17

等式约束最优化

可以写为:

引入拉格朗日乘子(\lambda _i\ne 0)把问题转换成拉格朗日函数


L(x,\lambda)=\left[  f(x) +\sum_{i=1}^{n}\lambda_ih_i(x) \right]

因为对于任何可行解,有h_i(x)=0,所以有

f(x)=L(x,\lambda ),也就是,minf(x)=minL(x,\lambda )

 

求解。对L(x,\lambda)分别求x,\lambda的偏导数为零,得到方程组求解极值点,然后从极值点挑出最值点。

 

x的偏导为零,使得目标函数和约束函数的法向量共线(梯度共线)。为什么梯度共线能求到极值?

 

绿线标出的是约束g(x,y)=c的点的轨迹。蓝线是f(x,y)的等高线,箭头是各个点的梯度。

 

从图上可以看到,蓝线(f(x,y)=d_1)与绿线相交,意味着肯定还存在其它的等高线(f(x,y)=d_2)在该条等高线的内部或者外部,使得新的等高线与目标函数的交点的值更大或者更小。所以当取到极值时,蓝线与绿线相切,而切点的梯度共线。

不等式约束最优化

可以写为:

引入拉格朗日乘子(\mu  _i\geq 0),定义上述问题的拉格朗日量(Lagrangian)如下

L(x,\mu )=\left[  f(x) +\sum_{i=1}^{m}\mu _ig_i(x) \right]

同时定义拉格朗日对偶函数(Lagrange dual function) 如下:

F(\mu )=inf_x L(x,\mu )=inf_x \left[ f(x)+\sum_{i=1}^{m}\mu _ig_i(x) \right]

一般情况下,L(x,\mu )是能取到最小值的,所以F(\mu )=inf_x L(x,\mu )=min_x L(x,\mu )

求解。当强对偶性成立时,通过KKT条件求解极值点,然后从极值点挑出最值点。

 

 

第一个条件使得目标函数和约束函数的法向量共线(梯度共线)。

 

 

最后一个条件称为互补松弛条件(Complementary Slackness Condition)。通过引入这个条件,增加了m个等式约束,使得等式的数量跟变量一样。


更一般地,我们把等式约束也加进来,优化问题可以写为:

 

KKT条件为

 

如果没有“不等式”约束条件,即 m=0,KKT条件就是拉格朗日乘数法中极值点满足的方程组。所以KKT条件是拉格朗日乘数法的推广,拉格朗日乘数法是KKT条件的特例。

注意到:

  1. KKT条件是强对偶性的必要条件,强对偶性下KKT条件才成立
  2. 一般仅用KKT条件来验证找到的解
  3. 当目标函数和约束都是线性时,优化问题为我们熟悉的线性规划(LP)
  4. 在线性规划里,\lambda ,\mu表示的是对应约束的影子价格

 

4,KKT与强对偶性

这里讨论只有不等式约束,并且强对偶性的情况

由强对偶性,有f(x)=max_\mu L(x,\mu ),也就是,min_xf(x)=min_x max_\mu  L(x,\mu )

 

原问题目标函数为f(x)=max_\mu L(x,\mu ),对应的对偶函数为F(\mu) =min_x L(x,\mu )

由强对偶性,我们有min_x f(x)= max_\mu F(\mu ),也就是min_x max_\mu  L(x,\mu ) = max_\mu min_x   L(x,\mu )

为什么强对偶下可以得到KKT条件?

首先看梯度共线。

x^*表示原问题取得最优值的解,也就是f(x^*)=min_x f(x)=min_x max_\mu  L(x,\mu )。由强对偶性,可得max_\mu min_x   L(x,\mu )=min_x max_\mu  L(x,\mu )=f(x^*)。也就是说,min_x   L(x,\mu )x=x^*处取得极值,也就是,偏导数为零。

然后看互补松弛条件。

 

x=x^*时,有max_\mu min_x   L(x,\mu )=max_\mu\left[  f(x^*) +\sum_{i=1}^{m}\mu _ig_i(x^*) \right]=f(x^*)+max_\mu\left[ \sum_{i=1}^{m}\mu _ig_i(x^*) \right]=f(x^*)

也就是,max_\mu\left[ \sum_{i=1}^{m}\mu _ig_i(x^*) \right]=0,也就是\mu _ig_i(x^*)=0,\forall i=1,...,m

5,拉格朗日函数与对偶性

对于不等式约束,

一般的,由\mu \geq 0,g_i(x)\leq 0,有f(x)\geq max_\mu L(x,\mu )。所以min_xf(x)\geq min_x max_\mu  L(x,\mu )

而根据拉格朗日对偶函数,有对偶问题为max_\mu F(\mu )=max_\mu min_x L(x,\mu )。由因为对偶问题是凸优化(Slater条件也满足),根据对偶问题的强对偶性,有max_\mu F(\mu )=max_\mu min_x L(x,\mu )= min_x max_\mu L(x,\mu )

 

所以,有min_xf(x)\geq min_x max_\mu  L(x,\mu )=max_\mu min_x L(x,\mu )=max_\mu F(\mu )。这就是原问题的对偶性。

当原问题有强对偶性时,由min_xf(x)=max_\mu F(\mu ),有min_xf(x)= min_x max_\mu  L(x,\mu )

 

6,参考

无约束最优化方法 - Orisun - 博客园 

拉格朗日乘子法和KKT条件 - Orisun - 博客园
【整理】深入理解拉格朗日乘子法(Lagrange Multiplier) 和KKT条件
KKT conditions深入理解拉格朗日乘子法(Lagrange Multiplier) 和KKT条件优化问题中的对偶性理论


http://chatgpt.dhexx.cn/article/arj7j9vB.shtml

相关文章

增广拉格朗日函数法

增广拉格朗日函数法 在二次罚函数法中,为了保证可行性,罚因子必须趋于正无穷。此时,子问题因条件数爆炸而难以求解。那么,是否可以通过对二次罚函数进行某种修正,使得对有限的罚因子,得到的毕竟最优解也是…

增广拉格朗日函数法(ALM)

增广拉格朗日函数法( Augmented Lagrangian method) 一、等式约束 考虑问题: min ⁡ x f ( x ) s . t . c i ( x ) 0 , i 1 , ⋯ , m . \begin{array}{ll} \min_x &f(x)\\ s.t. &c_i(x) 0, \quad i1,\cdots,m. \end{array} min…

广义拉格朗日函数的理解

1、拉格朗日函数: 求极值 求函数f(x,y,z)在条件φ(x,y,z)0下的极值 方法(步骤)是: 1.做拉格朗日函数Lf(x,y,z)λφ(x,y,z),λ称拉格朗日乘数 2.求L分别对x,y,z,λ求偏导,得方程组,求出驻点P(x,y,z) 如果这个实际问题的最大或…

各种拉格朗日函数

目录 一,拉格朗日函数 二,部分拉格朗日函数 三,增广拉格朗日函数 一,拉格朗日函数 以三元函数为例: 求f(x,y,z)的极值,s.t.g(x,y,z)0 拉格朗日函数L(x,y,z,a) f(x,y,z) a * g(x,y,z) 在极值点处一…

拉格朗日函数相关推导

优化问题(即高数中的求极值)可分为三类:无约束、等式约束、不等式约束。 对于无约束的优化问题:求导,令导数为零即可求解。 等式约束优化:(拉格朗日乘子法最开始就是求解等式约束优化的方法&…

HJ-浇花

链接:https://ac.nowcoder.com/acm/contest/322/M 来源:牛客网 时间限制:C/C 1秒,其他语言2秒 空间限制:C/C 32768K,其他语言65536K 64bit IO Format: %lld 题目描述 HJ养了很多花(9999999…

单片机——自动浇花系统

目录 1、图片 2、代码 1、图片 2、代码 #include <reg51.h> //调用单片机头文件 #define uchar unsigned char //无符号字符型 宏定义 变量范围0~255 #define uint unsigned int //无符号整型 宏定义 变量范围0~65535#include <intrins.h>sbi…

基于ESP32的开源定时浇花系统

基于ESP32的开源定时浇花系统 文章目录 基于ESP32的开源定时浇花系统前言一、软硬件环境二、模块连接图1.浇花功能说明2.Web界面展示 总结 前言 养了些许花花草草&#xff0c;需要按时浇灌&#xff0c;奈何总是要出差&#xff08;总想出去玩&#xff09;&#xff0c;又怕没人浇…

RISC-V MCU 自动浇花装置设计

第一部分 设计概述 1.设计目的 在家里养养盆花可以陶冶情操&#xff0c;丰富生活&#xff0c;因此&#xff0c;家庭盆栽如今被许多人喜爱。花草生长问题80%以上是由花儿浇灌问题引起的;好不容易种植几个月的花草&#xff0c;因为浇水不及时而缺水死亡。虽然市场上有卖盆花自…

【Proteus仿真】【51单片机】自动浇花灌溉系统设计

文章目录 一、功能简介二、软件设计三、实验现象联系作者 一、功能简介 本项目使用Proteus8仿真51单片机控制器&#xff0c;使用LCD1602液晶、按键、DS18B20、PCF8591 ADC、土壤湿度传感器、水位传感器、蜂鸣器模块等。 系统运行后&#xff0c;LCD1602显示传感器检测的温度、…

自动浇花系统的电路分析

目录 ​​​电源电路 复位电路和晶振 UCB转TTL ​​​电源电路 首先是电源电路&#xff0c;我们来看电源电路的原理图。 首先忽略D、D-&#xff0c;我们先来看电源本身&#xff0c;除了两头的接地和VCC&#xff08;这是USB接口的特性&#xff09;&#xff0c;我们观察到右…

51单片机wifi物联网的浇花控制系统设计

硬件设计 浇花控制系统采用51单片机与LCD液晶显示屏来实现&#xff0c;利用温度、湿度传感器及相应的显示、驱动执行机构、报警装置等实现温室作物生长环境控制器的设计。 硬件电路主要由51单片机最小系统lcd1602显示屏蜂鸣器报警模块设置按键微型水泵adc0832模数转换模…

毕业设计 远程智能浇花灌溉系统 - stm32 单片机 嵌入式 物联网

文章目录 0 前言1 简介2 主要器件3 实现效果4 设计原理5 部分关键代码5 最后 0 前言 &#x1f525; 这两年开始毕业设计和毕业答辩的要求和难度不断提升&#xff0c;传统的毕设题目缺少创新和亮点&#xff0c;往往达不到毕业答辩的要求&#xff0c;这两年不断有学弟学妹告诉学…

基于51单片机的智能浇花系统(可做毕设)

基于51单片机的智能浇花系统&#xff08;可做毕设&#xff09; 一、系统介绍二、仿真展示三、实物展示四、仿真过程五、代码1、ADC08322、LCD16023、按键4、水泵5、温湿度6、定时器7、main.c 五、完整工程 对LCD1602原理和操作掌握不好的可以看这篇&#xff1a; 快速掌握——LC…

毕设--自动浇花系统的设计

目录 毕设--自动浇花系统的设计1、作品实物图2、PCB原理图3、元器件清单4、土壤温湿度采集与显示5、硬件电路设计6、程序源码7、资料获取 毕设–自动浇花系统的设计 注&#xff1a;本毕设资源可在微信公众号&#xff1a;“Kevin的学习站” 中获取&#xff01; 1、作品实物图 2、…

单片机毕业设计 自动浇花灌溉系统设计

文章目录 1 简介2 主要器件3 实现效果4 设计原理5 关键代码 1 简介 Hi&#xff0c;大家好&#xff0c;今天向大家介绍一个学长做的单片机项目 基于单片机的自动浇花灌溉系统设计 大家可用于 课程设计 或 毕业设计 选题指导&#xff0c;项目分享&#xff1a; https://gitee…

基于android智能浇花装置,一种基于WiFi通讯远程控制的智能浇花装置的制作方法...

本实用新型涉及自动浇水装置技术领域&#xff0c;具体涉及一种基于WiFi 通讯远程控制的智能浇花装置。 背景技术&#xff1a; 花卉作为一种极具观赏性的植物&#xff0c;在室内和室外广泛种植。上班族工作紧张&#xff0c;空闲时间较少&#xff0c;长期不在家时需要花草被照料。…

基于51单片机的自动浇花系统

目录 一、项目需求 二、仿真 三、程序 四、资料清单 资料下载地址&#xff1a;基于51单片机的自动浇花系统 一、项目需求 1、自动检测土壤湿度、温度、光照强度&#xff1b; 2、土壤湿度过低驱动水泵进行浇花&#xff1b; 3、LCD1602显示当前土壤湿度、温度、光照强度…

【IoT开发】基于机智云物联网的智能浇花教程

摘要:随着近年来物联网技术的发展,相关的技术已经广泛应用于人们的生产和生活中。文章针对长期无人在家时花卉植物的浇水问题,设计了一套基于物联网的智能浇花系统。系统采用STM32与51增强型单片机作为控制器,esp8266物联网模块作为通信设备,底层采用MQTT协议,连接到物联网云平…

条件概率的链式法则

条件概率 条件概率是指事件A在事件B发生的条件下发生的概率。条件概率表示为&#xff1a;P&#xff08;A|B&#xff09;&#xff0c;读作“A在B发生的条件下发生的概率”。若只有两个事件A&#xff0c;B&#xff0c;那么&#xff0c; xx 事件发生时 yy 事件发生的概率: P(yy|x…