matpower 里的一些问题

sss3434

matpower搭建了个模型，修改case9文件中的 gen 数据，把全局变量windpower风电出力导入到gen(2,pg) ，建立风电接入模型。
function [compResults,PG,results] =  opfWind
& P' U, e& X) qglobal windPower;
6 m6 R1 ]' ?( t6 a- bwind = xlsread('wind.xls','I2:I8756');
compResults = zeros(8755,1);
3 d+ T) Y; m- @. gPG = zeros(8755,1);
% 循环进行潮流计算

1 {$ g) k$ U& _; D) {% W" qfor i = 1:length(wind)
$ p4 r3 d$ s7 G) T+ Y- H+ {4 o    windPower = wind(i,1)/100;
    loadcase('jycase9');
    results = runpf('jycase9');
    compResults(i,1) = results.branch(7,14);
. c3 w! M5 a" G; c9 e    PG(i,1) = results.branch(7,16);
2 v$ {/ a/ i% [2 [' Kend
# u: G) D# M  S  W( C# H
%*******************************核密度估计**********************************
2 M! e/ U2 n6 x3 \) ?% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
% S* g) y% @' T! h+ h. ]U1 = ksdensity(compResults,compResults,'function','cdf','width',9);
' s, ~  |+ i. H, L$ Q; g1 HV1 = ksdensity(PG,PG,'function','cdf');
  |) d. w3 d7 e' `% **********************绘制核密度估计图**********************
[Xsort,id] = sort(compResults);  % 为了作图的需要，对X进行排序
; T  v. {$ e% T2 U) L5 O( Ofigure;  % 新建一个图形窗口
, b2 W. Y8 {# ~hold on
plot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
- o6 A$ q. c$ a9 ?legend('核密度估计', 'Location','NorthWest'); % 加标注框
xlabel('节点8的电压幅值');  % 为X轴加标签
, e, n7 I4 K- ~& S6 b" l9 b' Eylabel('F(x)');  % 为Y轴加标签
* W( {/ j- T, H0 H: C  t但是输出的累积分布图（cdf)却没趋向于1，不知道哪里出了问题，是不是matpower里的程序的问题。