基于cvx求解储能调峰调频模型附matlab代码

✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，matlab项目合作可私信。

🍎个人主页：Matlab科研工作室

🍊个人信条：格物致知。

更多Matlab仿真内容点击👇

智能优化算法       神经网络预测       雷达通信      无线传感器        电力系统

信号处理              图像处理               路径规划       元胞自动机        无人机

⛄ 内容介绍

We consider using a battery storage system simultaneously for peak shaving and frequency regulation through a joint optimization framework which captures battery degradation, operational constraints and uncertainties in customer load and regulation signals. Under this framework, using real data we show the electricity bill of users can be reduced by up to 15%. Furthermore, we demonstrate that the saving from joint optimization is often larger than the sum of the optimal savings when the battery is used for the two individual applications. A simple threshold real-time algorithm is proposed and achieves this super-linear gain. Compared to prior works that focused on using battery storage systems for single applications, our results suggest that batteries can achieve much larger economic benefits than previously thought if they jointly provide multiple services.

⛄ 部分代码

function [total, elec, peak, bat, reg, regc, regp] = ...
    ps_only(fig_flag, lambda, battery, T, ts, tt, s)

cvx_begin quiet
    variables b(T,1);
    minimize (lambda.elec*ones(1,T)*(s-b)*ts + ...
        lambda.peak*max(s-b) + lambda.battery*norm(b,1)*ts);
    subject to
    ones(1,T)*b == 0;
    b >= -battery.power;       
    b <= battery.power;
    tril(ones(T))*b*ts <= ones(T,1)*...
        (battery.socini-battery.socmin)*battery.energy;
    tril(ones(T))*b*ts >= ones(T,1)*...
        (battery.socini-battery.socmax)*battery.energy;
cvx_end

% Bill
total = lambda.elec*ones(1,T)*(s-b)*ts + lambda.peak*max(s-b)...
    + lambda.battery*norm(b,1)*ts;
elec = lambda.elec*ones(1,T)*(s-b)*ts;
peak = lambda.peak*max(s-b);
bat = lambda.battery*norm(b,1)*ts;
reg = 0;
regc = 0;
regp = 0;

if fig_flag == 1
    % SoC
    SoC_3 = zeros(T,1);
    SoC_3(1) = battery.socini;
    for i = 2:T
        SoC_3(i) = (battery.socini*battery.energy-sum(b(1:i-1))*ts)/battery.energy;
    end

    % figure
    figure;
    subplot(2,1,1)
    hold all;
    plot(tt,s,'-.b','LineWidth',2);
    plot(tt,s-b,'--g','LineWidth',2);
    hold off;
    xlabel('time[s]');
    ylabel('power[MW]');
    grid on
    legend('grid consumption','after shaving');
    title(sprintf('Peak shaving only, shaved peak: %.2f', max(s)-max(s-b)));

    subplot(2,1,2)
    plot(tt,SoC_3,'b');
    xlabel('time[s]');
    ylabel('SoC[%]');
    grid on
    title('Battery only for peak shaving: SoC curve');
end
end

⛄ 运行结果

⛄ 参考文献

[1] Shi Y , Xu B , Wang D ,et al.Using Battery Storage for Peak Shaving and Frequency Regulation: Joint Optimization for Superlinear Gains[J].IEEE Transactions on Power Systems, 2017:1-1.DOI:10.1109/TPWRS.2017.2749512.

⛳️ 代码获取关注我

❤️部分理论引用网络文献，若有侵权联系博主删除

❤️ 关注我领取海量matlab电子书和数学建模资料

🍅 仿真咨询

1 各类智能优化算法改进及应用

生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化、公交排班优化、充电桩布局优化、车间布局优化、集装箱船配载优化、水泵组合优化、解医疗资源分配优化、设施布局优化、可视域基站和无人机选址优化

2 机器学习和深度学习方面

卷积神经网络(CNN)、LSTM、支持向量机(SVM)、最小二乘支持向量机(LSSVM)、极限学习机(ELM)、核极限学习机(KELM)、BP、RBF、宽度学习、DBN、RF、RBF、DELM、XGBOOST、TCN实现风电预测、光伏预测、电池寿命预测、辐射源识别、交通流预测、负荷预测、股价预测、PM2.5浓度预测、电池健康状态预测、水体光学参数反演、NLOS信号识别、地铁停车精准预测、变压器故障诊断

2.图像处理方面

图像识别、图像分割、图像检测、图像隐藏、图像配准、图像拼接、图像融合、图像增强、图像压缩感知

3 路径规划方面

旅行商问题(TSP)、车辆路径问题(VRP、MVRP、CVRP、VRPTW等)、无人机三维路径规划、无人机协同、无人机编队、机器人路径规划、栅格地图路径规划、多式联运运输问题、车辆协同无人机路径规划、天线线性阵列分布优化、车间布局优化

4 无人机应用方面

无人机路径规划、无人机控制、无人机编队、无人机协同、无人机任务分配

5 无线传感器定位及布局方面

传感器部署优化、通信协议优化、路由优化、目标定位优化、Dv-Hop定位优化、Leach协议优化、WSN覆盖优化、组播优化、RSSI定位优化

6 信号处理方面

信号识别、信号加密、信号去噪、信号增强、雷达信号处理、信号水印嵌入提取、肌电信号、脑电信号、信号配时优化

7 电力系统方面

微电网优化、无功优化、配电网重构、储能配置

8 元胞自动机方面

交通流 人群疏散 病毒扩散 晶体生长

9 雷达方面

卡尔曼滤波跟踪、航迹关联、航迹融合