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

@BoDRbIi · Регистрация: 21.07.2009

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

Пишу LU разложение матрицы
Почему то L и U равны, смотрел отладку L меняется с U! я чего то не замечаю, ткните плиз пальцем.

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
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
import random
 
def print_matrix(matrix):
   for i in matrix:
      print i
 
print "LU-decomposition"
 
size = int(raw_input("\nEnter size of matrix: "))
matrix = []
U = []
L = []
tmp_matrix = []
tmp_matrix2 = []
 
x,y = 0, 0
 
while x < int(size):
    y=0
    while y < int(size): 
        element = random.randint(0, 5)
        tmp_matrix.append(element)
        tmp_matrix2.append(None)
        y = y+1
   
    matrix.append(tmp_matrix)
    L.append(tmp_matrix2)
    U.append(tmp_matrix2)
 
    tmp_matrix = []
    tmp_matrix2 = []
 
    x = x+1
 
print_matrix(L)
print "\n\n"
print_matrix(U)
 
for i in range(size):
    for j in range(size):
        U[0][int(i)] = matrix[0][int(i)]
        L[i][0] = matrix[i][0] / U[0][0]
        sum = float(0)
        for k in range(i):
            sum += L[i][k] * U[k][j]
        U[i][j] = matrix[i][j] - sum
        if i > j:
            L[j][i] = 0
        else:
            sum = 0
            for k in range(i):
                sum += L[j][k] * U[k][j]
            L[j][i] = (matrix[j][i] - sum) / U[i][i]
 
print "\n"
print_matrix(matrix)
print "\n\n"
print_matrix(L)
print "\n\n"
print_matrix(U)
raw_input("\nPress any key to exit")

Добавлено через 54 секунды
Когда доходит до 43 строки, L и U одинаковые, такого быть не должно. магя какая-то

@freest · 19.12.2011, 01:02

учить матчасть надо!

Сообщение от BoDRbIi

L.append(tmp_matrix2)
U.append(tmp_matrix2)

тут передается адресс, а не копия
Следовательно L[len(L)] указывает тудаже куда и tmp_matrix2i i U[len(U)]

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
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
import random
 
def print_matrix(matrix):
   for i in matrix:
      print i
 
print "LU-decomposition"
 
size = int(raw_input("\nEnter size of matrix: "))
matrix = []
U = [] 
L = []
tmp_matrix = []
tmp_matrix2 = []
 
x,y = 0, 0
 
while x < int(size):
    y=0
    while y < int(size): 
        element = random.randint(0, 5)
        tmp_matrix.append(element)
        tmp_matrix2.append(None)
        y = y+1
   
    matrix.append(tmp_matrix)
    L.append(tmp_matrix2)
    U.append([i for i in tmp_matrix2])
    tmp_matrix = []
    tmp_matrix2 = []
 
    x = x+1
 
print_matrix(L)
print "\n\n"
print_matrix(U)
print "\n\n"
print_matrix(matrix)
 
for i in range(size):
    for j in range(size):
        U[0][i] = matrix[0][i]
        L[i][0] = matrix[i][0] / U[0][0]
        print_matrix(L)
        print "\n\n"
        print_matrix(U)
        print "\n\n"
        print_matrix(matrix)
        sum = float(0)
        for k in range(i):
            sum += L[i][k] * U[k][j]
        U[i][j] = matrix[i][j] - sum
        if i > j:
            L[j][i] = 0
        else:
            sum = 0
            for k in range(i):
                sum += L[j][k] * U[k][j]
            L[j][i] = (matrix[j][i] - sum) / U[i][i]

accept · 19.12.2011, 03:09

копия списка делается через lst[:]
если есть вложенные списки, то через copy.deepcopy(lst)

@BoDRbIi · 19.12.2011, 09:10 **[ТС]**

Спасибо

Новые блоги и статьи Все статьи Все блоги /
Установка 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 существует уже давно, и также давно существуют скрипты под нее. Тем не менее, прога живет, периодически что-то не спеша дополняется, улучшается. Что меня в первую очередь. . .	Отображение реквизитов в документе по условию и контроль их заполнения Maks 04.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеСпецтехники", разработанного в конфигурации КА2. Данный документ берёт данные из другого нетипового документа. . .	Фото всей Земли с борта корабля Orion миссии Artemis II kumehtar 04.04.2026 Это первое подобное фото сделанное человеком за 50 лет. Снимок называют новым вариантом легендарной фотографии «The Blue Marble» 1972 года, сделанной с борта корабля «Аполлон-17». Новое фото. . .
Вывод диалогового окна перед закрытием, если документ не проведён Maks 04.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: реализовать программный контроль на предмет проведения документа. . .	Программный контроль заполнения реквизитов табличной части документа Maks 02.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: 1. Реализовать контроль заполнения реквизита. . .	wmic не является внутренней или внешней командой Maks 02.04.2026 Решение: DISM / Online / Add-Capability / CapabilityName:WMIC~~~~ Отсюда: https:/ / winitpro. ru/ index. php/ 2025/ 02/ 14/ komanda-wmic-ne-naydena/	Программная установка даты и запрет ее изменения Maks 02.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: при создании документов установить период списания автоматически. . .

accept 4866 / 3287 / 468 Регистрация: 10.12.2008 Сообщений: 10,570
	19.12.2011, 03:09
	копия списка делается через lst[:] если есть вложенные списки, то через copy.deepcopy(lst) 0

@BoDRbIi 60 / 10 / 0 Регистрация: 21.07.2009 Сообщений: 248
	19.12.2011, 09:10 [ТС]
	Спасибо 0