非线性规划转化为线性规划问题(司守奎老师数模例题)-------数模学习笔记一
目录
- 一、非线性问题
- 1.问题描述与分析
- 2.问题解决思路
- 二、具体求解步骤
- 1.引入ui与vi变量
- 2.引入变量后新方程
- 3.matlab程序
- 总结
本文主要记录司守奎数模书中第一章所给例题
提示:以下是本篇文章正文内容,下面案例可供参考
一、非线性问题
1.问题描述与分析
1>、根据线性规划定义,目标函数和约束条件都要为线性方程的才算为LP问题。而上述所给目标函数中各个变量为绝对值,故为非线性。
2>、线性方程:因变量与自变量之间的关系为线性的关系,就是直角坐标系中一条直线;一个自变量对应一个因变量。而这里目标函数中有绝对值,所以很明显此方程肯定不是一个自变量对应一个因变量。
2.问题解决思路
1>、将非线性方程利用数学方法更改为线性方程,那么非线性规划问题就变为线性规划问题了。
2>、这里是由于绝对值导致方程非线性,故想方设法去掉绝对值。
二、具体求解步骤
1.引入ui与vi变量
构造xi与(ui,vi)的关系:xi = (ui,vi)
ui=(|xi|+xi)/2
vi=(|xi|-xi)/2
ui-vi=xi
ui+vi=|xi|
i = 1,2,3,4
2.引入变量后新方程
由于ui+vi=|xi| ,如:|x1| = u1+v1 2 |x2| = 2(u2+v2)
本来一个变量变为两个翻倍,故问题共有8个变量
目标函数minZ =(u1+v1)+2(u2+v2)+3(u3+v3)+4(u4+v4)=u1+2u2+3u3+4u4+v1+2v2+3v3+4v4 (整齐好看点)
s.t.约束条件 A1(Ui-Vi) = 0A2(Ui-Vi) = 1A3(Ui-Vi) = -0.5
3.matlab程序
%% 非线性转为线性问题
clc;clear;
c = [1,2,3,4,1,2,3,4]';
A1 = [1,-1,-1,1;1,-1,1,-3;1,-1,-2,3];
Aeq = [A1,-A1];
Beq = [0;1;-0.5];
Ib = zeros(8,1);
%三个等方式约束,不等式直接空矩阵
[x,f] = linprog(c,[],[],Aeq,Beq,Ib);
f
x = x(1:4)-x(5:8);
x
结果:
Optimal solution found.f =1.2500x =0.250000-0.2500
总结
难点应该是:由于变量增加一倍,故梳理起来肯定会更复杂,但是我们可以展开,就是把具体化表达,从u1—>u4,都在方程中列写出来,这样更加清晰。也要注意目标函数中变量的顺序,因为系数不同,不要搞错系数了,还有就是,有一个矩阵的拼接。
A1 = [1,-1,-1,1;
1,-1,1,-3;
1,-1,-2,3];
Aeq = [A1,-A1]; % -A1整个矩阵加在A1矩阵右侧
结果:
Aeq =1 -1 -1 1 -1 1 1 -11 -1 1 -3 -1 1 -1 31 -1 -2 3 -1 1 2 -3
非线性规划转化为线性规划问题(司守奎老师数模例题)-------数模学习笔记一相关推荐
- 司守奎《数学建模算法与应用》课后习题:线性规划
系列目录 司守奎<数学建模算法与应用>课后习题:线性规划 模拟退火算法解决旅行商问题详解 遗传算法解决旅行商问题 MATLAB实现 基于模拟退火优化的投影寻踪评价法求解供货能力评价模型 目 ...
- 数学建模算法与程序(司守奎):第一章 线性规划
第一章 线性规划 1 线性规划 1.1 线性规划的实例与定义 1.2 线性规划的MATLAB标准形式是 1.3线性规划问题解的概念 1.4 线性规划的图解法 1.5 求解线性规划的MATLAB解法 例 ...
- 【学习笔记】第六章 整数规划和非线性规划
目录 6.1 整数规划 6.1.1 整数规划问题与求解 6.2 非线性规划 6.2.1 非线性规划概念和理论 1. 非线性规划模型 2. 无约束非线性规划的求解 3. 有约束非线性规划的求解 6.2. ...
- 数学建模学习笔记(二):非线性规划模型例题与灵敏度分析
文章目录 前言 一.一个简单的非线性规划模型 二.问题的求解 1.模型的建立与求解 2.得出结论 三.灵敏度分析 总结 参考书目 前言 数学建模解决问题的过程一般分为五个步骤,称为五步方法,五个步骤如 ...
- 【数模学习笔记】【线性规划】巧用0-1变量
说明:本人是只挣扎在数模海洋里的小可怜,只是记录一下学习感悟.以下内容是在学习司守奎<数学建模算法与应用>(第三版) P21-23时的笔记. 1. 0-1变量的基本用法:0-1变量的特征 ...
- 【学习笔记】第五章 线性规划
目录 5.1 线性规划的概念和理论 1. 线性规划的一般模型 2. 线性规划解的概念及理论 3. 可转化为线性规划的问题 5.2 线性规划的Python求解 5.2.1 用scipy.optimize ...
- 陈宝林《最优化理论与算法》超详细学习笔记 (二)————补充知识(凸集) 第二章 线性规划的基本性质
陈宝林<最优化理论与算法>超详细学习笔记 (二)----补充知识 凸集 & 第二章 线性规划的基本性质 补充知识 凸集 方向与极方向 表示定理 择一定理 第一章 线性规划的基本性质 ...
- matlab 线性规划_从零开始的matlab学习笔记——(37)线性规划——后传
matlab应用--求极限,求导,求积分,解方程,概率统计,函数绘图,三维图像,拟合函数,动态图,傅里叶变换,随机数,优化问题....更多内容尽在个人专栏:matlab学习 呼,终于将matlab的优 ...
- 【四足机器人】学习笔记 欧拉角与旋转矩阵之间的转化关系
[四足机器人]学习笔记 欧拉角与旋转矩阵之间的转化 一.基础旋转矩阵 二.组合旋转矩阵 在研究四足机器人站立姿态的时候,需要将通过IMU模块采集到的欧拉角数据信息转变为旋转矩阵R,关于如何将获取的欧拉 ...
最新文章
- 硬盘重装Ubuntu12.04的感受
- SpringBoot底层注解-@ConfigurationProperties配置绑定
- 某道词典在线翻译JS解密,完全扣js代码解密 (小白基础篇)
- ApacheCN 大数据译文集 20211206 更新
- 反射(二) :整型数组能存字附---反射”暴力“超越数据类型的检测
- 企业的生存必须找到你的核心竞争力
- 一篇google developer 安全介绍的翻译
- 聚集索引与非聚集索引及其查询效率【转载】
- 锐起BSD5.0企业无盘下载免费试用
- 【Matlab图像隐写】DCT数字水印嵌入与提取【含GUI源码 943期】
- 患上糖尿病,为什么要做OCT眼底照相检查?
- numpy数组的拼接(扩维拼接和非扩维拼接)
- Opencv使用imread函数时,读到的数据为空
- docker 分析cpu占用过高
- (九)巴菲特与索罗斯的投资习惯:术业有专攻
- Web技术老矣,尚能Run否?| U4内核在Web开发平台的探索之路
- 全志 H6 芯片分析
- 四级英语测试软件,英语四级学习软件都有哪些好用啊
- Windows Server 2012 R2桌面化详细设置图解
- 【项目管理】工作分解结构