//牛拉法解潮流程序

#include
#include
#define?N?5??????????????//節(jié)點(diǎn)數(shù)
#define?n_PQ?3?????????//PQ節(jié)點(diǎn)數(shù)
#define?n_PV?1?????????//PV節(jié)點(diǎn)數(shù)
#define?n_br?5?????????//串聯(lián)支路數(shù)
void?main（）
{
void?disp_matrix（float?*disp_pint?disp_mint?disp_n）;???????//矩陣顯示函數(shù)
float?Us[2*N]={1.001.001.001.0501.050};?????????????????????//電壓初值
float?Ps[N]={-1.6-2-3.75};??????????????????????????????????????//有功初值
float?Qs[N]={-0.8-0-1-1.3};?????????????????????????????????????????????????//無功初值????????
float?G[N][N]B[N][N];?????????????????????????????????????????????????//各幾點(diǎn)電導(dǎo)電納
struct??????????????????????????????????????????????????????????????????????//阻抗參數(shù)
???{
????int?nl;?????????????//左節(jié)點(diǎn)
????int?nr;????????????//右節(jié)點(diǎn)
????float?R;??????????//串聯(lián)電阻值
????float?X;???????