基于Copula函数的风光联合场景生成与K-means聚类削减的MATLAB实现.pdf

基于Copula函数与K-means聚类的风光联合场景生成及削减研究（MATLAB实现）

最近在搞风光联合出力场景生成，发现个有意思的事儿——大伙儿都在说风光互补，但实操时候总把

风电光伏当两个独立变量处理。这不科学啊！风机叶片转得欢的时候，旁边光伏板可能正被同一片云挡着

光呢。咱今天就用MATLAB整点真实的，让风光出力牵着手出场。

先搞个高斯copula建立相关性模型。这玩意儿就像给风电和光伏套上隐形橡皮筋，既保持各自的出

力特性，又能体现空间上的联动效应。看代码：

```matlab

% 设置copula参数

rho = 0.7; % 相关系数

solar_scenes = norminv(u(:,2), solar_mu, solar_sigma);

```

这里用了正态逆变换，把copula生成的均匀分布变量映射到实际风光分布。注意rho参数控制着二

者的基情程度，调成0.3立马变成塑料兄弟情。

接着把这对好基友打包处理：

```matlab

joint_scenes = [wind_scenes, solar_scenes];

[cluster_idx, centroids] = kmeans(joint_scenes, 5, 'MaxIter',1000,'Replicates',10);

```

K-means这老伙计干活得哄着，多给几次初始值（Replicates）能避免掉进局部最优的坑。这里直接