电力系统随机潮流概率潮流计算MATLAB程序包含 蒙特卡洛模拟法、半不变量法+级数展开(Gram-Charlie,Cornish-Fisher); 考虑光伏不确定性(Beta分布)负荷不确定性(正态分布),以IEEE34节点为例,计算节点电压、支路潮流概率密度、累计概率并绘制曲线。 有注释,附带参考文献,直观清晰。
可分别做某节点电压幅值(某支路有功无功)概率密度或累计分布曲线,也可以将所有节点电压(所有支路有功无功)概率密度曲线或累计分布曲线。蒙特卡洛与半步变量法做对比图。
该程序每一部分都含有子程序,相应的计算过程及解释都包含在程序中,以及参考论文对应的页码。
%miu=0.150314263;
%sita=0.049758487;
M=textread('Guangzhou.txt');
maxM=max(M(:,1));
miuM=mean(M(:,1)./maxM); %最大值归一化
sitaM=sqrt(var(M(:,1)./maxM));
a=miuM^2*(1-miuM)/sitaM^2-miuM; %beta分布参数计算
b=miuM*(1-miuM)^2/sitaM^2-(1-miuM);
nbins=30; %设置直方图需要绘制的区间数
Max=max(M(:,1));
Min=min(M(:,1));
Length=(Max-Min)/nbins;
[nhist,c_points]=hist(M(:,1)/maxM,nbins); %获取所绘制的频数直方图各区间的中心以及各区间对应的元素的个数
figure
x=0:0.001:1;
y1=betapdf(x,a,b);
y2=betacdf(x,a,b);
subplot(1,2,1);
bar(c_points,nhist/size(M,1)/Length,'grouped','white');
hold on;
plot(x,y1,'LineWidth',2);
hold on;
title('Beta分布的概率密度函数','FontSize',16);
grid on;
box on;
subplot(1,2,2);
plot(x,y2,'LineWidth',2);
title('Beta分布的累积分布函数','FontSize',16);
grid on;
box on;
多图及多图部分如下:
代码可在下面(V:K25_Taurus ⑧这)链接获得