Метод преобразования многолучевой звезды в полный многоугольник (Лабораторная работа № 4), страница 9

C              массиве записаны номер узла, к которому сворачивается схема

C   KW – количество линейно независимых узлов схемы минус один

C   KH – количество ветвей схемы

C Выходные переменные и массивы подпрограммы

C   KNNE - количество ветвей с проводимостями Yij, которые  необходимо запомнить

C          при выполнении прямого хода метода Гаусса для того, чтобы выполнить

C          расчёт узловых напряжений при обратном ходе (колькасць

C          всех ненулевых элементов в верхней треугольной матрице коэфC          фициентов узловых уравнений после окончания прямого хода)

C   KHMPZWW(KW) - массив количества ветвей многолучевых звёзд удаляемых

C               узлов, которые были подключены к удаляемому узлу в момент

C               его выделения (используется при подсчёте напряжений)

C   YSMPZWW(KW) - массив сумм проводимостей ветвей многолучевых звёзд уда

C               ляемага узла, подключенных к удаляемогу узлу в моC               мент их удаления (диагональный коэффициент системы узлоC               вых уравнений для уравнения, которое соответствует номеру удаляC               емого узла)

C   JMPZWW(KW) - массив токов источников тока многолучевых звёзд, которые "наC               гнетали" ток в удаляемый узел в момент его удаления

C               (правая часть системы того узлового уравнения, которое

C              соответствует номеру удаляемого узла)

C   NWMPZWW(KNNE) - массив номеров вершин многоучевых звёзд удаляеC                   мых узлов, на которые опирается звезда в момент удаC                   ления этого узла (номера столбцов ненулевых элеменC                   тов в строках верхней треугольной матрицы коэффициенC                   тов системы узловых уравнений, которые соответствуют ноC                   меру удаляемого узла)

C   YHMPZWW(KNNE) - массив проводимостей ветвей многолучевх звёзд

C                  удаляемых узлов (ненулевые элементы тех строк верхней

C                   треугольной матрицы коэффициентов системы узловых

C  уравнений, соответствующие номеру удаляемого узла)

C Рабочие массивы в подпрограмме ZHORTKA

C KHDNW (KW) - Количество ветвей, присоединенных к неудаленным узлам, начиная

C от баллансирующего узла и заканчивая последним узлом в списке NW (KW)

C NWMPZ (KHMPZ) - Массив номеров узлов многолучевой звезды

C YHMPZ (KHMPZ) - массив проводимостей ветвей многолучевой звезды

C JHMPZ (KHMPZ) - массив тока источников тока в области многолучевой звезды                                  REAL JY,JMPZWW,JHMPZ

DIMENSION N1K(100),N2K(100),Y(100),JY(100),NW(70),

*KHMPZWW(70),YSMPZWW(70),JMPZWW(70),NWMPZWW(300),

*YHMPZWW(300),KHDNW(70),NWA(70),NWMPZ(30),YHMPZ(30),

*JHMPZ(30)

C Блок 101: Подготовка ко входу в основной цикл     

IF(NDRUK.EQ.0) GO TO 800

WRITE(3,*)' БЛОК 101',' KW=',KW,' NW(KW+1)',(NW(I),I=1,KW+1)

WRITE(3,*)' KH=',KH,' N1K(KH),N2K(KH),Y(KH),JY(KH)=',

*(N1K(I),N2K(I),Y(I),JY(I),I=1,KH)

800 KHB=KH

KNNE=0

DO 1 I=1,KW

JMPZWW(I)=0.

1 KHMPZWW(I)=0

C Начало цикла по удалению узлов из схемы; последовательно удаляются

C в порядке возрастания количества присоединенных к ним ветвей все узлы схемы,

C кроме последнего узла

DO 32 IUZ = 1, KW

C Блок 102: Определение количества ветвей, которые присоединены к неудаляемым узлам

C начиная с баллансирующего (удаляемого) узла и заканчивая последним узлом

C в массиве NW (KW)         

IF(NDRUK.EQ.0) GO TO 820

WRITE(3,*)' БЛОК 102',' IUZ=',IUZ,' KHB=',KHB

820 DO 7 J=IUZ,KW

KHDNW(J)=0

M=NW(J)

DO 2 K=1,KHB

IF(N1K(K).EQ.M.OR.N2K(K).EQ.M)KHDNW(J)=KHDNW(J)+1

2 CONTINUE

7 CONTINUE

IF(NDRUK.EQ.0) GO TO 830

WRITE(3,*)' БЛОК 102',' KW=',KW,' KHDNW(KW)=',(KHDNW(I),I=IUZ,KW)

C Блок 103: Определение наибольшего количества ветвей KHMAX, которые присоединены к                                                   C однму из неудаленных узлов

830 KHMAX=0

DO 3 J=IUZ,KW

IF(KHMAX.LT.KHDNW(J))KHMAX=KHDNW(J)

3 CONTINUE

IF(NDRUK.EQ.0) GO TO 840

WRITE(3,*)' БЛОК 103',' KHMAX=',KHMAX

C Блок 104: Перастановка неудаленных номеров узлов в массиве NW (KW) таким