function?[mpc]?=?bpa2mpc（bpafile）
%?BPA2MPC?Convert?BPA?format?file?to?Matpower?mpc?format?file
%?X類型母線（自動(dòng)控制投切電容電抗器）暫未考慮
%?buslist中母線名稱格式為:母線名+額定電壓等級(jí)因?yàn)閎pa中允許同一母線名對(duì)應(yīng)多個(gè)電壓等級(jí)
%?抽頭可調(diào)節(jié)變壓器的變比固定為設(shè)定變比
%?非對(duì)稱線路（兩側(cè)電納不相等）的情況按對(duì)稱線路考慮，兩側(cè)電納均取左側(cè)電納
%?mpc中沒有線路接地電導(dǎo)G和變壓器等效電導(dǎo)G
%?線路并聯(lián)高抗（L+支路）加到線路對(duì)地電納中b=-2*Mvar/U^2（只能認(rèn)為兩側(cè)高抗相等），也可以加到線路首末的母線上（第六列）感性為負(fù)。
%?如果提前用BPA計(jì)算了系統(tǒng)潮流，則兩端直流線路按其指定輸送功率等效為一正一負(fù)的兩個(gè)負(fù)荷，注意直流線路可以傳輸無功（P，Q），多端直流線路沒有考慮
%?mpc.branch中第14列存回路標(biāo)志，回路標(biāo)志是數(shù)字或字符的ascll碼
%?L卡中線路的額定電流不是線電流，而是一個(gè)等效的虛擬電流，其值為sqrt（3）倍線電流（待定）
%?case118.m?由118bpa.dat轉(zhuǎn)換的118mpc，用BPA計(jì)算的118bpadat?三者結(jié)果相同
%?By?Shijun?Tian?Xi‘a(chǎn)n?Jiaotong?University，tsjguoke@163.com?2017/10/20

%?Define?Default?Values
defaultVmax?=?1.052;
defaultVmin?=?0.95;
defaultQgmax?=?9999.0;
defaultQgmin?=?-9999.0;
defaultPgmin?=?0.0;
baseMVA?=?100;?
%?讀入BPA文件的所有內(nèi)容
bpa_info?=?cell（500001）;?%?默認(rèn).dat文件最大行數(shù)
bpa_info（:）?=?{blanks（80）};?%?數(shù)據(jù)行最大列數(shù)為80
disp（[‘讀入原始BPA潮流文件‘?bpafile?‘的內(nèi)容...‘]）;
%?打開BPA潮流文件
fid?=?fopen（bpafile?‘r‘‘n‘‘GBK‘）;
%?將所有內(nèi)容讀入到字符串集合中
ii?=?1;
while?（~feof（fid））
????newline?=?fgetl（fid）;
????bpa_info{ii}（1:length（newline））?=?newline;
????ii?=?ii?+?1;
end
fclose（fid）;%??關(guān)閉BPA潮流文件
bpa_info（ii:end）?=?[];?%釋放多余的空間
%?刪除注釋語句
temp=1;
for?ii?=?1?:?length（bpa_info）
????thisline?=?bpa_info{ii};
????if?（?strcmp（thisline（1）?‘.‘）?||?strcmp（thisline（1）?‘?‘）?）
????????deleteLine（temp）=ii;
????????temp=temp+1;
????end
end
bpa_info（deleteLine）=[];
%?統(tǒng)計(jì)某類型節(jié)點(diǎn)或支路數(shù)量
%?經(jīng)統(tǒng)計(jì)nwsgdata.dat中TP=0，BX=0，E=0LD=22LM=0+A=1
Num=0;
for?ii?=?1?:?length（bpa_info）
????thisline?=?bpa_info{ii};
????if?（strcmp（thisline（1:1）?‘+‘））
????????Num=Num+1;
????end
end
%?得到母線個(gè)數(shù)，支路個(gè)數(shù)，基準(zhǔn)容量
busNumA?=?0;?%?交流母線個(gè)數(shù)
branchNumA?=?0;?%?交流支路個(gè)數(shù)
busNumD?=?0;?%?直流母線個(gè)數(shù)
branchNumD?=?0;?%?直流支路個(gè)數(shù)
disp（‘獲得電網(wǎng)基本參數(shù)（母線個(gè)數(shù)，支路個(gè)數(shù)，基準(zhǔn)容量）...‘）;
for?ii?=?1?:?length（bpa_info）
????thisline?=?bpa_info{ii};
????%?獲得基準(zhǔn)容量
????if?（strcmp（thisline（1?:?9）?‘/MVA_base‘））
????????baseMVA?=?str2double（thisline（11?:?13））;
????end
????%?統(tǒng)計(jì)交流母線個(gè)數(shù)
????if?（strcmp（thisline（1）?‘B‘）?&&??~strcmp（thisline（2）?‘D‘）?&&?~strcmp（thisline（2）?‘M‘））
????????busNumA?=?busNumA?+?1;
????end
????%?統(tǒng)計(jì)交流支路個(gè)數(shù)
????if?（strcmp（thisline（1:2）?‘L?‘）?||?strcmp（thisline（1）?‘T‘）?||?strcmp（thisline（1）?‘E‘）?）
????????branchNumA?=?branchNumA?+?1;
????end
????%?統(tǒng)計(jì)直流母線個(gè)數(shù)
????if?（strcmp（thisline（1:2）?‘BD‘）?||??strcmp（thisline（1:2）?‘BM‘））
????????busNumD?=?busNumD?+?1;
????end
????%?統(tǒng)計(jì)直流支路個(gè)數(shù)
????if?（strcmp（thisline（1:2）?‘LD‘）?||?strcmp（thisline（1:2）?‘LM‘））
????????branchNumD?=?branchNumD?+?1;
????end
end
%?母線名稱可能出現(xiàn)漢字，為保證字符長度對(duì)應(yīng)bpa中的長度，要在漢字后面補(bǔ)充空格
for?ii?=?1?:?length（bpa_info）
????thisline?=?bpa_info{ii};
????if?（?~strcmp（thisline（1）?‘（‘）?&&?~strcmp（thisline（1）?‘/‘）?&&?~strcmp（thisline（1）?‘>‘）?）?%排除控制行，所有包含母線名稱的數(shù)據(jù)行都是7：14
????????temp?=?unicode2native（thisline（7:14）‘GBK‘）;
????????numChn?=?length（temp）-8;%?母線名中包含的中文字符數(shù)?0<=numChn<=4
????????thisline=[thisline（1:14-numChn）blanks（numChn）thisline（15-numChn:end）];
????end
????if?（strcmp（thislin

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-04-26?16:07??bpa2mpc\
?????文件???????94541??2018-04-26?16:07??bpa2mpc\118bpa.dat
?????文件???????21974??2018-04-26?16:07??bpa2mpc\bpa2mpc.m
?????文件???????33539??2018-04-26?16:07??bpa2mpc\case118.m
?????文件??????????77??2018-04-26?16:07??bpa2mpc\main.m