ftrbid_data.m

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

% FTR Market Bid data (FTR 拍卖基本数据数据)
function [b0,P_base,T2,T3]=FTRbid_data

%%用标幺值表示%%
%      (按ieee-30节点系统顺序排列)    %联络线路传输容量限制          
% b0=[1.3,1.3,0.65,1.3,1.3,0.65,0.9,0.7,1.3,0.32,0.65,0.32,0.65,0.65,0.65,0.65,0.32,0.32,0.32,0.16,0.16,0.16,0.16,0.32,0.32,0.32,0.32,0.32,0.32,0.16,0.16,0.16,0.16,0.16,0.16,0.65,0.16,0.16,0.16,0.32,0.32];                                      
b0=[130,130,65,130,130,65,90,70,130,32,65,32,65,65,65,65,32,32,32,16,16,16,16,32,32,32,32,32,32,16,16,16,16,16,16,65,16,16,16,32,32];

%%计算初始FTR分配矩阵%%
T1=zeros(30,6);           %FTR优先分配矩阵
T1(1,1)=1;   T1(3,1)=-0.95;    T1(20,1)=-0.05;
T1(2,2)=1;   T1(14,2)=-0.375;  T1(15,2)=-0.375;  T1(18,2)=-0.25;
T1(13,3)=1;  T1(21,3)=-1;
T1(22,4)=1;  T1(12,4)=-0.68;   T1(16,4)=-0.12;   T1(19,4)=-0.2;
T1(23,5)=1;  T1(7,5)=-0.6;     T1(10,5)=-0.4;    
T1(27,6)=1;  T1(8,6)=-0.25;    T1(10,6)=-0.75;
% FTR_base=[1,0.8,0.2,0.5,0.5,0.2];   %在该时段仍有效的FTR数量
FTR_base=[100,80,20,50,50,20];
P_base=T1*FTR_base';
  
T2=zeros(30,7);                     %FTR拍卖矩阵(最后两列为offer，对原有潮流是release)
T2(2,1)=1;   T2(25,1)=-0.2;   T2(27,1)=-0.4;   T2(30,1)=-0.4;
T2(8,2)=1;   T2(29,2)=-1;
T2(1,3)=0.5; T2(2,3)=0.5;     T2(5,3)=-0.5;    T2(23,3)=-0.5;
T2(3,4)=1;   T2(10,4)=-1; 
T2(14,5)=1;  T2(13,5)=-0.5;   T2(20,5)=-0.5;
T2(1,6)=-1;  T2(15,6)=1;          %offer数据
T2(7,7)=-1;  T2(22,7)=1;


%风力发电的节点分布数据(1节点)
T3=zeros(30,1);
T3(1,1)=1;       T3(2,1)=-0.3;    T3(3,1)=-0.3;     T3(4,1)=-0.15;  
T3(7,1)=-0.0625; T3(8,1)=-0.0625; T3(10,1)=-0.0625; T3(12,1)=-0.0625;
                              %风力发电的数据(2节点)
% T3=zeros(30,1);
% T3(2,1)=1;       T3(4,1)=-0.3;    T3(3,1)=-0.2;     T3(7,1)=-0.2;  
% T3(8,1)=-0.1;    T3(12,1)=-0.1;   T3(10,1)=-0.1; 

                               %风力发电的数据(11节点)
% T3=zeros(30,1);
% T3(11,1)=1;      T3(10,1)=-0.3;   T3(17,1)=-0.1;     T3(20,1)=-0.1;  
% T3(21,1)=-0.1;    T3(8,1)=-0.1;   T3(7,1)=-0.1;      T3(2,1)=-0.1;   T3(4,1)=-0.1; 

                              %风力发电的数据(23节点)
% T3=zeros(30,1);
% T3(23,1)=1;       T3(24,1)=-0.3;    T3(15,1)=-0.3;     T3(25,1)=-0.1;  
% T3(14,1)=-0.1;    T3(18,1)=-0.1;    T3(12,1)=-0.1; 

                              %风力发电的数据(22节点)
% T3=zeros(30,1);
% T3(22,1)=1;       T3(21,1)=-0.2;    T3(10,1)=-0.2;     T3(24,1)=-0.2;  
% T3(23,1)=-0.1; T3(25,1)=-0.1; T3(17,1)=-0.1;  T3(20,1)=-0.1; 

                              %风力发电的数据(27节点)
% T3=zeros(30,1);
% T3(27,1)=1;       T3(29,1)=-0.3;    T3(30,1)=-0.3;     T3(25,1)=-0.2;  
% T3(26,1)=-0.1;    T3(24,1)=-0.1;    

                               %风力发电的数据(6节点)
% T3=zeros(30,1);
% T3(6,1)=1;        T3(4,1)=-0.15;    T3(2,1)=-0.15;     T3(7,1)=-0.15;  T3(8,1)=-0.15;
% T3(10,1)=-0.15;   T3(3,1)=-0.0625;  T3(17,1)=-0.0625;  T3(20,1)=-0.0625;  T3(21,1)=-0.0625;

                               %风力发电的数据(13节点)
% T3=zeros(30,1);
% T3(13,1)=1;        T3(12,1)=-0.257;    T3(14,1)=-0.1225;     T3(15,1)=-0.1645;  T3(16,1)=-0.075;
% T3(17,1)=-0.105;   T3(18,1)=-0.065;    T3(19,1)=-0.07;       T3(23,1)=-0.0675;  T3(24,1)=-0.0735;

%风力发电的数据(27节点)
%T3=zeros(30,1);
%T3(27,1)=1;        T3(7,1)=-0.0122;    T3(8,1)=-0.2788;       T3(10,1)=-0.0083;  T3(17,1)=-0.0117;     T3(19,1)=-0.0083;                    
%T3(20,1)=-0.0029;  T3(21,1)=-0.0253;   T3(24,1)=-0.1321;      T3(26,1)=-0.1837;  T3(29,1)=-0.1223;     T3(30,1)=-0.2144; 

return