Определить токи в схеме, если известны по условию: матрица задающих токов – J, матрица

Определить токи в схеме, если известны по условию: матрица задающих токов – J, матрица коэффициентов распределения для дерева схемы – G0, блоки матрицы инциденции - Nα, Nβ, блоки матрицы сопротивлений ветвей - Zαα, Zββ. Матрицу токов следует вычислять по формуле: I=IαIβ=E-NαT∙Yk∙Nα∙Zαα-Yk∙Nα∙Zαα∙C0∙J где Yk – матрица контурных проводимостей, Yk=Zk-1, Zk – матрица контурных сопротивлений, Zk=N·ZB·NT, ZB – матрица сопротивлений ветвей, ZB=ZααZαβZβαZββ, Zαβ=Zβα=0, N – матрица инциденций, N=(Nα Nβ), Е - единичная матрица третьего порядка. J=123, C0=-1-100-1000-1, Nα=-1-11-1-10, Nβ=1001 Zαα=100020003, Zββ=4005

Запишем матрицу сопротивления ветвей:
ZB=1002000000000000300040005
Запишем матрицу инциденций:
N=-1-1110-1-1001
Транспонируем матрицу инциденций:
NT=-1-11-1-101001
Вычислим матрицу контурных сопротивлений:
ZK=N∙ZB∙NT=-1-1110-1-1001∙1002000000000000300040005∙-1-11-1-101001=-1-2340-1-2005∙-1-11-1-101001=10338
Её определитель ∆=80-9=71
Матрица контурных проводимостей: Yk=Zk-1=171∙10338=1071371371871
-Yk∙Nα∙Zαα=-171∙10338∙-1-11-1-10∙100020003=-171∙-10-3-10-310+0-3-8-3-83+0∙100020003=-171-13-1310-11-113∙100020003=-171-13+0+00-26+00+0+30-11+0+00-22+00+0+9=-171-13-2630-11-229=13712671-307111712271-971
E-NαT∙Yk∙Nα∙Zαα= =100010001--1-1-1-110∙171∙10338∙-1-11-1-10∙100020003=100010001-171-10-3-3-8-10-3-3-810+03+0∙-1-11-1-10∙100020003=100010001-171-13-11-13-11103∙-1-11-1-10∙100020003=100010001-17113+1113+11-13-013+1113+11-13+0-10-3-10-310+0∙100020003=100010001-1712424-132424-13-13-1310∙100020003=100010001-17124+0+00+48+00+0-3924+0+00+48+00+0-39-13+0+00-26+00+0+30=100010001-1712448-392448-39-13-2630=100010001-24714871-397124714871-3971-1371-26713071=4771-48713971-247123713971137126714171=17147-4839-242339132641
Матрица токов:
I=IαIβ=E-NαT∙Yk∙Nα∙Zαα-Yk∙Nα∙Zαα∙C0∙J==4771-48713971-24712371397113712671417113712671-307111712271-971∙-1-100-1000-1∙123=17147-4839-2423391359411326-301122-9∙-1-100-1000-1∙123=171-47+0+0-47+48+00+0-3924-0+024-23+00+0-39-13+0+0-13-59+00+0-41-13+0+0-13-26+00+0+30-11+0+0-11-22+00+0+9∙123=171-471-39241-39-13-72-41-13-3930-11-339∙123=171-47+2-11724+2-117-13-144-123-13-78+90-11-66+27=171-162-91-280-1-50=-16271-9171-28071-171-5071=-2,282-1,282-3,944-0,014-0,704
Блоки матрицы токов: Iα-2,282-1,282-3,944, Iβ=-0,014-0,704