Неточности вычислений

@you_fail_me · Регистрация: 16.02.2012

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

Ситуация такова, что начал замечать что в некоторых расчетах, в которых должен получаться аккуратно ноль, вместо нуля возникает малое но ненулевое значение порядка 10^-13—10^-16.
С чем это связано, это мой косяк или особенности пайтоновских вычислений?Как с этим бороться?(Занимаюсь такого рода расчетами что такие величины могут играть роль).
Для примера привожу код, в котором заметил такое:

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
def gauss1(array):
    for k in range(0, len(array)-1):
        for i in range(k+1, len(array)):
            lam = array[i][k]/array[k][k]
            for j in range(k, len(array[i])):
                array[i][j] = array[i][j]-lam*array[k][j]
    for k in range(len(array)-1, 0, -1):
        for i in range(k-1, -1, -1):
            if array[k][k]!=0:
                lam = array[i][k]/array[k][k]
                for j in range(len(array[i])):
                    array[i][j] = array[i][j]-lam*array[k][j]
rerurn array

Данная функция приводит матрицу к диагональному виду, вызвав ее для матрицы:
15 -3 8 12
2 5 -2 21
-5 16 2 2
Получаем результат вида:

Python
1
[[15.0, -4.440892098500626e-16, 8.881784197001252e-16, 45.64606741573033], [0.0, 5.4, 0.0, 8.261797752808986], [0.0, 0.0, 13.185185185185183, -47.888888888888886]]

Как видим, вместо нулей на месте второго и третьего членов первой строки имеем эти самые малые ненулевые значения.

Python 2.7.2, если что.

@Nameless One · 05.04.2012, 15:45

you_fail_me, это не проблема питона, это особенность представления чисел с плавающей точкой в компьютере (вкратце — произвольное вещественное число невозможно абсолютно точно отобразить в данном представлении). Так что можешь смело считать такие числа за ноль.

@you_fail_me · 05.04.2012, 15:55 **[ТС]**

А можно какуюто ссылку на подробное изложение данного нюанса? Каков «нулевой порог»?

@Nameless One · 05.04.2012, 15:56

Можно начать отсюда: http://ru.wikipedia.org/wiki/Ч... ей_запятой

@you_fail_me · 07.04.2012, 21:37 **[ТС]**

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

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
def gauss1(array):
    for k in range(0, len(array)-1):
        for i in range(k+1, len(array)):
            lam = array[i][k]/array[k][k]
            for j in range(k, len(array[i])):
                array[i][j] = array[i][j]-lam*array[k][j]
    for k in range(len(array)-1, 0, -1):
        for i in range(k-1, -1, -1):
            if array[k][k]!=0:
                lam = array[i][k]/array[k][k]
                for j in range(len(array[i])):
                    array[i][j] = array[i][j]-lam*array[k][j]
    for i in range(len(array)):
        for j in range(len(array[i])):
            array[i][j] = array[i][j]/array[i][i]
    return array

Но результат получился не таким, каким его хотелось бы видеть — почему-то не изменяет свободные члены, и после обработки этой функцией такой же исходной матрицы, какая приводилась в начале темы получается результат:

Python
1
[[1.0, -4.440892098500626e-16, 0.0, 4.851017441860465], [0.0, 1.0, 0.0, -1.2675872093023255], [-0.0, -0.0, 1.0, 4.716049382716049]]

Подскажите, пожалуйста, чем это вызвано, что я делаю не так. Буду очень признателен.

Добавлено через 18 минут
Понимаю, что можно сделать проще и чтоб работало, например, так:

Python
1
sol = list((array[i][-1]/array[i][i]) for i in range(len(array)))

но интересно, почему тот вариант не срабатывает.

Новые блоги и статьи Все статьи Все блоги /
Модель ЗдрввоСохранения 7: больше работников, больше ресурсов. anaschu 08.04.2026 работников и заданий может быть сколько угодно, но настроено всё так, что используется пока что только 20% kYBz3eJf3jQ	Дальние перспективы сервера - слоя сети с космологическим дизайном интефейса карты и логики. Hrethgir 07.04.2026 Дальнейшее ближайшее планирование вывело к размышлениям над дальними перспективами. И вот тут может быть даже будут нужны оценки специалистов, так как в дальних перспективах всё может очень сильно. . .	Горе от ума kumehtar 07.04.2026 Эта мне ментальная установка, что вот прямо сейчас, мол, мне для полного счастья не хватает (нужное вписать), и когда я этого достигну - тогда и полный кайф. Одна из самых сильных ловушек на пути. . . .	Использование значений реквизитов справочника в документе, с определенными условиями и правами Maks 07.04.2026 1. Контроль срока действия договора Алгоритм из решения ниже реализован на примере нетипового документа "ЗаявкаНаРаботу", разработанного в конфигурации КА2. Задача: уведомлять пользователя, если. . .
Доступность команды формы по условию Maks 07.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: сделать доступной кнопку (команда формы "ЗавершитьСписание") при. . .	Уведомление о неверно выбранном значении справочника Maks 06.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "НарядПутевка", разработанного в конфигурации КА2. Задача: уведомлять пользователя, если в документе выбран неверный склад. . .	Установка 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 существует уже давно, и также давно существуют скрипты под нее. Тем не менее, прога живет, периодически что-то не спеша дополняется, улучшается. Что меня в первую очередь. . .

@Nameless One 5828 / 3479 / 358 Регистрация: 08.02.2010 Сообщений: 7,448
	05.04.2012, 15:45
	you_fail_me, это не проблема питона, это особенность представления чисел с плавающей точкой в компьютере (вкратце — произвольное вещественное число невозможно абсолютно точно отобразить в данном представлении). Так что можешь смело считать такие числа за ноль. 1

@you_fail_me 0 / 0 / 0 Регистрация: 16.02.2012 Сообщений: 51
	05.04.2012, 15:55 [ТС]
	А можно какуюто ссылку на подробное изложение данного нюанса? Каков «нулевой порог»? 0

@Nameless One 5828 / 3479 / 358 Регистрация: 08.02.2010 Сообщений: 7,448
	05.04.2012, 15:56
	Можно начать отсюда: http://ru.wikipedia.org/wiki/Ч... ей_запятой 0