1 分布式能源的选址与定容
分布式电源（DG）作为发展智能电网的重要组成部分，具有减轻环境污染、降低网络损耗以及提
高电能质量等优点。按 DG 所使用的一次能源可以主要分为两类：可再生能源 DG 和不可再生能源
DG。可再生能源 DG 主要包括小型风力发电及太阳能光伏发电，其特点是出力受自然条件因素制约，不受调度人员控制；不可再生能源 DG 主要包括燃气轮机发电、燃料电池发电及天然气发电等，其特点是具有很好的调节能力，调度人员可对其发出的有功 P 和无功 Q 进行控制。研究表明，DG 并网对配电网有很大影响，包括电压水平、网损、可靠性等方面，且其影响程度与 DG 的安装位置及容量密切相关。合理的安装位置及容量可以有效改善配电网电压质量、减小有功损耗、提高系统负荷率。配置不合理将适得其反，甚至威胁电网的安全稳定运行。因此，为充分发挥DG 对配电网的积极影响，在规划阶段应对接入配电网的 DG 在最佳安装位置与装机容量两个方面进行重点研究。
2 运行结果
部分代码：
function ybus=ybuspg_ds33 line_dat=Ldat33; b_dat=Bdat33; fb=line_dat(:,2);tb=line_dat(:,3); Res=line_dat(:,4); Xe=line_dat(:,5); Bs=line_dat(:,6); % 接地导纳, B/2... a = line_dat(:,7); % 抽头设定值.. z = Res + 1i*Xe; % 阻抗矩阵... y = 1./z; % 节点导纳 b = 1i*Bs; % 使 B 为虚部... nb = max(max(fb),max(tb)); % 节点个数... ybus = zeros(nb,nb); % 初始化节点导纳矩阵.. br=max(line_dat(:,1)); % br= 支路数 bus=max(line_dat(:,3)); % bus= 节点个数 baseMva=max(b_dat(:,7)); % 基准功率 baseKv=12.62; % 基准电压 z=zeros(bus,bus); zpu=zeros(bus,bus); for ii=1:br x=line_dat(ii,2); y=line_dat(ii,3); z(x,y)=complex(line_dat(ii,4),line_dat(ii,5)); % 在 pu 中转换 a 和 y zpu(x,y)=z(x,y)*baseMva/(((baseKv)^2)*1000); end ybus=zeros(bus,bus); for ii=1:br x=line_dat(ii,2); y=line_dat(ii,3); ybus(x,y)= -(1/zpu(x,y))/a(x); ybus(y,x)= -(1/zpu(x,y))/a(x); ybus(x,x)= ((ybus(x,x)-ybus(x,y))/a(x))+(1i*(line_dat(ii,6)/2)); ybus(y,y)= (ybus(y,y)-ybus(y,x)/a(x)); end ybus;
function ybus=ybuspg_ds33 line_dat=Ldat33; b_dat=Bdat33; fb=line_dat(:,2);tb=line_dat(:,3); Res=line_dat(:,4); Xe=line_dat(:,5); Bs=line_dat(:,6); % 接地导纳, B/2... a = line_dat(:,7); % 抽头设定值.. z = Res + 1i*Xe; % 阻抗矩阵... y = 1./z; % 节点导纳 b = 1i*Bs; % 使 B 为虚部... nb = max(max(fb),max(tb)); % 节点个数... ybus = zeros(nb,nb); % 初始化节点导纳矩阵.. br=max(line_dat(:,1)); % br= 支路数 bus=max(line_dat(:,3)); % bus= 节点个数 baseMva=max(b_dat(:,7)); % 基准功率 baseKv=12.62; % 基准电压 z=zeros(bus,bus); zpu=zeros(bus,bus); for ii=1:br x=line_dat(ii,2); y=line_dat(ii,3); z(x,y)=complex(line_dat(ii,4),line_dat(ii,5)); % 在 pu 中转换 a 和 y zpu(x,y)=z(x,y)*baseMva/(((baseKv)^2)*1000); end ybus=zeros(bus,bus); for ii=1:br x=line_dat(ii,2); y=line_dat(ii,3); ybus(x,y)= -(1/zpu(x,y))/a(x); ybus(y,x)= -(1/zpu(x,y))/a(x); ybus(x,x)= ((ybus(x,x)-ybus(x,y))/a(x))+(1i*(line_dat(ii,6)/2)); ybus(y,y)= (ybus(y,y)-ybus(y,x)/a(x)); end ybus;
