知识分享|分段函数线性化及matlab测试

目录

[1 使用0-1变量将分段函数转换为线性约束](#1 使用0-1变量将分段函数转换为线性约束)

[2 连续函数采用分段线性化示例](#2 连续函数采用分段线性化示例)

[3 matlab程序测试](#3 matlab程序测试)

[4 matlab测试结果说明](#4 matlab测试结果说明)

[5 分段线性化应用](#5 分段线性化应用)


1 使用0-1变量将分段函数转换为线性约束

2 连续函数采用分段线性化示例

3 matlab程序测试

clc;clear all;gn=10;tn=1;x_pf=sdpvar(1, tn,'full');Pgone=8;gw1=sdpvar(gn+1,tn,'full');gz1=binvar(gn, tn,'full');gl1=10/gn;for i=1:1gl2(i,:)=0:gl1:10;endcon=[];con = [con, x_pf(1,:)==gl2(1,:).^2*gw1];con = [con, gw1(1,:)<=gz1(1,:)];for i=2:gn    con = [con, gw1(i,:)<=gz1(i-1,:)+gz1(i,:)];endcon = [con, gw1(gn+1,:)<=gz1(gn,:)];con = [con, sum(gw1)==ones(1,tn)];con = [con, sum(gz1)==ones(1,tn)];con = [con, Pgone(1,:)==gl2(1,:)*gw1];​con = [con, gw1>=0];f=1;ops=sdpsettings('solver','cplex');result=optimize(con,f,ops);

4 matlab测试结果说明

在上述测试程序中,Pgone是平方之前的变量,x_pf是平方项,为了测试运算效果,目标函数采用定值,这样就能验证约束是不是严格限制,上述程序将Pgone取值为8,最终得到x_pf得64,将Pgone改成5,则x_pf得25,可以看出,此分段线性化模型严格约束。

大家在测试过程中,可以调节gn的值,该值代表的是分段数量,调节该值能够发现,随着分段数量增多,分段线性化的精度也会越来越高,如,当Pgone=5,gn=3时,x_pf=27.78,随着gn增大,x_pf值也就更加精确。

5 分段线性化应用

在电力系统中,常规机组的成本涉及到机组功率的平方项,这时候就可以采用分段线性化,如日前日内多阶段多时间尺度源荷储协调调度(matlab代码)考虑源荷两侧不确定性的含风电电力系统低碳调度等程序中均有应用。

相关推荐
xrgs_shz2 小时前
MATLAB读入不同类型图像并显示图像和相关信息
图像处理·计算机视觉·matlab
火山口车神丶5 小时前
某车企ASW面试笔试题
c++·matlab
今天吃饺子6 小时前
2024年SCI一区最新改进优化算法——四参数自适应生长优化器,MATLAB代码免费获取...
开发语言·算法·matlab
Evand J18 小时前
集合卡尔曼滤波(Ensemble Kalman Filter),用于二维滤波(模拟平面上的目标跟踪),MATLAB代码
matlab·平面·目标跟踪
我是博博啦1 天前
matlab例题
人工智能·算法·matlab
2402_871321951 天前
MATLAB方程组
gpt·学习·线性代数·算法·matlab
Matlab程序猿小助手2 天前
【MATLAB源码-第218期】基于matlab的北方苍鹰优化算法(NGO)无人机三维路径规划,输出做短路径图和适应度曲线.
开发语言·嵌入式硬件·算法·matlab·机器人·无人机
Matlab精灵2 天前
使用MATLAB进行字符串处理
开发语言·matlab
机器学习之心2 天前
POD-Transformer多变量回归预测(Matlab)
matlab·回归·transformer·pod-transformer
WangYan20223 天前
高光谱遥感是什么?高光谱遥感数据如何处理?(基于Matlab和Python多案例解析)从小白到精通
matlab·高光谱遥感数据处理·混合像元分解