matpower 里的一些问题

sss3434

matpower搭建了个模型，修改case9文件中的 gen 数据，把全局变量windpower风电出力导入到gen(2,pg) ，建立风电接入模型。
function [compResults,PG,results] =  opfWind
global windPower;
9 \( P. h6 j$ E3 p/ c0 A# v1 Mwind = xlsread('wind.xls','I2:I8756');
compResults = zeros(8755,1);
+ {' i7 h0 _- W0 C/ t* BPG = zeros(8755,1);
( x+ G# |# ~. B7 Q  I% u% 循环进行潮流计算

: \- D; I* A* nfor i = 1:length(wind)
! t& u  ^$ H6 O4 b. ~  t0 C    windPower = wind(i,1)/100;
    loadcase('jycase9');
    results = runpf('jycase9');
    compResults(i,1) = results.branch(7,14);
    PG(i,1) = results.branch(7,16);
end
0 V) x$ @$ m, r0 J0 u+ H
: H3 m4 O, s: F( K5 N# r2 _% K/ P%*******************************核密度估计**********************************
% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
4 x  G3 T4 |3 Q2 B0 QU1 = ksdensity(compResults,compResults,'function','cdf','width',9);
V1 = ksdensity(PG,PG,'function','cdf');
% **********************绘制核密度估计图**********************
[Xsort,id] = sort(compResults);  % 为了作图的需要，对X进行排序
figure;  % 新建一个图形窗口
hold on
plot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
9 I3 w. V6 k! g' A+ m7 ?) |legend('核密度估计', 'Location','NorthWest'); % 加标注框
" E: U1 x  D+ d+ q* Txlabel('节点8的电压幅值');  % 为X轴加标签
. L: Q# {/ y  |* E1 jylabel('F(x)');  % 为Y轴加标签
( @+ y$ P( O4 l7 J. P但是输出的累积分布图（cdf)却没趋向于1，不知道哪里出了问题，是不是matpower里的程序的问题。