matpower 里的一些问题

sss3434

matpower搭建了个模型，修改case9文件中的 gen 数据，把全局变量windpower风电出力导入到gen(2,pg) ，建立风电接入模型。
9 ]1 p6 t9 ?3 yfunction [compResults,PG,results] =  opfWind
2 E3 x% n! l9 xglobal windPower;
; i) `9 n5 {! G+ |8 Iwind = xlsread('wind.xls','I2:I8756');
: k/ ]) u( s& v8 S# o$ ccompResults = zeros(8755,1);
3 ^4 _' G3 `1 ?& cPG = zeros(8755,1);
  F! }; \: p6 g% 循环进行潮流计算
/ c$ a: Y/ |+ Y+ j# c5 z8 E
3 h2 i& F9 _2 P3 ^# r8 F  E' g" \9 |/ X' bfor i = 1:length(wind)
    windPower = wind(i,1)/100;
    loadcase('jycase9');
; U+ l0 T! X+ @, U/ I5 o    results = runpf('jycase9');
    compResults(i,1) = results.branch(7,14);
7 U* ?) _) o: A) P8 N    PG(i,1) = results.branch(7,16);
- u4 n4 I9 {1 M- w  h% \4 mend
. Z7 o& @" K% b  `9 J7 Y
; [: |$ Y. i8 ]& ]%*******************************核密度估计**********************************
  i3 l+ a# P2 A  M- ~% B& l' C, Q% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
, E; P. k9 @' G# ^# _9 i2 RU1 = ksdensity(compResults,compResults,'function','cdf','width',9);
6 @9 q9 B: e7 xV1 = ksdensity(PG,PG,'function','cdf');
$ N7 ?) d9 H- g) `7 ~% **********************绘制核密度估计图**********************
( U! k& o2 |7 V5 x[Xsort,id] = sort(compResults);  % 为了作图的需要，对X进行排序
figure;  % 新建一个图形窗口
hold on
plot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
legend('核密度估计', 'Location','NorthWest'); % 加标注框
xlabel('节点8的电压幅值');  % 为X轴加标签
* V3 x$ W# k8 ^+ a9 }; {ylabel('F(x)');  % 为Y轴加标签
但是输出的累积分布图（cdf)却没趋向于1，不知道哪里出了问题，是不是matpower里的程序的问题。