PQ 分解法编程 版本一：C 语言 这里要提一下，本人使用的数据出自鞠平主编《电力工程》，PQ 分解法可能由于求逆程序的精确度问题而比matlab 迭代次数稍多，可将float 改为dobule float 一试，若有不当之处，还望包涵（另外，本人使用c 语言编译器为tc3.0） //PQ 分解法解潮流程序 #include #include #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; //形成导纳矩阵-------------------------------------------------------------------------------------------------- for(i=0;i
