节点电压法-自编的迷你spice资料整理

cir=[['R01',1,0,1.5e3],

     ['R02',7,0,300],

     ['R03',8,2,15e3],

     ['R04',8,9,15e3],

     ['R05',10,9,18e3],

     ['R06',10,11,18e3],

     ['R07',11,0,39e3],

     ['R11',1,2,25e3],

     ['R12',2,7,25e3],

     ['R21',1,3,25e3],

     ['R22',3,7,25e3],

     ['R31',1,4,25e3],

     ['R32',4,7,25e3],

     ['R41',1,5,25e3],

     ['R42',5,7,25e3],

     ['R51',1,6,25e3],

     ['R52',6,7,25e3],

     ['C01',3,8,33e-9],

     ['C02',4,9,10e-9],

     ['C03',5,10,4.7e-9],

     ['C04',6,11,1e-9],

     ['C05',8,0,4.7e-9],

     ['C06',9,0,6.8e-9],

     ['C07',10,0,2.2e-9],

     ['C08',11,0,1e-9],

     ['I01',1,0,1]]

##获得两个相同节点之间的所有元件，计算自电导用

def get_y1(cir,n1,n2):

    y=[]

    for e in cir:

        if(n1==e[1] or n2==e[2]):

            y.append(e)

    return y



##获得两个不同节点之间的所有元件，计算互电导用

def get_y2(cir,n1,n2):

    y=[]

    a=max(n1,n2)

    b=min(n1,n2)

    for e in cir:

        c=max(e[1],e[2])

        d=min(e[1],e[2])

        if(a==c and b==d):

            y.append(e)

    return y



##建立Y元件导纳矩阵子函数

def get_y(cir,n1,n2):

    if(n1==n2):

        return get_y1(cir,n1,n2)

    else:

        return get_y2(cir,n1,n2)



##建立Y元件导纳矩阵

def get_y_matrix(cir):

    Y=[]

    for i in range(1,12):

        YL=[]

        for j in range(1,12,1):

            YL.append(get_y(cir,i,j))

        Y.append(YL)

    return Y

##从元件获得复导纳值子函数，暂时只支持RLC，转化成数值导纳矩阵用

def get_y_value_from_element(e,f):

    s=1j*2*math.pi*f

    if(e[0].startswith('R')):

        y=1/e[3]

    elif(e[0].startswith('C')):

        y=(s*e[3])

    elif(e[0].startswith('L')):

        y=1/(s*e[3])

    else:

        y=0

    ##print(e[0],e[3],y)

    return y



##从一组元件求总导纳值，转化成数值导纳矩阵用   

def get_y_value_from_elements(elements,f):

    Y=0

    for e in elements:

        Y=Y+get_y_value_from_element(e,f)

    return Y



##从元件导纳矩阵获得数值导纳矩阵   

def get_y_value_matrix(Y,f):

    YV=[]

    for i in range(1,12):

        YVL=[]

        for j in range(1,12,1):

            #自电导为正，邱关源《电路》P69

            if(i==j):

                v=get_y_value_from_elements(Y[i-1][j-1],f)

            #互电导为负

            else:

                v=-1*get_y_value_from_elements(Y[i-1][j-1],f)

            YVL.append(v)

        YV.append(YVL)

    return YV

##king5555的5段均衡器数值传递函数

def get_hjw(cir,f):

    Y=get_y_matrix(cir)

    display_y_syms(Y)

    YM=numpy.mat(get_y_value_matrix(Y,f))

    I=numpy.mat([[1/(1.5e3)],

                 [0],

                 [0],

                 [0],

                 [0],

                 [0],

                 [0],

                 [0],

                 [0],

                 [0],

                 [0]])

    ##print(YM)

    V=YM.I*I

    return V

[[0.86937316]

 [0.40459092]

 [0.43641139]

 [0.44254197]

 [0.44322551]

 [0.44375548]

 [0.01829224]

 [0.35758804]

 [0.30048658]

 [0.22958955]

 [0.15727944]]

1 1 ['R01', 'R11', 'R21', 'R31', 'R41', 'R51', 'I01']

1 2 ['R11']

1 3 ['R21']

1 4 ['R31']

1 5 ['R41']

1 6 ['R51']

1 7 []

1 8 []

1 9 []

1 10 []

1 11 []

2 1 ['R11']

2 2 ['R03', 'R11', 'R12']

2 3 []

2 4 []

2 5 []

2 6 []

2 7 ['R12']

2 8 ['R03']

2 9 []

2 10 []

2 11 []

3 1 ['R21']

3 2 []

3 3 ['R21', 'R22', 'C01']

3 4 []

3 5 []

3 6 []

3 7 ['R22']

3 8 ['C01']

3 9 []

3 10 []

3 11 []

4 1 ['R31']

4 2 []

4 3 []

4 4 ['R31', 'R32', 'C02']

4 5 []

4 6 []

4 7 ['R32']

4 8 []

4 9 ['C02']

4 10 []

4 11 []

5 1 ['R41']

5 2 []

5 3 []

5 4 []

5 5 ['R41', 'R42', 'C03']

5 6 []

5 7 ['R42']

5 8 []

5 9 []

5 10 ['C03']

5 11 []

6 1 ['R51']

6 2 []

6 3 []

6 4 []

6 5 []

6 6 ['R51', 'R52', 'C04']

6 7 ['R52']

6 8 []

6 9 []

6 10 []

6 11 ['C04']

7 1 []

7 2 ['R12']

7 3 ['R22']

7 4 ['R32']

7 5 ['R42']

7 6 ['R52']

7 7 ['R02', 'R12', 'R22', 'R32', 'R42', 'R52']

7 8 []

7 9 []

7 10 []

7 11 []

8 1 []

8 2 ['R03']

8 3 ['C01']

8 4 []

8 5 []

8 6 []

8 7 []

8 8 ['R03', 'R04', 'C01', 'C05']

8 9 ['R04']

8 10 []

8 11 []

9 1 []

9 2 []

9 3 []

9 4 ['C02']

9 5 []

9 6 []

9 7 []

9 8 ['R04']

9 9 ['R04', 'R05', 'C02', 'C06']

9 10 ['R05']

9 11 []

10 1 []

10 2 []

10 3 []

10 4 []

10 5 ['C03']

10 6 []

10 7 []

10 8 []

10 9 ['R05']

10 10 ['R05', 'R06', 'C03', 'C07']

10 11 ['R06']

11 1 []

11 2 []

11 3 []

11 4 []

11 5 []

11 6 ['C04']

11 7 []

11 8 []

11 9 []

11 10 ['R06']

11 11 ['R06', 'R07', 'C04', 'C08']