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

个人主页:Matlab科研工作室

个人信条:格物致知。

更多Matlab仿真内容点击

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

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

内容介绍

将储能资源纳入电力现货市场有助于优化电力系统灵活性配置,促进可再生能源消纳,推动能源革命及低碳化转型.考虑到储能资源具有较好的调频性能,该文主要探讨包含独立储能的现货电能量与调频辅助服务市场出清协调机制.首先,总结当前国内外典型市场中现货电能量与调频辅助服务市场出清协调机制,分析各机制出清流程及定价方法等特点.其次,提出独立储能参与市场的特殊问题及解决思路,包括申报参数的选取,实时市场优化的处理,提供调频辅助服务机会成本的计算等,并给出包含独立储能的现货电能量与调频辅助服务市场顺次出清及联合出清两种机制的数学模型.最后,构建算例,分析对比各机制的实用性和有效性,讨论独立储能对市场电价及购电成本的影响,独立储能参与市场的效益及行为,独立储能申报及出清中相关参数对市场出清结果的影响等.

独立储能的现货电能量与调频辅助服务市场出清协调机制附Matlab代码_开发者

部分代码

function [fit_record,result]=Apso(option,data)

%% 开始计算

creat_x=option.creat_x;

fitness=option.aimFcn;

%% 初始种群生成 

% x=cell(option.M ,1);

% fit=zeros(option.M ,1);

% for i=1:option.M 

%     x{i}=creat_x(option,data);

%     [m,n]=size(x{i});

%     v{i}=2*rand(m,n)-1;

%     fit(i)=fitness(x{i},option,data);

% end

x=data.x;

fit=data.fit;

v=data.v;

%% 开始计算

%% 粒子群算法求解

[fit,position]=sort(fit);

% 记录种群位置

x=x(position);

fit_p=fit;

x_p=x;

best_fit_p=fit_p(1);

best_x_p=x_p(1);

temp_x=x;

fit_record(1)=best_fit_p;

% 迭代计算

for t=1:option.peiod

    t

    % 粒子更新

   %  option.c1=2-sin((t.*pi)/option.peiod);%这是学习因子的公式

   %  option.c2=1+sin((t.*pi)/option.peiod);

    %w = option.wmax-(t-1)*(option.wmax-option.wmin)/(option.peiod-1); %权重线性递减

    w = option.wmin+exp(-5*t/option.peiod)^2*(option.wmax-option.wmin);

    parfor i=1:option.M

        temp_v{i}= w*v{i} + option.c1*rand(size(x_p{i})).*(x_p{i}-x{i}) + option.c2*rand(size(x_p{i})).*(best_x_p{1} - x{i});%速度更新

        temp_x{i}=x{i} + temp_v{i};%位置更新

        temp_x{i}=test_x(temp_x{i},option,data);

        temp_fit(i)=fitness(temp_x{i},option,data);

    end

    [best_fit_p_1,position]=min(temp_fit);

    % 更新最优的粒子群位置

    if best_fit_p_1<best_fit_p

        best_fit_p=best_fit_p_1;

        best_x_p=temp_x(position);

    end

    for i=1:option.M

        if fit_p(i)>temp_fit(i)

            fit_p(i)=temp_fit(i);

            x_p(i)=temp_x(i);

            v{i}=temp_v{i};

        end

    end

    x=temp_x;

    v=temp_v;

    fit=temp_fit;

    fit_record(t+1)=best_fit_p;

    fit_record_mean(t+1)=mean(fit_p);

end

%%

result.best_x=best_x_p{1};

[result.best_fit,result.detail]=fitness(result.best_x,option,data);

end

运行结果

独立储能的现货电能量与调频辅助服务市场出清协调机制附Matlab代码_开发者_02

参考文献

[1]肖云鹏, 张兰, 张轩,等. 包含独立储能的现货电能量与调频辅助服务市场出清协调机制[J]. 中国电机工程学报, 2020.

[2]林阿竹, 柯清辉, 江岳文. 独立储能参与调频辅助服务市场机制设计[J]. 电力自动化设备, 2022, 42(12):9.

[3]王杰, 方日升, 温步瀛. 风储联合系统参与能量市场和调频辅助服务市场协同优化[J]. 低压电器, 2019.

️ 完整代码

️部分理论引用网络文献,若有侵权联系博主删除
️ 关注我领取海量matlab电子书和数学建模资料