matpower搭建三机九节点模型的问题

sss3434

matpower搭建了个模型，修改case9文件中的 gen 数据，把全局变量windpower风电出力导入到gen(2,pg) ，建立风电接入模型。
+ |+ Z) C) y& s2 ^  ^2 M2 jfunction [compResults,PG,results] =  opfWind
: O* J" X' @# }5 a2 i5 cglobal windPower;
0 E6 T( {3 Z4 {* _wind = xlsread('wind.xls','I2:I8756');
compResults = zeros(8755,1);
PG = zeros(8755,1);
- C$ q  m+ K1 }, ~& m% 循环进行潮流计算
. L5 {2 _' ^$ ^% H
9 Q4 u0 e2 t5 |! pfor i = 1:length(wind)
2 b. Y( G7 {- V    windPower = wind(i,1)/100;
    loadcase('jycase9');
    results = runpf('jycase9');
    compResults(i,1) = results.branch(7,14);
; H' v. H" P5 y1 N' j    PG(i,1) = results.branch(7,16);
+ {$ O, F2 K7 f+ L( f8 K/ S- Vend
2 L" E! P# s" Y& d
%*******************************核密度估计**********************************
% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
) B2 W( M& y% B. CU1 = ksdensity(compResults,compResults,'function','cdf','width',9);
V1 = ksdensity(PG,PG,'function','cdf');
- u! [1 A9 q2 l* g9 F$ C1 l+ v% **********************绘制核密度估计图**********************
4 B" T$ d9 p; K9 ?. {2 g[Xsort,id] = sort(compResults);  % 为了作图的需要，对X进行排序
* ?+ r4 Z4 @( z1 [8 Xfigure;  % 新建一个图形窗口
hold on
plot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
legend('核密度估计', 'Location','NorthWest'); % 加标注框
xlabel('节点8的电压幅值');  % 为X轴加标签
* M$ ^9 v+ Z9 ~ylabel('F(x)');  % 为Y轴加标签
+ i. r2 q# N2 E3 S' H7 k# v但是输出的累积分布图（cdf)却没趋向于1，不知道哪里出了问题，是不是matpower里的程序的问题。

itouka

我不太懂非参数估计的问题，我只是读到过一篇文章，里面说抽样次数少了就不能达到概率应达到的点，这个能不能推广到核函数呢？是不是你的概率分布处理出了问题？