电脑桌面
添加小米粒文库到电脑桌面
安装后可以在桌面快捷访问

牛顿-拉夫逊迭代法极坐标潮流计算C语言程序

牛顿-拉夫逊迭代法极坐标潮流计算C语言程序_第1页
1/21
牛顿-拉夫逊迭代法极坐标潮流计算C语言程序_第2页
2/21
牛顿-拉夫逊迭代法极坐标潮流计算C语言程序_第3页
3/21
/*最大矩阵阶数*//*迭代次数*//*循环控制变量*//*中间变量*//*节点数*//*支路数*//*PQ 节点数*//*PV 节点数*//*迭代精度*/tt;/*中间变量*//*复数运算函数的返回值*//*节点导纳矩阵中/*雅克比矩阵、不平/*利用牛顿—拉夫逊迭代法(极坐标形式),计算复杂电力系统潮流,具有收敛性好,收敛速度快等优点.所有参数应归算至标幺值下./*可计算最大节点数为 100,可计算 PQ,PV/*可计算非标准变比和平行支路*/#includevstdio。h〉#include〈math。h〉#include〈stdlib.h〉#defineM100#defineNl100inti,j,k,a,b,c;intt,l;doubleP,Q,H,J;intn,m,pq,pv;doubleeps;doubleaa[M],bb[M],cc[M],dd[M],max,rrdoublemo,c1,d1,c2,d2;doubleG[M][M],B[M][M],Y[M][M];的实部、虚部及其模方值*/doubleykb[M][M],D[M],d[M],dU[M];衡量矩阵*/inth,number;structjd{intnum,ty;型*/doublep,q,S,U,zkj,dp,dq,du,dj;模值,电压模值,阻抗角平衡量、电压修正量*/}jd[M];structzl{intnumb;intp1,p2;doublekx;doubler,x;}zl[M];FILE*fp1,*fp2;voiddata(){/*节点结构体*//*num 为节点号,ty 为节点类/*节点有功、无功功率,功率牛顿——拉夫逊中功率不/*支路结构体*//*numb 为支路号*//*支路的两个节点*//*非标准变比*//*支路的电阻与电抗*//*读取数据*/fpl=fopen("input。txt","r");fscanf(fp1,”%d,%d,%d,%d,%lf\n",&n,&m,&pq,&pv,&eps);/*输入节点数,支路数,PQ 节点数,PV 节点数和迭代精度*/for(i=l;i<=n;i++)/*输入节点编号、类型、输入功率和电压初值*/{fscanf(fpl,”%d,%d”,&number,&h);if(h==1)/*类型 h=1是 PQ 节点*/{fscanf(fp1,"%lf,%lf,%lf,%lf\n",&jd[i].p,&jd[i].q,&jd[i]。U,&jd[i].zkj);jd[i]。num=number;jd[i]。ty=h;}if(h==2)/*类型 h=2是 pv 节点*/{fscanf(fp1,”,%lf,%lf,%lf\n”,&jd[i].p,&jd[i]。U,&jd[i].zkj);jd[i].num=number;jd[i].ty=h;jd[i].q=-1。567;}if(h==3)/*类型 h=3是平衡节点*/{fscanf(fp1,",%lf,%lf\n”,&jd[i]。U,&jd[i]。zkj);jd[i].num=number;jd[i]。ty=h;}}for(i=1;i〈=m;i++)/*输入支路阻抗*/fscanf(fp1,”%d,%lf,%d,%d,%lf,%lf\n”,&zl[i]。numb,&zl[i]。kx,&zl[i].p1,&zl[i].p2,&zl[i].r,&zl[i]。x);fclose(fp1);if((fp2=fopen(”output.txt",”w"))==NULL){print...

1、当您付费下载文档后,您只拥有了使用权限,并不意味着购买了版权,文档只能用于自身使用,不得用于其他商业用途(如 [转卖]进行直接盈利或[编辑后售卖]进行间接盈利)。
2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。
3、如文档内容存在违规,或者侵犯商业秘密、侵犯著作权等,请点击“违规举报”。

碎片内容

牛顿-拉夫逊迭代法极坐标潮流计算C语言程序

确认删除?
VIP
微信客服
  • 扫码咨询
会员Q群
  • 会员专属群点击这里加入QQ群
客服邮箱
回到顶部