Метод Гаусса прямым ходом

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

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

Проблема с получением нижней треугольной матрицы. С верхней треугольной матрицей вроде нет проблем. Вот что получается в ответе:
U
1.0 5e-05 3.3333333333333335e-05 2.5e-05 2e-05 1.6666666666666667e-05
0.0 1.0 3.333000016665e-05 2.49975001249875e-05 1.9998000099990004e-05 1.6665000083325e-05
0.0 0.0 1.0 2.499625037495417e-05 1.9997000299963337e-05 1.6664166916636114e-05
0.0 0.0 0.0 1.0 1.9996333824926536e-05 1.6663611520772113e-05
0.0 0.0 0.0 0.0 1.0 1.6663194999904353e-05
0.0 0.0 0.0 0.0 0.0 1.0

L
9.99999985501611 0.0 0.0 0.0 2.710302164310409e-20 0.0
0.0019998100132212946 19.99999990500661 0.0 0.0 5.4208753466420714e-20 0.0
0.002999815007499694 0.001499907503749847 29.999999938335833 0.0 -2.032804537729561e-24 0.0
0.003999853335999951 0.0019999266679999756 0.001333284445333317 39.999999963334 1.0842021724855044e-19 0.0
0.004999916666666666 0.002499958333333333 0.001666638888888889 0.0012499791666666666 49.99999998333333 0.0
0.006 0.003 0.002 0.0015 0.0012 60.
код:

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
import numpy as np
n=6
r=1
p=2
q=1
A=np.zeros((n,n))
LU=np.zeros((n,n))
#заполнение массива А
for i in range(n): 
    for j in range(n): 
        if i==j:
            A[i,j]=10*(i+1)**(p/2)
        else:
            A[i,j]=10**(-3)*((i+1)/((j+1)**r))**(1/q)   
B=np.zeros((n))
 #заполнение массива В
for i in range(n):
    B[i]=9*i**(p/2)
 
for i in range(len(A)): 
    for j in range(len(A[i])):
        print(A[i,j], end = ' ')
    print()
D=np.copy(A)
# прямой ход
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
for i in range(len(U)): 
    for j in range(len(U[i])):
        print(U[i,j], end = ' ')
    print()
# матрица L
for i in range(5,0,-1):
    for j in range (i-1,-1,-1):
        D[j,:]-=D[i,:]*D[j,i]/D[i,i]
L=D
for i in range(len(L)): 
    for j in range(len(L[i])):
        print(L[i,j], end = ' ')
    print()

@Gdez · 29.11.2020, 11:36

Gromjk, в чем именно проблема?

Добавлено через 1 час 11 минут
Посмотри тут, может поможет
https://orion1401.gitbooks.io/... henie.html

Новые блоги и статьи Все статьи Все блоги /
Валидация и контроль данных табличной части документа перед записью Maks 22.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа, разработанного в КА2. Задача: контроль и валидация данных табличной части документа перед записью с учетом регламента компании. . .	Отчёт о затраченных материалах за определенный период с макетом печатной формы Maks 21.04.2026 Отчёт из решения ниже размещён в конфигурации КА2. Задача: разработка отчёта по затраченным материалам за определённый период, с возможностью вывода печатной формы отчёта с шапкой и подвалом. В. . .	Отчёт о спецтехнике находящейся в ремонте Maks 20.04.2026 Отчёт из решения ниже размещен в конфигурации КА2. Задача: отобразить спецтехнику, которая на данный момент находится в ремонте. Есть нетиповой документ "Заявка на ремонт спецтехники" который. . .	Памятка для бота и "визитка" для читателей "Semantic Universe Layer (Слой семантической вселенной)" Hrethgir 19.04.2026 Сгенерировано для краткого описания по случаю сборки и компиляции скелета серверного приложения. И пусть после этого скажут, что статьи сгенерированные AI - туфта и не интересно. И это не реклама -. . .
Запрет удаления строк ТЧ документа при определённом условии Maks 19.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "Аккумуляторы", разработанного в конфигурации КА2. У данного документа есть ТЧ, в которой в зависимости от прав доступа. . .	Модель заражения группы наркоманов alhaos 17.04.2026 Условия задачи сформулированы тут Суть: - Группа наркоманов из 10 человек. - Только один инфицирован ВИЧ. - Колются одной иглой. - Колются раз в день. - Колются последовательно через. . .	Мысли в слух. Про "навсегда". kumehtar 16.04.2026 Подумалось тут, что наверное очень глупо использовать во всяких своих установках понятие "навсегда". Это очень сильное понятие, и я только начинаю понимать край его смысла, не смотря на то что давно. . .	My Business CRM MaGz GoLd 16.04.2026 Всем привет, недавно возникла потребность создать CRM, для личных нужд. Собственно программа предоставляет из себя базу данных клиентов, в которой можно фиксировать звонки, стадии сделки, а также. . .

@Gdez 8851 / 4502 / 1864 Регистрация: 27.03.2020 Сообщений: 7,318
	29.11.2020, 11:36
	Gromjk, в чем именно проблема? Добавлено через 1 час 11 минут Посмотри тут, может поможет https://orion1401.gitbooks.io/... henie.html 0

Опции темы