电力系统短路程序

各种短路程序，有vb还有c语言还有matlab的程序

% NS=input（‘请输入短路点的数目：NS=‘）;
    NS=1;
for nsn=1:NS
%     n=input（‘请输入独立节点数：n=‘）;
%     l=input（‘请输入支路数：l=‘）;
%     B=input（‘请输入支路参数矩阵：B=‘）;
%     V0=input（‘请输入由各节点的初始电压标幺值形成的列阵：V0=‘）;
%     D=input（‘请输入由短路号、短路点阻抗组成的行阵：D=‘）;

n=4;
l=5;
B=[0 1 0.2i 0 1 0;0 2 4i 0 1 0;1 3 0.51i 0 1 0;2 3 0.59i 0 1 0;3 4 1.43i 0 1 0];
V0=[1;1;1;1];
D=[4;0];
    m=0;
    Z=zeros（n）;%阻抗矩阵
    V=zeros（n1）;%电压矩阵
    I=zeros（l1）;%电流矩阵
    for kk=1:l
        p=B（kk1）;
        q=B（kk2）;
        if B（kk6）==0
            k=1./B（kk5）;
        else
            k=B（kk5）;
        end
        if p==0
            if q>m
                Z（qq）=B（kk3）;
                m=m+1;
            else
                for i=1:m
                    Z（im+1）=-Z（iq）;
                    Z（m+1i）=-Z（qi）;
                end
                Z（m+1m+1）=Z（qq）+B（kk3）;
                for i=1:m
                    for j=1:m
                        Z（ij）=Z（ij）-Z（im+1）*Z（m+1j）./Z（m+1m+1）;
                    end
                    Z（im+1）=0;
                end
                for i=1:m+1
                    Z（m+1i）=0;
                end
            end
        else
            if q>m
                for i=1:m
                    Z（iq）=Z（ip）*k;
                    Z（qi）=Z（pi）*k;
                end
                Z（qq）=k^2*Z（pp）+k^2*B（kk3）;
                m=m+1;
            else
                for i=1:m
                    Z（im+1）=k*Z（ip）-Z（iq）;
                    Z（m+1i）=k*Z（pi）-Z（qi）;
                end
                Z（m+1m+1）=k^2*Z（pp）+Z（qq）-2*k*Z（pq）+k^2*B（kk3）;
                for i=1:m
                    for j=1:m
                        Z（ij）=Z（ij）-Z（im+1）*Z（m+1j）./Z（m+1m+1）;
                    end
                    Z（im+1）=0;
                end
                for i=1:m+1
                    Z（m+1i）=0;
                end
            end
        end
    end
end
disp（‘阻抗矩阵Z=‘）;
disp（Z）;
I（D（11）1）=V0（D（11）1）./（Z（D（11）D（11））+D（21））;
disp（‘短路点电流标幺值If=‘）;
disp（I（D（11）1））;
for i=1:n
    V（i1）=V0（i1）-I（D（11）1）*Z（iD（11））;
end
for i=1:l
    if B（i6）==0
        k=B（i5）;
    else
        k=1./（B（i5））;
    end
    p=B（i1）;
    q=B（i2）;
    if p==0
        e=0;
        b=B（i3）;
        I（i1）=（e-V（q1）./k）./b;
    else
        I（i1）=（V（p1）-V（q1）./k）./B（i3）;
    end
end
disp（‘各节点电压标幺值V=‘）;
disp（V）;
disp（‘各支路电流标幺值I=‘）;
disp（I）;
  

        
        
        
        
        
        
        
        
        
        
        
        
        
        

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2866  2012-07-03 20:33  短路程序\BASIC短路电流计算小程序.txt
     文件        2749  2012-07-08 21:01  短路程序\duanlu.m
     文件       11990  2012-07-08 21:01  短路程序\gezhongduanlujisuan11.m
     文件        1356  2012-07-03 20:23  短路程序\shortcircuit.m
     文件       22223  2012-07-03 20:29  短路程序\短路计算程序.cpp
     目录           0  2013-04-18 16:48  短路程序\