|
马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
matpower搭建了个模型,修改case9文件中的 gen 数据,把全局变量windpower风电出力导入到gen(2,pg) ,建立风电接入模型。
1 o- P3 _, |+ Ofunction [compResults,PG,results] = opfWind" \2 c9 v2 W4 R' n
global windPower;
, z( z% Z4 a# j( A8 ?5 jwind = xlsread('wind.xls','I2:I8756');8 n$ G! U: W( Y# T' Y
compResults = zeros(8755,1);& I4 a T, j% I4 S' d
PG = zeros(8755,1);
9 c( n! b8 l# {- I7 D1 h% 循环进行潮流计算
' `1 ~& V8 F+ I0 A# e0 [5 U! [3 i0 P8 f8 r/ Q1 w- ^/ b) w
for i = 1:length(wind)
1 `: V& [) p' T windPower = wind(i,1)/100;7 O. o/ w3 d9 d& G
loadcase('jycase9');
) W4 P3 n8 \# R$ F7 ]* j results = runpf('jycase9');9 D+ ?( e* f& a; O# H1 {; k
compResults(i,1) = results.branch(7,14);/ C8 ?7 x* e1 R8 |- ?1 t# P
PG(i,1) = results.branch(7,16);. r: O+ l' U7 h7 H2 ?8 A
end
! j! f- O4 X, q& D n1 V1 C i3 Q2 b6 ]: B
%*******************************核密度估计**********************************3 \- M5 M9 b$ d5 e2 C$ c! j
% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
% _* B9 T; G( U" h, ?6 QU1 = ksdensity(compResults,compResults,'function','cdf','width',9);' T. u& U6 u$ o
V1 = ksdensity(PG,PG,'function','cdf');
' _% u3 o* Z) x% **********************绘制核密度估计图**********************
& t8 M8 r h6 p[Xsort,id] = sort(compResults); % 为了作图的需要,对X进行排序8 E+ J; t4 V( p# [" `
figure; % 新建一个图形窗口" o Z# R. J9 S3 U5 K- h2 ^
hold on
4 g5 B/ C, u/ `' z3 g5 rplot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
7 b- a8 g# Q* B* Elegend('核密度估计', 'Location','NorthWest'); % 加标注框* }9 _5 B& O! Z1 c! ^: a! i
xlabel('节点8的电压幅值'); % 为X轴加标签6 A6 d' G5 p1 T0 a
ylabel('F(x)'); % 为Y轴加标签
s: F9 y- J2 H+ t但是输出的累积分布图(cdf)却没趋向于1,不知道哪里出了问题,是不是matpower里的程序的问题。 |
|