潮流计算源程序

《潮流计算源程序》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、vl.O/*矩阵阶数*//*迭代次数勺inti,j,k,a,b,c;intt,l,g;intn,m,pq，PV,duidi;/*循环控制变量*/严中间变量*//*节点数*//*支路数*//*PQ节点数勺/*PV节点数*/严对地支路数*/doubleeps;doubleaa[M],bb[M],cc[M],dd[M],max,temp,rr,tt;doublemo,c1,d1,c2,d2,fushuqiujiao;doubleG[M][M],B[M][M],Y[M][M];虚部及其模方值*/doubleykb[M][M],D[M],dU[M];阵勺/*精度*//*中间变量*//*复

2、数运算函数的返回值*/代节点导纳矩阵中的实部、/咅雅克比矩阵、不平衡量矩structjd{intnum,s;类型*/doublep,q,S,e,f,U,zkj,dp,dq,du,de,df;电压纵、横分量，电压模值，阻抗角，/*节点结构体*//*num为节点号，s为节点/*节点有功、无功功率，功率模值,牛顿•■拉夫逊中功率不平衡量、电压不平衡量*/}jd[M];structzhl/*支路结构体*/潮流计算源程序/***********************************************************文件名称：潮流计算程序.cpp作者：版本:说明：修

3、改记录：•丄“•土•丄■•丄■•丄■•土•丄■«£••丄■•丄■•丄■*£••丄■«£••丄■*£••丄■•卜•丄■«£••丄■•卜•丄■•丄■•丄■•£*•丄■•丄■•丄■•£*•丄■/^7^^7w^7^^7^叫.#T%#T%#T%#T%#T%#T^#T%^T%吓・#T%吓・#T%吓・#T%吓・#T%吓・*■#7^#T>f#include#include#include#defineM100#defineN100/*numb为支路号*//*支路的两个

4、节点*//*支路的电阻与电抗*/{intnumb;intpl,p2;doubler,x；}zhlLMJ;FILE*fpl,*fp2;函数名称：voidReadData()函数功能：读取数据入口参数：出口参数：备注：voidReadData(){inth,number;fp1=fopen(,,input.txt,,,,,r");if(fp1==NULL)则跳出勺{printf(ucannotopenfile!u);exit(O);}fscanf(fp1,"%d,%d,%d,%d,%d,%lfH,&n,&m,&pq,&pv,&duidi,&eps);支路数,PQ节点数,PV

5、节点数对地支路数和精度*//*如果输入为空/*输入节点数,j=l;k=pq+l;for(i=l;i<=n;i4-+)编号、类型、输入功率和电压初值勺/*输入节点fscanf(fp1,“％d,%d”,&number,&h);if(h==l)/*类型h=l是PQ节点*/{fscanf(fpl,”,％lf,%lf,%lf,%lf”,&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=2fscanf(fpl,“,％lf,%lf”,&jd[k].p,&jd

6、[k].U);jd[k].num=number;jd[k].s=h;jd[k].q=O;k++;/*类型h=3/*输入支路if(h==3)是平衡节点*/{fscanf(fpl,”,％lf,%lf”,&jd[n]・e,&jd[n].f);jd[n].num=number;jd[n].s=h;}}fbr(i=l;i<=m;i++)阻抗勺fscanf(lp1；,%d,%d,%d,%lf,%lfn,&zhl[i].numb,&zhl[i].p1，&zhl[i].p2,&zhl[i],r,&zhl[i].x);fclose(fpl);if((fp2=fopen("output.t

7、xt",HwH))==NULL){printf(Hcannotopenfile!n);exit(O);fprintf(fp2/200801000705M);fprintf(fp2；*fprintf(fp2,u**********原始数据E);潮流上机实习华北电力大学电气化0807段春姝fprintf(fp2,n节点数:％d支路数:％dPQ节点数:％dPV节点数:％d对地支路数:％d精度:%fn,n,m,pq,pv,duidi,eps);fprintf(fp2,"n);