潮流计算求解导纳矩阵的通用程序

给定任何网络拓扑结构的导纳矩阵都能被求解、、包括带变压器

clear;clc;
format compact
% n=input（‘请输入节点数：n=‘）;
% pr=input（‘请输入精度误差：pr=‘）;
% Z=input（‘请输入由各节点间阻抗形成的距阵：Z=‘）;
% B=input（‘请输入由各节点直接对地的导纳距阵：B=‘）;
% K=input（‘请输入由节点间变压器变比值：K=‘）;
% N=input（‘请输入当变压器两端节点对地等效导纳值确定条件：N=‘）;
%数据输入环节%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n=5;
pr=0.00001;
Z=[0 0.04+0.25i 0.1+0.35i 0 0;
   0.04+0.25i 0 0.08+0.30i 0.015i 0;
   0.1+0.35i 0.08+0.30i 0 0 0.03i;
   0 0.015i 0 0 0;
   0 0 0.03i 0 0];
B=[0.25i 0.5i 0.25i 0 0];
K=[1 1 1 1 1;
   1 1 1 1.05 1;
   1 1 1 1 1.05;
   1 1.05 1 1 1;
   1 1 1.05 1 1;];
N=[0 1 1 0 0];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Y=zeros（n）;C=zeros（n）%（经变压器变换后对地的等效导纳）;
A=zeros（n）;D=zeros（1n）%（经变压器变换后总的对地导纳）;
H=zeros（1n）%（经变压器变换后节点间互导纳）;
for i=1:n
    for j=1:n
        if Z（ij）~=0
         A（ij）=1/Z（ij）
        end
        if Z（ij）==0
          A（ij）=0
        end
    end
end
for i=1:n
    for j=1:n
        H（ij）=A（ij）/K（ij）
        if K（ij）==1
           C（ij）=0
        end
        if K（ij）~=1
           if N（i）==1
              C（ij）=（1-K（ij））*H（ij）/K（ij）
           end
            if N（i）==0
                C（ij）=（K（ij）-1）*H（ij）
           end
        end
        Q=C（i:） %（取C矩阵的某一行元素）
        D（i）=sum（sum（Q））+B（i）
    end
end
for i=1:n
    for j=1:n
        if i~=j
            Y（ij）=-H（ij）
        end
        if i==j
            E=H（i:） %（取H矩阵的某一行元素）
            Y（ij）=sum（sum（E））+D（i）
        end
    end
end       %求导纳距阵
disp（‘导距Y为:‘）;
disp（Y）

            
 
       

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1715  2010-04-08 23:22  daonaijuzhentongyongchengxu.m

----------- ---------  ---------- -----  ----

                 1715                    1