最近写小论文做数值仿真时遇到了问题,模型的目标函数中包含上下限为决策变量的正态分布概率密度函数以及分布函数,使用yalmip调cplex求解发现结果有点问题,网上求助无果,最后求助于matlab自带的fmincon函数。由于模型比较复杂,因此下文把核心问题列出来进行说明。

上式为模型目标函数中的一部分,其中为模型的决策变量, 、均为正态分布概率密度函数,但服从不同的均值和方差。

首先,对积分项进行表达,通过matlabFunction将int输出的syms转化为句柄函数,然后使用integral对外层进行积分,保证输出为数值型。由于fmincon要求输入为句柄函数,因此最外层利用@将整个式子转化为句柄函数。完整代码如下:

Matlab代码

syms x;
f = @(v)pr*(integral(matlabFunction(normpdf(x,mu_qn,sigma_qn)*(int(@(y)y.*normpdf(y,mu_qr,sigma_qr),0,v(2)+v(1)-x)

未考虑清楚的问题:

1、最里层使用integral会报错,提示上下限有问题;

2、内部乘法运算必须写为点乘,否则报错维度不一致;

补充:

使用fmincon求解时,由于模型复杂度较高,因此算法采用默认的interior - point无法在受限的迭代次数中求得最优解,此时可考虑将算法改为精确度更高的sqp,代码如下:

options = optimoptions('fmincon','Display','iter','Algorithm','sqp');

Matlab Fmincon 解决带积分的二元非线性规划问题相关推荐

  1. Matlab+YALMIP+CPLEX解决带储能的微电网调度优化

    Matlab+YALMIP+CPLEX解决带储能的微电网优化调度问题,以微电网运行成本最小为目标函数,能量设备包含风电.光伏以及蓄电池,同时考虑到电价和与电网的交互,约束条件包含蓄电池SOC约束.交换 ...

  2. Matlab+YALMIP+CPLEX解决带储能的微电网优化调度问题(完整代码分享)

    2Matlab+YALMIP+CPLEX解决带储能的微电网优化调度问题 请先确保YALMIP工具箱和CPLEX正确安装,MATLAB导入对应文件,否则无法运行程序!! 如果对你有用,请给博主点歌赞!! ...

  3. Matlab中解决出现的错误使用 svmtrain (line 234) Y must be a vector or a character array.问题

    Matlab中解决出现的错误使用 svmtrain (line 234) Y must be a vector or a character array.问题 目录 解决问题 解决思路 解决方法 解决 ...

  4. 基带信号及其眼图MATLAB仿真实现,Matlab通信仿真——带限系统下的基带信号

    Matlab通信仿真--带限系统下的基带信号 Matlab通信仿真--带限系统下的基带信号 1 余弦滚降特性 满足消除码间串扰条件的H(f)有很多种,容易想到的一种极限情况,就是H(f)为理想低通型. ...

  5. matlab 眼图 值,Matlab通信仿真——带限系统下的基带信号

    Matlab通信仿真--带限系统下的基带信号 Matlab通信仿真--带限系统下的基带信号 1 余弦滚降特性 满足消除码间串扰条件的H(f)有很多种,容易想到的一种极限情况,就是H(f)为理想低通型. ...

  6. matlab图像导数求积分_第二讲matlab求微分方程导数积分

    第二讲matlab求微分方程导数积分 第二讲 导数与微分方程 一.实验内容 1.实际引例 (牛顿冷却模型)警察上午9点钟发现一被谋杀者,并测得尸体温 度为32.4℃,一小时以后,尸体的温度变为31.7 ...

  7. matlab收敛半径,运用MATLAB语言解决级数及其相关问题 李娟娟

    运用MATLAB语言解决级数及其相关问题 李娟娟 MATLAB 语言课程论文运用 MATLAB 语言解决级数及其相关问题姓名:李娟娟学号:12010245220专业:电子信息工程班级:2010 级电子 ...

  8. Matlab:解决“内存不足”错误

    Matlab:解决"内存不足"错误 问题 可能的解决方案 在不再需要时清除变量 使用适当的数据存储 使用相应的数值类. 减少存储数据时的开销量. 尽可能使数组稀疏. 使用适当的 M ...

  9. 2021-05-11 Matlab遗传算法工具箱的使用及实例(非线性规划)

    Matlab遗传算法工具箱的使用及实例(非线性规划) 本文将介绍MATLAB遗传算法工具箱求解非线性规划问题.在阅读本文之前,建议读者阅读上一期"MATLAB遗传算法工具箱求解线性规划问题& ...

最新文章

  1. osg中实现HUD(OSG初级篇1)
  2. JavaWeb三大组件(Servlet,Filter,Listener 自己整理,初学者可以借鉴一下)
  3. Windows Server 2012正式版RDS系列⑥
  4. junit 静态方法_预期异常规则和模拟静态方法– JUnit
  5. 计算机三种引用方式,单元格的引用方式有哪几种
  6. HTML5 FileReader API 测试(一)
  7. Linq to Xml (1) 创建和查询包含命名空间的XML文档
  8. 防止数据中心停机需要采取什么措施
  9. hiho 第1周 最长回文子串
  10. pyplot设置刻度字体大小以及标签字体大小
  11. Python单下划线与双下划线
  12. amazon sqs java_Amazon SQS消息多传递
  13. kvaser二次开发
  14. 吴翰清:云计算安全是“皇帝的新装”么?
  15. 支付宝沙箱钱包的下载与登录
  16. SpringBoot整合redis,启动提示错误信息: An exception was thrown by org.redisson.misc.RedissonPromise$$Lambda$682
  17. 当远程连接MySQL数据库的时候显示Can't connect to MySQL server (10060
  18. unity中多屏显示
  19. UM5202EEDFSOT-143工作电压5V用于高速线路保护的2线ESD保护二极管阵列管UM5202EEDF USB端口和以太网端口保护TVS二极管阵列,低泄漏电流和箝位电压25V结电容1PF
  20. 吉林大学计算机专业,考生吉林大学计算机专业,心里闷闷不乐怎么办?

热门文章

  1. Bootloader(启动引导程序)--->u-boot
  2. 晶体生长计算机实验报告,食盐晶体生长实验报告
  3. 2019年中国社交电商行业市场分析:社交电商巨头已浮出水面
  4. excel查找空值快捷键_excel定位选取:再谈快捷键Ctrl+G的妙用
  5. 军品研制过程所需文件-进阶版
  6. 双向长短期记忆网络模型_基于深度双向长短期记忆网络的空气质量预测方法与流程...
  7. 概率论-4.3随机变量序列的两种收敛性(待补充)
  8. 工具型产品如何建立自己的商业化模式?
  9. MPLAB 安装 mmc 报错需要安装web相关插件
  10. MPLAB常见问题及解决方法