windpower_weibull.m

考虑风电的FTR拍卖模型研究论文基于蒙特卡罗模拟的遗传算法源程序

%产生风机随机功率的函数
function Pw=windpower_weibull(n)

V=14*weibrnd(2,8.5,n,1);             %生成n个风速随机样本
Pr=36;                               %风机额定容量(600KW*60)可修改
Vci=3;                               %切入风速
Vco=25;                              %切出风速
Vr=14;                               %额定风速
Pw1=zeros(n,1);
for i=1:n                            %分情况产生风机随机功率
    if V(i)<Vci|V(i)>=Vco
        Pw1(i)=0;
    end
    if V(i)>=Vci & V(i)<Vr
        Pw1(i)=(Pr/(Vr^3-Vci^3))*(V(i)^3-Vci^3);
    end
    if V(i)>=Vr & V(i)<Vco
        Pw1(i)=Pr;
    end
end

%%用标幺值计算%%
Pw= Pw1;
return