Методы вычисления. Лекции 1-2
.pdfМашинная арифметика, стандарт IEEE-754 Стандарт IEEE-754
Точность представления числа в стандарте IEEE-754
В ЭВМ представимы лишь конечный набор рациональных чисел специального вида.
Эти числа образуют представимое множество вычислительной машины.
Для всех остальных чисел x возможно лишь их приближенное представление с ошибкой, которую принято называть ошибкой представления (или ошибкой округления).
Пример непредставимого числа 0:110 = 0:000110011001100110011 : : :.
Кафедра ТМ (СГАУ) |
Методы вычислений |
19 февраля 2012 г. |
41 / 50 |
Машинная арифметика, стандарт IEEE-754 Стандарт IEEE-754
Точность представления числа в стандарте IEEE-754
Абсолютная максимальная ошибка для числа в формате IEEE-754 равна в пределе половине шага чисел.
Шаг чисел удваивается с увеличением экспоненты двоичного числа на единицу.
Чем дальше от нуля, тем шире шаг чисел в формате IEEE754 по числовой оси.
Кафедра ТМ (СГАУ) |
Методы вычислений |
19 февраля 2012 г. |
42 / 50 |
Машинная арифметика, стандарт IEEE-754 Стандарт IEEE-754
Точность представления числа в стандарте IEEE-754
Предел максимальной абсолютной ошибки будет равен 1/2 шага числа:
single : A(x ) = 2Es 22 127=2 = 2(Es 150)
double : A(x ) = 2Es 51 1023=2 = 2(Es 1075)
Кафедра ТМ (СГАУ) |
Методы вычислений |
19 февраля 2012 г. |
43 / 50 |
Машинная арифметика, стандарт IEEE-754 Стандарт IEEE-754
Точность представления числа в стандарте IEEE-754
Относительная погрешность нормализованного числа
(x ) = |
2E 150 |
|
= |
1 |
|
2E 127 1 + |
M |
|
223 + M |
||
|
223 |
|
Относительная погрешность денормализованного числа
(x ) = |
2E 150 |
= |
1 |
|||
2E 126 |
M |
2M |
||||
|
|
|||||
|
223 |
|
|
Кафедра ТМ (СГАУ) |
Методы вычислений |
19 февраля 2012 г. |
44 / 50 |
Машинная арифметика, стандарт IEEE-754 Стандарт IEEE-754
Погрешность чисел одинарной точности
x |
a (x ) |
|
(x ), % |
|
2 149 1:401298 10 45 |
2 150 |
0:700649 10 45 |
50 |
|
2 148 2:802597 10 45 |
2 150 |
0:700649 10 45 |
25 |
|
1:0 |
2 23 |
1:192 10 7 |
11:9209 10 6 |
|
100 |
2 17 |
7:6294 10 6 |
7:62939 10 6 |
|
1:0 1010 |
210 1:024 103 |
10:24 10 6 |
Кафедра ТМ (СГАУ) |
Методы вычислений |
19 февраля 2012 г. |
45 / 50 |
Машинная арифметика, стандарт IEEE-754 Стандарт IEEE-754
Погрешность чисел двойной точности
x |
a (x ) |
|
(x ), % |
|
2 1074 4:940656 10 324 |
2 1075 |
2:470328 10 324 |
50 |
|
2 1073 9:881313 10 324 |
2 1075 |
2:470328 10 324 |
25 |
|
1:0 |
2 52 |
2:220446 10 16 |
11:9209 10 6 |
|
100 |
2 17 |
7:6294 10 6 |
2:220446 10 14 |
|
1:0 1010 |
2280 1:942669 1084 |
1:699641 10 14 |
Кафедра ТМ (СГАУ) |
Методы вычислений |
19 февраля 2012 г. |
46 / 50 |
Машинная арифметика, стандарт IEEE-754 Ошибки
Сложение и вычитание чисел
Значащих цифр в мантиссе двоичного числа в формате Single не более 24.
Если числа отличаются более чем в 223 (для single) и 252 (для double), то операции сложения и вычитания между этими числами невозможны.
300 =
1.00101100000000000000000 28
0:00001 =
0.00000000000000000000000010100111110001011011 28
300 + 0:00001 = 300
Кафедра ТМ (СГАУ) |
Методы вычислений |
19 февраля 2012 г. |
47 / 50 |
Задания
Задание 1.
Вычислите приближенное значение функции и оцените абсолютную погрешность результата. Запишите результат в форме f A(f ).
|
№ |
f (x; y) |
x |
A(x ) |
y |
A(y ) |
|
|
|
|
|
|||||||||||||
1 |
x + 2xy |
0.13 |
0.002 |
1.121 |
0.001 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
2 |
2x + xy |
0.13 |
0.002 |
1.121 |
0.001 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
3 |
x+2y |
0.13 |
0.002 |
1.121 |
0.001 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
x |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
4 |
x+2y |
0.13 |
0.002 |
1.121 |
0.001 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
2yx x |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
5 |
0.13 |
0.002 |
1.121 |
0.001 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6 |
2x+y |
0.13 |
0.002 |
1.121 |
0.001 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
2x |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
7 |
x + 2xy |
0.203 |
0.002 |
1.506 |
0.001 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
8 |
2x + xy |
0.205 |
0.002 |
1.505 |
0.001 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
9 |
x+2y |
0.051 |
0.002 |
0.911 |
0.001 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
x |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
10 |
x+2y |
0.063 |
0.002 |
1.121 |
0.001 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
x |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Кафедра ТМ (СГАУ) |
|
|
Методы вычислений |
|
|
|
|
19 февраля 2012 г. |
48 / 50 |
Задания
№ f (x; y) |
x |
A(x ) |
y |
A(y ) |
||
11 |
x+2y |
0.325 |
0.002 |
1.221 |
0.001 |
|
y |
|
|||||
12 |
2x+y |
0.235 |
0.002 |
2.121 |
0.001 |
|
2x |
2 |
|||||
13 |
x+y |
0.051 |
0.002 |
0.911 |
0.001 |
|
x |
2 |
|||||
14 |
1+y |
0.063 |
0.002 |
1.121 |
0.001 |
|
x |
|
|||||
15 |
x2+y |
0.325 |
0.002 |
1.221 |
0.001 |
|
y |
|
|||||
16 |
x+y |
|
0.235 |
0.002 |
2.121 |
0.001 |
1+x |
|
|||||
17 |
2x2 + y |
0.203 |
0.002 |
1.506 |
0.001 |
|
18 |
2x + y2 |
0.205 |
0.002 |
1.505 |
0.001 |
|
19 |
2x2 + y |
0.203 |
0.003 |
1.506 |
0.002 |
|
20 |
2x + y2 |
0.205 |
0.004 |
1.505 |
0.002 |
Кафедра ТМ (СГАУ) |
Методы вычислений |
19 февраля 2012 г. |
49 / 50 |
Задания
Задание 2.
Запишите десятичное число n10 в формате single2 в соответствии с правилами стандарта IEEE-754.
|
|
№ |
n10 |
№ |
n10 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
12:15 11 0:251 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
2 |
129:1 12 1:253 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
3 |
135:6 13 9:445 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
4 |
158:6 |
14 |
125:3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
5 |
0:153 15 6:525 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
6 |
253:2 |
16 |
3:214 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
7 |
31:45 |
17 |
2:125 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
8 |
21:21 |
18 |
1:562 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
9 |
15:12 |
19 |
1:521 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
10 |
12:81 20 0:093 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
Кафедра ТМ (СГАУ) |
|
Методы вычислений |
|
|
19 февраля 2012 г. 50 / 50 |