|
|
马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
matpower搭建了个模型,修改case9文件中的 gen 数据,把全局变量windpower风电出力导入到gen(2,pg) ,建立风电接入模型。7 }! j: }8 k, \9 B9 q
function [compResults,PG,results] = opfWind
0 L% i) t* {7 y8 ~/ l9 C4 z% _global windPower;, T9 {3 W; p3 v$ W# K3 m% t
wind = xlsread('wind.xls','I2:I8756');
+ G1 y% ?7 i! @3 Z: g6 pcompResults = zeros(8755,1);5 c- _, {2 T% X9 F' O) q, _
PG = zeros(8755,1);
% {% ?' I2 C5 a- P( {6 F' Q3 A% 循环进行潮流计算' Z' s2 o+ D/ L
# ~+ B8 G. A* N, G9 e3 ffor i = 1:length(wind)+ X) ~; J9 q! Q4 M. Q3 C: ]
windPower = wind(i,1)/100;6 A: }! n- Q( E, J$ x
loadcase('jycase9'); N7 L+ J5 @" ^5 M
results = runpf('jycase9');& O) V! M2 K! G ^3 d& D2 `
compResults(i,1) = results.branch(7,14);
+ u% f4 N& Y" f, m& O PG(i,1) = results.branch(7,16);
9 X' Z& J# ~6 o# d& f( T% Dend5 a: K' ~$ A* J$ K: E3 B% c) [
, Q4 Z5 y" {! k) N) o%*******************************核密度估计**********************************% I' ^; J* B$ h! I% q/ t
% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
E9 o- J; p0 ]1 ?( k9 f+ E! {+ W9 iU1 = ksdensity(compResults,compResults,'function','cdf','width',9);
' v2 b% ^1 K% xV1 = ksdensity(PG,PG,'function','cdf');
& M3 V/ f* }& ` c4 S4 o/ [% **********************绘制核密度估计图**********************4 j+ ^" C" R" |9 x
[Xsort,id] = sort(compResults); % 为了作图的需要,对X进行排序
+ i% E0 B; T& ~/ J% pfigure; % 新建一个图形窗口
: r2 j' A( z7 L, {hold on
' o& e2 ] d/ S. h, Q8 C. {plot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
G# o. ] U; Glegend('核密度估计', 'Location','NorthWest'); % 加标注框; b9 v D! t; B( h/ c
xlabel('节点8的电压幅值'); % 为X轴加标签3 [$ m6 {+ V9 V4 c: j" O
ylabel('F(x)'); % 为Y轴加标签
; g* T' Z1 r) C但是输出的累积分布图(cdf)却没趋向于1,不知道哪里出了问题,是不是matpower里的程序的问题。 |
|