利用matlab求解线性优化问题【基于matlab的动力学模型学习笔记_11】

最新推荐文章于 2023-11-02 20:33:43 发布
最新推荐文章于 2023-11-02 20:33:43 发布
阅读量1.5k 收藏 15
点赞数 1
分类专栏： 基于matlab的动力学模型研究 文章标签： matlab 学习 开发语言
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/L1234X/article/details/126052235
版权

题目：

 （1）画出可行域范围

question1.m

%% 标出函数
r=4;
a1=0;
a2=0;
theta=0:pi/20:2*pi;
x_1=a1+r*cos(theta);
x_2=a2+r*sin(theta);
plot(x_1,x_2);
hold on
text(2,4,'16-(x_1)^2-(x_2)^2=0','color','b'); %在坐标点(6.8,4)显示x1=7这个函数线

L2=[-2,-4;5,3]; 
plot(L2(:,1),L2(:,2));hold on    %x2最大值为3
text(3,1,'2-x_1-x_2=0','color','b'); %从点L2(:,1)到点L2(:,2)

L3=[-5,0;5 0]; 
plot(L3(:,1),L3(:,2));hold on 
text(3,0,'x_1=0','color','b') 

L4=[0,-5;0,5]; 
plot(L4(:,1),L4(:,2)); 
text(0,3,'x_2=0','color','b') 

grid on

%% 填充
[X1,X2]=meshgrid(0:0.01:5,0:0.01:5);%画出区域
idX1=(X1.*X1+X2.*X2<=16)&(-X2+X1<=2)&(X1>=0)&(X2>=0);
X1=X1(idX1);
X2=X2(idX1);
k=convhull(X1,X2); %计算面积
h=fill(X1(k),X2(k),'g');  %绿色填充        
set(h,'edgealpha',0,'facealpha',0.3)  %边界，透明度

 （2）利用fmincon分别求解无约束、约束最优解

Matlab中的fmincon函数可以用于求解带约束的非线性多变量函数的最小值，这里用来求解此处的最优解。

fmincon函数参考：

Matlab求解非线性规划，fmincon函数的用法总结_小朱同学~的博客-CSDN博客_fmincon函数用法

官方帮助文档：https://ww2.mathworks.cn/help/optim/ug/fmincon.html? searchHighlight=fmincon&s_tid=srchtitle_fmincon_1

(2.1)无约束最优解

目标函数fun1.m：

%目标函数
function f=fun1(x)
f=(x(1)-2).^2+(x(2)-5).^2
end

无约束最优解question2_1.m：

%% 标出函数
r=4;
a1=0;
a2=0;
theta=0:pi/20:2*pi;
x_1=a1+r*cos(theta);
x_2=a2+r*sin(theta);
plot(x_1,x_2);
hold on
text(2,4,'16-(x_1)^2-(x_2)^2=0','color','b'); %在坐标点(6.8,4)显示x1=7这个函数线

L2=[-2,-4;5,3]; 
plot(L2(:,1),L2(:,2));hold on    %x2最大值为3
text(3,1,'2-x_1-x_2=0','color','b'); %从点L2(:,1)到点L2(:,2)

L3=[-5,0;5 0]; 
plot(L3(:,1),L3(:,2));hold on 
text(3,0,'x_1=0','color','b') 

L4=[0,-5;0,5]; 
plot(L4(:,1),L4(:,2)); 
text(0,3,'x_2=0','color','b') 

grid on

%% 填充
[X1,X2]=meshgrid(0:0.01:5,0:0.01:5);%画出区域
idX1=(X1.*X1+X2.*X2<=16)&(-X2+X1<=2)&(X1>=0)&(X2>=0);
X1=X1(idX1);
X2=X2(idX1);
k=convhull(X1,X2); %计算面积
h=fill(X1(k),X2(k),'g');  %绿色填充        
set(h,'edgealpha',0,'facealpha',0.3)  %边界，透明度


%问题2.1主函数
options=optimset;
x0=[0;0];%给定初值
lb=[0;0];%函数下限
ub=[5;5];%函数上限
[x,y]=fmincon('fun1',x0,[],[],[],[],lb,ub)

%加标注
text(-3,2,'X*(1)=2.0000') 
text(-2.1,1.6,'4.9994') 
text(-3,1.2,'f(X*(1))=4.1847e-07')

 (2.2)约束最优解

目标函数fun1.m：

%目标函数
function f=fun1(x)
f=(x(1)-2).^2+(x(2)-5).^2
end

非线性约束条件函数fun2.m：

%非线性约束条件函数
function[g,h]=fun2(x)
%matlab中默认g<=0，若不对应需取反
g(1)=-16+x(2).^2+x(1).^2;
g(2)=-2+x(1)+x(2);
h=[];%没有等式约束的时候用空值代替
end

无约束最优解question2_2.m：

%% 标出函数
r=4;
a1=0;
a2=0;
theta=0:pi/20:2*pi;
x_1=a1+r*cos(theta);
x_2=a2+r*sin(theta);
plot(x_1,x_2);
hold on
text(2,4,'16-(x_1)^2-(x_2)^2=0','color','b'); %在坐标点(6.8,4)显示x1=7这个函数线

L2=[-2,-4;5,3]; 
plot(L2(:,1),L2(:,2));hold on    %x2最大值为3
text(3,1,'2-x_1-x_2=0','color','b'); %从点L2(:,1)到点L2(:,2)

L3=[-5,0;5 0]; 
plot(L3(:,1),L3(:,2));hold on 
text(3,0,'x_1=0','color','b') 

L4=[0,-5;0,5]; 
plot(L4(:,1),L4(:,2)); 
text(0,3,'x_2=0','color','b') 

grid on

%% 填充
[X1,X2]=meshgrid(0:0.01:5,0:0.01:5);%画出区域
idX1=(X1.*X1+X2.*X2<=16)&(-X2+X1<=2)&(X1>=0)&(X2>=0);
X1=X1(idX1);
X2=X2(idX1);
k=convhull(X1,X2); %计算面积
h=fill(X1(k),X2(k),'g');  %绿色填充        
set(h,'edgealpha',0,'facealpha',0.3)  %边界，透明度

%问题2.2主函数
options=optimset;
x0=[0;0];%给定初值
lb=[0;0];%函数下限
ub=[5;5];%函数上限
[x,y]=fmincon('fun1',x0,[],[],[],[],lb,ub,'fun2')

%加标注
text(-3,2,'X*(2)=0.0000') 
text(-2.1,1.6,'2.0000') 
text(-3,1.2,'f(X*(2))=13')

 （3）线性约束最优解

目标函数fun1.m：

%目标函数
function f=fun1(x)
f=(x(1)-2).^2+(x(2)-5).^2
end

线性约束条件函数fun3.m：

%线性约束条件函数
function[g,h]=fun3(x)
%matlab中默认g<=0，若不对应需取反
g(1)=-16+x(2).^2+x(1).^2;
g(2)=-2+x(1)+x(2);

%线性约束条件
h=x(1)-x(2);
end

线性约束最优解question3.m：

%% 标出函数
r=4;
a1=0;
a2=0;
theta=0:pi/20:2*pi;
x_1=a1+r*cos(theta);
x_2=a2+r*sin(theta);
plot(x_1,x_2);
hold on
text(2,4,'16-(x_1)^2-(x_2)^2=0','color','b'); %在坐标点(6.8,4)显示x1=7这个函数线

L2=[-2,-4;5,3]; 
plot(L2(:,1),L2(:,2));hold on    %x2最大值为3
text(3,1,'2-x_1-x_2=0','color','b'); %从点L2(:,1)到点L2(:,2)

L3=[-5,0;5 0]; 
plot(L3(:,1),L3(:,2));hold on 
text(3,0,'x_1=0','color','b') 

L4=[0,-5;0,5]; 
plot(L4(:,1),L4(:,2)); 
text(0,3,'x_2=0','color','b') 

grid on

%% 填充
[X1,X2]=meshgrid(0:0.01:5,0:0.01:5);%画出区域
idX1=(X1.*X1+X2.*X2<=16)&(-X2+X1<=2)&(X1>=0)&(X2>=0);
X1=X1(idX1);
X2=X2(idX1);
k=convhull(X1,X2); %计算面积
h=fill(X1(k),X2(k),'g');  %绿色填充        
set(h,'edgealpha',0,'facealpha',0.3)  %边界，透明度
%问题3主函数
options=optimset;
x0=[0;0];%给定初值
lb=[0;0];%函数下限
ub=[5;5];%函数上限
[x,y]=fmincon('fun1',x0,[],[],[],[],lb,ub,'fun3')

%加标注
text(-3,2,'X*(3)=1.0000') 
text(-2.1,1.6,'1.0000') 
text(-3,1.2,'f(X*(3))=17.0000')

 

歪卜巴比
关注
  • 1
    点赞
  • 踩
  • 15
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
Matlab线性/非线性规划优化算法（7）
01-06
介绍如何使用quadprog and mpcqpsolver. 依然是解决一类典型的二次规划问题： 通过看matlab文档中的例子可以基本了解如何使用。下面给出几个例子： quadprog: %% the use of quadprog H = [1 -1; -1 2]; f = [-2; -6]; A = [1 1; -1 2; 2 1]; b = [2; 2; 3]; [x,fval,exitflag,output,lambda] = ... quadprog(H,f,A,b) %% H = [1 -1; -1 2]; f = [-2; -6]; Aeq = [1 1]; beq
线性优化 matlab,matlab中线性规划优化计算方法和实例？
weixin_34489135的博客
03-16 1080
1、我们首先需要知道，matlab中用于线性规划优化计算的是linprog()函数，公式是［x,fval,exitflag,output,lambda]=linprog(f,A,b,Aeq,beq,lb,ub);，其中各个参数的意思可以看下面的注释，如下图所示：2、线性规划优化计算求最优解的方法很多，有单纯形法，大M法，内点法等，linprog函数集中了这几种线性规划算法，如果感兴趣的话，可以去了...
【优化配置】基于matlab两阶段鲁棒优化算法求解微网多电源容量配置优化问题【含Matlab源码 2257期】
订阅付费专栏（初级版），可赠送该付费专栏任意代码1份；
12-08 1863
两阶段鲁棒优化算法求解微网多电源容量配置优化问题 完整的代码，方可运行；可提供运行操作视频！适合小白！
Matlab优化算法以及应用案例分析
最新发布
vipfanxu的博客
11-02 934
无论是无约束优化、约束优化还是全局优化问题，Matlab都提供了相应的函数接口和工具箱来支持。通过灵活使用Matlab的优化算法，可以帮助我们解决实际问题，提高问题的求解效率。无约束优化问题是指在没有约束条件的情况下，找到使得目标函数达到最小或最大值的变量取值。A,b是线性约束条件，Aeq,beq是线性等式条件，lb,ub是变量的取值范围，非线性约束条件通过nlcon函数定义。在上述代码中，fun是目标函数的定义，gs是全局优化求解器，problem是优化问题的设置。通过run函数来求解全局最优解。
matlab求解最优化问题(数学建模)
热门推荐
wmx1117的博客
06-08 6万+
matlab求解最优化问题（数学建模） 1.线性规划 matlab中线性规划优化计算方法和实例 在matlab中用于线性规划优化计算的是linprog()函数。 公式:[x,fval,exitflag,output,lambda]=linprog(c,A,b,Aeq,beq,lb,ub,x0); %x:表示最优解，fval:表示目标函数最优值，exitflag:表示求解的结果是成功还是失败，1代表成功。 %output优化过程中的各种输出信息，lambda：结构体，包含最优解处的拉格朗日乘子。 %c:目标
matlab中线性规划优化计算方法和实例
phymat.nico的专栏
09-20 516
https://jingyan.baidu.com/article/8065f87fbc11682330249843.html matlab解决一般的一元/多元函数的非线性优化问题_Konata的博客-CSDN博客
matlab线性规划搜索法优化,Matlab线性规划优化算法（1）
weixin_39719729的博客
03-18 555
在Matlab中解形如下式的线性规划问题：其中包括优化对象 f’ * x, 不等式约束，等式约束，以及约束变量的上下界。在Matlab中提供了linprog函数进行线性优化的求解：eg:[x,fval,exitflag,output,lambda] = linprog(f,A,b,Aeq,beq,lb,ub,options)函数的输入f, 即为优化对象f, A,b, Aeq, beq, lb,u...
Matlab优化算法-线性规划
dele
08-19 1606
优化算法-线性规划
Matlab的最优化实现方法之线性规划：linprog
lianghao11的博客
07-31 3601
一、概述 所谓最优化方法，即指采用数学方法研究系统的优化途径及方案，为决策者提供科学决策的依据，诸如低成本、高回报、长寿命等问题。 最优化问题的数学描述： 1、优化模型的三要素：变量、目标函数、约束条件； 变量（Variables）：设计变量之间要求线星独立，不存在确定的数学关系；依据设计变量的数量，将优化问题分为一维问题、N维问题，依据设计变量的离散性，分为离散优化和连续优化。 目标函数（obje...
matlab求解最优化线性规划问题代码以及思维导图
03-25
思维导图，包含了最优化问题中所有的算法记起分类和相关实现图片。附华南理工大学matlab实现最优化问题求解的代码课件以及相关代码示例，详细版。
Matlab有约束规划，无约束规划
01-19
数学建模必须要掌握的模型就是线性规划，这个rar包括乐线性规划与非线性约束规划。
MATLAB最优化计算20例
04-05
黄金搜索法求解无约束最优化问题 二次插值法求解无约束最优化问题 Nelder-Mead算法求解无约束最优化问题 最速下降法求解无约束最优化问题 牛顿法求解无约束最优化问题 无约束最优化问题求解综合 遗传算法求解无约束最优化问题 拉格朗日乘子法求解约束最优化问题 惩罚函数法求解约束最优化问题 无约束最优化函数应用实例 约束最优化函数应用实例 线性规划函数应用实例 最大利润问题综合实例 最大容积问题综合实例 最优生产决策问题综合实例 投资问题综合实例 最小费用问题综合实例 最佳定位问题综合实例
基于MATLAB求解气泡动力学方程
06-14
基于MATLAB求解RP方程以及一二阶气泡动力学方程。
【Matlab算法】粒子群算法求解二维线性优化问题（附MATLAB代码）
AlbertDS的博客
03-26 2835
二维线性优化问题指的是在二维空间中，对于一个由线性函数构成的目标函数，通过限制自变量的范围或满足特定的约束条件，寻找一个最优解（最小值或最大值）。
MATLAB关于线性优化和整数线性优化
weixin_43326436的博客
10-22 309
Aeq，beq是等式约束中的矩阵Aeq与向量beq，A，b是不等式约束Ax<=b中的矩阵A与向量b，ub，lb是决策变量上下界。2.混合整数线性优化函数。
利用Lingo和MATLAB求解线性规划问题模型D的案例
05-28
好的，我将简要介绍如何使用Lingo和MATLAB求解线性规划问题模型D的案例。 首先，让我们看一下模型D的数学表示： 最大化 Z = 5x1 + 4x2 约束条件： 2x1 + x2 ≤ 10 x1 + 2x2 ≤ 8 x1, x2 ≥ 0 接下来，我们将...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值