ieee33节点仿真及潮流计算.zip

  • 风过无痕1
    了解作者
  • matlab
    开发工具
  • 28.1KB
    文件大小
  • zip
    文件格式
  • 0
    收藏次数
  • 1 积分
    下载积分
  • 1
    下载次数
  • 2021-12-13 15:34
    上传日期
IEEE33节点潮流计算程序可以参考下。
ieee33节点仿真及潮流计算.zip
  • ieee33节点仿真及潮流计算
  • mieee33.m
    2.3KB
  • IEEE33.slx
    27.1KB
内容介绍
clear all; Bus=[1,0,0 ; 2,100, 60; 3,90,40; 4,120,80; 5,60,30; 6,60,20; 7, 200, 100 ; 8,200, 100 ; 9, 60, 20; 10,60, 20; 11,45,30; 12,60, 35; 13,60, 35; 14,120,80; 15,60, 10; 16,60, 20; 17,60, 20; 18,90, 40; 19,90, 40; 20,90, 40; 21,90, 40; 22,90, 40; 23,90, 40; 24,420, 200; 25,420, 200 ; 26,60, 25; 27,60 ,25; 28,60, 20; 29,120, 70; 30,200,600 ; 31,150, 70; 32,210, 100 ; 33,60,40; ]; Branch=[1 ,1 ,2 ,0.0922,0.0407; 2 ,2 ,3 ,0.4930,0.2511; 3 ,3 ,4 ,0.3660,0.1864; 4 ,4 ,5 ,0.3811,0.1941; 5 ,5 ,6 ,0.8190,0.7070; 6 ,6 ,7 ,0.1872,0.6188; 7 ,7 ,8 ,0.7144,0.2351; 8 ,8 ,9 ,1.0300,0.7400; 9 ,9 ,10,1.0440,0.7400; 10,10,11,0.1966,0.065; 11,11,12,0.3744,0.1238; 12,12,13,1.4680,1.1550; 13,13,14,0.5416,0.7129; 14,14,15,0.5910,0.5260; 15,15,16,0.7463,0.5450; 16,16,17,1.2890,1.7210; 17,17,18,0.7320,0.5740; 18,2, 19,0.1640,0.1565; 19,19,20,1.5042,1.3554; 20,20,21,0.4095,0.4784; 21,21,22,0.7089,0.9373; 22,3, 23,0.4512,0.3083; 23,23,24,0.8980,0.7091; 24,24,25,0.8960,0.7011; 25,6, 26,0.2030,0.1034; 26,26,27,0.2842,0.1447; 27,27,28,1.0590,0.9337; 28,28,29,0.8042,0.7006; 29,29,30,0.5075,0.2585; 30,30,31,0.9744,0.9630; 31,31,32,0.3105,0.3619; 32,32,33,0.3410,0.5302; ]; [busnum,row]=size(Bus); [branchnum,row]=size(Branch); soubus=Branch(:,2); mobus=Branch(:,3); Vbus=ones(busnum,1); Vbus(:,1)=12.66; Vbus1=Vbus; Ploss=zeros(busnum,1); Qloss=zeros(busnum,1); e=1; k=0; Branch1=Branch; n=1; while ~isempty(Branch1) m=1; [s,row]=size(Branch1); while s>0 t=find(Branch1(:,2)== Branch1(s,3)); if isempty(t) T1(n,:)= Branch1(s,:); n=n+1; else T2(m,:)= Branch1(s,:); m=m+1; end; s=s-1; end; Branch1=T2; T2=[]; end; while e>1.0e-05 P=zeros(busnum,1); Q=zeros(busnum,1); for s=1:branchnum i=T1(s,2); j=T1(s,3); R=T1(s,4); X=T1(s,5); Pload=Bus(j,2); Qload=Bus(j,3); II=((Pload+P(j))^2+(Qload+Q(j))^2)/(Vbus(j)^2*1000); Ploss(i,j)=II*R; Qloss(i,j)=II*X; P(i,j)=Pload+Ploss(i,j)+P(j); Q(i,j)=Qload+Ploss(i,j)+Q(j); P(i)=P(i)+P(i,j); Q(i)=Q(i)+Q(i,j); end; for s=branchnum:-1:1 i=T1(s,2); j=T1(s,3); R=T1(s,4); X=T1(s,5); Vbus(j)=(Vbus(i)-(P(i,j)*R+Q(i,j)*X)/(Vbus(i)*1000))^2+((P(i,j)*X-Q(i,j)*R)/(Vbus(i)*1000))^2; Vbus(j)=sqrt(Vbus(j)); end; e=max(abs(Vbus1-Vbus)); Vbus1=Vbus; k=k+1; end
评论
    相关推荐
    • 潮流计算程序.rar
      潮流计算程序,带数据文件,未进行稀疏技术优化
    • 潮流计算程序.zip
      这是基于matlab编程实现的潮流计算,包括PQ计算、雅克比计算
    • 潮流计算程序9.zip
      九节点的潮流计算matlab程序,可以用来进行电力系统的潮流计算
    • MATLAB潮流计算程序
      此软件包是基于MATLAB2010平台编写的潮流计算软件,用于计算电力系统潮流,针对3机9节点,及经典39节点进行验证结果正确。压缩包里有代码的详细介绍。
    • matlab潮流计算程序
      matlab潮流计算程序matlab潮流计算程序潮流计算程序
    • 潮流计算程序.rar
      用matlab平台编写潮流计算程序的具体过程,通过实例证明基于牛顿-拉夫逊法的潮流计算程序具有收敛速度快、占用内存小的优点
    • matlab潮流计算程序
      这是matlab编写的33节点潮流计算软件。
    • 电力系统潮流计算程序
      潮流计算的软件,可以计算节点电流电压,用标幺值表示。算法有些复杂
    • 风机潮流计算程序
      这是个通用的快速解耦法的潮流计算程序,对于MATLAB编程的初学者特别有用。
    • 开式网络潮流计算程序
      程序测试使用方法:在E盘根目录下建立输入文件: 输入文件名:input.txt;...程序功能说明:本程序可以计算任意长度线型开始网络潮流; 支持多电压等级下的计算; 可在除供电节点外任意节点引出负载;