psse2psat

一个较好的MATLAB潮流程序

		$ltc_r[$nltc] = $ltc_r[$nltc]/$h;
	    }    	    
	}	
    } else {
	# three-winding transformer
	$ntw++;
	$bustw1[$ntw] = int($data1[0]);
	$bustw2[$ntw] = int($data1[1]);
	$bustw3[$ntw] = int($data1[2]);
	$x12[$ntw] = $data2[1]*$pbas/$data2[2];
	$x23[$ntw] = $data2[4]*$pbas/$data2[5];
	$x13[$ntw] = $data2[7]*$pbas/$data2[8];
	if ($data1[5] == 3) {
	    $j = $kvbase{$bustw1[$ntw]};
	    $h = $data2[2]*$data2[2]/3/$j/$j;
	    $r12[$ntw] = $data2[0]/$h;
	    $j = $kvbase{$bustw2[$ntw]};
	    $h = $data2[5]*$data2[5]/3/$j/$j;
	    $r23[$ntw] = $data2[3]/$h;
	    $j = $kvbase{$bustw3[$ntw]};
	    $h = $data2[8]*$data2[8]/3/$j/$j;
	    $r13[$ntw] = $data2[6]/$h;
	} else {
	    $r12[$ntw] = $data2[0]*$pbas/$data2[2];
	    $r23[$ntw] = $data2[3]*$pbas/$data2[5];
	    $r13[$ntw] = $data2[6]*$pbas/$data2[8];
	}
	$rt1a[$ntw] = $data3[3]/$pbas;
	$rt1b[$ntw] = $data3[4]/$pbas;
	$rt1c[$ntw] = $data3[5]/$pbas;
	$rt2a[$ntw] = $data4[3]/$pbas;
	$rt2b[$ntw] = $data4[4]/$pbas;
	$rt2c[$ntw] = $data4[5]/$pbas;
	$rt3a[$ntw] = $data5[3]/$pbas;
	$rt3b[$ntw] = $data5[4]/$pbas;
	$rt3c[$ntw] = $data5[5]/$pbas;
    }
}

# area data
while (<IN>) {
    @data = mysplit($_);
    if ($data[0] == 0) { last; }        
}

# two terminal DC data
while (<IN>) {
    @data = mysplit($_);
    if ($data[0] == 0) { last; }        
}

# VSC DC line data
while (<IN>) {
    @data = mysplit($_);
    if ($data[0] == 0) { last; }        
}

# switched shunt data
my ($nsvc,@svcbus,@svcbma,@svcbmi,@svcref);
$nsvc = -1;
while (<IN>) {
    @data = mysplit($_);
    if ($data[0] == 0) { last; }        
    #if (isidx(int($data[0]))) {next;}
    
    #$j = pvidx(int($data[0]));
    
    # add SVC only if there is no PV at the same bus
    #if ($j) {
    #printf "SVC at PV bus " . $data[0] . " will be ignored.\n" ;
    #} elsif ($data[1] == 0) { # fixed SVCs are treated as a shunt
    
    $nsh++;
    $idxsh[$nsh] = int($data[0]);
    $psh[$nsh] = 0;	
    $qsh[$nsh] = $data[6]/$pbas;		
    
    #} else {
    #$nsvc++;
    #$svcbus[$nsvc] = int($data[0]);
    #if ($data[1] < 3) {
    #    $svcref[$nsvc] = 0.5*($data[2]+$data[3]);
    #} else {
    #    $svcref[$nsvc] = 1;	
    #}
    #$svcbma[$nsvc] = 0;
    #for ($i = @data-1; $i > 6; $i-=2) {
    #    $svcbma[$nsvc] += $data[$i]*$data[$i-1]; 
    #}
    #$svcbma[$nsvc] = abs($svcbma[$nsvc])/$pbas; 
    #$svcbmi[$nsvc] = -$svcbma[$nsvc]; 
    #if ($svcbma[$nsvc] == 0) {$svcbma[$nsvc] = 1;}
    #if ($svcbmi[$nsvc] == 0) {$svcbmi[$nsvc] = -1;}
    
    ## add PV generator to initialize SVC component
    #$npv++;
    #$pvsb[$npv] = $pbas;
    #$pvbus[$npv] = int($data[0]);
    #$pvpg[$npv] = 0;
    #if ($data[2] && $data[1] > 2) {
    #    $pvqt[$npv] = $data[2];
    #} else {
    #    $pvqt[$npv] = 99;
    #}
    #if ($data[3] && $data[1] > 2) {
    #    $pvqt[$npv] = $data[3];
    #} else {
    #    $pvqt[$npv] = -99;
    #}
    #$pvqb[$npv] = -99;
    #$pvvs[$npv] = $svcref[$nsvc];
    #$pvpt[$npv] = 0;
    #$pvpb[$npv] = 0;
    #}
}


# -----------------------------------------------------------------------
# close input data file
# -----------------------------------------------------------------------
close(IN) || die "cannot close $ARGV[0]: $!\n";

# -----------------------------------------------------------------------
# open output data file
# -----------------------------------------------------------------------
print "Writing PSAT file \"$ARGV[1]\"...\n";
carret();
open(OUT,">$ARGV[1]") || die "cannot open $ARGV[1]: $!\n";

# -----------------------------------------------------------------------
# write output data file
# -----------------------------------------------------------------------
print OUT "% File generated by PSAT from PSS/E data file.\n";
print OUT "% "."-" x 78 . "\n";
print OUT "% Author:   Federico Milano\n";
print OUT "% E-mail:   fmilano\@thunderbox.uwaterloo.ca\n";
print OUT "% Web-site: http://thunderbox.uwaterloo.ca/~fmilano\n";
print OUT "% "."-" x 78 . "\n";
print OUT "% $title1\n";
print OUT "% $title2\n\n";

# -----------------------------------------------------------------------
# write Bus.con
# -----------------------------------------------------------------------
$format = "%4d %8.2f %8.5f %8.5f %2d %2d";
$nbus >= 0 && print OUT "Bus.con = [ ...\n";
for ($i = 0; $i <= $nbus-1; $i++) {
    printf OUT "$format;\n",$busidx[$i],$kvb[$i],$ein[$i],$ang[$i],
    $kae[$i],$kzo[$i];
}
$i = $nbus;
$nbus >= 0 && printf OUT "$format];\n\n", $busidx[$i], $kvb[$i], 
    $ein[$i],$ang[$i],$kae[$i],$kzo[$i];

# -----------------------------------------------------------------------
# write SW.con
# -----------------------------------------------------------------------
if ($nsw >= 0) {
    $format = "%4d %7.2f %6.2f %8.5f 0.00 %8.5f %8.5f 1.1 0.9 %8.5f 1";
    printf OUT "SW.con = [$format];\n\n", $swbus,$swsb,$kvbase{$swbus},
    $swvs,$swqt,$swqb,$swpg;
}

# -----------------------------------------------------------------------
# write PV.con
# -----------------------------------------------------------------------
$npv >= 0 && printf OUT "PV.con = [ ...\n";
$format = "%4d %7.2f %6.2f " . "%8.5f " x 4 . "1.1 0.9 1";
for ($i = 0; $i <= $npv-1; $i++) {
    #$h = $idxpv[$i];
    printf OUT "$format;\n", $pvbus[$i],$pvsb[$i],$kvbase{$pvbus[$i]},
    $pvpg[$i],$pvvs[$i],$pvqt[$i],$pvqb[$i];
}
#$h = $idxpv[$npv];
$i = $npv;
$npv >= 0 && printf OUT "$format];\n\n",$pvbus[$i],$pvsb[$i],
    $kvbase{$pvbus[$i]},$pvpg[$i],$pvvs[$i],$pvqt[$i],$pvqb[$i];

# -----------------------------------------------------------------------
# write PQ.con
# -----------------------------------------------------------------------
$npq >= 0 && printf OUT "PQ.con = [ ...\n";
$format = "%4d $pbas %8.2f %8.5f %8.5f 1.1 0.9 1";
for ($i = 0; $i <= $npq-1; $i++) {
    $h = $pqbus[$i];
    printf OUT "$format;\n", $h,$kvbase{$h},$pqpl[$i],$pqql[$i];
}
$h = $pqbus[$npq];
$i = $npq;
$npq >= 0 && printf OUT "$format];\n\n",$h,$kvbase{$h},$pqpl[$i],
    $pqql[$i];

# -----------------------------------------------------------------------
# write Pl.con
# -----------------------------------------------------------------------
$npl >= 0 && print OUT "Pl.con = [ ...\n";
$format = "%4d $pbas %8.2f 60 %8.5f %8.5f 0 %8.5f %8.5f 0 0";
for ($i = 0; $i <= $npl-1; $i++) {
    printf OUT "$format;\n", $plbus[$i],$kvbase{$plbus[$i]},$plyp[$i],
    $plip[$i],$plyq[$i],$pliq[$i];
}
$h = $plbus[$npl];
$i = $npl;
$npl >= 0 && printf OUT "$format];\n\n",$plbus[$i],$kvbase{$plbus[$i]},
    $plyp[$i],$plip[$i],$plyq[$i],$pliq[$i];

# -----------------------------------------------------------------------
# write Shunt.con
# -----------------------------------------------------------------------
$nsh >= 0 && print OUT "Shunt.con = [ ...\n";
$format = "%4d $pbas %8.2f 60 %8.5f %8.5f";
for ($i = 0; $i <= $nsh-1; $i++) {
    $h = abs($idxsh[$i]);
    printf OUT "$format;\n", $h,$kvbase{$h},$psh[$i],$qsh[$i];
}
$h = $idxsh[$nsh];
$i = $nsh;
$nsh >= 0 && printf OUT "$format];\n\n",$h,$kvbase{$h},$psh[$i],$qsh[$i];

# -----------------------------------------------------------------------
# write Line.con
# -----------------------------------------------------------------------
$nline >= 0 && print OUT "Line.con = [ ...\n";
#$format = "%4d %4d %8.2f %8.2f 60 0 %5.2f " . "%8.5f " x 3 . 
#    " 0 %8.5f %8.3f %8.3f %8.3f";
$format = "%4d %4d %8.2f %8.2f 60 0 %5.2f " . "%8.5f " x 5 . 
    " %8.3f %8.3f %8.3f";
for ($i = 0; $i <= $nline-1; $i++) {
    $k = $kvbase{$busfr[$i]}/$kvbase{$bustt[$i]};
    if ($k == 1) {$k = 0;}
    printf OUT "$format;\n", $busfr[$i],$bustt[$i],$lbas[$i],
    $kvbase{$busfr[$i]},$k,$rest[$i],$reat[$i],$susc[$i],$ratio[$i],$phas[$i],
    $ratea[$i],$rateb[$i],$ratec[$i];
}
$i = $nline;
if ($nline >= 0) {
    $k = $kvbase{$busfr[-1]}/$kvbase{$bustt[-1]};
    if ($k == 1) {$k = 0;}
    printf OUT "$format];\n\n",$busfr[$i],$bustt[$i],$lbas[$i],
    $kvbase{$busfr[$i]},$k,$rest[$i],$reat[$i],$susc[$i],$ratio[$i],$phas[$i],
    $ratea[$i],$rateb[$i],$ratec[$i];
}

# -----------------------------------------------------------------------
# write Ltc.con
# -----------------------------------------------------------------------
$nltc >= 0 && print OUT "Ltc.con = [ ...\n";
$format = "%4d %4d %8.2f %8.2f 60 1 0.001 0.1 " . 
    "%8.5f " x 6 . " %4d %4d";
for ($i = 0; $i <= $nltc-1; $i++) {
    printf OUT "$format;\n", $ltc_busfr[$i],$ltc_busto[$i],$ltc_bas[$i],
    $kvbase{$ltc_busfr[$i]},$ltc_vma[$i],$ltc_vmi[$i],$ltc_dm[$i],
    $ltc_ref[$i],$ltc_x[$i],$ltc_r[$i],$ltc_busct[$i],$ltc_ctr[$i];
}
$i = $nltc;
if ($nltc >= 0) {
    printf OUT "$format];\n\n", $ltc_busfr[$i],$ltc_busto[$i],$ltc_bas[$i],
    $kvbase{$ltc_busfr[$i]},$ltc_vma[$i],$ltc_vmi[$i],$ltc_dm[$i],
    $ltc_ref[$i],$ltc_x[$i],$ltc_r[$i],$ltc_busct[$i],$ltc_ctr[$i];
}

# -----------------------------------------------------------------------
# write Phs.con
# -----------------------------------------------------------------------
$nphs >= 0 && print OUT "Phs.con = [ ...\n";
$format = "%4d %4d %8.2f %8.2f %8.2f 60 0.001 10 50 " . "%8.5f " x 5;
for ($i = 0; $i <= $nphs-1; $i++) {
    printf OUT "$format;\n",$phs_busfr[$i],$phs_busto[$i],$phs_bas[$i],
    $kvbase{$ltc_busfr[$i]},$kvbase{$ltc_busto[$i]},$phs_ref[$i],
    $phs_r[$i],$phs_x[$i],$phs_vma[$i],$phs_vmi[$i];
}
$i = $nphs;
if ($nphs >= 0) {
    printf OUT "$format];\n\n",$phs_busfr[$i],$phs_busto[$i],
    $phs_bas[$i],$kvbase{$ltc_busfr[$i]},$kvbase{$ltc_busto[$i]},
    $phs_ref[$i],$phs_r[$i],$phs_x[$i],$phs_vma[$i],$phs_vmi[$i];
}

# -----------------------------------------------------------------------
# write Twt.con
# -----------------------------------------------------------------------
$ntw >= 0 && print OUT "Twt.con = [ ...\n";
$format = "%4d %4d %4d $pbas 60 %5.2f %5.2f %5.2f " . "%8.5f " x 6 . 
    " 1 " . "%7.3f " x 9;
for ($i = 0; $i <= $ntw-1; $i++) {
    printf OUT "$format;\n", $bustw1[$i],$bustw2[$i],$bustw3[$i],
    $kvbase{$bustw1[$i]},$kvbase{$bustw2[$i]},$kvbase{$bustw3[$i]},
    $r12[$i],$r13[$i],$r23[$i],$x12[$i],$x13[$i],$x23[$i],
    $rt1a[$i],$rt1b[$i],$rt1c[$i],$rt2a[$i],$rt2b[$i],$rt2c[$i],
    $rt3a[$i],$rt3b[$i],$rt3c[$i];
}

$i = $ntw;
if ($ntw >= 0) {
    printf OUT "$format];\n\n", $bustw1[$i],$bustw2[$i],$bustw3[$i],
    $kvbase{$bustw1[$i]},$kvbase{$bustw2[$i]},$kvbase{$bustw3[$i]},
    $r12[$i],$r13[$i],$r23[$i],$x12[$i],$x13[$i],$x23[$i],
    $rt1a[$i],$rt1b[$i],$rt1c[$i],$rt2a[$i],$rt2b[$i],$rt2c[$i],
    $rt3a[$i],$rt3b[$i],$rt3c[$i];
}

#my ($nsvc,@svcbus,@svcbma,@svcbmi,@svcref);
# -----------------------------------------------------------------------
# write Svc.con
# -----------------------------------------------------------------------
$nsvc >= 0 && print OUT "Svc.con = [ ...\n";
$format = "%4d $pbas %8.2f 60 1 1 10 " . "%8.5f " x 3;
for ($i = 0; $i <= $nsvc-1; $i++) {
    printf OUT "$format;\n",$svcbus[$i],$kvbase{$svcbus[$i]},
    $svcref[$i],$svcbma[$i],$svcbmi[$i];
}
$i = $nsvc;
if ($nsvc >= 0) {
    printf OUT "$format];\n\n",$svcbus[$i],$kvbase{$svcbus[$i]},
    $svcref[$i],$svcbma[$i],$svcbmi[$i];
}


# -----------------------------------------------------------------------
# write Varname.bus
# -----------------------------------------------------------------------
for ($i = 0; $i <= $nbus; $i++) {
    $busname[$i] =~ s/\'//g;
    $busname[$i] = deblank ($busname[$i]);    
}
$nbus >= 0 && print OUT "Varname.bus = { ...\n";
$h = ($nbus+1) % 5;
if ($h == 0) {$h = 5;}
if (($nbus+1) > 5) {
    for ($i = 0; $i <= $nbus-$h; $i+=5) {
	print OUT "  '$busname[$i]'; '$busname[$i+1]'; " . 
	    "'$busname[$i+2]'; '$busname[$i+3]'; '$busname[$i+4]';\n";
    }
}
print OUT "  ";
for ($i = $nbus-$h+1; $i <= $nbus-1; $i++) {
    print OUT "'$busname[$i]'; ";
}
print OUT "'$busname[$nbus]'};\n\n";

# -----------------------------------------------------------------------
# close output data file
# -----------------------------------------------------------------------
close(OUT) || die "cannot close $ARGV[1]: $!\n";
print "Conversion completed.\n";

#########################################################################
#########################################################################

# -----------------------------------------------------------------------
# function for formatting data and/or assigning default values
# -----------------------------------------------------------------------
sub assign {
    my $param;
    $param = substr($_[0],$_[1],$_[2]);
    $param =~ s/\s*//g;
    unless ($param) {$param = "0";}
    #$param =~ s/^0*\.0*//;
    unless ($param =~ /[1-9]/) {$param = $_[3];}
    return $param;
}

# -----------------------------------------------------------------------
# function for formatting data and/or assigning default values
# -----------------------------------------------------------------------
sub deblank {
    my $name = $_[0];
    $name =~ s/( *)$//;
    return $name;
}

# -----------------------------------------------------------------------
# function for writing a separator
# -----------------------------------------------------------------------
sub carret {
    $verbose && print "-" x 50 . "\n";
}

# -----------------------------------------------------------------------
# function for reading comma separated data
# -----------------------------------------------------------------------
sub mysplit {
    my $string = $_[0];
    $string =~ s/\/(.*)//;
    my @mydata = split /\s*,\s*/, $string;
    if ($mydata[0] eq "") { shift(@mydata); } 
    return @mydata;
}

# -----------------------------------------------------------------------
# function for finding multiple PQ loads at the same bus
# -----------------------------------------------------------------------
sub pqidx {
    my $mypqbus = $_[0];
    my $check = 0;

    for ($i = 0; $i <= $npq; $i++) {
	if ($pqbus[$i] == $mypqbus) {
	    $check = $i; 
	    last;
	}	    
    }	
    return $check;
}

# -----------------------------------------------------------------------
# function for finding multiple PV generators at the same bus
# -----------------------------------------------------------------------
sub pvidx {
    my $mypvbus = $_[0];
    my $check = 0;

    for ($i = 0; $i <= $npv; $i++) {
	if ($pvbus[$i] == $mypvbus) {
	    $check = $i; 
	    last;
	}	    
    }	
    return $check;
}

# -----------------------------------------------------------------------
# function for finding multiple PV generators at the same bus
# -----------------------------------------------------------------------
sub isidx {
    my $myisbus = $_[0];
    my $check = 0;

    for ($i = 0; $i <= $nisl; $i++) {
	if ($islbus[$i] == $myisbus) {
	    $check = 1; 
	    last;
	}	    
    }	
    return $check;
}