matpower 里的一些问题

sss3434

matpower搭建了个模型，修改case9文件中的 gen 数据，把全局变量windpower风电出力导入到gen(2,pg) ，建立风电接入模型。
/ C3 y- e' K+ _function [compResults,PG,results] =  opfWind
# L( m& d- M" F3 n; N; oglobal windPower;
wind = xlsread('wind.xls','I2:I8756');
compResults = zeros(8755,1);
PG = zeros(8755,1);
9 ?! U. I  P  i1 N" P% 循环进行潮流计算

$ p% J; t: U' B5 `$ C( c3 M! lfor i = 1:length(wind)
/ V" P" X, Q( E, Q5 x7 d( \/ }    windPower = wind(i,1)/100;
    loadcase('jycase9');
: l; i1 ~; O# k, }& N; n  ?+ K! s    results = runpf('jycase9');
" h/ P  {( ^0 J: R6 J7 A    compResults(i,1) = results.branch(7,14);
& |8 y3 v( d1 l    PG(i,1) = results.branch(7,16);
end

%*******************************核密度估计**********************************
% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
. C9 S! e. ]) }9 [/ |  YU1 = ksdensity(compResults,compResults,'function','cdf','width',9);
: o' n/ u4 v1 t1 e+ O6 F! ZV1 = ksdensity(PG,PG,'function','cdf');
6 ^. e+ B3 Q7 s, b( X% **********************绘制核密度估计图**********************
[Xsort,id] = sort(compResults);  % 为了作图的需要，对X进行排序
2 I- s7 s4 B. \1 d( Z! g! T; n$ Efigure;  % 新建一个图形窗口
hold on
: O* N7 H0 ?- B+ `plot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
$ z  _* C* T! v7 j3 m4 H. Alegend('核密度估计', 'Location','NorthWest'); % 加标注框
1 z/ A1 a; Y( y) O  u* W/ O9 k5 ^3 y4 Gxlabel('节点8的电压幅值');  % 为X轴加标签
ylabel('F(x)');  % 为Y轴加标签
; T, m6 F: M+ z: L. s但是输出的累积分布图（cdf)却没趋向于1，不知道哪里出了问题，是不是matpower里的程序的问题。