机器学习之拉格朗日乘子法和 KKT

有约束的最优化问题

最优化问题一般是指对于某一个函数而言,求解在其指定作用域上的全局最小值问题,一般分为以下三种情况(备注:以下几种方式求出来的解都有可能是局部极小值,只有当函数是凸函数的时候,才可以得到全局最小值):

1、无约束问题:求解方式一般为梯度下降法、牛顿法、坐标轴下降法等;

                                                                                         \min f(x)

2、等式约束条件:求解方式一般为拉格朗日乘子法

                                                                              \begin{center} \min f(x) \\s.t.\ \ h_k(x)=0 \ \ \ k=1, 2, \dots, p \\F(x, \alpha )=f(x)+\sum_{k=1}^p\alpha_kh_k(x)\end{center}

3、不等式约束条件:求解方式一般为KKT条件

                                                       \begin{center} \min f(x) \\ s.t.\ \ \ \ \ h_k(x)=0 \ \ \ \ \ \ k=1, 2, \dots, p \\\ \ \ \ \ \ \ \ \ g_j(x) \le 0 \ \ \ \ \ \ j=1, 2, \dots, q \\L(x, \alpha, \beta)=f(x)+\sum_{k=1}^p\alpha_kh_k(x)+\sum_{j=1}^q\beta_jg_j(x) \end{center}

拉格朗日乘子法

拉格朗日乘子法就是当我们的优化函数存在等值约束的情况下的一种最优化求解方式;

其中参数α被称为拉格朗日乘子,要求α不等于0

\begin{center} \min f(x) \\s.t.\ \ h_k(x)=0, \ \ \ k=1, 2, \dots, p \\F(x, \alpha )=f(x)+\sum_{k=1}^{p}\alpha_kh_k(x),\ \ \ \alpha_k \neq 0\end{center}

假设现在有一个二维的优化问题

画出图像加深理解

数学证明可参考链接:https://wenku.baidu.com/view/ac56710e2e3f5727a5e962a7.html

对偶问题

\underset{\beta }{min}\ \underset{x}{min}\ L(x,\beta)=\underset{x}{min}\ \underset{\beta }{max}\ L(x,\beta)

在优化问题中,目标函数f(x)存在多种形式,如果目标函数和约束条件都为变量x的线性函数,则称问题为线性规划;如果目标函数为二次函数,则称最优化问题为二次规划;如果目标函数或者约束条件为非线性函数,则称最优化问题为非线性优化。每个线性规划问题都有一个对应的对偶问题。对偶问题具有以下几个特性:

  1. 对偶问题的对偶是原问题;
  2. 无论原始问题是否是凸的,对偶问题都是凸优化问题;
  3. 对偶问题可以给出原始问题的一个下界;
  4. 当满足一定条件的时候,原始问题和对偶问题的解是完美等价的。

KKT条件

KKT条件是泛拉格朗日乘子法的一种形式;主要应用在当我们的优化函数存在不等值约束的情况下的一种最优化求解方式;KKT条件即满足不等式约束情况下的条件。

\begin{center} \min f(x) \\ s.t.\ \ \ \ \ h_k(x)=0 \ \ \ \ \ \ k=1, 2, \dots, p \\\ \ \ \ \ \ \ \ \ g_j(x) \le 0 \ \ \ \ \ \ j=1, 2, \dots, q \\L(x, \alpha, \beta)=f(x)+\sum_{k=1}^p\alpha_kh_k(x)+\sum_{j=1}^q\beta_j g_j(x) ,\ a_k\neq 0,\beta_j \geq 0\end{center}

可行解必须在约束区域g(x)之内,由图可知可行解x只能在g(x)<0和g(x)=0的区域取得;

 

当可行解x在g(x)<0的区域中的时候,此时直接极小化f(x)即可得到;

当可行解x在g(x)=0的区域中的时候,此时直接等价于等式约束问题的求解。

 

KKT条件理解

当可行解在约束内部区域的时候,令β=0即可消去约束。

对于参数β的取值而言,在等值约束中,约束函数和目标函数的梯度只要满足平行即可,而在不等式约束中,若β≠0,则说明可行解在约束区域的边界上,这个时候可行解应该尽可能的靠近无约束情况下的解,所以在约束边界上,目标函数的负梯度方向应该远离约束区域朝无约束区域时的解,此时约束函数的梯度方向与目标函数的负梯度方向应相同;从而可以得出β>0。

 

 

对偶问题的直观理解:最小的里面的那个最大的要比最大的那个里面的最小的大;从而就可以为原问题引入一个下界。

KKT 案例

 

这里利用该KKT条件满足对偶条件:

对偶问题的直观理解:最小的里面的那个最大的要比最大的那个里面的最小的大;从而就可以为原问题引入一个下界

KKT条件总结

KKT条件为下列五个

  1. 拉格朗日取得可行解的充要条件;
  2. 将不等式约束转换后的一个约束,称为松弛互补条件;
  3. 初始的约束条件;
  4. 初始的约束条件;
  5. 不等式约束需要满足的条件。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/454343.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

定长顺序串的实现

string.h #define MAXSTRLEN 255#define ERROR 0#define OK 1 typedef int Status;typedef char String[MAXSTRLEN 1]; //初始化字符串Status StrAssign(String T, char e); //有串S复制得串TStatus StrCopy(String T,String S); //比较两个串的大小Status StrCompare(String …

pmp思维导图 第六版_PMP考试技巧攻略(上)

PMP考试需要有保证足够的时间投入&#xff1a;获得PMP 考试并拿到5A 成绩&#xff0c;并且还需要理解性记忆&#xff1a;PMP 指定教材PMBOK第六版&#xff08;教材为必看三遍以上&#xff09;&#xff0c;学习起来是有趣的&#xff0c;同时也是痛苦的。因为看书时字面的字我们认…

程序员应该具备的素质(来自csdn)

程序员是一种技术工作&#xff0c;在IT的发展中有相当重要的地位&#xff0c;从底层硬件通讯协议的建立&#xff0c; 到数据传输层的处理&#xff0c;到操作系统的建设&#xff0c;到数据库平台的建设&#xff0c;一直到应用层上各种数 据营销平台的搭建&#xff0c;程序员在里…

linux的du使用方法

该命令的各个选项含义如下&#xff1a; -s 对每个Names参数只给出占用的数据块总数。 -a 递归地显示指定目录中各文件及子孙目录中各文件占用的数据块数。若既不指定-s&#xff0c;也不指定-a&#xff0c;则只显示Names中的每一个目录及其中的各子目录所占的磁盘块数-b 以字节为…

浅谈MVC MVP MVVM

复杂的软件必须有清晰合理的架构&#xff0c;否则无法开发和维护。 MVC&#xff08;Model-View-Controller&#xff09;是最常见的软件架构之一&#xff0c;业界有着广泛应用。 它本身很容易理解&#xff0c;但是要讲清楚&#xff0c;它与衍生的 MVP 和 MVVM 架构的区别就不容易…

机器学习接口和代码之 线性回归

线性回归sklearn 接口和代码 官网api&#xff1a;https://scikit-learn.org/stable/modules/linear_model.html#ordinary-least-squares LinearRegression class sklearn.linear_model.LinearRegression(fit_interceptTrue, normalizeFalse, copy_XTrue, n_jobs1)参数说明&a…

中国智慧VS西方智慧-看中国IT风云与IT产业怪状

为什么国外没有一家互联网公司在中国取得成功&#xff0c;为什么他们都水土不服&#xff0c;为什么他们都在中国都混不下去&#xff0c;YAHOO, EBAY等等这样享誉全球的互联网公司都在中国无法取得成功&#xff01;为什么连让IT巨无霸微软都觉得发抖&#xff0c;让比尔盖茨夜夜做…

商务搜索引擎_外贸研修 | 世界各国常用搜索引擎,开发客户必备!

我们平时生活中也好&#xff0c;开发客户也好&#xff0c;搜索引擎是我们离不开的工具。最佳没有之一的当属谷歌了。谷歌网址&#xff1a;www.google.com谷歌高级搜索&#xff1a;https://www.google.com/advanced_search (通过设置/排除一些字词缩小精确搜索范围)作为普通使用…

HaProxy+Keepalived+Mycat高可用群集配置

概述 本章节主要介绍配置HaProxyKeepalived高可用群集&#xff0c;Mycat的配置就不在这里做介绍&#xff0c;可以参考我前面写的几篇关于Mycat的文章。 部署图&#xff1a; 配置 HaProxy安装 181和179两台服务器安装haproxy的步骤一致 --创建haproxy用户 useradd haproxy--…

奇怪的bug,不懂Atom在添加markdown-themeable-pdf,在配置好phantomjs的情况下报错

本来打算用一下atom但是导出pdf报错&#xff0c;可是在预览的情况下就没有问题&#xff0c;顺便吐槽一下谷歌浏览器自己的markdown在线预览插件无法适配&#xff0c;用搜狗搭载谷歌的插件才能导出pdf&#xff0c;一下感觉逼格少了很多&#xff0c;等忙完这阵再来看一下。先贴出…

机器学习之 sklearn.preprocessing 模块

sklearn.preprocessing.PolynomialFeatures 多项式扩展。 它是使用多项式的方法来进行的&#xff0c;如果有a&#xff0c;b两个特征&#xff0c;那么它的2次多项式为&#xff08;1,a,b,a^2,ab, b^2&#xff09;&#xff0c;这个多项式的形式是使用poly的效果。 api class s…

Python 面试题

Python面试315道题第一部 Python面试题基础篇&#xff08;80道&#xff09;1、为什么学习Python&#xff1f;2、通过什么途径学习的Python&#xff1f;3、Python和Java、PHP、C、C#、C等其他语言的对比&#xff1f;PHPjavacc#c4、简述解释型和编译型编程语言&#xff1f;编译型…

周鸿祎,高司令

还是感到有必要将自己的一些想法快速记下来。 首先是对周鸿祎新员工演讲的看法。 就说实话这一点来说&#xff0c;周鸿祎比很多人强。所以我比较喜欢引用他的话&#xff0c;确实比较实在&#xff0c;不装逼。 至于一个公司招人的风格&#xff0c;是公司自己定的&#xff0c;别人…

JDBC与JNDI应用比较

JNDI用了多年但是一直没去弄懂其和JDBC的区别&#xff0c;今天在网上搜了下&#xff0c;发下些资料说明的还不错记录下。 JNDI是 Java 命名与目录接口&#xff08;Java Naming and Directory Interface&#xff09;&#xff0c;在J2EE规范中是重要的规范之一&#xff0c;不少专…

bzoj1038500AC!

序列dp 先开始想了一个类似区间dp的东西...少了一维 然后发现似乎不太对&#xff0c;因为女生的最大差和男生的最大差并不相等 dp[i][j][x][y]表示当前有i个人&#xff0c;j个男生&#xff0c;男生和女生的后缀最大差是x&#xff0c;女生和男生最大差是y&#xff0c;x,y>0,转…

机器学习接口代码之 Ridge、Lasso、Elasitc Net

目录 Ridge Regression &#xff08;岭回归&#xff09; Lasso Regression Elasitc Net&#xff08;弹性网络&#xff09; 案例&#xff1a;葡萄酒质量预测 官网地址https://scikit-learn.org/stable/modules/linear_model.html Ridge Regression &#xff08;岭回归&…

公司技术管理角度看C++游戏程序员发展

公司技术管理角度看C游戏程序员发展 H3D 这是我多年来招聘培训游戏程序员的一点想法。一直想汇总一下。主要目的是为了更好的对公司新进C程序员进行培训&#xff0c;并且建立起游戏程序员培训&#xff0c;发展&#xff0c;成才&#xff0c;成为核心骨干&#xff0c;管理层&am…

android生命周期_Android开发 View的生命周期结合代码详解

咱们以TextView控件为例&#xff1a;/*** Created by SunshineBoy on 2020/9/23.*/public class TestTextView extends android.support.v7.widget.AppCompatTextView {public TestTextView(Context context) {super(context);Log.e("TestTextView","TestTextVi…

salt

安装服务端和客户端服务端(marster)yum install salt-master -y客户端(slave)yum install salt-minion -ymarster192.168.11.17/etc/init.d/salt-master start配置文件: vi /etc/salt/mastercat master|egrep -v ;|#|^$auto_accept:True #设置自动接受日志: /var/log/salt/mas…

python | 查看pip支持的文件名和版本

python | 查看pip支持的文件名和版本win下查询大哥推荐已经解决win下查询 import pip._internalprint(pip._internal.pep425tags.get_supported())64位的需要在pip后面加个_internal 如果不行试试下边的 大哥推荐已经解决 import wheel.pep425tags as w print(w.get_suppor…