最优化方法复习——线性规划之对偶问题

一、线性规划对偶问题定义

原问题:

\begin{gathered} (LP)\mathrm{Max}\quad z=c^\top x \\ \quad \quad s.t. \quad Ax \leq b \\ \quad \quad \quad \quad x \geq 0 \\ \end{gathered}

对偶问题:

\begin{gathered}(DP)\mathrm{Min}\quad f=b^\top y \\ \quad \quad s.t. \quad A^\top y\geq c \\ \quad \quad \quad \quad y\geq0 \\ \end{gathered}

(1)若一个模型为目标求 “极大”,约束为“小于等于” 的不等式,则它的对偶模型为目标求“极小”,约束是“大于等于”的不等式。即“Max, ≤”和“Min,≥”相对应。

(2)从约束系数矩阵看:一 个模型中为A,则另一个模型中为AT 。一个模型是m个约束、 n个变量,则它的对偶模型为 n个约束、m个变量。

(3)从数据b、c的位置看:在两个规划模型中,b和c的位置对换。

(4)两个规划模型中的变量皆非负。

一般称不具有对称形式的一对线性规划为非对称形式的对偶规划。 对于非对称形式的规划,可以按照下面的对应关系直接给出其对偶规划。

(1)将模型统一为“Max,≤”或“Min,≥” 的形式,对 于其中的等式约束按下面(2)、(3)中的方法处理。

(2)若原规划的某个约束条件为等式约束,则在对偶规划中与此约束对应的那个变量取值没有非负限制。

(3)若原规划的某个变量的值没有非负限制,则在对偶问题中 与此变量对应的那个约束为等式。

二、对偶定理与影子价格

  • 定理(弱对偶定理) 若 x,y 分别为(LP) 和(DP)的可行解, 那么c^Tx \leq b^Ty。 推论:设(LP)有可行解,那么若(LP) 无有界最优解,则(DP)无可行解。
  • 定理 (最优性准则定理) 若x^* ,y^*分别(LP),(DP)的可行解,且c^Tx^* =b^Ty^*,那么x^* ,y^*分别为(LP)和(DP) 的最优解。
  • 定理(主对偶定理) 若(LP)有最优解,则(DP)也有最优解。反之也成立,且最优值相等。

影子价格是一个向量,它的分量表示最优目标值随相应资源数量变化的变化率。

影子价格反映了不同的局部或个体的增量可以获得不同的整体经济效益。如果为了扩大生产能力,考虑增加设备,就应该从影子价格高的设备入手。这样可以用较少的局部努力,获得较大的整体效 益。

三、对偶单纯形法

1、基本思想

从原规划的一个基本解出发,此基本解不一定可行,但它对应着一个对偶可行解(检验数非正),所以也可以说是从 一个对偶可行解出发;然后检验原规划的基本解是否可行,即是否有负的分量,如果有小于零的分量,则进行迭代,求另一个基本解,此基本解对应着另一个对偶可行解(检验数非正)。

如果得到的基本解的分量皆非负,则该基本解为最优解。也就是说,对偶单纯形法在迭代过程中始终保持对偶解的可行性(即检验数非正),使原规划的基本解由不可行逐步变为可行,当同时得到对偶规划与原规划的可行解时,便得到原规划的最优解

什么情况下使用对偶单纯性法呢?

在初始计算中$ \sigma^{T}=(-c^{T},0^{T})\leq0$,即对偶可行,但是由于$-b\leq0$,所以\begin{cases}x_{s}=-b\leq0\\x=0\end{cases},所以原问题不可行。

应用前提:有一个基,其对应的基满足:

(1) 单纯形表的检验数行全部非正(对偶可行);

(2) 变量取值可有负数(非可行解)。

注:通过矩阵行变换运算,使所有相应变量取值均为非负数即得到最优单纯形表。

2、步骤

(1)建立初始对偶单纯形表,对应一个基本解,所有检验数均非正,转 (2)。

(2)若b\geq 0,则得到最优解,停止;否则,若有b_k<0则选k行的基变量为出基变量,转(3)。

(3)若所有a_{kj} \geq 0(j = 1,2,…,n),则原问题无可行解,停止(因为以任何a_{kj}为主元做主元消去时,都不可能使b_k变为正数); 否则,若有a_{kj}<0则选\theta=\min\{\sigma_{j}/a_{kj}|a_{kj}<0\}=\sigma_{r}/a_{kr}, 那么x_r为进基变量,转(4)。

(4)以a_{kr}为主元,作矩阵行变换使其变为1,该列其他元变为0,利用主元消去计算A,b,\sigma,转(2)。

$\theta$的选取保证新的检验数$\sigma_j^{\prime}\leq0$,因为

$ \sigma_j^{\prime}=\sigma_j-\frac{a_{kj}}{a_{kr}}\sigma_r=a_{kj}(\frac{\sigma_j}{a_{kj}}-\frac{\sigma_r}{a_{kr}}) $

主元消去运算后,对偶问题的目标函数值减小(至少不增大),因为


$ -(c_B^Tb)^{\prime}=-c_B^Tb-\frac{\sigma_r}{a_{kr}}b_k $

由于\frac{\sigma_r}{a_{kr}}b_k\leq0,故$-(c_B^Tb)'\geq-c_B^Tb$,即(c_B^Tb)'\leq c_B^Tb

注:由于主元消去前a_{kj}b_k同为负数,因此主元消去后右端列第k个分量变成正数,这有利于基本解向着满足可行性的方向转化。

3、例子

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

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

相关文章

2024年甘肃省职业院校技能大赛信息安全管理与评估三阶段理论样题一

2024年甘肃省职业院校技能大赛高职学生组电子与信息大类信息安全管理与评估赛项样题一 第六部分 理论技能与职业素养&#xff08;100 分&#xff09; 【注意事项】 1.该部分答题时长包含在第三阶段竞赛时长内&#xff0c;请在临近竞赛结束前提交。 2.参赛团队可根据自身情况…

数据库系统概论复习资料

数据库系统概论考试需知 一、分值分布 1、判断题&#xff08;10分&#xff09; 1分一个 2、填空题&#xff08;20分&#xff09; 2分一个 3、选择题&#xff08;20分&#xff09; 2分一个 4、分析题&#xff08;30分&#xff09; 第一题10分&#xff0c;第二题…

如何设置echart图表在vue页面屏幕比例缩放自适应问题

问题&#xff1a;页面的echart图表在浏览器缩放屏幕比例时无法随着屏幕的比例自动改变大小 解决方式&#xff1a; 可以通过监听窗口的 resize 事件&#xff0c;并在事件回调函数中重新调整图表的大小。 <template><div ref"chartContainer" style"w…

Enterprise Architect 12版本使用教程

Enterprise Architect 12版本使用教程 1.下载安装Enterprise Architect 122.Enterprise Architect原始DDL模板配置及存在的问题1.DDL Column Definition原始模板&#xff08;没有default值&#xff1a;可忽略&#xff09;2.DDL Data Type原始模板&#xff08;timestamp等时间字…

Apollo新版本Beta自动驾驶技术沙龙参会体验有感—百度自动驾驶开源框架

在繁忙的都市生活中&#xff0c;我们时常对未来的科技发展充满了好奇和期待。而近日&#xff0c;我有幸参加了一场引领科技潮流的线下技术沙龙&#xff0c;主题便是探索自动驾驶的魅力——一个让我们身临其境感受创新、了解技术巨擘的机会。 在12月2日我有幸参加了Apollo新版本…

智能优化算法应用:基于沙猫群算法无线传感器网络(WSN)覆盖优化 - 附代码

智能优化算法应用&#xff1a;基于沙猫群算法无线传感器网络(WSN)覆盖优化 - 附代码 文章目录 智能优化算法应用&#xff1a;基于沙猫群算法无线传感器网络(WSN)覆盖优化 - 附代码1.无线传感网络节点模型2.覆盖数学模型及分析3.沙猫群算法4.实验参数设定5.算法结果6.参考文献7.…

kotlin协程反编译java学习原理

前情提要 GlobalScope.launch(Dispatchers.Main) { // 默认是Default异步 // 1.从当前协程体作用域Dispatchers.Main 协程依附的线程 到 挂起点 挂起函数 到 执行 请求耗时操作 而 切换到 IO异步线程 // 2.IO异步线程执行完成后&#xff0c;开始恢复&#xff0c;当前作…

modbus转profinet网关解决plc插槽号不够用的情况

PLC作为常用的控制设备之一&#xff0c;其插槽号有时会限制外部设备的连接数量。然而&#xff0c;通过使用modbus转profinet网关&#xff0c;可以解决这一问题。这种设备能够将modbus协议转换为profinet协议&#xff0c;实现PLC与更多外部设备的连接。 modbus转profinet网关还具…

游戏盾的防御原理以及为什么程序类型更适合接入游戏盾。

游戏盾是一种专门用于游戏服务器的安全防护服务&#xff0c;旨在抵御各种网络攻击。它的原理主要包括以下几个方面&#xff1a; 流量清洗和过滤&#xff1a;游戏盾会对进入游戏服务器的流量进行实时监测、分析和过滤。它通过识别恶意流量和攻击流量&#xff0c;过滤掉其中的攻击…

浏览器渲染页面的过程以及原理

什么是浏览器渲染 简单来说&#xff0c;就是将HTML字符串 —> 像素信息 渲染时间点 浏览器什么时候开始渲染&#xff1f; 网络线程发送请求&#xff0c;取回HTML封装为渲染任务并将其传递给渲染主线程的消息队列。 问题&#xff1a;只取回HTML吗&#xff1f;那CSS和JS呢&am…

面试经典150题(1-2)

leetcode 150道题 计划花两个月时候刷完&#xff0c;今天完成了两道(1-2)150&#xff1a; &#xff08;88. 合并两个有序数组&#xff09;题目描述&#xff1a; 给你两个按 非递减顺序 排列的整数数组 nums1 和 nums2&#xff0c;另有两个整数 m 和 n &#xff0c;分别表示 n…

元宇宙vr党建云上实景展馆扩大党的影响力

随着科技的飞速发展&#xff0c;VR虚拟现实技术已经逐渐融入我们的日常生活&#xff0c;尤其在党建领域&#xff0c;VR数字党建展馆更是成为引领红色教育新风尚的重要载体。今天&#xff0c;就让我们一起探讨VR数字党建展馆如何提供沉浸式体验&#xff0c;助力党建工作创新升级…

十年OpenCV开发以后发布的作品 - OpenCV实验大师

OpenCV介绍 OpenCV是知名的计算机视觉框架&#xff0c;支持数十个不同的视觉处理模块&#xff0c;提供了超过2000多个传统算法&#xff0c;其核心功能支持图像处理、图像分析、特征提取、对象检测、深度学习模型推理等。当前支持C、Python、JS、C#等多种语言SDK&#xff0c;支…

智能优化算法应用:基于袋獾算法无线传感器网络(WSN)覆盖优化 - 附代码

智能优化算法应用&#xff1a;基于袋獾算法无线传感器网络(WSN)覆盖优化 - 附代码 文章目录 智能优化算法应用&#xff1a;基于袋獾算法无线传感器网络(WSN)覆盖优化 - 附代码1.无线传感网络节点模型2.覆盖数学模型及分析3.袋獾算法4.实验参数设定5.算法结果6.参考文献7.MATLAB…

雷达点云数据.pcd格式转.bin格式

雷达点云数据.pcd格式转.bin格式 注意&#xff0c;方法1原则上可行&#xff0c;但是本人没整好pypcd的环境 方法2是绝对可以的。 方法1 1 源码如下&#xff1a; def pcb2bin1(): # save as bin formatimport os# import pypcdfrom pypcd import pypcdimport numpy as np…

python pandas dataframe常用数据处理总结

最近一直在做数据处理相关的工作&#xff0c;有几点经常遇到的情况总结如下&#xff1a; 数据中存在为空数据如何处理 处理方式1&#xff1a;丢弃数据行 # 实现方式1 data data.dropna(subset[id]) # 若id列中某行数值为空&#xff0c;丢弃整行数据 # 实现方式2 data df[df…

Ant Design Vue 年选择器

文章目录 参考文档效果展示实现过程 参考文档 提示&#xff1a;这里可以添加本文要记录的大概内容&#xff1a; DatePicker 日期选择框 大佬&#xff1a;搬砖小匠&#xff08;Ant Design vue 只选择年&#xff09; 提示&#xff1a;以下是本篇文章正文内容&#xff0c;下面案…

光标虽小,但是能让其凭空消失的原因不少,需要仔细检查

在Windows 10中遇到鼠标光标有问题吗?我们已经为你提供了所需的修复程序。 光标消失的原因 光标不断消失可能是由各种原因引起的,因此有可能找到各种各样的解决方案。光标可能根本无法工作,或者在特定情况下可能会消失。鼠标按钮甚至可能在光标隐藏时工作。 以下是用户注…

如何验证一个URL是否合法

在JavaScript中&#xff0c;可以使用正则表达式&#xff08;RegExp&#xff09;或使用内置的URL对象来校验一个URL。下面是一些常用的方法以及对应的代码示例&#xff1a; 使用正则表达式进行校验&#xff1a; function validateURL(url) {const pattern /^(https?:\/\/)?…

使用Caliper对Fabric地basic链码进行性能测试

如果你需要对fabric网络中地合约进行吞吐量、延迟等性能进行评估&#xff0c;可以使用Caliper来实现&#xff0c;会返回给你一份网页版的直观测试报告。下面是对test-network网络地basic链码地测试过程。 目录 1. 建立caliper-workspace文件夹2. 安装npm等3. calipe安装4. 创建…