Машинная арифметика в рациональных числах

Содержание

Слайд 2

Устранение причины возможной ошибки Задача. Составить программу решения квадратного уравнения a2

Устранение причины возможной ошибки

Задача.
Составить программу решения квадратного уравнения
a2 + bx +

c = 0

Формулы корней квадратного уравнения достаточно сложны по числу и набору действий:

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

Слайд 3

Из анализа формул корней, видно, что в одной из формул в

Из анализа формул корней, видно, что в одной из формул в

числителе обязательно будет разность чисел.
Пусть b>0, тогда при ac<

Здесь можно полностью устранить причину возможной ошибки за счет тождественного преобразования для х2

Получили другую формулу для одного из корней квадратного уравнения

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

Слайд 4

Если b Очевидно, что новые формулы для корней квадратного уравнения пригодны

Если b<0, то разность близких чисел возможна в формуле для х1.

Соответствующими преобразованиями можно и ее привести к виду:

Очевидно, что новые формулы для корней квадратного уравнения пригодны и для «нормальных случаев», когда разность близких чисел не возникает. Т.о. проводить анализ близости соответствующих чисел нет необходимости. Достаточно, определив знак b, выбрать для вычисления каждого корня ту формулу, в которой разности чисел нет.

Слайд 5

Примеры вычислительно неустойчивых алгоритмов Рассмотрим задачу вычисления функции ex . Известно,

Примеры вычислительно неустойчивых алгоритмов

Рассмотрим задачу вычисления функции ex . Известно, что

эта задача хорошо обусловлена.

при x<0

Пример.
Найти значение функции ex при x= -15.
Верное значение e-15 =1 / e15 ≈ 0.000000305902

1. Традиционные вычисления
После выполнения 82 итераций было получено: e-15 ≈ 0.000000256502
Относительная погрешность
составила 19,2%.

2. Решение
e-x = 1/ ex

Слайд 6

Применение высокоточных вычислений для численного решения уравнения теплопроводности с разномасштабными коэффициентами

Применение высокоточных вычислений для численного решения уравнения теплопроводности с разномасштабными коэффициентами

Нестационарное

уравнение теплопроводности:

.

Уравнение теплопроводности описывает распространение тепла вдоль стенок труб
теплообменных устройств, поверхностей нагрева котлов, и других объектов.

- коэффициент температуропроводности

При

максимальная относительная погрешность решения

системы методом Гаусса-Зейделя составила около 100% при числе итераций порядка 5000.
Физическая интерпретация этих параметров может быть такой: ищется распределение температуры магистрального трубопровода длиной 1 км через каждые 1 метр, каждые 3 часа.
Применение высокоточных вычислений позволило получить решение системы с требуемой точностью

Система конечно-разностных линейных уравнений:

Слайд 7

Подход к решению проблемы высокоточных вычислений на основе модулярной арифметики К

Подход к решению проблемы высокоточных вычислений на основе модулярной арифметики

К

настоящему времени модулярная арифметика использовалась как средство повышения быстродействия в криптографии, нейронных сетях, цифровой обработке сигналов и др.
Проведенные исследования показали качественно новые возможности применения модулярной арифметики в повышении точности вычислений и ослаблении зависимости времени вычислений от точности, для некоторых частных задач:
решение дифференциальных уравнений методами Рунге-Кутта,
нахождение скалярного произведения векторов,
решения систем линейных уравнений методами Гаусса-Зейделя,
релаксации,
дискретном преобразовании Фурье .
Слайд 8

Схема вычислений с исключением ошибок округления Рост числителей и знаменателей!

Схема вычислений с исключением ошибок округления

Рост числителей и знаменателей!

Слайд 9

ВОЗМОЖНЫЕ ПРИЛОЖЕНИЯ ВЫЧИСЛЕНИЙ С ИСКЛЮЧЕНИЕМ ОШИБОК ОКРУГЛЕНИЯ 1. Точное вычисление обобщенных

ВОЗМОЖНЫЕ ПРИЛОЖЕНИЯ ВЫЧИСЛЕНИЙ С ИСКЛЮЧЕНИЕМ ОШИБОК ОКРУГЛЕНИЯ

1. Точное вычисление обобщенных обратных

матриц. Например, таких как, g-обратная матрица Мура-Пенроуза. Многие алгоритмы требуют умения распознавать значение численного ранга, а это является трудной задачей при наличии ошибок округления.
2. Целочисленное решение систем линейных уравнений. Примером могут служить построение оптимальных решений в задачах целочисленного программирования.
3. Точное вычисление характеристического многочлена матрицы. Вследствие ошибок округления будут получены приближенные значения коэффициентов. Если многочлен плохо обусловлен, то корни "приближенного» характеристического уравнения могут быть плохими приближениями к корням истинного уравнения.
4. Обращение матриц Гильберта, Адамара и др. особо чувствительных к ошибкам округления.
5. Для решения промежуточных между классами корректных и некорректных задач.
Класс задач, изменяющих корректность при решении. Это расчет устойчивости систем управления, выч. собств. знач. систем лин. одн. урав. и др.
Слайд 10

Схема вычислений с исключением ошибок округления

Схема вычислений с исключением ошибок округления

Слайд 11

Высокоточные вычисления с иррациональными числами Модулярная система счисления позволяет представлять: -

Высокоточные вычисления с иррациональными числами

Модулярная система счисления позволяет представлять:
- целые числа,

рациональные числа
комплексные числа с рациональной мнимой и вещественной частью

Иррациональные числа представлять в виде цепных дробей.

Точное значение 10/7 это [1,2,3].
Приближенное значение [1,2]

Вычисления с иррациональными числами с контролируемой точностью.

Слайд 12

Способы нахождения обратного по модулю Подбор С помощью малой теоремы Ферма С помощью расширенного алгоритма Евклида

Способы нахождения обратного по модулю

Подбор
С помощью малой теоремы Ферма
С помощью расширенного

алгоритма Евклида
Слайд 13

Слайд 14

Алгоритм Евклида

Алгоритм Евклида

Слайд 15

Слайд 16

Слайд 17

Алгоритм Евклида

Алгоритм Евклида

Слайд 18

Слайд 19

Слайд 20

Алгоритм нахождения обратного

Алгоритм нахождения обратного

Слайд 21

Алгоритм обратного преобразования чисел в дроби Фарея

Алгоритм обратного преобразования чисел в дроби Фарея

Слайд 22

Алгоритм обратного преобразования чисел в дроби Фарея

Алгоритм обратного преобразования чисел в дроби Фарея

Слайд 23

Преимущество модулярной арифметики - отсутствие переносов при сложении, вычитании и умножении

Преимущество модулярной арифметики - отсутствие переносов при сложении, вычитании и умножении чисел -

единое представление чисел различной природы (целых, комплексных)

Задачи на модулярную арифметику Найти последние 2 цифры числа 4^30.

Слайд 24

Табличная реализация модулярная арифметики Таблица умножения по модулю m = 11

Табличная реализация модулярная арифметики

Таблица умножения по модулю m = 11

Слайд 25

Китайская теорема об остатках и обратное преобразование Китайская теорема об остатках

Китайская теорема об остатках и обратное преобразование

Китайская теорема об остатках –

существует только одно число, имеющее остатки по модулям в диапазоне до произведения модулей минус один