牛顿拉普逊源程序.docVIP

牛顿拉普逊源程序 #includestdio.h #includemath.h #define M 100 /*矩阵阶数*/ #define N 100 /*迭代次数*/ int i,j,k,a,b,c; /*循环控制变量*/ int t,l,g; /*中间变量*/ int n, /* 节点数 */ m, /* 支路数 */ pq, /* PQ节点数 */ pv; /* PV节点数 */ duidi; /*对地支路数*/ K; /*变比*/ double eps; /* 精度 */ double AA[M],BB[M],CC[M],DD[M],max,temp, rr,tt; /*中间变量*/ double mo,c1,d1,c2,d2,jiao; /*复数运算函数的返回值*/ double G[M][M],B[M][M],Y[M][M]; /*节点导纳矩阵中的实部、虚部及其模方值*/ double ykb[M][M],D[M],dU[M]; /*雅克比矩阵、不平衡量矩阵*/ struct jd /* 节点结构体 */ { int num,s; /* num为节点号，s为节点类型*/ double p,q,S,e,f,U,zkj,dp,dq,du,de,df; /*节点有功、无功功率，功率模值，电压纵、横分量，电压模值,阻抗角 牛顿--拉夫逊中功率不平衡量、电压不平衡量*/ } jd[M]; struct zhl /* 支路结构体 */ { int numb; /*numb为支路号*/ int p1,p2; /*支路的两个节点*/ double kx; /*变比*/ double r,x; /*支路的电阻与电抗*/ } zhl[M]; FILE *fp1,*fp2; void data() /* 读取数据 */ { int h,number; fp1=fopen(input.txt,r);if(fp1==NULL) { printf( can not open file !\n); exit(0); } fscanf(fp1,%d,%d,%d,%d,%d,%lf\n,n,m,pq,pv,duidi,eps); /*输入节点数,支路数,PQ节点数,PV节点数和精度*/ j=1;k=pq+1; for(i=1;i=n;i++) /*输入节点编号、类型、输入功率和电压初值*/ { fscanf(fp1,%d,%d,nu