matpower 里的一些问题

sss3434

matpower搭建了个模型，修改case9文件中的 gen 数据，把全局变量windpower风电出力导入到gen(2,pg) ，建立风电接入模型。
function [compResults,PG,results] =  opfWind
: {6 ]5 Q" q0 O1 m0 H; n9 w4 cglobal windPower;
wind = xlsread('wind.xls','I2:I8756');
1 n$ B5 o! o( F/ x; A0 n$ McompResults = zeros(8755,1);
PG = zeros(8755,1);
5 _2 I; [+ ^3 R; `& Z% 循环进行潮流计算
) U' U- Z# |5 t% o) _
$ p7 z$ m1 h' B' r# i, ]$ Ffor i = 1:length(wind)
    windPower = wind(i,1)/100;
7 J0 b- C0 P0 f& v    loadcase('jycase9');
2 z3 g; \9 c) q  t    results = runpf('jycase9');
# G- Y$ Y6 s0 C( d* n    compResults(i,1) = results.branch(7,14);
: P! ~4 |  Z  B( f$ x    PG(i,1) = results.branch(7,16);
end
( x* Q1 s8 t8 e9 s( C
%*******************************核密度估计**********************************
% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
U1 = ksdensity(compResults,compResults,'function','cdf','width',9);
V1 = ksdensity(PG,PG,'function','cdf');
% **********************绘制核密度估计图**********************
[Xsort,id] = sort(compResults);  % 为了作图的需要，对X进行排序
7 m6 R! c: r6 i" ufigure;  % 新建一个图形窗口
$ ]! M" f" \5 X8 S( B" k1 ohold on
plot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
  V  U0 Q! {7 Q7 }0 D" B" r' [3 blegend('核密度估计', 'Location','NorthWest'); % 加标注框
. E1 ?5 N  Q1 Y' ^xlabel('节点8的电压幅值');  % 为X轴加标签
ylabel('F(x)');  % 为Y轴加标签
7 s, V1 {6 z5 U$ @5 {但是输出的累积分布图（cdf)却没趋向于1，不知道哪里出了问题，是不是matpower里的程序的问题。