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

sss3434

matpower搭建了个模型，修改case9文件中的 gen 数据，把全局变量windpower风电出力导入到gen(2,pg) ，建立风电接入模型。
function [compResults,PG,results] =  opfWind
! y0 h, k+ h) J9 h% l: p, X4 [global windPower;
wind = xlsread('wind.xls','I2:I8756');
* o5 K  m, G0 K, B$ F. l0 QcompResults = zeros(8755,1);
PG = zeros(8755,1);
6 ^3 k# @2 J1 N# v% T% 循环进行潮流计算

) a2 S5 p8 A6 i/ c% ?  Q1 zfor i = 1:length(wind)
/ n, G1 P. f0 f( Y6 B    windPower = wind(i,1)/100;
; J& ~  k8 B2 X4 W9 r    loadcase('jycase9');
$ Z" g6 P& }: l( g( ?* T8 B( x    results = runpf('jycase9');
    compResults(i,1) = results.branch(7,14);
- A0 i# ]/ j% a( j, j, O) t    PG(i,1) = results.branch(7,16);
# `+ z4 o- D+ ~- u- Q% {end
. [3 D0 T% V, z1 W2 T4 v
%*******************************核密度估计**********************************
% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
7 k0 V, s" I- ?& N" c% Z* dU1 = ksdensity(compResults,compResults,'function','cdf','width',9);
$ Z, t& e) b5 O3 q2 yV1 = ksdensity(PG,PG,'function','cdf');
' S3 ?5 ?. O7 C+ D5 x9 x  v$ a% **********************绘制核密度估计图**********************
[Xsort,id] = sort(compResults);  % 为了作图的需要，对X进行排序
  f& U( z% Z5 f3 t" e5 dfigure;  % 新建一个图形窗口
( O+ O9 T4 {: \7 ~* P) yhold on
( F0 ~; y. y; y# F8 e  fplot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
/ _: h' O& J! c5 B& Olegend('核密度估计', 'Location','NorthWest'); % 加标注框
# k" [# q6 {) j* _2 C7 v/ dxlabel('节点8的电压幅值');  % 为X轴加标签
ylabel('F(x)');  % 为Y轴加标签
但是输出的累积分布图（cdf)却没趋向于1，不知道哪里出了问题，是不是matpower里的程序的问题。

itouka

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