matpower 里的一些问题

sss3434

matpower搭建了个模型，修改case9文件中的 gen 数据，把全局变量windpower风电出力导入到gen(2,pg) ，建立风电接入模型。
function [compResults,PG,results] =  opfWind
global windPower;
wind = xlsread('wind.xls','I2:I8756');
; [2 x) q$ p4 G& p& DcompResults = zeros(8755,1);
/ }9 A/ ^9 ?. g3 Z* ePG = zeros(8755,1);
% 循环进行潮流计算

, T7 u3 O, ]$ d) _for i = 1:length(wind)
! Q# b3 R* U. J. p* ^5 \    windPower = wind(i,1)/100;
9 d% W# ^  W6 C" @3 a! l, c    loadcase('jycase9');
7 K# Y0 h! I: R* j8 F( C: R1 J    results = runpf('jycase9');
0 h, j0 t+ m8 n7 o, w4 I( E2 M    compResults(i,1) = results.branch(7,14);
    PG(i,1) = results.branch(7,16);
end

%*******************************核密度估计**********************************
% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
# p# s5 Q- p' G( G, f* X# jU1 = ksdensity(compResults,compResults,'function','cdf','width',9);
V1 = ksdensity(PG,PG,'function','cdf');
% **********************绘制核密度估计图**********************
0 h: x+ I; n' I. h8 `[Xsort,id] = sort(compResults);  % 为了作图的需要，对X进行排序
1 q3 L0 X+ _! x' ifigure;  % 新建一个图形窗口
hold on
6 ^7 X- Q, v( K0 [# L  dplot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
legend('核密度估计', 'Location','NorthWest'); % 加标注框
xlabel('节点8的电压幅值');  % 为X轴加标签
ylabel('F(x)');  % 为Y轴加标签
但是输出的累积分布图（cdf)却没趋向于1，不知道哪里出了问题，是不是matpower里的程序的问题。