• 伊万拉基蒂奇
    了解作者
  • matlab
    开发工具
  • 74KB
    文件大小
  • rar
    文件格式
  • 0
    收藏次数
  • 10 积分
    下载积分
  • 15
    下载次数
  • 2017-06-18 15:06
    上传日期
牛顿拉夫逊潮流计算的函数文件,可以计算系统的潮流,考虑到了有分接头的变压器。
牛顿拉夫逊法计算潮流.rar
  • 牛顿拉夫逊法计算潮流
  • Initial_e_f.m
    98B
  • Node.xls
    18.5KB
  • Node33.xls
    20.5KB
  • initial_e_f8.xls
    18KB
  • Ratio.xls
    18KB
  • Power_flow_computation.m
    3.2KB
  • ReadData.m
    1.3KB
  • Stopping_criterion.xls
    17.5KB
  • Stopping.m
    106B
  • Sensibility.m
    2.4KB
  • initial_e_f33.xls
    19.5KB
  • Branch8.xls
    18.5KB
  • Branch.xls
    18KB
  • Admittance_matrix.m
    3.8KB
  • initial_e_f.xls
    18KB
  • Branch33.xls
    21.5KB
  • Ratio33.xls
    18KB
  • Node8.xls
    18.5KB
  • Per_unit.m
    301B
内容介绍
function [AdmittanceMatrix,Admittance_of_earth]=Admittance_matrix(S_base,U_base) % Node = % 节点类型 节点编号 节点电压实部 节点电压虚部 电压有效值 节点注入有功 节点注入无功 % 0.0000 3.0000 NaN NaN NaN -1.6000 -0.8000 % 0.0000 1.0000 NaN NaN NaN -2.0000 -1.0000 % 0.0000 2.0000 NaN NaN NaN -3.7000 -1.3000 % 1.0000 4.0000 NaN NaN 1.0500 5.0000 NaN % 2.0000 5.0000 1.0500 0 NaN NaN NaN % Branch = % 支路类型 支路编号 起始节点 终止节点 支路电阻 支路电抗 % 0.0000 1.0000 3.0000 1.0000 0.0400 0.2500 % 0.0000 2.0000 3.0000 2.0000 0.1000 0.3500 % 0.0000 3.0000 1.0000 2.0000 0.0800 0.3000 % 1.0000 4.0000 1.0000 4.0000 0 0.0150 % 1.0000 5.0000 2.0000 5.0000 0 0.0300 % Ratio= % 变压器支路编号 起始节点 终止节点 变压器变比(终节点电压/始节点电压) % 4 2 4 0.95 % 5 3 5 1.05 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %代码要求支路上的起始节点编号必须小于终止节点编号 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% [Node,Branch,Ratio]=Per_unit(S_base,U_base); Node=sortrows(Node,1); Branch=sortrows(Branch,1); [node_number,~]=size(Node); [branch_number,~]=size(Branch); AdmittanceMatrix=zeros(node_number,node_number); Admittance_of_earth=zeros(node_number,node_number); Chain_index=zeros(node_number,node_number); for ii=1:node_number for jj=1:node_number for kk=1:branch_number if (Node(ii,2)==Branch(kk,3))&&(Node(jj,2)==Branch(kk,4)) Chain_index(ii,jj)=kk; Chain_index(jj,ii)=Chain_index(ii,jj); AdmittanceMatrix(ii,jj)=-(Branch(Chain_index(ii,jj),5)+Branch(Chain_index(ii,jj),6).*1i).^(-1); AdmittanceMatrix(jj,ii)=-(Branch(Chain_index(ii,jj),5)+Branch(Chain_index(ii,jj),6).*1i).^(-1); end end end end for ii=1:node_number for jj=1:node_number if Chain_index(ii,jj)==0 continue; end if Branch(Chain_index(ii,jj),1)~=1 AdmittanceMatrix(ii,ii)=AdmittanceMatrix(ii,ii)+(Branch(Chain_index(ii,jj),5)+Branch(Chain_index(ii,jj),6).*1i).^(-1); elseif ii<jj kk=find(Ratio(:,1)==Branch(Chain_index(ii,jj),2)); AdmittanceMatrix(ii,ii)=AdmittanceMatrix(ii,ii)+(Ratio(kk,4).*(Branch(Chain_index(ii,jj),5)+Branch(Chain_index(ii,jj),6).*1i)).^(-1)+(Ratio(kk,4)-1)./(Ratio(kk,4).*(Branch(Chain_index(ii,jj),5)+Branch(Chain_index(ii,jj),6).*1i)); Admittance_of_earth(ii,jj)=(Ratio(kk,4)-1)./(Ratio(kk,4).*(Branch(Chain_index(ii,jj),5)+Branch(Chain_index(ii,jj),6).*1i)); else ll=find(Ratio(:,1)==Branch(Chain_index(ii,jj),2)); AdmittanceMatrix(ii,ii)=AdmittanceMatrix(ii,ii)+(Ratio(ll,4).*(Branch(Chain_index(ii,jj),5)+Branch(Chain_index(ii,jj),6).*1i)).^(-1)+(Ratio(ll,4)-1)./((Ratio(ll,4).^2).*(Branch(Chain_index(ii,jj),5)+Branch(Chain_index(ii,jj),6).*1i)); Admittance_of_earth(ii,jj)=(Ratio(ll,4)-1)./((Ratio(ll,4).^2).*(Branch(Chain_index(ii,jj),5)+Branch(Chain_index(ii,jj),6).*1i)); end end end % Admittance_of_earth=zeros(5,5); % AdmittanceMatrix=[1.378742-6.291665i -0.624025+3.900156i -0.754717+2.641509i 0 0 % -0.624025+3.900156i 1.453901-66.980821i -0.829876+3.112033i 63.492063i 0 % -0.754717+2.641509i -0.829876+3.112033i 1.584593-35.737858i 0 31.746032i % 0 63.492063i 0 -66.666667i 0 % 0 0 31.746032i 0 -33.333333i];
评论
    相关推荐