Две последние программы при использовании ПМК с входным языком ЯМК52 можно несколько упростить, но для л>-5 уравнений системы для ее решения приходится использовать полуавтоматический режим вычислений с записью промежуточных результатов в вычислительный бланк.

Программа 85. Решение системы из л12 линейных уравнений по методу Жордана с использованием вычислительного бланка

1 Сх ИП1 -i- ПО 1 ПС .-

С/П ИПД ИПС 1 + ПС - х>0 26 <-

кипе ИПО X - БП 10 -н ИПО БП 02

Инструкция. Для вычисления элементов р-й части бланка (табл. 3.2), в нулевой части которого записаны столбцы коэффициентов, столбец свободных членов qi = Wi,n+i и контрольный столбец, элементы которого равны суммам элементов соответствующих строк, выполнить: л = РД, =Р1. ms" =

> = Рл (л < 12); w\P~l\ = РХ В/О С/П РХ = 0; поочередно

= Р2, ,,,

вводить остальные элементы (р+1)-го и следующих столбцов (р-1)-й части бланка в регистр РХ, нажимая после каждого ввода только клавишу С/П и

Таблица 3.2. Решение системы уравнений методом Жордана с циклической перестановкой строк вычислительного бланка

регистрируя результат на одну строку того же столбца выше, чем в предыдущей части бланка записан введенный элемент (после ввода первого элемента следующего столбца результат записывается в последней строке предыдущего столбца); после ввода последнего элемента контрольного столбца и выполнения программы записать результат в предпоследней строке контрольного столбца, еще раз нажать клавишу С/П и записать результат в последней строке контрольного столбца; после выполнения п шагов с заполнением всех п частей бланка в п-м столбце п-й части бланка будут записаны значения искомых решений (корней) системы уравнения xt; вычисляемые значения элементов контрольного столбца должны быть равными сумме элементов соответствующих строк бланка, включая р единиц, записываемых в р нижних строках р-го столбца р-й части бланка - нарушение этого равенства свидетельствует об ошибке в вычислениях.

Для проверки правильности ввода и выполнения программы, а также порядка заполнения бланка можно воспользоваться данными табл. 3.2. С целью уменьшения операционных погрешностей целесообразно иеред началом решения уравнения системы переставить так, чтобы на главной диагонали матрицы коэффициентов находились наибольшие по абсолютному значению числа.

Прн определении функций безынерционных цепей по выражениям, подобным формулам (3.4) - (3.7), требуется вычислить определители, включая алгебраические дополнения, которые могут лишь знаком отличаться от миноров, равных определителям матриц с вычеркнутыми строками н столбцами.

Программа 86. Вычисление определителя Д матрицы второго порядка

-.- X П8 X ИП8 - С/П

Инструкция. ш„ = РТ, tiy,j=PZ, a)2, = PY, Ш22 = РХ В/О С/П РХ=Д. Пример

1 2 3 4

= -2.

Программа 87. Вычисление определителя Д матрицы третьего порядка

П7 П8 «• ИП7 -4 П9 2 ПО 7 П1 С/П ПА ИП8 X - КП1 ИПА ИП9

X - КП1 L0 И ИПб X ИП4 ИПб Х - ИП7 X С/П

Инструкция. syi, = PZ, tt),2=PY, ш,з = РХ В/О С/П tti2i = PZ, Ш22 = РУ, «2з=РХ СП (03i = PZ, ft)32=PY, Шзз = РХ СП РХ=Д. Пример.

3 4 1

4 10 6 1 2 3

= 28.

Программа 88. Вычисление определителя Д матрицы порядка

Инструкция. гг = РО, ш„ = РХ В/О С/П ш,2 = РХ С/П ... ш,п = РХ С/П а)2, = РХ С/П Ш22 = РХ С/П ... а)2п = РХ С/П ... а)„, = РХ С/П ш„2=РХ С/П... а)„„ = РХ С/П РХ=А.

10 4 5 6

= 40(<=к23 с);

5 6 -1 4

3 3

= 69,000001 (<~ 72 с);

= -60 (t~ 140 с);

= 360(/sk260 с).

Для вычисления определителей матриц порядка иЗ можно исполыовать полуавтоматический режим с записью промежуточных результатов в вычислительный бланк. В этом случае целесообразно использовать схему единственного деления Гаусса [14], так как достаточно выполнить вычисления лишь на ее прямом ходе с относительно небольшим числом строк вычислительного бланка.

Программа 89. Вычисление определителя А матрицы порядка иЮ по схеме единственного деления Гаусса с использованием вычислительного бланка

ПВ ИПС X ПС ИПД 1 - ПД ПО ИПС С/П ИПВ ~ КПО ИПО 1 - х=0 10 С/П ПА ИПД ПО С/П КИПО ИПА у - ИПО 1 - х=0 36 БП 19 -.- БП 23

Инструкция. Записать в нулевую часть бланка столбцы коэффициентов матрицы и в контрольный столбец суммы элементов каждой строки матрицы (табл. 3.3); гг-(-2=РД; 1=РС; для каждой р-й части бланка aj/fjl*, =РХ В/О С/П РХ=Д, остальные элементы первого столбца (р-1)-й части ввести