题目
两节点潮流计算题目如下,这个题目用Newton-Raphson(牛拉法)求解,MATLAB编程
分析
- 写出节点导纳矩阵Y
- 计算节点2的注入功率和方程
- 求出雅克比(Jacobi)矩阵
代码
Matlab源代码如下
%Newton-Raphson iteration program to solve two-bus system
clc;
clear;
tic;
%set initial values:V2---x1 theta2---x2 V1---balance node voltage amplitude
V1=1.05;
x1=1;x2=0;
P=1;Q=0.5;
xk=[x1;x2];
xk_1=[0;0];
xList=zeros(20,2);%存放每次迭代后的x1 x2
xList(1,:)=xk';
delta_x=xk_1-xk;
k=0;%记录迭代次数
while max(abs(delta_x))>10e-5 %判断是否收敛,误差为10e-5
%计算当前迭代次数k的f(x)
f0=[-P+5*V1*x1*cos(x2)-15*V1*x1*sin(x2)-5*x1^2;