Математическая теория обработки результатов экспериментов (На примере машиностроения )

получаемых решений.

1. Графический метод обработки результатов

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

исследуемыми величинами с последующим определением уравнения зависимости

между ними.

Графики строят прежде всего в равномерных шкалах. Если характер связи

между исследуемыми величинами неизвестен, то сначала проверяют совпадение

экспериментальных точек с заданной кривой. Если предварительные сведения о

характере уравнения отсутствуют, то первым этапом обработки данных является

нахождение кривой, совпадающей с опытными точками. Эта задача решается

методом подбора. Можно использовать эталон ( кальку с предварительно

вычерченным на ней семейством кривых с различными параметрами. Естественно,

что масштаб кальки и эмпирической кривой должен быть одинаков.

Построенный по опытным данным отрезок кривой может совпадать с

большим количеством различных кривых, проходящих достаточно близко к

опытным точкам. В этом случае выбирают кривую с наиболее простым и удобным

в использовании уравнением. Иногда эмпирическая кривая может иметь перегибы

или состоять из отдельных ярко выраженных участков. Однако при этом

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

Уравнение зависимости между исследуемыми величинами при графическом

методе просто определяется тогда, когда эмпирические точки достаточно

хорошо совпадают с прямой линией, т.е. описываются уравнением y = ax + b,

где a, b ( коэффициенты, подлежащие определению.

Определение коэффициентов при графическом методе основано на (способе

натянутой нити(. Нанеся результаты эксперимента на график (лучше, если он

выполнен на миллиметровке), подбираем графическую прямую, ближе всего

подходящую к нанесенным точкам. Выбрав положение прямой, определяем две

произвольные точки на этой прямой (не обязательно являющиеся точками

эксперимента), определяем их координаты (x1; y1), (х2; y2). И для

определения коэффициентов а и b получаем два простых уравнения

ах1 + b = y1;

ах2 + b = y2.

На рис. 10 приведена иллюстрация этого метода. Точки ( результаты,

полученные в эксперименте. Прямая проведена на глаз как можно ближе к

экспериментальным точкам. На прямой выбраны точки М (2; 4) и N (13; 10).

Коэффициент а характеризует угол наклона прямой.

Поэтому

[pic]

[pic].

Таким образом y = 0,55х + 2,9.

Рис. 10. Графический метод интерполяции

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

характер, то графическим способом в системе координат с равномерными

шкалами определить коэффициенты кривой затруднительно. Но достаточно

большой класс нелинейных зависимостей путем замены переменных и

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

и далее использовать способ натянутой нити.

2. Функциональные шкалы и их применение

Пусть функция y = ((х) непрерывна и монотонна на некотором промежутке

( a; b (. Возьмем ось ОМ, на которой будет строиться шкала, выберем на ней

точку начала отсчета О и установим масштаб (. Функциональная шкала строится

следующим образом.

Разбив интервал ( а; b ( на равные части, вычисляем значение функции

((х) в каждой из точек деления и отложим на оси ОМ для каждой точки отрезок

(((х). Получающаяся при этом точка снабжается отметкой х, т.е.

откладывается в выбранном масштабе значение функции, а надписывается

значение аргумента.

Иногда начало шкалы помещают в первую точку отсчета, т.е. точку с

надписью а совмещают с 0. Тогда точка х будет находиться в конце отрезка (

( ((х) ( ((а) (. Полученная шкала позволяет судить о поведении функции на

рассматриваемом участке: большие промежутки между отметками укажут, что

функция изменяется быстрее, чем там, где эти промежутки малы.

Выбор масштаба ( определяет длину шкалы. Чаще поступают наоборот:

задаются длиной шкалы l и определяют масштаб.

[pic] ( ( = [pic].

Пример. Построим функциональную шкалу для функции y = x2 на участке (

1; 2 (. Зададимся длиной шкалы l = 12 см. Тогда ( =

[pic] см. Разобьем отрезок ( 1; 2 ( на десять равных частей и вычислим

значения функции во всех точках деления. Совместим начало шкалы с точкой

отсчета х = 1. Результаты расчета сведены в табл. 2, а функциональная шкала

приведена на рис. 11.

Таблица 2

Расчет функциональной шкалы y = x2

| х |1,0 | 1,1 | 1,2 | 1,3 | 1,4 | 1,5 | 1,6 | 1,7| 1,8| 1,9 | 2,0|

| х2 |1,0 |1,21 |1,44 |1,69 |1,96 |2,25 |2,56 |2,89|3,24|3,61 |4,00 |

|х2(1 | 0 |0,21 |0,44 |0,69 |0,96 |1,25 |1,56 |1,89|2,24|2,26 |3,00 |

|4(х2(1)|0 |0,84 |1,76 |2,76 |3,84 |5,00 |6,24 |7,56|8,94|10,44|12,0 |

1,0 1,1 1,2 1,3 1,4 1,5 1,6 1,7 1,8 1,9 2,0

Рис. 11. Функциональная шкала y = x2

С помощью функциональных шкал графики многих функций могут быть

преобразованы к прямолинейному виду.

Например, уравнение параболы y = x2. Если на оси OY нанести

равномерную шкалу, а на оси OX1 шкалу квадратов х1 = х2, то получится

сетка, где уравнение параболы имеет изображение прямой линии ( y = x1 ),

проходящей через начало координат.

Особенно часто используются различные логарифмические функции, с

помощью которых можно (выпрямлять( графики степенных и показательных

функций. Например, y = aebx; lg y = (b lg е) х + lg a. Полагая lg y = y1,

lg a = A, b lg e = B запишем исходное уравнение в виде y1 = А + Вх, откуда

видно, что оставив равномерной шкалу х и построив логарифмическую шкалу y1,

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

сетка называется полулогарифмической.

Очевидно, что такого рода преобразования возможны и в более общем

случае. Всякая неявная функция, заданная соотношением вида

а((х) + b((y) + с = 0,

где a, b, с ( постоянные, будет изображаться прямой линией на

функциональной сетке, где на оси ОХ построена шкала ((х), а на оси OY (

шкала функции ((y). Естественно, что функции ((х) и ((y) должны

удовлетворять условиям непрерывности и монотонности. В табл. 3 приведены

преобразования для некоторых функций.

Таблица 3

Линеаризация некоторых функций

|y=[pic] | ( |[pic] | |

|y=[pic] |[pic] |[pic] | |

| | |[pic] | |

|y=[pic] |[pic] | | |

Из сказанного ясна роль функциональных сеток при обработке

результатов эксперимента. Если результаты эксперимента располагаются вблизи

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

типа функцией ее лучше всего приближать. Переведя полученные

экспериментальные данные на функциональные сетки можно оценить на какой из

них эти данные ближе всего подходят к прямой и, следовательно, какой

функцией лучше всего описываются.

3. Аналитические методы обработки результатов

Графический метод обработки результатов обладает наглядностью,

относительной простотой, однако его результаты содержат определенную

субъективность и относительно низкую точность.

Аналитические методы лишены в какой ( то степени указанных

недостатков и позволяют получить результат для более широкого класса

функций с большей точностью, чем графический метод.

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

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

Рассмотрим некоторые из существующих способов.

1. Способ средней

Допустим, что имеется n сочетаний xi, yi, полученных при

эксперименте. Даже в том случае, если между х и y теоретически установлена

функциональная связь ( в данном случае предположим, что линейная ), то

наблюдаемые значения yi будут отличаться от ахi + b вследствие наличия

экспериментальных ошибок. Обозначим через (i соответствующую ошибку

(i = yi ( axi ( b (i = 1, 2, ..., n)

Если выбирать параметры а и b так, чтобы для всех n наблюдений ошибки

уравновешивались, т.е. [pic], то это привело бы нас к одному уравнению,

тогда как для нахождения двух коэффициентов (а, b) их требуется два.

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

произведенных наблюдений в целом, но и для каждой группы, содержащей

половину ( или почти половину ) всех наблюдений в отдельности.

В этом случае можно прийти к системе уравнений

[pic] ,

где m ( число наблюдений в первой группе.

Данную систему уравнений запишем теперь в виде

[pic].

Изложенное показывает, что метод средних (уравновешивает(

положительные и отрицательные отклонения теоретической кривой от

экспериментальных значений.

Пример.Используя данные рис. 10 определим коэффициенты а, b методом

средней. Для этого семь измерений разделим на две группы m = 3 первых

значений, n ( m = 4 последующих

[pic]; [pic];

[pic] ; [pic].

Получаем систему

[pic]

Решая систему находим

[pic];

b = [pic]

Таким образом способ средней дает прямую

y = 0,55х + 3,11.

В сравнении с графическим способом коэффициенты а совпадают и имеется

различие в коэффициенте b.

3.3.2. Метод наименьших квадратов

В методе средних при определении коэффициентов уравнения

использовалось условие равенства нулю алгебраической суммы отклонений

результатов эксперимента от теоретической кривой ( в частном случае прямой

). Очевидно, что при этом (i могут быть значительной величины. Имеет

значение только (уравновешивание( положительных и отрицательных отклонений.

Поставим теперь задачу нахождения по результатам наблюдений наиболее

вероятные значения неизвестных коэффициентов.

Предположим, что искомая зависимость y = ((х) существует. Тогда

параметры этой линии необходимо выбрать таким образом, чтобы точки yi

располагались по обе стороны кривой y = ((х) как можно ближе к последней.

Предположим, что разброс точек yi относительно y = ((х) подчиняется закону

нормального распределения. Тогда мерой разброса является дисперсия (2 или

ее приближенное выражение ( средний квадрат отклонений.

[pic].

И требование минимального разброса будет удовлетворено, если минимизировать

выражение ( (yi )2. Как известно, необходимым условием того, что функция

приобретает минимальное значение, является то, что ее первая производная (

или частные производные для функции многих переменных ) равна нулю.

Применение метода наименьших квадратов имеет смысл, если число

экспериментальных точек n больше числа определяемых коэффициентов.

Рассмотрим реализацию метода наименьших квадратов применительно к

уравнению вида y = ax + b.

Для нахождения коэффициентов а, b искомой прямой необходимо

минимизировать сумму квадратов расстояний (yi по ординате от точки (хi; yi)

до прямой ( см. рис. 12 ). Расстояния (yi определятся

(yi = yi ( axi ( b.

Рис. 12. К способу наименьших квадратов

Для минимизации [pic] приравниваем к нулю производные этой суммы по

параметрам а, b:

[pic];

[pic].

Преобразуем эту систему

[pic]

Получим систему нормальных уравнений метода наименьших квадратов.

Решая ее относительно а, b получаем:

[pic] ; [pic].

Вычисляя из n опытов необходимые суммы и производя указанные

действия, получаем величину коэффициентов а, b.

Как видно, способ наименьших квадратов достаточно громоздок и при его

применении широко используется вычислительная техника. Метод наименьших

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

если определяются параметры квадратичной зависимости:

y = ах2 + bx + с,

то

[pic].

Дифференцируя это соотношение по а, b, с получаем систему нормальных

уравнений:

[pic][pic]

Из этой системы можно определить параметры а, b, с.

При использовании метода наименьших квадратов при других

нелинейностях, удобнее будет линеаризовать исходные зависимости.

В табл. 4 приведены системы нормальных уравнений для некоторых

исходных уравнений.

Таблица 4

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

|Исходное |Система нормальных уравнений |

|уравнение | |

| |[pic] |

|y=axb | |

| |[pic] |

|y=a(lgx+b | |

| |[pic] |

|y=eax+b | |

| |[pic] |

|y=aebx[pic] | |

| |[pic] |

|y=[pic] | |

| |[pic] |

|y=[pic] | |

| |[pic] |

|y=[pic] | |

Примечания: 1. Величины х, y обозначают значения величин хi, yi в i(ом

опыте;

2. Знак ( обозначают сумму величин от i = 1 до i = n, где n

( число равноточных измерений.

3. Интерполирование функций

Известно, что под интерполированием понимают отыскание значений

функции, соответствующих промежуточным значениям аргумента, отсутствующим в

таблице логарифмов, тригонометрических и др. функций.

В общем смысле можно сказать, что задача интерполирования обратна

задаче табулирования функций. При интерполировании по таблице значений

функции строится ее аналитическое выражение, т.е. по значениям функции yo,

y1, ..., yn при значениях аргумента хо, х1, ..., хn определяется выражение

неизвестной функции.

Понятно, что через данные точки ( даже большого числа ) можно

провести множество различных кривых. Поэтому существует интерполирование в

различных функциях F (х). Чаще всего требуют, чтобы функция F(х) была

многочленом степени на единицу меньшей, чем число известных значений.

Таким образом, задачу интерполирования функций можно сформулировать

следующим образом.

Для данных значений х ( хо, х1, ..., хn и y ( yo, y1, ..., yn найти

многочлен y = F (х) степени n, удовлетворяющий условиям F (хо) ( yo, F (х1)

( y1, ..., F (хn) ( yn. Точки хо, х1, ..., хn называют узлами интерполяции.

Многочлен F (х) ( интерполяционным многочленом , а формулы его построения (

интерполяционными формулами.

Как видно из описания сущности интерполирования, в отличии от

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

метода наименьших квадратов ), интерполяционный многочлен опишет кривую,

проходящую точно через заданные точки.

4. Параболическое интерполирование

При параболическом интерполировании в качестве интерполяционного

многочлена F (х) принимают многочлен n ( ой степени вида

F (х) ( ао + а1х + а2х2 + ... + аnxn.

Используя свойство прохождения функции F (х) через заданные точки

для неизвестных коэффициентов аi можно составить n + 1 уравнений с n + 1

неизвестным:

ао + а1хо + а2хо2 + ... + аnхоn ( yo;

ао + а1х1 + а2х12 + ... + аnх12 ( y1;

....................................................

ао + а1хn + а2хn2 + ... + аnхn2 ( yn.

Эта система имеет единственное решение, если значения хi отличны друг

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

системы. Перед рассмотрением общего способа решения, рассмотрим простой

пример.

Дано: хо = 0, х1 = 1, х2 = 2, yо = 1, y1 = 1, y2 = 3. Определить

многочлен F (х).

Записывая многочлен F (х) в виде

F (х) = ао + а1х + а2х2

составим систему уравнений

[pic] или [pic]

откуда ао = 1, а1 = ( 1, а2 = 1 и интерполирующий многочлен имеет вид

F (х) = 1 ( х + х2.

Теперь рассмотрим общий подход к отысканию интерполяционного

многочлена F (х), не решением системы, а непосредственной записью.

Определим выражение для многочлена, принимающего в точке х = хо

значение yо = 1, а в точках х = х1, х2, ..., хn ( значения y1 = y2 = ... =

yn = 0. Очевидно, что многочлен будет иметь вид

[pic].

Здесь при х = хо числитель и знаменатель равны, а при х = х1, х2,

..., хn ( числитель равен нулю.

Теперь построим многочлен Fо (х), принимающий в точке хо значение yо

и обращающийся в нуль для значений х = х1, х2, ..., хn. Учитывая предыдущее

построение можно записать

[pic].

Теперь можно записать многочлен F (х) для произвольного значения хi (

i = 0, 1, 2, ..., n ) принимающего значения F (хi) = yi, а во всех

остальных точках х ( хi значение, равное нулю

[pic].

Как видно из записи, числитель не будет содержать выражения (х ( хi),

а знаменатель ( (хi ( хi), т.е. выражений, обращающих числитель и

знаменатель в нуль.

Искомый многочлен будет равен сумме

[pic],

т.е. снова в каждой точке хi одно из слагаемых принимает нужное значение

yi, а все остальные обращаются в нуль.

В развернутом виде

[pic]

=[pic]

... + [pic].

Полученная формула называется интерполяционной формулой Лагранжа.

Используя формулу Лагранжа запишем многочлен F (х) для разобранного

выше примера.

[pic]

= [pic]

=[pic].

Получили тоже самое выражение, что и ранее.

Контрольные вопросы

1. Назначение графического метода обработки результатов;

2. Сущность графического метода обработки результатов;

3. Понятие и назначение функциональной шкалы;

4. Выбор масштаба функциональной шкалы;

Страницы: 1, 2, 3, 4, 5

МЕНЮ

Математическая теория обработки результатов экспериментов (На примере машиностроения )

ИНТЕРЕСНОЕ