基于C++语言的牛顿-拉夫逊法潮流计算

《基于C++语言的牛顿-拉夫逊法潮流计算》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、电力系统课程设计实验报告姓名：李贺贺学号：0912002053班级：电093电093李贺贺1图中节点0,1,2,3为PQ节点，4为平衡节点要求：使用牛顿法潮流计算、计算机语言为C++使用软件：VisualStudio2010过程：流程图如下：开始计算导纳矩阵计算节点功率计算节点功率失配量计算雅可比矩阵根据修正方程算出节点电压幅值和相角N最大功率残差小于?Y结束电093李贺贺21.计算导纳矩阵因为有5个节点，所以导纳矩阵为5阶。1⎧−(?≠?)⎪r??+???????=1⎨�(?=?)⎪r??+????⎩??

2、?其中???表示所有与节点m相连的支路。设电压初始值为U?=1(?=0,1,2,3,4)，?i=0(?=0,1,2,3,4)。2.计算节点功率2∙?+?�?(?cos?+?sin?)⎧??=??????????????⎪???⎨?=−?2∙?+?�?(?sin?−?cos?)⎪??????????????⎩???其中???表示与节点?邻接的所有节点；???=?i−??为节点电压相角差。3.计算节点功率失配量∆P?=???−??�∆??=???−??4.计算雅可比矩阵???−?????−?????−??H=�

3、�，N=��，J=(???−?????−??)，L=???−?????−?????−?????−??−?∙?2−?(?=?)???????=�??∙??�???sin???−???cos???�(?≠?)?∙?2+?(?=?)???????=�??∙??�???cos???+???sin???�(?≠?)−?∙?2+?(?=?)????J??=�−??∙??�???cos???+???sin???�(?≠?)−?∙?2+?(?=?)???????=�??∙??�???sin???−???cos???�(?≠

4、?)5.修正方程∆??????−?????−?????−??∆???∆???�???−?????−?????−??�∙�??�=�∆???�???−?????−?????−??∆???��∆??求出∆θ和∆U。6.修正电压幅值和相角U?=??+∆U�??=??+∆θ7.判断收敛max��∆??，∆??��

5、ude#include#include#defineM10usingnamespacestd;intpower(intm,intn){intsum=1;if(n==0&&m!=0)return1;else{for(inti=0;i

6、{A=a[0][0]*a[1][1]-a[0][1]*a[1][0];returnA;}elseif(n==3){A=a[0][0]*a[1][1]*a[2][2]+a[0][1]*a[1][2]*a[2][0]+a[0][2]*a[1][0]*a[2][1]-a[0][2]*a[1][1]*a[2][0]-a[0][1]*a[1][0]*a[2][2]-a[0][0]*a[1][2]*a[2][1];returnA;}else{A=0;doublec[M][M];for(inti=0;i

7、i]=a[i][0]*power(-1,i+1+1);for(inti=0;i

8、;i++)for(intj=0;j