33_IEEE33节点_ieee33_节点_电力系统潮流计算_源码_ieee33节点,ieee33节点系统数据资源-CSDN下载

共1个文件

m：1个

版权申诉

IEEE33节点

ieee33

电力系统潮流计算

5星 · 超过95%的资源 47 浏览量 2021-10-02 01:47:38 上传评论 15 收藏 1KB RAR 举报

在电力系统领域，潮流计算是一项基础且至关重要的任务，它涉及到电力网络中电压、电流、功率等关键参数的计算。本资源"33_IEEE33节点_ieee33_节点_电力系统潮流计算_源码"聚焦于利用牛顿-拉夫逊法对IEEE33节点电力系统的潮流进行计算，其特点是算法收敛速度快，能够高效解决实际问题。我们需要了解IEEE33节点系统。这是一个被广泛用于测试和研究电力系统潮流问题的标准模型。它包含了各种不同类型的元件，如发电机、变压器、线路、负荷等，可以模拟复杂的城市或工业区域电力网。通过分析这个系统，我们可以更好地理解和优化电力系统的运行状态。接着，我们来探讨牛顿-拉夫逊法。这是一种非线性方程求解的迭代方法，常用于电力系统的潮流计算。在电力系统中，潮流计算的目标是找到一组电压和电流值，使得所有设备（发电机、变压器、线路等）的功率平衡和电气约束都得以满足。牛顿-拉夫逊法通过迭代更新系统状态变量（如电压幅值和相角），直到达到一个平衡点，即所有的功率差和电压差接近于零，从而获得系统的稳态运行情况。该算法的基本步骤包括： 1. 初始化：设定初始的电压和功率估计值。 2. 建立雅可比矩阵：根据当前状态计算出系统中所有节点和支路的功率和电压的偏导数，形成雅可比矩阵。 3. 计算残差：比较实际功率与计算功率的差值，以及实际电压与计算电压的差值，得到残差向量。 4. 更新状态：利用雅可比矩阵和残差向量，计算下一迭代状态的修正值。 5. 检查收敛性：判断修正值是否足够小，若满足预设的收敛条件，则结束迭代，否则返回步骤2。牛顿-拉夫逊法的优势在于其快速的收敛特性，尤其对于大型电力系统，相比其他算法，它能在较少的迭代次数内找到解，节省了计算时间和资源。在压缩包中的源码，很可能是用一种编程语言（如Python、C++或Matlab）实现的牛顿-拉夫逊法，用于解决33节点电力系统的潮流问题。通过阅读和理解这段代码，不仅可以加深对算法的理解，还可以直接应用到实际项目中，或者作为进一步研究的基础。 "33_IEEE33节点_ieee33_节点_电力系统潮流计算_源码"资源提供了学习和实践电力系统潮流计算的宝贵材料，尤其是对牛顿-拉夫逊法的运用，有助于电力工程专业的学生和从业人员提升技能，解决实际问题。

资源详情

资源评论

资源推荐

收起资源包目录

33.rar （1个子文件）

jiedian_33.m 5KB

clear; clc; n=33; m=32; Isb=33; pr=0.00001; a=1; k=1; B1=[1 33 0.0922+0.0470i 0.000; 2 1 0.4930+0.2511i 0.000; 3 2 0.3660+0.1864i 0.000; 4 3 0.3811+0.1941i 0.000; 5 4 0.8190+0.7070i 0.000; 6 5 0.1872+0.6188i 0.000; 7 6 0.7114+0.2351i 0.000; 8 7 1.0300+0.7400i 0.000; 9 8 1.0440+0.7400i 0.000; 10 9 0.1966+0.0650i 0.000; 11 10 0.3744+0.1238i 0.000; 12 11 1.4680+1.1550i 0.000; 13 12 0.5416+0.7129i 0.000; 14 13 0.5910+0.5260i 0.000; 15 14 0.7463+0.5450i 0.000; 16 15 1.2890+1.7210i 0.000; 17 16 0.7320+0.5740i 0.000; 18 1 0.1640+0.1565i 0.000; 19 18 1.5042+1.3554i 0.000; 20 19 0.4095+0.4784i 0.000; 21 20 0.7089+0.9373i 0.000; 22 1 0.4512+0.3083i 0.000;% 23 22 0.8980+0.7091i 0.000; 24 23 0.8960+0.7011i 0.000; 25 5 0.2030+0.1034i 0.000; 26 25 0.2842+0.1447i 0.000; 27 26 1.0590+0.9337i 0.000; 28 27 0.8042+0.7006i 0.000; 29 28 0.5075+0.2585i 0.000; 30 29 0.9744+0.9630i 0.000; 31 30 0.3105+0.3619i 0.000; 32 31 0.3410+0.5302i 0.000; 7 20 2+2i 0.000; 8 14 2+2i 0.000; 11 21 2+2i 0.000; 17 32 0.5+0.5i 0.000; 24 28 0.5+0.5i 0.000 ];%支路阻抗 Y=zeros(n); for i=1:37 p=B1(i,1);%首段节点 q=B1(i,2);%末端节点 Y(p,q)=Y(p,q)-1/(B1(i,3)); Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/(B1(i,3)); Y(q,q)=Y(q,q)+1/(B1(i,3));%导纳矩阵 end G=real(Y); B=imag(Y); B2=[ 1 0 0 0 -0.0010 -0.0006 1 0; 2 0 0 0 -0.0009 -0.0004 1 0; 3 0 0 0 -0.0012 -0.0008 1 0; 4 0 0 0 -0.0006 -0.0003 1 0; 5 0 0 0 -0.0006 -0.0002 1 0; 6 0 0 0 -0.0020 -0.0010 1 0; 7 0 0 0 -0.0020 -0.0010 1 0; 8 0 0 0 -0.0006 -0.0002 1 0; 9 0 0 0 -0.0006 -0.0002 1 0; 10 0 0 0 -0.00045 -0.0003 1 0; 11 0 0 0 -0.0006 -0.00035 1 0; 12 0 0 0 -0.0006 -0.00035 1 0; 13 0 0 0 -0.0012 -0.0008 1 0; 14 0 0 0 -0.0006 -0.0001 1 0; 15 0 0 0 -0.0006 -0.0002 1 0; 16 0 0 0 -0.0006 -0.0002 1 0; 17 0 0 0 -0.0009 -0.0004 1 0; 18 0 0 0 -0.0009 -0.0004 1 0; 19 0 0 0 -0.0009 -0.0004 1 0; 20 0 0 0 -0.0009 -0.0004 1 0; 21 0 0 0 -0.0009 -0.0004 1 0; 22 0 0 0 -0.0009 -0.0005 1 0; 23 0 0 0 -0.0042 -0.0020 1 0; 24 0 0 0 -0.0042 -0.0020 1 0; 25 0 0 0 -0.0006 -0.00025 1 0; 26 0 0 0 -0.0006 -0.00025 1 0; 27 0 0 0 -0.0006 -0.0002 1 0; 28 0 0 0 -0.0012 -0.0007 1 0; 29 0 0 0 -0.0020 -0.0060 1 0; 30 0 0 0 -0.0015 -0.0007 1 0; 31 0 0 0 -0.0021 -0.0010 1 0; 32 0 0 0 0.00266 0.0009 1 0; 33 -1 0 0 0 0 1 0]; e=zeros(n,1); f=zeros(n,1); dP=zeros(n-1,1); dQ=zeros(n-1,1); Ps=zeros(n,1); Qs=zeros(n,1);%定义矩阵 for i=1:n e(i,1)=B2(i,7); f(i,1)=B2(i,8); Ps(i)=B2(i,5); Qs(i)=B2(i,6); end while a>pr Pi=zeros(n-1,1); Qi=zeros(n-1,1); E=zeros(n-1,1); F=zeros(n-1,1); for i=1:n-1 for j=1:n Pi(i)=Pi(i)+e(i)*(G(i,j)*e(j)-B(i,j)*f(j))+f(i)*(G(i,j)*f(j)+B(i,j)*e(j)); Qi(i)=Qi(i)+f(i)*(G(i,j)*e(j)-B(i,j)*f(j))-e(i)*(G(i,j)*f(j)+B(i,j)*e(j)); end dP(i)=Ps(i)-Pi(i); dQ(i)=Qs(i)-Qi(i); end jacobi=zeros(64,64); A=zeros(n-1,1); S=zeros(n-1,1); for i=1:n-1 for j=1:n A(i)=A(i)+G(i,j)*e(j)-B(i,j)*f(j); S(i)=S(i)+G(i,j)*f(j)+B(i,j)*e(j); end end for i=1:n-1 for j=1:n-1 if i~=j jacobi(2*i-1,2*j-1)=0-G(i,j)*e(i,1)-B(i,j)*f(i,1); jacobi(2*i-1,2*j)=B(i,j)*e(i,1)-G(i,j)*f(i,1); jacobi(2*i,2*j-1)=B(i,j)*e(i,1)-G(i,j)*f(i,1); jacobi(2*i,2*j)=G(i,j)*e(i,1)+B(i,j)*f(i,1); else jacobi(2*i-1,2*j-1)=0-A(i)-G(i,i)*e(i)-B(i,i)*f(i); jacobi(2*i-1,2*j)=0-S(i)+B(i,i)*e(i)-G(i,i)*f(i); jacobi(2*i,2*j-1)=S(i)+B(i,i)*e(i)-G(i,i)*f(i); jacobi(2*i,2*j)=0-A(i)+G(i,i)*e(i)+B(i,i)*f(i); end end end dW=zeros(2*n-2,1); for i=1:n-1 dW(2*i-1)=dP(i); dW(2*i)=dQ(i); end dV=-inv(jacobi)*dW; de=zeros(n-1,1); df=zeros(n-1,1); for i=1:n-1 de(i)=dV(2*i-1); df(i)=dV(2*i); end for i=1:n-1 e(i)=e(i)+de(i); f(i)=f(i)+df(i); end a=max(max(abs(dP)),max(abs(dQ))); k=k+1; end V=zeros(n,1); for i=1:n V(i)=e(i)+sqrt(-1)*f(i); end disp('节点电压V='); disp(V); disp(real(V)); I=zeros(37,1); U=zeros(37,1); R1=zeros(37,1); for i=1:37 p1=B1(i,1); q1=B1(i,2); U(i)=V(q1)-V(p1); R1(i,1)=real(B1(i,3)); I(i)=abs(U(i)/R1(i,1)); end disp(I);