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

sss3434

matpower搭建了个模型，修改case9文件中的 gen 数据，把全局变量windpower风电出力导入到gen(2,pg) ，建立风电接入模型。
3 z% o4 L5 @% ^* ?2 ]4 ufunction [compResults,PG,results] =  opfWind
6 W- m/ X3 `1 jglobal windPower;
7 z) v) M8 Z! d# m6 k$ {" h* |wind = xlsread('wind.xls','I2:I8756');
" y0 |% C9 ?0 I& K7 g- Z0 A. FcompResults = zeros(8755,1);
PG = zeros(8755,1);
5 O% `4 R. E8 C0 N" S# f1 x) @% 循环进行潮流计算

# b' C7 E: M0 u, G5 ]for i = 1:length(wind)
3 M2 j. w$ ^; z4 a6 f( @2 f+ l, g    windPower = wind(i,1)/100;
    loadcase('jycase9');
    results = runpf('jycase9');
! l+ {5 r8 M" \. s    compResults(i,1) = results.branch(7,14);
/ y5 \$ I5 _' g# e) H    PG(i,1) = results.branch(7,16);
end

%*******************************核密度估计**********************************
% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
# J1 C" `6 l1 t, H% u* WU1 = ksdensity(compResults,compResults,'function','cdf','width',9);
& }; S4 `& h. J- u3 S, q, A5 SV1 = ksdensity(PG,PG,'function','cdf');
% **********************绘制核密度估计图**********************
[Xsort,id] = sort(compResults);  % 为了作图的需要，对X进行排序
figure;  % 新建一个图形窗口
hold on
plot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
legend('核密度估计', 'Location','NorthWest'); % 加标注框
8 s, E& Y) ~, }" T9 q4 sxlabel('节点8的电压幅值');  % 为X轴加标签
3 W: P6 U7 A% Jylabel('F(x)');  % 为Y轴加标签
5 r7 M7 ~; Y* u* r6 i& T但是输出的累积分布图（cdf)却没趋向于1，不知道哪里出了问题，是不是matpower里的程序的问题。

itouka

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