|
|
马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
matpower搭建了个模型,修改case9文件中的 gen 数据,把全局变量windpower风电出力导入到gen(2,pg) ,建立风电接入模型。
7 H a2 C Z. r% w% \function [compResults,PG,results] = opfWind
/ n# f9 p+ W& m. w9 J( }global windPower;4 f: O2 m, r# F- `
wind = xlsread('wind.xls','I2:I8756');
" \0 ~4 c2 \% u% \9 Z; YcompResults = zeros(8755,1);) n, j" v5 C0 Z" V
PG = zeros(8755,1);
- w4 `- f" v8 I% 循环进行潮流计算
; c6 ?4 i H( m8 t/ x
: l) W9 v1 A8 o6 i/ R" W% C" a% Jfor i = 1:length(wind)/ y: i6 f9 v6 w8 O: \
windPower = wind(i,1)/100;
) g* A4 `) ^$ c loadcase('jycase9');
0 @$ X7 j/ Q; ]5 Z results = runpf('jycase9');
" t: V$ c! X. r3 F p# R9 L compResults(i,1) = results.branch(7,14);4 R$ \: Q5 ^" ~% S2 t& R
PG(i,1) = results.branch(7,16);
& x8 V* h% q1 C {end* V, j& |9 T9 K! d
# U W4 T% M; n' b% ^
%*******************************核密度估计**********************************
8 i. Z' U* s# h1 H% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
3 q+ b- b6 c/ B2 CU1 = ksdensity(compResults,compResults,'function','cdf','width',9);( ^! P5 k- x* M9 A& G8 w% D
V1 = ksdensity(PG,PG,'function','cdf');
3 W E$ t: ]5 U! J% **********************绘制核密度估计图**********************6 M3 S0 j; N5 z; g Q3 p! X' p
[Xsort,id] = sort(compResults); % 为了作图的需要,对X进行排序
6 Y( H% F- P! b$ j5 D/ A. t/ Xfigure; % 新建一个图形窗口
$ f+ o) r; b' j- ^- Whold on
. e7 b7 a6 m$ c* L3 ? Gplot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图% D2 W9 ~% k* C
legend('核密度估计', 'Location','NorthWest'); % 加标注框
6 h* m% Q! Q9 pxlabel('节点8的电压幅值'); % 为X轴加标签" b F* m! `4 o" e
ylabel('F(x)'); % 为Y轴加标签
( V" C; H+ q, S1 y; `但是输出的累积分布图(cdf)却没趋向于1,不知道哪里出了问题,是不是matpower里的程序的问题。 |
|