【硕士论文复现】可再生能源发电与电动汽车的协同调度策略研究(Matlab代码实现)

2024-01-21

???????????????? 欢迎来到本博客 ❤️❤️????????

????博主优势: ???????????? 博客内容尽量做到思维缜密,逻辑清晰,为了方便读者。

⛳️ 座右铭: 行百里者,半于九十。

???????????? 本文目录如下: ????????????

目录

????1 概述

1.1 上层问题

1.2 下层问题

????2 运行结果

????3 参考文献

????4 Matlab代码、数据


????1 概述

文献来源:

1.1 上层问题

在第3章所建立的上层模型中,目标函数为二次凸函数,且约束条件均为线性,属于典型的二次规划问题,二次规划问题是最基本的非线性规划问题,其主要特点就是目标函数为二次型且约束均为线性约束,一般形式可以表示如下:

1.2 下层问题

下层对各充电站有功无功充放量以及各可再生能源发电站无功出力的优化过程,实际上就是调度对象仅包含电动汽车储能和可再生能源的配电网的有功无功协调优化过程,属于最优潮流问题。最优潮流的求解方法主要分为智能算法和数值分析法两大类。智能算法对问题的形式没有任何要求,具有极强的通用性,并且从理论上具有全局寻优的能力。但在实际应用中,该类算法容易陷入局部最优,而且存在计算量大,寻优速度慢,求解结果不稳定等弊端[42]。
本文提出的调度策略是实时进行的,具有一定的计算速度的要求,从而也限制了对智能算法的使用。数值分析法大都只具备局部寻优能力,只有凸规划方法具有全局寻优性。此外,下层模型中计及充电桩及可再生能源并网逆变器的无功响应,配电网的运行优化过程涉及有功及无功的协调优化,该模型中的最优潮流问题是一个非凸非线性问题,基于凸优化理论的凸规划能够通过在算法上进行松弛为非凸问题提供一个有效的解决途径。
目前,凸优化算法中的半正定规划(semi-definite programming, SDP)和二阶锥规划(second order conic programming, SOCP)在电力系统中的应用较为广泛,其中SOCR作为锥规划中的一个特例,使用和求解较其他锥规划更加容易,因此对SOCR 的研究和应用也最为广泛。
二阶锥规划作为一种凸优化,拥有凸规划的所有优势,因此在二阶锥规划问题中,所有的局部最优解都是全局最优的,且求解效率及准确性相较于非线性规划有了明显提高。但它要求优化模型的目标函数必须为线性43,因此在下层模型的目标函数中,使用支路电流描述系统的有功网损,即可将原目标函数转换为线性函数:

详细讲解见第4部分。

???? 2 运行结果

部分代码:

%% 1.设参
mpc = IEEE33BW;
pload = mpc.Pload(:,t)*(pload1(t)/sum(mpc.Pload(:,t)))/100;%节点有功负荷
qload = mpc.Qload(:,1);%节点无功负荷
branch = mpc.branch_CG(1:32,:);
branch(:,3) = branch(:,3)*100/(12.66^2);%求阻抗标幺值
r=real(branch(:,3));
x=imag(branch(:,3));
r=r(1:32);
x=x(1:32);
T = 1;%时段数为1小时
nb = 33;%节点数,根节点为33
nl = 32;%支路数
nc = 5;%联络开关数

upstream=zeros(nb,nl);%代表流入节点支路
dnstream=zeros(nb,nl);%代表流出节点支路
for i=1:32
upstream(i,i)=1;
end
% upstream(20,33)=1;%支路33为20-7支路,流入节点20
% upstream(14,34)=1;%支路34为14-8支路,流入节点14
% upstream(21,35)=1;%支路35为21-11支路,流入节点21
% upstream(32,36)=1;%支路36为32-17支路,流入节点32
% upstream(28,37)=1;%支路37为28-24支路,流入节点28

????3 参考文献

文章中一些内容引自网络,会注明出处或引用为参考文献,难免有未尽之处,如有不妥,请随时联系删除。

[1]赵浩然. 可再生能源发电与电动汽车的协同调度策略研究[D].山东大学,2020.DOI:10.27272/d.cnki.gshdu.2020.001215.

???? 4 Matlab代码、数据

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

【硕士论文复现】可再生能源发电与电动汽车的协同调度策略研究(Matlab代码实现) 的相关文章

  • Matlab:掩码/创建一个知道其原点且具有一定半径的圆形 roi

    只是一个简单的问题 我有一张图像 并且提取了某个点 特征 我知道每个帧中该点的坐标 说 x1 和 y1 我需要一个圆形 ROI 形式 该点在图像上具有我选择的半径 我尝试了 impoly 和 roipoly 当我知道图像中的要点时 不知道如
  • 计算数组中接下来的 n 个元素的乘积

    我想计算下一个的乘积n矩阵的相邻元素 号码n要相乘的元素数应在函数的输入中给出 例如 对于此输入 我应该从第一个开始计算每 3 个连续元素的乘积 p ind max product 1 2 2 1 3 1 3 这给出了 1 2 2 2 2
  • 检查Matlab中脚本需要使用的函数

    我有一个别人写的代码包 我正在运行一个脚本 它调用一些函数 这些函数又调用更多函数 等等 我想获取不是 MATLAB 内置函数但属于包的一部分的函数列表 我尝试使用matlab codetools requiredFilesAndProdu
  • 扩展 MATLAB 函数名称的最大长度

    我编写了一个 MATLAB 程序 可以动态创建自定义 MATLAB 函数 并使用以下命令在其他 MATLAB 实例中启动它们unix命令 我使用这个程序来自动化 fMRI 神经影像分析 使用 SPM8 for MATLAB 一切正常 但是
  • 两个 y 轴与相同的 x 轴[重复]

    这个问题在这里已经有答案了 可能的重复 在单个图中绘制 4 条曲线 具有 3 个 y 轴 https stackoverflow com questions 1719048 plotting 4 curves in a single plo
  • 如何使用matlab生成不同频率的正弦波?

    对于我的项目 我需要使用 matlab 生成一个正弦波 它有 100 000 个样本 并且频率在每 10 000 个样本后随机变化 采样率和频率可以根据方便而定 matlab中有没有函数可以生成这个 好的另一个例子 生成 5 个随机频率 r
  • 在另一列中添加具有特定条件的一列,如 excel 的 sumif

    我有一个像这样的矩阵 A 1 2 2 3 3 4 4 5 5 6 6 8 7 9 8 5 9 4 现在我想添加第二列 条件是如果 limit 0 interval 3 且 limit limit interval 或者换句话说 当第 1 列
  • 在 matlab 中求 3d 峰的体积

    现在我有一个带有峰值的 3D 散点图 我需要找到其体积 我的数据来自图像 因此 x 和 y 值表示 xy 平面上的像素位置 z 值是每个像素的像素值 这是我的散点图 scatter3 x y z 20 z filled 我试图找到数据峰值的
  • 句柄类和值类的区别

    我有一些 C 背景 想使用 Matlab 中的类 句柄和值类有什么区别 我知道如果我想定义一个带有重载运算符 例如 和 的矩阵类 我会使用值类 然而 有时 当我选择一个手柄类时 事情似乎只对我有用 MathWorks 提供了一些有关其用途的
  • Matlab 字段名索引[重复]

    这个问题在这里已经有答案了 所以我有一个包含多个表的元胞数组 我试图访问表的第一个列名称 c table1 table2 table3 以下两行都给了我错误 fieldnames c 1 1 fieldnames c 1 1 Error i
  • Matlab没有优化以下内容吗?

    我有一个很长的向量 1xrv 和一个很长的向量w1xs 和一个矩阵Arxs 它是稀疏的 但维度非常大 我期望 Matlab 对以下内容进行优化 这样我就不会遇到内存问题 A v w 但看起来 Matlab 实际上是在尝试生成完整的v w矩阵
  • MATLAB 特征函数

    我很好奇哪里可以找到完整的描述FEATURE功能 它接受哪些论点 没有找到文档 我只听说过memstats and getpid 还要别的吗 gt gt which feature built in undocumented 注意 更完整的
  • MATLAB 教程中的 SIFT 实现

    我正在寻找 MATLAB 中的一些基本 SIFT 实现 我需要从第一原则来写它 另外 我正在寻找一些可以解释程序中发生的事情的内容 Vedali 的代码和 David Lowe 的代码超出了我的理解范围 如果您是 Matlab 用户 您一定
  • 如何在 Matlab 中将数组打印到 .txt 文件?

    我才刚刚开始学习Matlab 所以这个问题可能非常基本 我有一个变量 a 2 3 3 422 6 121 9 4 55 我希望将值输出到 txt 文件 如下所示 2 3 3 422 6 121 9 4 55 我怎样才能做到这一点 fid f
  • Matlab:如何更改矩阵的存储方式?从 1x1x3 到 1x3?

    我目前有 val 1 0 7216 val 2 0 7216 val 3 0 7216 但我想要 0 7216 0 716 0 721 我可以做什么样的操作来做到这一点 The reshape函数将在这里解决问题 Arrange the e
  • 如何告诉 mex 链接到 /usr/lib 中的 libstdc++.so.6 而不是 MATLAB 目录中的 libstdc++.so.6?

    现在 MATLAB 2012a 中的 mex 仅正式支持 gcc 4 4 6 但我想使用 gcc 4 7 风险自负 现在如果我直接用 mex 编译一些东西 它会抱怨 usr lib gcc i686 linux gnu 4 7 cc1plu
  • matlab中的排列函数是如何工作的

    这是一个有点愚蠢的问题 但我似乎无法弄清楚排列在 matlab 中是如何工作的 以文档为例 A 1 2 3 4 permute A 2 1 ans 1 3 2 4 到底是怎么回事 这如何告诉 matlab 3 和 2 需要交换 哇 这是我迄
  • 如何更改Plotyy第二轴的颜色和字体大小?

    我使用 MATLAB 的plotyy 函数绘制了两条曲线 AX H1 H2 plotyy voltage span amplitude voltage span Ca SR The problem is that I cannot chan
  • Matlab 的 imresize 函数中用于插值的算法是什么?

    我正在使用 Matlab Octaveimresize 对给定的二维数组重新采样的函数 我想了解如何使用特定的插值算法imresize works 我在Windows上使用八度 e g A 1 2 3 4 是一个二维数组 然后我使用命令 b
  • 检测数据集中线性行为的算法

    我已经发布了一个关于对数据集的一部分进行多项式拟合的算法 https stackoverflow com q 17595932 2320757前一段时间收到一些建议去做我想做的事 但我现在面临另一个问题 我尝试应用答案中建议的想法 我的目标

随机推荐