matpower 里的一些问题

sss3434

matpower搭建了个模型，修改case9文件中的 gen 数据，把全局变量windpower风电出力导入到gen(2,pg) ，建立风电接入模型。
function [compResults,PG,results] =  opfWind
global windPower;
wind = xlsread('wind.xls','I2:I8756');
compResults = zeros(8755,1);
" n) W# v) @# n& `PG = zeros(8755,1);
% 循环进行潮流计算
* j) ~$ k: t% ]9 v/ Q7 k+ X
for i = 1:length(wind)
    windPower = wind(i,1)/100;
* }  T7 N; U+ Y, a* Y, s    loadcase('jycase9');
    results = runpf('jycase9');
) x  i! B9 a3 N/ q7 J3 e    compResults(i,1) = results.branch(7,14);
" |$ V- ?4 c. f! }9 n: G    PG(i,1) = results.branch(7,16);
' d$ n4 I6 e* g% D- }) Iend

%*******************************核密度估计**********************************
  V' {2 k$ V2 g: e$ C  f; M% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
  y6 ^0 |1 z  T2 C' EU1 = ksdensity(compResults,compResults,'function','cdf','width',9);
2 j9 I: `4 \5 Y: L$ W5 tV1 = ksdensity(PG,PG,'function','cdf');
! S: W4 |  Y  W# U5 m* K2 P% **********************绘制核密度估计图**********************
. E( y" p7 Z0 s" r$ [3 v$ Z2 Z# c9 @[Xsort,id] = sort(compResults);  % 为了作图的需要，对X进行排序
figure;  % 新建一个图形窗口
hold on
plot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
- L- N! |3 L% {legend('核密度估计', 'Location','NorthWest'); % 加标注框
7 J3 B/ u( [2 |& w" Jxlabel('节点8的电压幅值');  % 为X轴加标签
4 I( \+ q; z. `& |/ f# a3 }ylabel('F(x)');  % 为Y轴加标签
7 Q" p8 A/ K4 b但是输出的累积分布图（cdf)却没趋向于1，不知道哪里出了问题，是不是matpower里的程序的问题。