Метод Ньютона выводит nan

@JewishBoy · Регистрация: 08.04.2022

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

Здравствуйте, делаю курсовую работу. Суть задания написать код для метода Ньютона и с помощью программы расчитать силу тока. Код скорее всего правильный, так как ввожу данные другого варианта, и ответы получаются +- одинаковыми, погрешность пару процентов. Что я выяснил, так это то, что на 4 круг цикла значение h становится настолько маленьким(-5^28) что в следущий раз оно превращается в inf и далее элементы делятся на этот inf. Как итог, в ответе только nan.

Кликните здесь для просмотра всего текста

https://sun9-80.userapi.com/im... type=album

Вот сам код, все исходные данные сразу ввёл, чтобы по 100 раз не вводить тоже самое.

Кликните здесь для просмотра всего текста

Code
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
#include <math.h>
#include <iostream>
using namespace std;
int main()
{
    int i, j, N, k, m, s, r;
    float A[5][5], B[5], X[5], E1, E2, E5, R1, R2, R3, R4, R5, Ri1, Ri2, Ri5, e, f[5], ltx, h, bb[5], aa[5][5];
    E1 = 22; E2 = 12; E5 = 15; R1 = 12; R2 = 20; R3 = 8; R4 = 5; R5 = 10; Ri1 = 0.3; Ri2 = 0.1; Ri5 = 0.1; N = 5;
    for (i = 0; i < N; i++)
    {
        B[i] = 0;
        for (j = 0; j < N; j++)
            A[i][j] = 0;
    }
 
    A[0][0] = 0;          A[0][1] = 0;          A[0][2] = 1; A[0][3] = -1; A[0][4] = -1;      B[0] = 0;
    A[1][0] = -1;         A[1][1] = 1;          A[1][2] = -1; A[1][3] = 0;  A[1][4] = 0;       B[1] = 0;
    A[2][0] = -(R1 + Ri1); A[2][1] = -(R2 + Ri2); A[2][2] = 0; A[2][3] = 0;  A[2][4] = 0;       B[2] = -E1 - E2;
    A[3][0] = 0;          A[3][1] = R2 + Ri2;   A[3][2] = R3; A[3][3] = R4; A[3][4] = 0;       B[3] = E2;
    A[4][0] = 0;          A[4][1] = 0;          A[4][2] = 0; A[4][3] = -R4; A[4][4] = R5 + Ri5; B[4] = E5;
 
    m = 100;
    e = 0.0001;
    s = 0;
    for (i = 0; i < N; i++)
    {
        X[i] = 0.1;
    }
Nachalo:    for (i = 0; i < N; i++)
{
    f[i] = 0;
    for (j = 0; j < N; j++)
    {
        f[i] += A[i][j] * X[j];
    }
    f[i] = f[i] - B[i];
}                //МАТРИЦА BB
for (i = 0; i < N; i++)
{
    bb[i] = -f[i];
}
for (j = 0; j < N; j++)
{
    ltx = X[j];
    h = e * fabs(ltx);
    X[j] = ltx + h;
    f[j] = 0;
    for (i = 0; i < N; i++)
    {
        f[j] += A[i][j] * X[i];
    }
    f[j] = f[j] - B[j];
    for (i = 0; i < N; i++)
    {
        aa[i][j] = (f[i] + bb[i]) / h;
    }
    X[j] = ltx;
}
s = s + 1;
if (s == m + 1) goto Konec;
else
{
    for (i = 0; i < N - 1; i++)
    {
        for (j = i + 1; j < N; j++)
        {
            aa[j][i] = -aa[j][i] / aa[i][i];               
            for (k = i + 1; k < N; k++)
            {
                aa[j][k] = aa[j][k] + aa[j][i] * aa[i][k];
            }
            bb[j] = bb[j] + aa[j][i] * bb[i];
        }
    }
    f[N - 1] = bb[N - 1] / aa[N - 1][N - 1];
    for (i = N - 2; i > 0; i--)
    {
        h = bb[i];
        for (j = i + 1; j < N; j++)
        {
            h = h - f[i] * aa[i][j];
        }
        f[i] = h / aa[i][i];
    }
    r = 0;
    for (i = 0; i < N; i++)
    {
        X[i] = X[i] + f[i];
        if (fabs(f[i] / X[i]) > e) r = 1;
    }
    if (r == 1) goto Nachalo;
    else
    {
    Konec:for (i = 0; i < N; i++)
        cout << "X[" << i << "]=" << X[i] << endl;
    }
}
}

@PencilTorch · 08.04.2022, 22:56

Если немного причесать и убрать олдскульное goto, то начинает как-то работать:

C++
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
#include <iostream>
#include <math.h>
 
using namespace std;
 
int main() {
    const int N = 5;
    bool r = true;
    int m = 100, s = 0;
    double A[5][5]{ 0 }, B[5]{ 0 }, X[5], E1, E2, E5, R1, R2, R3, R4, R5, Ri1, Ri2, Ri5, e = 0.0001, f[5]{ 0 }, ltx, h, bb[5], aa[5][5];
    E1 = 22; E2 = 12; E5 = 15; R1 = 12; R2 = 20; R3 = 8; R4 = 5; R5 = 10; Ri1 = 0.3; Ri2 = 0.1; Ri5 = 0.1;
 
    A[0][0] = 0;          A[0][1] = 0;          A[0][2] = 1; A[0][3] = -1; A[0][4] = -1;      B[0] = 0;
    A[1][0] = -1;         A[1][1] = 1;          A[1][2] = -1; A[1][3] = 0;  A[1][4] = 0;       B[1] = 0;
    A[2][0] = -(R1 + Ri1); A[2][1] = -(R2 + Ri2); A[2][2] = 0; A[2][3] = 0;  A[2][4] = 0;       B[2] = -E1 - E2;
    A[3][0] = 0;          A[3][1] = R2 + Ri2;   A[3][2] = R3; A[3][3] = R4; A[3][4] = 0;       B[3] = E2;
    A[4][0] = 0;          A[4][1] = 0;          A[4][2] = 0; A[4][3] = -R4; A[4][4] = R5 + Ri5; B[4] = E5;
 
    for (int i = 0; i < N; i++) {
        X[i] = 0.1;
    }
    
    while (r) {
        for (int i = 0; i < N; i++) {
            f[i] = 0;
            for (int j = 0; j < N; j++)
                f[i] += A[i][j] * X[j];
            f[i] -= B[i];
        }
        //МАТРИЦА BB
        for (int i = 0; i < N; i++)
            bb[i] = -f[i];
        for (int j = 0; j < N; j++) {
            ltx = X[j];
            h = e * fabs(ltx);
            X[j] = ltx + h;
            f[j] = 0;
            for (int i = 0; i < N; i++)
                f[j] = f[j] + A[i][j] * X[i];
            f[j] -= B[j];
            for (int i = 0; i < N; i++)
                aa[i][j] = (f[i] + bb[i]) / h;
            X[j] = ltx;
        }
        s++;
        if (s == m + 1)
            break;
        else {
            for (int i = 0; i < N - 1; i++) {
                for (int j = i + 1; j < N; j++) {
                    aa[j][i] = -aa[j][i] / aa[i][i];
                    for (int k = i + 1; k < N; k++)
                        aa[j][k] = aa[j][k] + aa[j][i] * aa[i][k];
                    bb[j] = bb[j] + aa[j][i] * bb[i];
                }
            }
            f[N - 1] = bb[N - 1] / aa[N - 1][N - 1];
            for (int i = N - 2; i > 0; i--) {
                h = bb[i];
                for (int j = i + 1; j < N; j++)
                    h = h - f[i] * aa[i][j];
                f[i] = h / aa[i][i];
            }
            r = 0;
            for (int i = 0; i < N; i++) {
                X[i] = X[i] + f[i];
                if (fabs(f[i] / X[i]) <= e) {
                    r = false;
                    break;
                }
            }
        }
    }
    for (int i = 0; i < N; i++)
        cout << "X[" << i << "]=" << X[i] << endl;
 
return 0;
}

@JewishBoy · 08.04.2022, 23:13 **[ТС]**

Спасибо, хоть ответы не те, но может что-то получится если оттолкнуться от вашего кода

@PencilTorch · 09.04.2022, 09:05

JewishBoy, с утра глянул на свежую голову - там я похоже лишнее оставил на 64 строке ( r = 0

. Без неё получше.

Новые блоги и статьи Все статьи Все блоги /
PhpStorm 2025.3: WSL Terminal всегда стартует в ~ and_y87 14.12.2025 PhpStorm 2025. 3: WSL Terminal всегда стартует в ~ (home), игнорируя директорию проекта Симптом: После обновления до PhpStorm 2025. 3 встроенный терминал WSL открывается в домашней директории. . .	Как объединить две одинаковые БД Access с разными данными VikBal 11.12.2025 Помогите пожалуйста !! Как объединить 2 одинаковые БД Access с разными данными.	Новый ноутбук volvo 07.12.2025 Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro	Музыка, написанная Искусственным Интеллектом volvo 04.12.2025 Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .	От async/await к виртуальным потокам в Python IndentationError 23.11.2025 Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .
Поиск "дружественных имён" СОМ портов Argus19 22.11.2025 Поиск "дружественных имён" СОМ портов На странице: https:/ / norseev. ru/ 2018/ 01/ 04/ comportlist_windows/ нашёл схожую тему. Там приведён код на С++, который показывает только имена СОМ портов, типа,. . .	Сколько Государство потратило денег на меня, обеспечивая инсулином. Programma_Boinc 20.11.2025 Сколько Государство потратило денег на меня, обеспечивая инсулином. Вот решила сделать интересный приблизительный подсчет, сколько государство потратило на меня денег на покупку инсулинов. . . .	Ломающие изменения в C#.NStar Alpha Etyuhibosecyu 20.11.2025 Уже можно не только тестировать, но и пользоваться C#. NStar - писать оконные приложения, содержащие надписи, кнопки, текстовые поля и даже изображения, например, моя игра "Три в ряд" написана на этом. . .	Мысли в слух kumehtar 18.11.2025 Кстати, совсем недавно имел разговор на тему медитаций с людьми. И обнаружил, что они вообще не понимают что такое медитация и зачем она нужна. Самые базовые вещи. Для них это - когда просто люди. . .	Создание Single Page Application на фреймах krapotkin 16.11.2025 Статья исключительно для начинающих. Подходы оригинальностью не блещут. В век Веб все очень привыкли к дизайну Single-Page-Application . Быстренько разберем подход "на фреймах". Мы делаем одну. . .

@JewishBoy 0 / 0 / 0 Регистрация: 08.04.2022 Сообщений: 2
	08.04.2022, 23:13 [ТС]
	Спасибо, хоть ответы не те, но может что-то получится если оттолкнуться от вашего кода 0

@PencilTorch 90 / 60 / 32 Регистрация: 06.08.2020 Сообщений: 220
	09.04.2022, 09:05
	JewishBoy, с утра глянул на свежую голову - там я похоже лишнее оставил на 64 строке ( r = 0. Без неё получше. 1