matpower 里的一些问题

sss3434

matpower搭建了个模型，修改case9文件中的 gen 数据，把全局变量windpower风电出力导入到gen(2,pg) ，建立风电接入模型。
9 p; G3 W" \+ q0 w$ c$ Ofunction [compResults,PG,results] =  opfWind
* E6 I: z, }. H) y, M% zglobal windPower;
" W, I- {1 z! x+ B7 W6 Q3 ~wind = xlsread('wind.xls','I2:I8756');
8 o8 }1 V3 J) `0 r% g+ ~$ O& t! V! s1 RcompResults = zeros(8755,1);
. q- b' K2 G( c/ q6 f: nPG = zeros(8755,1);
4 T) J7 {: C6 o: C: ^7 e% 循环进行潮流计算
8 t- p% |4 O; ]8 ^
  C5 I0 B/ ]. i& Zfor i = 1:length(wind)
    windPower = wind(i,1)/100;
) D# [$ d7 ]# V5 v* y6 w9 d    loadcase('jycase9');
6 p* o, t) N7 K+ _$ O! Q    results = runpf('jycase9');
+ A, ?0 B0 r/ V6 c& c* t, N! O    compResults(i,1) = results.branch(7,14);
    PG(i,1) = results.branch(7,16);
end

9 z! i+ k9 z4 s" G3 \2 D' `%*******************************核密度估计**********************************
% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
6 L9 Z2 a1 N3 q+ i4 h; I! W  aU1 = ksdensity(compResults,compResults,'function','cdf','width',9);
9 q$ A! t  c' ^V1 = ksdensity(PG,PG,'function','cdf');
; O. }* C# J* U* i) f% **********************绘制核密度估计图**********************
7 t' a. M* x7 l2 b2 V0 K. p6 w! b[Xsort,id] = sort(compResults);  % 为了作图的需要，对X进行排序
: l0 N4 Y) M$ k! A6 R, A( I( Ufigure;  % 新建一个图形窗口
hold on
6 [5 T  W; x' ?5 G; w) ^plot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
2 S9 s) K, d7 u) S" p" _legend('核密度估计', 'Location','NorthWest'); % 加标注框
7 i9 Y! P( u; F  s* u7 kxlabel('节点8的电压幅值');  % 为X轴加标签
# @* w. u- ~$ ]# Z" ^ylabel('F(x)');  % 为Y轴加标签
/ `; H: A6 g9 q但是输出的累积分布图（cdf)却没趋向于1，不知道哪里出了问题，是不是matpower里的程序的问题。