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

sss3434

matpower搭建了个模型，修改case9文件中的 gen 数据，把全局变量windpower风电出力导入到gen(2,pg) ，建立风电接入模型。
function [compResults,PG,results] =  opfWind
0 L% i) t* {7 y8 ~/ l9 C4 z% _global windPower;
wind = xlsread('wind.xls','I2:I8756');
+ G1 y% ?7 i! @3 Z: g6 pcompResults = zeros(8755,1);
PG = zeros(8755,1);
% {% ?' I2 C5 a- P( {6 F' Q3 A% 循环进行潮流计算

# ~+ B8 G. A* N, G9 e3 ffor i = 1:length(wind)
    windPower = wind(i,1)/100;
    loadcase('jycase9');
    results = runpf('jycase9');
    compResults(i,1) = results.branch(7,14);
+ u% f4 N& Y" f, m& O    PG(i,1) = results.branch(7,16);
9 X' Z& J# ~6 o# d& f( T% Dend

, Q4 Z5 y" {! k) N) o%*******************************核密度估计**********************************
% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
  E9 o- J; p0 ]1 ?( k9 f+ E! {+ W9 iU1 = ksdensity(compResults,compResults,'function','cdf','width',9);
' v2 b% ^1 K% xV1 = ksdensity(PG,PG,'function','cdf');
& M3 V/ f* }& `  c4 S4 o/ [% **********************绘制核密度估计图**********************
[Xsort,id] = sort(compResults);  % 为了作图的需要，对X进行排序
+ i% E0 B; T& ~/ J% pfigure;  % 新建一个图形窗口
: r2 j' A( z7 L, {hold on
' o& e2 ]  d/ S. h, Q8 C. {plot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
  G# o. ]  U; Glegend('核密度估计', 'Location','NorthWest'); % 加标注框
xlabel('节点8的电压幅值');  % 为X轴加标签
ylabel('F(x)');  % 为Y轴加标签
; g* T' Z1 r) C但是输出的累积分布图（cdf)却没趋向于1，不知道哪里出了问题，是不是matpower里的程序的问题。

itouka

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