PV matlab编程

matlab编程，牛拉法计算潮流，选取了34节点的配电网网络模型，通过对单个以及多个分布式电源的接入位置以及容量的不同情况对34节点配电网的网损以及节点电压状况进行了分析。

function [ref pv pq] = bustypes（bus gen）
%创建参考节点，pq节点和pv节点
[PQ PV REF NONE BUS_I BUS_TYPE PD QD GS BS BUS_AREA VM ...
    VA base_KV ZONE VMAX VMIN LAM_P LAM_Q MU_VMAX MU_VMIN] = idx_bus;
[GEN_BUS PG QG QMAX QMIN VG Mbase GEN_STATUS PMAX PMIN ...
    MU_PMAX MU_PMIN MU_QMAX MU_QMIN PC1 PC2 QC1MIN QC1MAX ...
    QC2MIN QC2MAX RAMP_AGC RAMP_10 RAMP_30 RAMP_Q APF] = idx_gen;
nb = size（bus 1）;
ng = size（gen 1）;
Cg = sparse（gen（: GEN_BUS） （1:ng）‘ gen（: GEN_STATUS） > 0 nb ng）;  
                                        
bus_gen_status = Cg * ones（ng 1）;     

ref = find（bus（: BUS_TYPE） == REF & bus_gen_status）;   %% 参考节点索引
pv  = find（bus（: BUS_TYPE） == PV  & bus_gen_status）;   %% PV 节点索引
pq  = find（bus（: BUS_TYPE） == PQ | ~bus_gen_status）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2278  2016-05-22 17:15  idx_gen.m

     文件        319  2016-05-22 17:15  isload.m

     文件       7698  2016-05-22 17:15  loadcase.m

     文件        988  2016-05-22 17:15  makeSbus.m

     文件       2094  2016-05-22 17:15  makeYbus.m

     文件      10884  2016-05-22 17:15  mpoption.m

     文件       2581  2016-05-22 17:15  newtonpf.m

     文件       3222  2016-05-22 17:15  pfsoln.m

     文件      37711  2016-05-22 17:15  printpf.m

     文件       2504  2016-05-22 17:15  runpf.asv

     文件       2592  2016-05-22 17:15  runpf.m

     文件       1004  2016-05-22 17:15  bustypes.m

     文件       1366  2016-05-22 17:15  case9.m

     文件       5385  2016-05-22 17:15  case34.m

     文件       1492  2016-05-22 17:15  dSbus_dV.m

     文件       1728  2016-05-22 17:15  idx_brch.m

     文件       1519  2016-05-22 17:15  idx_bus.m

----------- ---------  ---------- -----  ----

                85365                    17