matpower 里的一些问题

sss3434

matpower搭建了个模型，修改case9文件中的 gen 数据，把全局变量windpower风电出力导入到gen(2,pg) ，建立风电接入模型。
$ p$ Y9 l6 Y0 Yfunction [compResults,PG,results] =  opfWind
: }. W, i  L  ?global windPower;
9 q" }" i3 p2 f8 s2 qwind = xlsread('wind.xls','I2:I8756');
compResults = zeros(8755,1);
PG = zeros(8755,1);
% B9 Z/ g) E' @, S9 E. R( ]% 循环进行潮流计算
) `* T$ N, l( q2 n/ x1 f! L  A) }
for i = 1:length(wind)
$ q) I# v9 l" f7 e3 T6 k/ B    windPower = wind(i,1)/100;
    loadcase('jycase9');
    results = runpf('jycase9');
. n! H* X! T( M! x% Q/ k8 S1 O6 X    compResults(i,1) = results.branch(7,14);
    PG(i,1) = results.branch(7,16);
end

%*******************************核密度估计**********************************
% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
U1 = ksdensity(compResults,compResults,'function','cdf','width',9);
V1 = ksdensity(PG,PG,'function','cdf');
' X$ ^6 {  f: k0 E% **********************绘制核密度估计图**********************
2 [5 Z" Z! m& Q4 O+ {" H: V) V* W4 O) L# D[Xsort,id] = sort(compResults);  % 为了作图的需要，对X进行排序
# P& ?1 A/ W- b2 X4 Dfigure;  % 新建一个图形窗口
hold on
plot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
" S6 [8 t. E) n: nlegend('核密度估计', 'Location','NorthWest'); % 加标注框
xlabel('节点8的电压幅值');  % 为X轴加标签
ylabel('F(x)');  % 为Y轴加标签
/ ?& l! T! C9 `  \5 W+ H) U但是输出的累积分布图（cdf)却没趋向于1，不知道哪里出了问题，是不是matpower里的程序的问题。