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

sss3434

matpower搭建了个模型，修改case9文件中的 gen 数据，把全局变量windpower风电出力导入到gen(2,pg) ，建立风电接入模型。
function [compResults,PG,results] =  opfWind
global windPower;
) T0 I) E& p- o8 Q4 f. a! pwind = xlsread('wind.xls','I2:I8756');
compResults = zeros(8755,1);
9 q; f; @- l. `9 XPG = zeros(8755,1);
% 循环进行潮流计算

  K9 u9 |; }$ k2 C0 Nfor i = 1:length(wind)
" G4 L1 s$ D  y& _# A: ], w    windPower = wind(i,1)/100;
    loadcase('jycase9');
$ J; B4 e- T" x( Z- ?1 k8 x    results = runpf('jycase9');
* l! d  ^- b* x9 B0 G( O- S    compResults(i,1) = results.branch(7,14);
    PG(i,1) = results.branch(7,16);
end
/ F# }3 x7 u& G3 H+ a' A. G6 w
6 O  |; l5 N7 r! G%*******************************核密度估计**********************************
7 p) n) T) d/ d4 h% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
6 b0 ?; ^: @* j: `8 ?2 Q' QU1 = ksdensity(compResults,compResults,'function','cdf','width',9);
/ M( S  a9 K  L8 aV1 = ksdensity(PG,PG,'function','cdf');
5 Q- T" q4 x8 L( q# F  Y$ E% **********************绘制核密度估计图**********************
( ^& n1 ]( q, q[Xsort,id] = sort(compResults);  % 为了作图的需要，对X进行排序
( x9 Q0 I& f9 g& }. ~+ |2 Nfigure;  % 新建一个图形窗口
hold on
plot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
legend('核密度估计', 'Location','NorthWest'); % 加标注框
xlabel('节点8的电压幅值');  % 为X轴加标签
ylabel('F(x)');  % 为Y轴加标签
但是输出的累积分布图（cdf)却没趋向于1，不知道哪里出了问题，是不是matpower里的程序的问题。

itouka

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