|
|
马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
matpower搭建了个模型,修改case9文件中的 gen 数据,把全局变量windpower风电出力导入到gen(2,pg) ,建立风电接入模型。
U% Y2 u* g' ufunction [compResults,PG,results] = opfWind2 Y5 d2 b( t- k4 q3 w$ [
global windPower;
n' k8 j* j# g |: A5 o* uwind = xlsread('wind.xls','I2:I8756');
: f0 J* U1 ?4 m9 H4 hcompResults = zeros(8755,1);& f8 B) F9 Z& D! r; N1 Z" I
PG = zeros(8755,1);
( Y, U: A; P# d! ?. m% 循环进行潮流计算: d O+ ~2 e! s2 G
: i! Q2 H0 f3 h$ x( x: N+ b
for i = 1:length(wind)
# a, l7 N% O9 l8 O4 R" V windPower = wind(i,1)/100;
4 W: X, A% b5 c2 s$ i1 k loadcase('jycase9');$ H2 A: c& b# s N7 Z' G8 E! ?
results = runpf('jycase9');+ `6 T8 K- q( x" @4 p* L4 ?
compResults(i,1) = results.branch(7,14);
+ W+ y/ A0 P8 C2 K# ]# b8 k PG(i,1) = results.branch(7,16);
% v, q2 ?" l' Q# j# y2 |end
* d. j: q! r. v" l* D3 Y. A, q, ^4 G k. V
%*******************************核密度估计**********************************7 S, H0 m% e6 Z, Q t
% 调用ksdensity函数分别计算原始样本X和Y处的核分布估计值, A d* P& w. s3 l
U1 = ksdensity(compResults,compResults,'function','cdf','width',9);, J9 W# V+ i( G/ V
V1 = ksdensity(PG,PG,'function','cdf');
* G7 B1 t" p1 v+ p% a: t5 ~1 I% **********************绘制核密度估计图**********************
* f+ H) q1 `/ y: O8 X: t[Xsort,id] = sort(compResults); % 为了作图的需要,对X进行排序
( k9 a, Z ]* N& s6 _1 qfigure; % 新建一个图形窗口
, l$ A* K+ o. B, k! L: {1 yhold on) ^; F. f' f0 Y% ^4 M& }& h
plot(Xsort,U1(id),'k-.','LineWidth',2); % 绘制核分布估计图
' T: `8 G% o& w& R" r3 z6 _4 vlegend('核密度估计', 'Location','NorthWest'); % 加标注框
* k5 @0 R% b( y6 R" p. v5 nxlabel('节点8的电压幅值'); % 为X轴加标签
$ k3 R$ F5 U4 } qylabel('F(x)'); % 为Y轴加标签
* Q1 x) |; d& o8 c3 k% P! N但是输出的累积分布图(cdf)却没趋向于1,不知道哪里出了问题,是不是matpower里的程序的问题。 |
|