Оптимизировать алгоритм Гаусса для разреженной матрицы

25.05.2023, 23:01. Показов 616. Ответов 0

Метки нет (Все метки)

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

вместо звездочек число = 4
кроме диагонали главной, на ней: np.random.rand()*20+4

нужно оптимизировать алгоритм гаусса LU разложения без перестановки строк, чтобы учитывалось как можно меньше нулевых элементов. Нельзя использовать никакие библиотеки кроме numpy.

вот код для заполнения матрицы:

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
import numpy as np
 
def stupenki_sleva(start,finish,matrica):
    i=0
    for j in range(finish, start,-1):
        matrica[i,:j] = 4
        i+=1
    return matrica
 
matrica =np.eye((20), dtype=float)*4
matrica[8:12,:] = 4
matrica = stupenki_sleva(0,6,matrica)
i=0
for j in range(14, 20):
        matrica[i,j::] = 4
        i+=1
i=14
for j in range(1, 7):
        matrica[i,:j] = 4
        i+=1
i=14
for j in range(19,13,-1):
        matrica[i,j::] = 4
        i+=1
matrica[:,9:11]=4
for k in range(20):
    x=np.random.rand()*20+24
    matrica[k][k]=x
for i in range(20):
    print('\n')
    for j in range(20):
        print(matrica[i][j], end=' ')

вот структура матрицы на рисунке:

Миниатюры

Оптимизировать алгоритм Гаусса для разреженной матрицы

IT_Exp

Эксперт

34794 / 4073 / 2104

Регистрация: 17.06.2006

Сообщений: 32,602

Блог

25.05.2023, 23:01

Ответы с готовыми решениями:

Для заданной сильно разреженной упакованной матрицы А найти след матрицы АА'
Для заданной сильно разреженной упакованной матрицы А(m,n) найти след матрицы АА' (следом назовем сумму элементов главной диагонали).С...

Как перегрузить оператор [][] для разреженной матрицы?
как перегрузить оператор для разреженной матрицы?знаю, что такого оператора нет, есть только , но все равно не понимаю как его перегрузить

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

Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.

BasicMan

Эксперт

29316 / 5623 / 2384

Регистрация: 17.02.2009

Сообщений: 30,364

Блог

25.05.2023, 23:01

Помогаю со студенческими работами здесь

Написать алгоритм решения матрицы методом Гаусса, чтобы он подходил к любому размеру матрицы
Всем привет, Нужно в маткаде написать алгоритм решения матрицы методом гауса, чтобы он подходил к любому размеру матрицы(10х10,30х30,20х10...

Массив: Модифицировать класс разреженной матрицы для работы с тремя измерениями – x,y,z.
Задание : модифицировать класс разреженной матрицы (проект SparseMatrix) для работы с тремя измерениями – x,y,z. Вывод элементов в методе...

Реализовать алгоритм нахождения обратной матрицы методом Гаусса
Прошу помощи, так как сам уже даже и не знаю. Нужно реализовать алгоритм нахождения обратной матрицы методом Гаусса Вот мой класс: ...

Решение разреженной матрици методом Гаусса единственного деления и методом простой итерации
Помогите написать программу решения разреженной матрицы методом Гаусса единственного деления и методом простой итерации. И посчитать...

Создание разреженной матрицы
Здравствуйте, помогите с заданием Нужно создать разреженную матрицу с целыми числами, где все нулевые элементы размещены выше главной...

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:

Новые блоги и статьи Все статьи Все блоги /
Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .	SDL3 для Web (WebAssembly): Сборка библиотек SDL3 и Box2D из исходников с помощью CMake и Emscripten 8Observer8 27.02.2026 Недавно вышла версия SDL 3. 4. 2 библиотеки SDL3. На странице официальной релиза доступны исходники, готовые DLL (для x86, x64, arm64), а также библиотеки для разработки под Android, MinGW и Visual. . .	SDL3 для Web (WebAssembly): Реализация движения на Box2D v3 - трение и коллизии с повёрнутыми стенами 8Observer8 20.02.2026 Содержание блога Box2D позволяет легко создать главного героя, который не проходит сквозь стены и перемещается с заданным трением о препятствия, которые можно располагать под углом, как верхнее. . .	Конвертировать закладки radiotray-ng в m3u-плейлист damix 19.02.2026 Это можно сделать скриптом для PowerShell. Использование . \СonvertRadiotrayToM3U. ps1 <path_to_bookmarks. json> Рядом с файлом bookmarks. json появится файл bookmarks. m3u с результатом. # Check if. . .
Семь CDC на одном интерфейсе: 5 U[S]ARTов, 1 CAN и 1 SSI Eddy_Em 18.02.2026 Постепенно допиливаю свою "многоинтерфейсную плату". Выглядит вот так: https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11617&stc=1&d=1771445347 Основана на STM32F303RBT6. На борту пять. . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .	И ясному Солнцу zbw 12.02.2026 И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.	«Знание-Сила» zbw 12.02.2026 «Знание-Сила» «Время-Деньги» «Деньги -Пуля»

Наверх