matpower 里的一些问题

sss3434

matpower搭建了个模型，修改case9文件中的 gen 数据，把全局变量windpower风电出力导入到gen(2,pg) ，建立风电接入模型。
function [compResults,PG,results] =  opfWind
1 T( |8 T& |% u9 S4 t6 s8 |global windPower;
! g$ ]0 P2 e, Qwind = xlsread('wind.xls','I2:I8756');
  \+ ?  ~  t$ J5 M; ncompResults = zeros(8755,1);
PG = zeros(8755,1);
% T2 A4 q* L8 _7 A* y0 s8 E% 循环进行潮流计算

for i = 1:length(wind)
    windPower = wind(i,1)/100;
    loadcase('jycase9');
    results = runpf('jycase9');
; B/ m: C0 }* j, D2 q- z" ?    compResults(i,1) = results.branch(7,14);
* c) N% e4 j" D( m# G    PG(i,1) = results.branch(7,16);
- u  H; B! Z+ c- D- hend

1 m! H( V' ?; ~; |2 @8 i%*******************************核密度估计**********************************
% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
U1 = ksdensity(compResults,compResults,'function','cdf','width',9);
V1 = ksdensity(PG,PG,'function','cdf');
% **********************绘制核密度估计图**********************
3 F& u7 Q1 _1 ~, q6 F% b3 g: g1 p) ^[Xsort,id] = sort(compResults);  % 为了作图的需要，对X进行排序
figure;  % 新建一个图形窗口
hold on
/ S4 b2 J/ n  f( [plot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
3 u6 l; H, F6 x, C) Xlegend('核密度估计', 'Location','NorthWest'); % 加标注框
xlabel('节点8的电压幅值');  % 为X轴加标签
ylabel('F(x)');  % 为Y轴加标签
$ H, Y3 l& x1 @6 ?% }( u! ~, O" E但是输出的累积分布图（cdf)却没趋向于1，不知道哪里出了问题，是不是matpower里的程序的问题。