pq分解法程序[借鉴].pdfVIP

PQ 分解法编程 版本一： C 语言 这里要提一下，本人使用的数据出自鞠平主编《电力工程》， PQ 分解法可能由于求逆程序 的精确度问题而比 matlab 迭代次数稍多， 可将 float 改为 dobule float 一试，若有不当之处， 还望包涵（另外，本人使用 c 语言编译器为 tc3.0 ） //PQ 分解法解潮流程序 #includestdio.h #includemath.h #define N 4 //节点数 #define n_PQ 2 //PQ 节点数 #define n_PV 1 //PV 节点数 #define n_br 5 // 串联支路数 #define PI 3.1415 void main() { void disp_matrix(float *disp_p,int disp_m,int disp_n); //矩阵显示函数 float Us[2*N]={1.0,0,1.0,0,1.05,0,1.05,0}; // 电压初值 float Ps[N]={0,-0.5,0.2}; // 有功初值 float Qs[N]={0,-0.3}; //无功初值 float G[N][N],B[N][N],B1[N-1][N-1],B2[n_PQ][n_PQ],invB1[N-1][N-1],invB2[n_PQ][n_PQ]; struct // 阻抗参数 { int nl; //左节点 int nr; //右节点 float R; // 串联电阻值 float X; // 串联电抗值 float Bl; // 左节点并联电纳 float Br; // 右节点并联电纳 }ydata[n_br]={ {1,2,0,0.1880,-0.6815,0.6040}, {1,3,0.1302,0.2479,0.0129,0.0129}, {1,4,0.1736,0.3306,0.0172,0.0172}, {3,4,0.2603,0.4959,0.0259,0.0259}, {2,2,0,0.05,0,0} }; float Z2; //Z^2=R^2+X^2 各串联阻抗值的平方 float U_amp[N],U_ang[N]; //e,f 存储电压的 x 轴分量和 y 轴分量 ,dfe 存储电压修正值 float dU_amp[N],dU_ang[N]; // 电压幅值和相角修正值 float mid1[N],mid2[N],dP[N-1],dQ[n_PQ]; //mid1 、 mid2 存储计算雅 克比行列式对角线元素的中间值 ,dS 存储 PQU 的不平衡量 float dPQU=1.0; //PQU 不平衡量最大值 int kk=0; // 迭代次数 int i,j,k; float t; // 形成导纳矩阵 fo