Лабораторная работа № 1. Определение абсолютной и относительной погрешностей приближенных чисел. Оценка погрешностей результата

При проведении расчетов необходимо уметь практически оценивать погрешности ре-

зультата. Разумная оценка погрешности позволяет удерживать оптимальное число знаков при вычислениях, оптимизируя трудоемкость расчетов. Абсолютная и относительная погрешности вычисляются непосредственно по формулам (1.2.1) и (1.2.2). Число верных знаков определяется в соответствии с теоремами 1-3 параграфа 1.2.

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

Если, например, то Тогда Аналогично, если то Здесь при округлении числа использовано обычное правило округления по дополнению.

Погрешности арифметических операций вычисляются по формулам параграфа 1.3. При этом, например, погрешность суммы большого числа слагаемых, вычисленная по формуле (1.3.1), оказывается сильно завышенной. Поэтому часто применяют статистическую оценку погрешности суммы:

(1.6.1)

где - число слагаемых в сумме , то есть - номер десятичного разряда, до которого произведено округление. Формулу (1.6.1) используют при

Относительная погрешность разности двух положительных чисел, определяемая по формуле (1.3.3), может быть очень большой, особенно когда эти числа близки между собой. Чтобы избежать потери точности, в этом случае часто используют преобразования типа преобразований суммы в произведение.

Пример. Найтиалгебраическую сумму приближенных чисел и указать ее погрешность: (все знаки верные).

Начнем со слагаемых . Будем считать, что абсолютная погрешность каж-
дого из них не превосходит половины единицы младшего оставленного разряда,
то есть Тогда По формуле (1.3.1) , то есть абсолютная погрешность суммы не превосходит абсолютной погрешности наименее точного из слагаемых. Сохраним в расчетах только один запасной знак и будем округлять все слагаемые до 0.01. Тогда Окончательный результат округляем до 0.1. При этом к вычисленной ранее абсолютной погрешности 0.05 добавляется погрешность округления, также равная 0.05. Таким образом, окончательно Отсюда

Для оценки относительных погрешностей результата при умножении и делении формулы (1.3.4) используются редко. Обычно на практике выполняются условия и используется приближенное равенство

. (1.6.2)

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

Пример. Высота и радиус основания цилиндра измерены с точностью до 0.5% и оказались равными 2.06 и 6.5 см. Какова предельная относительная и абсолютная погрешности при вычислении объема цилиндра?

. Число может быть записано с произвольно высокой точностью. Для и имеем Таким образом, результат содержит не более двух верных знаков. При вычислении сохраним один запасной знак, округляя все числа до трех знаков: Тогда

Погрешности вычисления значений функций подсчитываются по формулам (1.4.1) и (1.4.2). Если значения функций положительны, то для относительной погрешности можно использовать формулу

. (1.6.3)

Для функции одной переменной легко можно определить допустимую погрешность аргумента по допустимой погрешности функции

. (1.6.4)

Для функции нескольких аргументов эта задача решается при введении дополнительного предположения – так называемого принципа равных влияний. При этом полагают, что в формуле (1.4.1) все слагаемые равны между собой, тогда

(1.6.5)

Пример. С каким числом верных знаков должен быть известен свободный член уравнения чтобы получить корни с четырьмя верными знаками?

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

По условию По формуле (1.4.2)

Если бы оба коэффициента влияли на точность корней, то можно было бы воспользоваться для вычисления принципом равных влияний. Однако в нашем случае известен точно, то есть и Вычислим

Тогда

Таким образом, для того чтобы получить корни с четырьмя верными знаками, необходимо в исходном уравнении задать не менее чем с тремя верными знаками.

Пример. Для определения модуля Юнга по прогибу стержня прямоугольного сечения применяется формула где - длина стержня, - основание и высота поперечного сечения, - стрела прогиба, - нагрузка. С какой точностью следует измерить длину и стрелу , чтобы погрешность не превышала 5.5% при условии, что известна с точностью до 0.1%, величины известны с точностью до 1%,

Применение формулы (1.4.2) для относительной погрешности дает . По условию задачи известно, что Тогда По принципу равных влияний на долю и приходится
по 0.7%. Тогда равно 0.7%, то есть надо измерять с точностью 0.7%=0.007, необходимо измерять точнее, Отсюда то есть

Задание № 1. Величина подъемной силы крыла самолета оценивается по формуле где - площадь проекции крыла на горизонтальную плоскость, - скорость натекания воздуха на крыло, - плотность атмосферы на заданной высоте, - угол атаки, отсчитываемый от направления нулевой подъемной силы, - коэффициент, зависящий от формы крыла. Требуется вычислить при заданных значениях и заданных абсолютных или относительных значениях этих величин. (В таблице 1 для нечетных вариантов задано , для четных .

Таблица 1

Вари- ант или или  
  0.005 0.001   0.01 0.99 0.01       0.001  
  0.004 0.001   0.01 0.98 0.01   0.01   0.04  
  0.005 0.001   0.01 0.99 0.01       0.001  
  0.004 0.001   0.01 0.98 0.01   0.01   0.02  
  0.005 0.001   0.01 0.99 0.01       0.001  
  0.004 0.001   0.01 0.98 0.01   0.01   0.04  
  0.005 0.001   0.01 0.99 0.01       0.001  
  0.004 0.001   0.01 0.98 0.01   0.01   0.02  
  0.005 0.001   0.01 0.99 0.01       0.001  
  0.004 0.001   0.01 0.96 0.01   0.01   0.04  
  0.005 0.001   0.01 0.95 0.01       0.001  
  0.004 0.001   0.01 0.95 0.01   0.01   0.04  
  0.005 0.001   0.01 0.95 0.01       0.001  
  0.004 0.001   0.01 0.95 0.01   0.01   0.02  
  0.005 0.001   0.01 0.95 0.01       0.001  
  0.004 0.001   0.01 0.95 0.01   0.01   0.04  
  0.005 0.001   0.01 0.95 0.01       0.001  
  0.004 0.001   0.01 0.95 0.01   0.01   0.02  
  0.005 0.001   0.01 0.95 0.01       0.001  
  0.004 0.001   0.01 0.95 0.01   0.01   0.04  
  0.005 0.001   0.01 0.89 0.01       0.001  
  0.004 0.001   0.01 0.88 0.01   0.01   0.04  
  0.005 0.001   0.01 0.89 0.01       0.001  
  0.004 0.001   0.01 0.88 0.01   0.01   0.02  
  0.005 0.001   0.01 0.89 0.01       0.001  
  0.004 0.001   0.01 0.88 0.01   0.01   0.04
  0.005 0.001   0.01 0.89 0.01       0.001
  0.004 0.001   0.01 0.88 0.01   0.01   0.02
  0.005 0.001   0.01 0.89 0.01       0.001
  0.004 0.001   0.01 0.86 0.01   0.01   0.04

Образец выполнения задания № 1

Пусть Относительная погрешность функции равна , абсолютная

Так как относительная погрешность велика, то значение функции следует вычислять не более чем с двумя – тремя знаками (с запасом), то есть (а не 59073 и не 59000, так как тогда было бы ). Отсюда Поскольку то результат целесообразно округлить до двух знаков. Окончательно,

Задание № 2. Найти абсолютную погрешность вычисления функции при заданных значениях аргументов.

1.

       
17.25¸0.01 2.1¸0.05 6.089¸0.003
34.725¸0.005 18.4¸0.03 0.121¸0.001
18.3¸0.01 2.06¸0.005 167.18¸0.01

2.

       
21.7¸0.02 3.45¸0.01 0.956¸0.004
9.272¸0.001 17.654¸0.0005 5.322¸0.0005
11.01¸0.005 2.19¸0.02 1.315¸0.0004

3.

       
1.5¸0.005 0.21¸0.03 3.78¸0.002
6.28¸0.02 11.615¸0.001 139.96¸0.006
5.137¸0.0003 11.09¸0.005 125.3¸0.01

4.

       
3.456¸0.001 2.5¸0.03 86.2¸0.05
7.12¸0.004 0.854¸0.0005 44.831¸0.001
9.7¸0.1 1.16¸0.004 67.18¸0.006

5.

       
49.83¸0.01 0.85¸0.003 4.523¸0.003
2.35¸0.001 0.536¸0.002 0.014¸0.0005
11.7¸0.1 3.7¸0.02 1.76¸0.002
3.147¸0.0001 0.7071¸0.00005 54¸0.5
1.78¸0.001 0.0465¸0.0001 3.14¸0.005

6.

       
2.435¸0.005 72¸0.3 1.84¸0.006
1.27¸0.01 3.374¸0.002 0.8345¸0.0004

7.

       
3.889¸0.001 145.5¸0.08 27.3¸0.04
0.8454¸0.0002 28.6¸0.1 0.93¸0.001
0.643¸0.0005 0.1736¸0.00005 1.73¸0.03

8.

       
1.245¸0.001 15.8¸0.03 3.45¸0.001
0.3746¸0.0001 0.9848¸0.0003 1.8387¸0.00005

9.

       
9.14¸0.005 3.6¸0.04 23.76¸0.003
1.2¸0.05 0.2518¸0.0001 0.85¸0.006
0.56¸0.01 0.34¸0.005 1.8¸0.01

10.

       
5.93¸0.005 2.786¸0.0006 0.75¸0.003
3.1415¸0.0001 0.734¸0.004 0.0761¸0.0001

11.

       
52.34¸0.01 0.0545¸0.00005 0.65¸0.02
2.0435¸0.0001 0.82¸0.01 6.3¸0.02

12.

       
0.285¸0.0002 0.2731¸0.0002 5.843¸0.001
0.64¸0.004 10.8¸0.02 4.17¸0.001

13.

       
8.23¸0.005 0.5324¸0.0001 2.5¸0.03
0.281¸0.001 0.032¸0.003 1.7¸0.01
11.7¸0.05 1.81¸0.005 31.456¸0.002

14.

       
23.76¸0.003 3.236¸0.002 17.3¸0.03
5.8¸0.02 1.315¸0.0004 4.431¸0.001
3.45¸0.001 1.57¸0.006 0.28¸0.006

15.

       
7.751¸0.0005 29.72¸0.01 3.14¸0.005
0.12¸0.01 78¸0.3 0.6384¸0.0002
1.0¸0.05 0.512¸0.0005 2.1951¸0.0001

16.

       
2.4¸0.02 0.121¸0.0002 10.536¸0.001
3.35¸0.005 5.8¸0.02 7.18¸0.006
8.235¸0.001 21.1¸0.1 0.8¸0.1
0.258¸0.0002 0.1756¸0.0001 0.4834¸0.0001

17.

       
0.3575¸0.0002 5.12¸0.01 3.6¸0.04
28.6¸0.1 1.27¸0.002 7.221¸0.003

18.

       
0.184¸0.006 0.357¸0.0005 0.63¸0.02
0.3575¸0.0002 0.25¸0.003 0.0562¸0.0001

19.

       
81.26¸0.01 341.17¸0.03 4.2¸0.05
0.142¸0.0003 0.642¸0.0005 0.28¸0.006

20.

       
4.218¸0.001 10.5¸0.1 1.0397¸0.00005
0.15¸0.005 0.0545¸0.0001 81.26¸0.02

21.

       
3.0567¸0.0001 0.678¸0.003 5.72¸0.02
28.3¸0.02 7.45¸0.01 3.2817¸0.0001
6.3¸0.02 0.82¸0.01 4.2¸0.05

22.

       
4.16¸0.005 7.27¸0.01 12.163¸0.002
5.205¸0.001 6.3¸0.04 3.244¸0.001
2.39¸0.01 4.198¸0.0005 9.3¸0.1

23.

       
4.22¸0.004 3.4¸0.1 5.82¸0.005
7.5¸0.05 8.13¸0.002 8.632¸0.003
5.205¸0.001 1.315¸0.0005 0.14¸0.01

24.

       
0.28¸0.006 0.156¸0.0005 0.08¸0.01
10.536¸0.002 7.18¸0.01 31.110¸0.0005

25.

       
0.235¸0.0003 0.318¸0.0002 1.8¸0.01
48.39¸0.001 0.46¸0.002 3.108¸0.0003

26.

       
0.44¸0.03 0.772¸0.0005 1.49¸0.01
9.542¸0.001 4.218¸0.001 29.725¸0.0005

27.

       
3.35¸0.01 6.17¸0.02 5.441¸0.0005
0.383¸0.0002 0.891¸0.0005 1.3¸0.01

28.

       
17.29¸0.005 28.3¸0.1 32.14¸0.02
72¸0.3 42.08¸0.005 116.929¸0.0005

29.

       
1.29¸0.005 11.3461¸0.0001 1.84¸0.006
2.316¸0.001 6.83¸0.03 0.8345¸0.0004
3.274¸0.0005 9.1¸0.1 32.7¸0.04

30.

       
2.45¸0.01 7.834¸0.001 0.15¸0.01
0.6294¸0.0001 0.1113¸0.00005 1.29¸0.005
17.3¸0.03 11.99¸0.005 2.4¸0.02

Образец выполнения задания № 2

Дано

   
1.1753¸0.0003
0.318¸0.0005

Абсолютные погрешности аргументов Относительные погрешности будут равны: . По формуле (1.4.2) получим где Само значение функции равно Следует округлить значение функции до трех значащих цифр. Итак,

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

По ходу выполнения лабораторных работ по курсу “Численные методы” научимся быстро и легко решать в среде Mathcad простейшие математические задачи. Подробная характеристика и описание пакета даны в специальной литературе справочного характера. Рассмотрим конфигурацию Mathcad 7 Professional, входящего в математическое обеспечение вычислительной лаборатории кафедры вычислительной и прикладной математики БГТУ.

Пакет запускается стандартным образом – двойным щелчком левой кнопки мыши по пиктограмме Mathcad. Происходит автоматическая загрузка всех программ, необходимых для работы пакета, затем открывается окно приложения, имеющее стандартную Windows–структуру.

 
 

В первой строке приведено имя приложения – Mathcad Professional, затем имя рабочего файла. Справа в строке находятся три функциональные кнопки: свернуть, развернуть на полный экран, закрыть окно приложения. Вторая строка – это строка меню, состоящая как из стандартных для Windows–приложений пунктов: Файл (File), Редактировать (Edit), Просмотр (View), Формат (Format), Окно (Window), Помощь (Help), так и специфических, созданных специально для пакета Mathcad: Вставка (Insert), Математика (Math), Символьные операции (Symbolics). Щелчок по любому пункту меню открывает подменю со списком доступных операций.

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

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

Общие арифметические операторы;

Знаки отношения;

Различные двух -и трехмерные графики;

Матричные и векторные операции;

Производные, интегралы и пределы. Ряды и

произведения;

Программные структуры;

Греческие буквы;

Ключевые слова символьных вычислений.

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

Под панелями инструментов находится окно рабочего документа Mathcad. Содержимое этого окна можно редактировать, печатать, сохранять в файлах на диске и так далее. Для просмотра всего пространства служат две полосы прокрутки: вертикальная и горизонтальная.

Содержание первой лабораторной работы таково, что Mathcad можно использовать как простой калькулятор. Опишем пример подобных вычислений, определив из образца выполнения задания № 1.

После щелчка в любом месте рабочего документа появляется небольшой крестик. С этого места в рабочем документе будет помещаться весь ввод с клавиатуры. Введем символы в следующей последовательности: . После набора
знака = Mathcad вычисляет выражение и выводит результат. На экране отобразится формула в точности в том виде, как ее пишут на доске: .

Для удаления любой цифры в выражении надо щелкнуть мышью справа внизу от нее и нажать клавишу - (Backspace). Для удаления всего выражения щелкаем мышью по любому месту в выражении и нажимаем клавишу -(Space) до тех пор, пока все оно не будет выделено угловой синей рамкой. После этого нажимаем клавиши (Backspace) и -(Delete). Выражение исчезнет с экрана.

Чтобы определить любую переменную, нужно напечатать ее имя, затем набрать символ . В рабочем документе этот символ отобразится, как знак присваивания . После этого нужно набрать значение, присваиваемое набранной переменной, например, , и либо щелкнуть мышкой по свободному месту вне поля ввода, либо нажать клавишу - (Enter).

Аналогично вычисляется значение функции. Для этого нужно начать с определения функции, например, набрав . На экране отобразятся следующие выражения: . Нижние индексы вводятся нажатием клавиши . Что бы ни печаталось после [, все остается в нижнем индексе, пока не будет нажата клавиша , чтобы выйти оттуда.

Теперь можно использовать определенную функцию , чтобы вычислить ее значение для произвольных значений аргументов. Чтобы вычислить функцию для , достаточно напечатать . Mathcad возвратит соответствующее значение . Этого же результата можно добиться, подставляя числовые значения аргументов в функцию: .

Определим теперь дискретные значения аргументов . Это делается следующим образом. Сначала печатается имя дискретного аргумента. Это должно быть простое имя. Никакие нижние индексы или функциональные определения недопустимы. Затем набирается двоеточие и диапазон значений, например, . На экране отобразится: . Это самое простое определение дискретного аргумента. Здесь величина шага равна единице. Mathcad допускает дискретные аргументы со значениями, расположенными от любого числа до любого другого значения с произвольным постоянным шагом. Например, введем .

Здесь 0 - первое значение, принимаемое аргументом , 3 - второе значение аргумента в диапазоне, 6 - последнее значение. Величина шага определяется автоматически: . На экране появится выражение . Далее вводим . На экране отобразится:

Mathcad содержит восемь переменных, значения которых определяются сразу после запуска программы. Это предопределенные или встроенные переменные. Одна из таких переменных - - определяет индекс первого элемента массива. По умолчанию в Mathcad координаты векторов и строки матриц нумеруются с нуля (). Эти восемь переменных можно переопределять в любом месте рабочего документа. Как и в электронных таблицах, любое изменение содержимого рабочего документа Mathcad вызывает обновление всех зависимых результатов и перерисовку графиков.

 
 

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

Наиболее часто употребляемые типы графиков:

X-Y Plot – графики функций одной переменной в декартовых координатах;

Polar Plot - графики функций одной переменной в полярных координатах;

Surface Plot – графики поверхностей (функций двух переменных) в декартовых

координатах;

Contour Plot – графики линий уровня функций двух переменных в декартовых

координатах;

3D Scatter Plot – изображения точек в трехмерном пространстве в декартовых

координатах;

3D Bar Chart – трехмерные гистограммы.

Все перечисленные графики, как это принято во всех Windows–приложениях, можно вызвать из главного меню выбором соответствующих операций: Вставка (Insert) – Графика (Graph).

Построим график функции из образца выполнения задания № 2. Это функция , ее график называется циссоидой. Полярное уравнение этой кривой имеет вид . Для построения графика в декартовой и полярной системе необходима явная зависимость или . Определим . Для двух ветвей графика необходимо определить две разные функции. Введем с клавиатуры . Теперь выберем пункт Декартов график из меню Графика. Появляется пустой график с полями ввода для выражений, отображаемых по осям графика. Множество точек, из которых состоит график, определяется дискретными аргументами. В поле ввода под осью абсцисс нужно ввести имя аргумента , а в поле ввода напротив середины оси ординат ввести и через запятую . Оставшиеся поля предназначены для ввода границ на полях (максимального и минимального значений и ). Если оставить эти поля пустыми, Mathcad автоматически заполнит их при создании графика. После щелчка вне графика Mathcad вычисляет и строит точки обоих графиков. На экране появится

.

Двойной щелчок в любом месте графика вызывает панель форматирования. Здесь можно установить все доступные характеристики графика. Те же действия выполняются при построении графика в полярной системе координат. С учетом периодичности и нечетности синуса и тангенса для отображения обеих ветвей графика надо набрать , затем щелкнуть по пиктограмме полярного графика; по оси аргумента в пустом графике ввести , по оси функции . На экране отразится . Число входит в восемь предопределенных переменных. В численных расчетах Mathcad использует значение с учетом 15 значащих цифр. Чтобы ввести , надо нажать .

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

.

Далее, нажав на кнопку графика поверхности в панели графиков, введем

в первом графике в поле ввода матрицу М, во втором - матрицу М1. На экране монитора получим следующую картину:

 
 

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

,
а ее линии уровня определяются уравнением и получаются сечением поверхности плоскостями, параллельными координатной плоскости . Итак, введем

На экране отобразится:

 
 

С помощью панели форматирования первый график повернут на , и на нем убраны невидимые линии, оси второго графика оцифрованы по в пределах от 0 до 30.


Понравилась статья? Добавь ее в закладку (CTRL+D) и не забудь поделиться с друзьями:  




Подборка статей по вашей теме: