matpower 里的一些问题

sss3434

matpower搭建了个模型，修改case9文件中的 gen 数据，把全局变量windpower风电出力导入到gen(2,pg) ，建立风电接入模型。
function [compResults,PG,results] =  opfWind
* M2 I, A* @& k7 F# r7 xglobal windPower;
9 N: ~4 ~+ f0 l7 q6 I4 d4 zwind = xlsread('wind.xls','I2:I8756');
! B& e4 o, E: W6 e: D( O6 J5 McompResults = zeros(8755,1);
PG = zeros(8755,1);
% 循环进行潮流计算
& {* R: N9 S$ c* y
5 y8 q4 d. z' y4 @# s" H7 Yfor i = 1:length(wind)
    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

" W, `& I5 |! J0 g8 W) a" [/ E1 o%*******************************核密度估计**********************************
% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
U1 = ksdensity(compResults,compResults,'function','cdf','width',9);
: _7 R. k, j8 I% i$ PV1 = ksdensity(PG,PG,'function','cdf');
( \/ s& }& I8 `6 i6 N; t% **********************绘制核密度估计图**********************
[Xsort,id] = sort(compResults);  % 为了作图的需要，对X进行排序
figure;  % 新建一个图形窗口
hold on
0 ?6 s  J' G7 k8 P" F  l. pplot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
! \2 t2 y! l3 e2 ]- T4 g4 [legend('核密度估计', 'Location','NorthWest'); % 加标注框
xlabel('节点8的电压幅值');  % 为X轴加标签
ylabel('F(x)');  % 为Y轴加标签
/ O/ u& B2 q& q5 \6 u* H9 {但是输出的累积分布图（cdf)却没趋向于1，不知道哪里出了问题，是不是matpower里的程序的问题。