matpower 里的一些问题

sss3434

matpower搭建了个模型，修改case9文件中的 gen 数据，把全局变量windpower风电出力导入到gen(2,pg) ，建立风电接入模型。
function [compResults,PG,results] =  opfWind
global windPower;
wind = xlsread('wind.xls','I2:I8756');
) o) q2 L' j2 ]; S' @- \$ vcompResults = zeros(8755,1);
PG = zeros(8755,1);
' z+ B9 u; l  S9 V% 循环进行潮流计算
  P+ Z7 r* e8 X/ \
for i = 1:length(wind)
    windPower = wind(i,1)/100;
    loadcase('jycase9');
    results = runpf('jycase9');
    compResults(i,1) = results.branch(7,14);
( [* H* _0 ]" I6 y  F, i    PG(i,1) = results.branch(7,16);
end
( o) f. ~0 ^& C0 m( T, G6 j. M- q
( H% _' A) w1 ]8 A%*******************************核密度估计**********************************
% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
U1 = ksdensity(compResults,compResults,'function','cdf','width',9);
9 |" v7 A* `+ x  W1 R5 N; v- I" GV1 = ksdensity(PG,PG,'function','cdf');
! x6 @! G. c2 }; L% **********************绘制核密度估计图**********************
( `& K& a  q  P, m( G[Xsort,id] = sort(compResults);  % 为了作图的需要，对X进行排序
figure;  % 新建一个图形窗口
! l8 n% _- H: t8 n4 R0 X" whold on
plot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
! N( O8 w2 `- r4 j( u8 |legend('核密度估计', 'Location','NorthWest'); % 加标注框
/ N5 s( L9 x7 M3 Z5 x( c# Axlabel('节点8的电压幅值');  % 为X轴加标签
/ Y4 |. H) c: {+ x0 w4 w  r& y5 cylabel('F(x)');  % 为Y轴加标签
- t6 J1 `+ B. o7 l但是输出的累积分布图（cdf)却没趋向于1，不知道哪里出了问题，是不是matpower里的程序的问题。