Я прекрасно понял, то что ты хотел сказать, но это неподходит, т.к:
1. Мой элемент размещен на диагонали случайно и почти не бывает последним;
2. Вся модификация в твоем примере будет сводиться к установке точки выхода из Ж-Г при расчете требуемого элемента, ведь если мой элемент является первым, то обращение будет проведено полностью, что совсем не то, что мне хотелась!
Дело в том, что данную проблему я решаю применительно к расчету сопротивлений электрической схемы:
диагональная матрица сопротивлений ветвей - Z;
матрица содинений схемы замещения - М;
матрица узловых проводимостей - Yu=M*1/Z*MT;
! матрица узловых сопротивлений - Zu=Revers(Yu);//азъ есим - обращение матрицы
Но мне нужно посчитать собственное сопротивление только одного узла, т.е. Zu[x,x].
Сначала я думал, что Zu[i,i]=1/Yu[i,i], но это утверждение работало для некоторых матриц, а для некоторых нет.
Из твоих слов я понял, что это справедливо только для последнего элемента в треугольной матрице, получаемой при прямом ходе в методе Гаусса. Тогда может необходимо делать какие-нибудь перестановки перед этим, что бы мой элемент оказывался в конце, например в самой матрице соединений М. А вот этого лучше не делать!
Поэтому я и ищу модификацию именного самого алгоритма обращения. |
, но отдельно было добавлено, что