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

sss3434

matpower搭建了个模型，修改case9文件中的 gen 数据，把全局变量windpower风电出力导入到gen(2,pg) ，建立风电接入模型。
. w8 y8 z2 T# A+ i/ ofunction [compResults,PG,results] =  opfWind
7 e/ F8 H6 E4 fglobal windPower;
+ w% B7 u. V8 ^9 ^/ g! awind = xlsread('wind.xls','I2:I8756');
) ]. {/ I+ D! ?8 j4 o4 |9 S2 ^( \! mcompResults = zeros(8755,1);
PG = zeros(8755,1);
% 循环进行潮流计算
5 G8 K! T( x& l" ]4 w" Z, K
for i = 1:length(wind)
" U( ?4 }; w% _; {  Z3 ^    windPower = wind(i,1)/100;
: ]( {3 O# {% M" g" A! `" Z    loadcase('jycase9');
    results = runpf('jycase9');
: D& c: _" x' c1 S    compResults(i,1) = results.branch(7,14);
8 G1 F6 |' K5 a" k( S    PG(i,1) = results.branch(7,16);
# b* `3 R: A" H, j3 xend
4 X" {9 C  ~( n$ H
%*******************************核密度估计**********************************
% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
. T5 `# f$ k* o0 E, B: W2 ]U1 = ksdensity(compResults,compResults,'function','cdf','width',9);
V1 = ksdensity(PG,PG,'function','cdf');
% **********************绘制核密度估计图**********************
[Xsort,id] = sort(compResults);  % 为了作图的需要，对X进行排序
figure;  % 新建一个图形窗口
: q. r/ u1 B& r! z$ v7 ihold on
plot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
legend('核密度估计', 'Location','NorthWest'); % 加标注框
xlabel('节点8的电压幅值');  % 为X轴加标签
# x8 b5 M# I% a) f& x: E) I+ S) iylabel('F(x)');  % 为Y轴加标签
8 s: L( ?7 g6 P5 w* o9 c5 q6 L9 {但是输出的累积分布图（cdf)却没趋向于1，不知道哪里出了问题，是不是matpower里的程序的问题。

itouka

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