matpower 里的一些问题

sss3434

matpower搭建了个模型，修改case9文件中的 gen 数据，把全局变量windpower风电出力导入到gen(2,pg) ，建立风电接入模型。
6 ]0 l+ c- |, ?' M. Pfunction [compResults,PG,results] =  opfWind
global windPower;
wind = xlsread('wind.xls','I2:I8756');
compResults = zeros(8755,1);
2 b9 F: d% a* X* o" X( [$ APG = zeros(8755,1);
/ r1 H1 F3 D/ P  U3 U9 ]% 循环进行潮流计算

for i = 1:length(wind)
# K3 C9 C% x/ |+ R    windPower = wind(i,1)/100;
6 B' A) s* i, [; y    loadcase('jycase9');
9 G: n5 [% ?( o$ f7 X* h3 O    results = runpf('jycase9');
    compResults(i,1) = results.branch(7,14);
    PG(i,1) = results.branch(7,16);
end

%*******************************核密度估计**********************************
% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
' j* d! V! v3 O8 a* a3 E! vU1 = ksdensity(compResults,compResults,'function','cdf','width',9);
V1 = ksdensity(PG,PG,'function','cdf');
3 m; f' C7 e2 d8 @% **********************绘制核密度估计图**********************
[Xsort,id] = sort(compResults);  % 为了作图的需要，对X进行排序
* w& m4 ~6 Q7 F& f# vfigure;  % 新建一个图形窗口
; u7 K: E) M0 `/ x, Ghold on
plot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
' |+ A5 c& T2 rlegend('核密度估计', 'Location','NorthWest'); % 加标注框
, X; {! R: v, y+ m$ `; p7 G' nxlabel('节点8的电压幅值');  % 为X轴加标签
2 }0 I0 p6 B9 Vylabel('F(x)');  % 为Y轴加标签
但是输出的累积分布图（cdf)却没趋向于1，不知道哪里出了问题，是不是matpower里的程序的问题。