基于二元Frank-Copula函数的风光出力场景生成方法附Matlab代码

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

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

🍊个人信条：格物致知。

⛄ 内容介绍

新能源的随机性,波动性及间歇性为电力系统规划带来困扰,对风,光出力的变化规律进行合理刻画,生成典型出力场景是新能源规划的常用方法.

⛄ 部分代码

function [sys,x0,str,ts] = spacemodel(t,x,u,flag)

switch flag,

case 0,

    [sys,x0,str,ts]=mdlInitializeSizes;

case 3,

    sys=mdlOutputs(t,x,u);

case {2,4,9}

    sys=[];

otherwise

    error(['Unhandled flag = ',num2str(flag)]);

end

function [sys,x0,str,ts]=mdlInitializeSizes

sizes = simsizes;

sizes.NumContStates  = 0;

sizes.NumDiscStates  = 0;

sizes.NumOutputs     = 4;

sizes.NumInputs      = 0;

sizes.DirFeedthrough = 1;

sizes.NumSampleTimes = 1;

sys = simsizes(sizes);

x0  = [];

str = [];

ts  = [0 0];

function sys=mdlOutputs(t,x,u)

q1d=sin(3*t);

dq1d=3*cos(3*t);

q2d=cos(3*t);

dq2d=-3*sin(3*t);

sys(1)=q1d;

sys(2)=dq1d;

sys(3)=q2d;

sys(4)=dq2d;

⛄ 运行结果

⛄ 参考文献

[1]段偲默等. "一种基于Copula函数的风光联合出力建模方法.", CN109038648A. 2018.

[2]别朝红, 潘超琼, 王灿,等. 一种基于Copula函数的时空相关风光序列模拟方法:. 

⛳️ 完整代码

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

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