基于牛顿-拉夫逊法的潮流计算例题编程报告

=-((inv(J))*x2); % x1 为所求△x 的列向量% 求节点电压新值,准备下一次迭代 for i=1:4 oa(i)=x1(i); oU(i)=x1(i+4)*U(i); end for i=1:4 a(i)=a(i)+oa(i); U(i)=U(i)+oU(i); end k=k+1; end k,U,a % 求节点注入功率 i=5; for j=1:5 P(i)=U(i)*U(j)*(G(i,j)*cos(a(i)-a(j))+B(i,j)*sin(a(i)-a(j)))+P(i); Q(i)=U(i)*U(j)*(G(i,j)*sin(a(i)-a(j))-B(i,j)*cos(a(i)-a(j)))+Q(i); end S(5)=P(5)+Q(5)*sqrt(-1); S% 求节点注入电流 I=Y*U' 四运行结果 1节点导纳矩阵 2经过五次迭代后的雅克比矩阵 3迭代次数以及节点电压的幅值和相角(弧度数) 4节点注入功率和电流五结果分析在这次学习和实际操作过程里:首先,对电力系统分析中潮流计算的部分特别是潮流计算的计算机算法中的牛顿- 拉夫逊法进行深入的研读,弄明白了其原理、计算过程、公式推导以及设计流程。牛顿- 拉夫逊法是求解非线性方程的迭代过程,其计算公式为 F J X ? ??,式中 J 为所求函数的雅可比矩阵; X?为需要求的修正值; F?为不平衡的列向量。利用 x (*)=x (k+1) +X?(k+1) 进行多次迭代,通过迭代判据得到所需要的精度值即准确值 x (*)。六结论通过这个任务,自己在 matlab 编程,潮流计算,word 文档的编辑功能等方面均有提高,但也暴漏出一些问题:理论知识储备不足,对 matlab 的性能和特点还不能有一个全面的把握,对 word 软件也不是很熟练,相信通过以后的学习能弥补这些不足,达到一个新的层次。