Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Методы вычисления. Лекции 1-2

.pdf
Скачиваний:
33
Добавлен:
16.03.2015
Размер:
693.29 Кб
Скачать

Машинная арифметика, стандарт 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