matpower搭建三机九节点模型的问题
matpower搭建了个模型,修改case9文件中的 gen 数据,把全局变量windpower风电出力导入到gen(2,pg) ,建立风电接入模型。function =opfWind
global windPower;
wind = xlsread('wind.xls','I2:I8756');
compResults = zeros(8755,1);
PG = zeros(8755,1);
% 循环进行潮流计算
for i = 1:length(wind)
windPower = wind(i,1)/100;
loadcase('jycase9');
results = runpf('jycase9');
compResults(i,1) = results.branch(7,14);
PG(i,1) = results.branch(7,16);
end
%*******************************核密度估计**********************************
% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值
U1 = ksdensity(compResults,compResults,'function','cdf','width',9);
V1 = ksdensity(PG,PG,'function','cdf');
% **********************绘制核密度估计图**********************
= sort(compResults);% 为了作图的需要,对X进行排序
figure;% 新建一个图形窗口
hold on
plot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
legend('核密度估计', 'Location','NorthWest'); % 加标注框
xlabel('节点8的电压幅值');% 为X轴加标签
ylabel('F(x)');% 为Y轴加标签
但是输出的累积分布图(cdf)却没趋向于1,不知道哪里出了问题,是不是matpower里的程序的问题。 我不太懂非参数估计的问题,我只是读到过一篇文章,里面说抽样次数少了就不能达到概率应达到的点,这个能不能推广到核函数呢?是不是你的概率分布处理出了问题?
页:
[1]