LU - разложение матрицы

@Gromjk · Регистрация: 22.12.2019

Студворк — интернет-сервис помощи студентам

Проблема с получением матрицы L, не совсем понятно как прийти к занулению коэффициентов над главной диагональю

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
A=[[10.0,0.0005,0.0003333333333333333,0.00025,0.0002],[0.002,20.0,0.0006666666666666666,0.0005,0.0004],[0.003,0.0015,30.0,0.00075,0.0006],[0.004,0.002,0.0013333333333333333,40.0,0.0008],[0.005,0.0025,0.0016666666666666668,0.00125,50.0],[0.006,0.003,0.002,0.0015,0.0012,60.0]]
D=np.copy(A)
# прямой ход (получение верхней треугольной матрицы U)
for i in range(5):
    A[i,:]/=A[i,i]
    for j in range (i+1,6):
        A[j,:]+=-A[i,:]*A[j,i]/A[i,i]
A[5,:]/=A[5,5]
U=A
#получение нижней треугольной матрицы L
for i in range(5,0,-1):
    for j in range (i-1,-1,-1):
        D[j,:]+=-D[i,:]*D[j,i]
L=D

@dondublon · 02.11.2020, 11:52

Что конкретно не получается?
Такие вещи хорошо на numpy делать.

@Gromjk · 02.11.2020, 11:59 **[ТС]**

numpy это хорошо, просто задание - написать разложение без него. Я не совсем понимаю какой надо взять цикл, для того чтобы он считал значения выше главной диагонали. Он у меня берет каждое значение и производит с ним операцию

@dondublon · 02.11.2020, 12:05

А у вас там уже numpy

Фомула есть?

Новые блоги и статьи Все статьи Все блоги /
Модель ЗдрввоСохранения 7: больше работников, больше ресурсов. anaschu 08.04.2026 работников и заданий может быть сколько угодно, но настроено всё так, что используется пока что только 20% kYBz3eJf3jQ	Дальние перспективы сервера - слоя сети с космологическим дизайном интефейса карты и логики. Hrethgir 07.04.2026 Дальнейшее ближайшее планирование вывело к размышлениям над дальними перспективами. И вот тут может быть даже будут нужны оценки специалистов, так как в дальних перспективах всё может очень сильно. . .	Горе от ума kumehtar 07.04.2026 Эта мне ментальная установка, что вот прямо сейчас, мол, мне для полного счастья не хватает (нужное вписать), и когда я этого достигну - тогда и полный кайф. Одна из самых сильных ловушек на пути. . . .	Использование значений реквизитов справочника в документе, с определенными условиями и правами Maks 07.04.2026 1. Контроль срока действия договора Алгоритм из решения ниже реализован на примере нетипового документа "ЗаявкаНаРаботу", разработанного в конфигурации КА2. Задача: уведомлять пользователя, если. . .
Доступность команды формы по условию Maks 07.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: сделать доступной кнопку (команда формы "ЗавершитьСписание") при. . .	Уведомление о неверно выбранном значении справочника Maks 06.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "НарядПутевка", разработанного в конфигурации КА2. Задача: уведомлять пользователя, если в документе выбран неверный склад. . .	Установка Qt Creator для C и C++: ставим среду, CMake и MinGW без фреймворка Qt 8Observer8 05.04.2026 Среду разработки Qt Creator можно установить без фреймворка Qt. Есть отдельный репозиторий для этой среды: https:/ / github. com/ qt-creator/ qt-creator, где можно скачать установщик, на вкладке Releases:. . .	AkelPad-скрипты, структуры, и немного лирики.. testuser2 05.04.2026 Такая программа, как AkelPad существует уже давно, и также давно существуют скрипты под нее. Тем не менее, прога живет, периодически что-то не спеша дополняется, улучшается. Что меня в первую очередь. . .

@dondublon 4653 / 2073 / 366 Регистрация: 17.03.2012 Сообщений: 10,183 Записей в блоге: 6
	02.11.2020, 11:52
	Что конкретно не получается? Такие вещи хорошо на numpy делать. 0

@Gromjk 0 / 0 / 0 Регистрация: 22.12.2019 Сообщений: 10
	02.11.2020, 11:59 [ТС]
	numpy это хорошо, просто задание - написать разложение без него. Я не совсем понимаю какой надо взять цикл, для того чтобы он считал значения выше главной диагонали. Он у меня берет каждое значение и производит с ним операцию 0

@dondublon 4653 / 2073 / 366 Регистрация: 17.03.2012 Сообщений: 10,183 Записей в блоге: 6
	02.11.2020, 12:05
	А у вас там уже numpy Фомула есть? 0