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

sss3434

matpower搭建了个模型，修改case9文件中的 gen 数据，把全局变量windpower风电出力导入到gen(2,pg) ，建立风电接入模型。
4 a/ x3 e7 p+ ]function [compResults,PG,results] =  opfWind
1 V/ Y( e- ?7 P! q7 w4 `% \global windPower;
! F1 J" C8 e7 g; Ywind = xlsread('wind.xls','I2:I8756');
compResults = zeros(8755,1);
PG = zeros(8755,1);
5 Q+ t; x$ d. {! N8 q  T% 循环进行潮流计算
* M0 M. h/ Z: c% O4 m/ l
- p4 \5 Z" n9 Y+ i2 |$ Jfor i = 1:length(wind)
& @6 C. R& o% L2 X- K0 @4 y4 b& w    windPower = wind(i,1)/100;
8 ]$ v! W; G0 h2 [, k9 o) q    loadcase('jycase9');
/ _9 r. h0 C. D. K    results = runpf('jycase9');
6 h" s; X9 U0 V    compResults(i,1) = results.branch(7,14);
    PG(i,1) = results.branch(7,16);
$ B! X8 N! f  k# D6 Oend
& {  m1 K' c  A2 W3 J; f3 @
, |' {1 F% t; N+ E* E%*******************************核密度估计**********************************
0 H4 B: t4 m; e  w. q% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
U1 = ksdensity(compResults,compResults,'function','cdf','width',9);
V1 = ksdensity(PG,PG,'function','cdf');
% **********************绘制核密度估计图**********************
[Xsort,id] = sort(compResults);  % 为了作图的需要，对X进行排序
figure;  % 新建一个图形窗口
9 A) ?# a# C3 Zhold on
. Q4 o* h8 G- M6 ]" J4 wplot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
8 [' j! M" w  s& A* S3 ulegend('核密度估计', 'Location','NorthWest'); % 加标注框
xlabel('节点8的电压幅值');  % 为X轴加标签
ylabel('F(x)');  % 为Y轴加标签
- w- `8 F* t" b+ i但是输出的累积分布图（cdf)却没趋向于1，不知道哪里出了问题，是不是matpower里的程序的问题。

itouka

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