Система линейных уравнений. Решение методом Гаусса

@Krajben · Регистрация: 24.05.2014

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

Доброго времени суток.
Сегодня столкнулся с проблемой. Лазил на форуме, находил решения, но все равно не верно получается... Вполне возможно упустил что.

В чем суть:
Имеется код, все считает отлично, но если в элемент (1,1) подставить значение 0, то все летит к чертям. Учусь в строительном, поэтому с программированием плоховато дружу.

Visual Basic
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
Dim a(5, 5), b(5) As Double
        Dim x(5), f, k As Double
        Label1.Text = " "
        Label2.Text = " "
        Label3.Text = " "
'Вывод матрицы на лэйбл
        For i = 1 To 5
            For j = 1 To 5
                Label1.Text = Label1.Text & a(i, j) & " "
            Next
            Label1.Text = Label1.Text & vbCrLf
        Next
'Вычисление элементов
        For u = 1 To 5
            For i = u + 1 To 5
                If a(i, u) = 0 Then
                    i = i
                Else
                    k = -a(u, u) / a(i, u)
                    If k = 0 Then
                        For j = u To 5
                            a(i, j) = a(i, j) * 0 + a(u, j)
                        Next j
                    Else
                        For j = u To 5
                            a(i, j) = a(i, j) * k + a(u, j)
                        Next j
                        b(i) = b(i) * k + b(u)
                    End If
                End If
            Next i
        Next u
'Обратный ход, нахождение Х
        For i = 5 To 1 Step -1
            f = 0
            For j = i + 1 To 5
                f = f + a(i, j) * x(j)
            Next
            x(i) = (b(i) - f) / a(i, i)
        Next
'Вывод результата
        For i = 1 To 5
            For j = 1 To 5
                Label2.Text = Label2.Text & a(i, j) & " "
            Next
            Label2.Text = Label2.Text & vbCrLf
        Next
        For i = 1 To 5
            label3.text = label3.text & x(i) & vbCrLf
        Next

Вот и не могу сообразить что не так... До элемента (4,3) считает все правильно, а дальше переплет в дебаге вижу...

Проверяю по данному варианту:

Система линейных уравнений. Решение методом Гаусса

Надеюсь, изложил достаточно информации и надеюсь на Вашу помощь

@Антихакер32 · 24.05.2014, 10:25

«Решение систем линейных уравнений. Метод Гаусса»

@Krajben · 24.05.2014, 10:33 **[ТС]**

FelixMacintosh, благодарен! Попробую разобраться, на такое не натыкался.

@Антихакер32 · 24.05.2014, 10:52

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

Надеюсь, изложил достаточно информации и надеюсь на Вашу помощь

во всяком случае, там более подробнее, и даже если не разберёшся,
то найдуться охотники разобрать это...

Добавлено через 2 минуты
я бы и сдесь выложил, переписав в своём хакерском стиле..
только меня бесит движок этого сайта, и здешний редактор формул, вобщем беда...

@Krajben · 24.05.2014, 11:02 **[ТС]**

FelixMacintosh, в любом случае большое спасибо! Я посмотрел - есть непонятные моменты, разберусь, это не сложно

Добавлено через 3 минуты
FelixMacintosh, мне вообще дали наводку на то, чтобы просветить коэф. k на равенство 0... Видимо, ошибся где-то и не могу догнать.

@Krajben · 04.06.2014, 18:02 **[ТС]**

UPD
Для тех, кто любит готовенькое...
Добил этот код до конца, прошу использовать грамотно. Выкладываю полностью программу.

Visual Basic
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
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
 Dim a(5, 5), b(10) As Double
 
        Dim x(5), f, k, swap As Double
 
        a(1, 1) = Val(TextBox1.Text)
        a(2, 1) = Val(TextBox2.Text)
        a(3, 1) = Val(TextBox3.Text)
        a(4, 1) = Val(TextBox4.Text)
        a(5, 1) = Val(TextBox5.Text)
        a(1, 2) = Val(TextBox10.Text)
        a(2, 2) = Val(TextBox9.Text)
        a(3, 2) = Val(TextBox8.Text)
        a(4, 2) = Val(TextBox7.Text)
        a(5, 2) = Val(TextBox6.Text)
        a(1, 3) = Val(TextBox15.Text)
        a(2, 3) = Val(TextBox14.Text)
        a(3, 3) = Val(TextBox13.Text)
        a(4, 3) = Val(TextBox12.Text)
        a(5, 3) = Val(TextBox11.Text)
        a(1, 4) = Val(TextBox20.Text)
        a(2, 4) = Val(TextBox19.Text)
        a(3, 4) = Val(TextBox18.Text)
        a(4, 4) = Val(TextBox17.Text)
        a(5, 4) = Val(TextBox16.Text)
        a(1, 5) = Val(TextBox25.Text)
        a(2, 5) = Val(TextBox24.Text)
        a(3, 5) = Val(TextBox23.Text)
        a(4, 5) = Val(TextBox22.Text)
        a(5, 5) = Val(TextBox21.Text)
        b(1) = Val(TextBox30.Text)
        b(2) = Val(TextBox29.Text)
        b(3) = Val(TextBox28.Text)
        b(4) = Val(TextBox27.Text)
        b(5) = Val(TextBox26.Text)
 
        Label1.Text = " "
        Label2.Text = " "
        Label3.Text = " "
 
        For i = 1 To 5
            For j = 1 To 5
                Label1.Text = Label1.Text & a(i, j) & " "
            Next
            Label1.Text = Label1.Text & vbCrLf
        Next
 
 
        For u = 1 To 5
            For i = u + 1 To 5
                If a(u, u) = 0 Then 'Если элемент главной диагонали равен 0, то нужно поменять строки местами
                    For j = u To 5
                        swap = a(i, j)
                        a(i, j) = a(u, j)
                        a(u, j) = swap
                    Next
                    swap = b(i)
                    b(i) = b(u)
                    b(u) = swap
                    If a(i, u) = 0 Then
                        i = i
                    Else
                        k = -a(u, u) / a(i, u)
                        For j = u To 5
                            a(i, j) = a(i, j) * k + a(u, j)
                        Next j
                        b(i) = b(i) * k + b(u)
                    End If
                Else
                    If a(i, u) = 0 Then
                        i = i
                    Else
                        k = -a(u, u) / a(i, u)
                        For j = u To 5
                            a(i, j) = a(i, j) * k + a(u, j)
                        Next j
                        b(i) = b(i) * k + b(u)
                    End If
                End If
            Next i
        Next u
        'Обратный ход, нахождение Х
        For i = 5 To 1 Step -1
            f = 0
            For j = i + 1 To 5
                f = f + a(i, j) * x(j)
            Next
            x(i) = (b(i) - f) / a(i, i)
        Next
 
        For i = 1 To 5
            For j = 1 To 5
                Label2.Text = Label2.Text & a(i, j) & " "
            Next
            Label2.Text = Label2.Text & vbCrLf
        Next
 
        For i = 1 To 5
            Label3.Text = Label3.Text & x(i) & vbCrLf
        Next

SoftIce · 04.06.2014, 21:27

А если матрица будет 100х100 Вы 10000 строк напишите?

А(99,99)=val(textbox10000.text)

@Антихакер32 · 05.06.2014, 09:45

10000 строк писать вручную !!!

Добавлено через 6 минут

Visual Basic
1
2
3
        For f = 1 To 5
            a(f, 1) = CallByName(Me, "TextBox" & f, VbGet)
        Next

Добавлено через 7 минут

Visual Basic
1
2
3
        For f = 1 To 5
            a(f, 1) = CallByName(Me, "TextBox" & f, VbGet).Text
        Next

Новые блоги и статьи Все статьи Все блоги /
Валидация и контроль данных табличной части документа перед записью 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, для личных нужд. Собственно программа предоставляет из себя базу данных клиентов, в которой можно фиксировать звонки, стадии сделки, а также. . .

@Антихакер32 1201 / 473 / 46 Регистрация: 06.01.2014 Сообщений: 1,797 Записей в блоге: 19
	24.05.2014, 10:25
	«Решение систем линейных уравнений. Метод Гаусса» 0

@Krajben 0 / 0 / 1 Регистрация: 24.05.2014 Сообщений: 4
	24.05.2014, 10:33 [ТС]
	FelixMacintosh, благодарен! Попробую разобраться, на такое не натыкался. 0

@Krajben 0 / 0 / 1 Регистрация: 24.05.2014 Сообщений: 4
	24.05.2014, 11:02 [ТС]
	FelixMacintosh, в любом случае большое спасибо! Я посмотрел - есть непонятные моменты, разберусь, это не сложно Добавлено через 3 минуты FelixMacintosh, мне вообще дали наводку на то, чтобы просветить коэф. k на равенство 0... Видимо, ошибся где-то и не могу догнать. 0

SoftIce es geht mir gut 11274 / 4760 / 1183 Регистрация: 27.07.2011 Сообщений: 11,439
	04.06.2014, 21:27
	А если матрица будет 100х100 Вы 10000 строк напишите? А(99,99)=val(textbox10000.text) 0

Опции темы