潮流计算源程序.docxVIP

vl.O/*矩阵阶数*//*迭代次数勺int i,j,k,a,b,c;int t,l,g;int n,m,pq，PV, duidi; vl.O /*矩阵阶数*/ /*迭代次数勺 int i,j,k,a,b,c; int t,l,g; int n, m, pq， PV, duidi; /*循环控制变量*/ 严中间变量*/ /*节点数*/ /*支路数*/ /*PQ节点数勺 /* PV节点数*/ 严对地支路数*/ double eps; double aa[M],bb[M],cc[M],dd[M],max,temp, rr,tt; double mo,c 1 ,d 1 ,c2,d2,fushuqiujiao; double G[M][M],B[M][M],Y[M][M]; 虚部及其模方值*/ double ykb[M][M],D[M],dU[M]; 阵勺 /*精度*/ /*中间变量*/ /*复数运算函数的返回值*/ 代节点导纳矩阵中的实部、 /咅雅克比矩阵、不平衡量矩 struct jd { int num,s; 类型*/ double p,q,S,e,f,U,zkj,dp,dq,du,de,df; 电压纵、横分量，电压模值，阻抗角， /*节点结构体*/ /* num为节点号，s为节点 /*节点有功、无功功率，功率模值, 牛顿?■拉夫逊中功率不平 衡量、电压不平衡量*/ }jd[M]; struct zhl /*支路结构体*/ 潮流计算源程序 /* ********************************************************** 文件名称：潮流计算程序.cpp 作者： 版本: 说明： 修改记录： ?丄“ ?土 ?丄■ ?丄■ ?丄■ ?土 ?丄■ ?￡? ?丄■ ?丄■ ?丄■ *￡? ?丄■ ?￡? ?丄■ *￡? ?丄■ ?卜?丄■ ?￡? ?丄■ ?卜?丄■ ?丄■ ?丄■ ?￡* ?丄■ ?丄■ ?丄■ ?￡* ?丄■ / ^7^ ^7w Tw ^7w Tw ^7^ 叫.?卜 #i ^7^ T ^7^ 叫.Tw Tw T T T #T% #T% #T% #T% #T% #T^ #T% ^T% 吓? #T% 吓? #T% 吓? #T% 吓? #T% 吓? *■ #7^ #T f #includestdio.h #includemath.h #includestdlib.h #define M 100 #define N 100 /*numb /*numb为支路号*/ /*支路的两个节点*/ /*支路的电阻与电抗*/ { int numb; intpl,p2; double r,x； } zhlLMJ; FILE *fpl,*fp2; 函数名称：void ReadData() 函数功能：读取数据 入口参数： 出口参数： 备注： void ReadData() { int h,number; fp 1 =fopen(,,input.txt,,,,,r);if(fp 1 ==NULL) 则跳出勺 { printf(u can not open file !\nu); exit(O); } fscanf(fp 1, %d,%d,%d,%d,%d,%lf\nH,n,m, pq,pv,duidi,eps); 支路数,PQ节点数,PV节点数对地支路数和精度*/ /*如果输入为空 /*输入节点数, j=l;k=pq+l; for(i=l;i=n;i4-+) 编号、类型、输入功率和电压初值勺 /*输入节点 fscanf(fp 1 ,“％d,%d”,number,h); if(h==l) /*类型 h=l 是PQ节点*/ { fscanf(fpl,”,％lf,%lf,%lf,%lf\n”,jd[j].p,jd[j]?q,jd[j].e,jd[j].f); jd[j].num=number; jdfj].s=h; j++； if(h==2) 是pv节点*/ /*类型h=2 fscanf(fpl,“,％lf,%lf\n”,jd[k].p,jd[k].U); jd[k].num=number; jd[k].s=h; jd[k].q=O; k++; /*类型h=3 /*类型h=3 /*输入支路 是平衡节点*/ { fscanf(fpl,”,％lf,%lf\n”,jd[n]?e,jd[n].f); jd[n].num=number; jd[n].s=h; } } fbr(i=l;i=m;i++) 阻抗勺 fscanf(lp 1 ；,%d,%d,%d,%lf,%lf\nn,zhl[i].numb,zhl [i] .p 1，zhl [i] .p2,zhl [i] ,r,zhl[i].x); fclose(fpl); if((