资源描述：
《c 编写潮流上机程序（c writing flow program）》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、c++编写潮流上机程序（C++writingflowprogram）#includesstdio.h><#includesmath.h><#includesstdlib.h><#definem矩阵阶数100/**/#definen100迭代次数/**/inti,j,k,a,b,c;/*循环控制变量*/intt,l,g中间变量;/**/intn;/*节点数*/m,支路数/**/fp,pq节点数/**/minutes,pv节点数/**/duidi对地支路数;/**/double精度eps;/**/doubleaa[m],bb[m]
2、,cc[m],dd[m],max,time,rr,中间变量);/**/doublemo,c1,d1,c2,d2,复数运算函数的返回值jiao;/**/doubleg[m][m],b[m][m],y[m][m];/*节点导纳矩阵中的实部、虚部及其模方值*/doubleykb[m][m],d[m],du[m];/*雅克比矩阵、不平衡量矩阵*/struct节点结构体jd/**/(intnum,d;/*num为节点号,s为节点类型*/doublep,q,s,e,f,u,zkj,dp,dq,du,de,df;/*节点有功、无功功率,功率
3、模值,电压纵、横分量,电压模值,阻抗角牛顿-拉夫逊中功率不平衡量、电压不平衡量*/}jd[m].struct支路结构体zhl/**/{intnumbnumb为支路号;/**/intp1,p2支路的两个节点;/**/doubler,x;/*支路的电阻与电抗*/}zhl[m].fp1file*,*fp2;voiddate(读取数据)/**/{inth,number;fp1=fopen("input.txt","r");if(fp1==null){printf("cannotopenfile!");exit(0);}fscanf
4、(fp1,"%d%d%d%d%d%lf",&n&m,&fp&minutes,&duidi,&eps);/*输入节点数,支路数,pq节点数,pv节点数和精度*/j=1;k=pq+1;for(i=1,the<=n;i++)/*输入节点编号、类型、输入功率和电压初值*/{fscanf(fp1,"%d%d"&number&h).if(b==1)/*类型h=1是pq节点*/{fscanf(fp1,,"%lf,lf,lf%,%,%lf",&jd[j].p,&jd[j]-he'ssaying,&jd[j]and[j].f&jd);j
5、d[j].num=number;jd[j].s=h;j++.}if(b==2)/*类型h=2是pv节点*/{fscanf(fp1,,"%lf%lf",&jd[k].p,&jd[k].u);jd[k].num=number;jd[k].s=h;jd[k]-he'ssaying=0;k++;}if(a==3)=3是平衡节点类型h/**/{fscanf(fp1,"%%lf,lfin,&jd.jd&[n],[n].f);jd=[n].numnumber;jd.s[n]=h;}}for(i=1;i<=m;i++)输入支路阻抗/**/
6、fscanf(fp1"%d%d%d%,%,lf,lfin&[i]",zhl.numb,&zhl[i]&[i].p1,zhl.p2,&zhl[i]&[i].r,zhl.x);fclose(fp1);if((fp2=fopen("output.txt","w"))==null){printf("cannotopenfile!in');exit(0);}fprintf(fp2"潮流上机实习onn,n');fprintf(fp2,"班级电力实07inn:");fprintf(fp2,"姓名刘波磊inn:");fprintf(fp2,"
7、学号200709030114inn:");fprintf(fp2"原始数据onn,n');fprintf(fp2,"n");fprintf(fp2,"节点数支路数:%d%d%d:pq节点数:pv节点数:%d:%d:%对地支路数精度fin","n,m,pq,pv,duidi,eps);fprintf(fp2,"n");for(i=1;i<=pq;i++)fprintf(fp2,"节点pq节点:%d%d]=[p,q][%d%d%di=n","jd[i][].num.num,jd(jd),[],[],.num.p,jd(jd),[]
8、.q);for(i=pq+1;i<=pq+pv;i++)fprintf(fp2,"节点pv节点:%d%d%p[f]=[u]=%d%d%d%初值q[f]=in,jd[i][].num.num,jd(jd),[],[],.num.p,jd(jd),[],[],.num.u,jd(
