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

книги / Теория инженерного эксперимента

..pdf
Скачиваний:
18
Добавлен:
12.11.2023
Размер:
13.91 Mб
Скачать

9 - 3 . »

Напряжение на сетке

лампы

 

(X ), в

 

0

Анодное напряжение в момент

 

начала разряда (F ),

в

0

1

Т ю

 

о

50 100

— 15 — 20 — 25 — 30

125

180

250

315

9.4.Рассматривая отклонения точек от прямой, полу­ ченной в задачах 9.1, 9 . 2 и 9.3 методом наименьших квад­ ратов, вычислите среднее квадратическое отклонение для переменной F. (Заметим, что в задаче 9.1 переменная Y известна точно, но предполагается, что недостаточный контроль над экспериментом оказывает влияние только на переменную F.) В задачах 9.2 и 9.3 рассмотрите ве­ личину отрезков, отсекаемых прямой на координатных осях, и сравните величину отклонения прямой от начала координат со средним квадратическим отклонением.

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

Сопротивление, ом

5 2 ,5

5 3 ,5

5 4 ,5

5 5 ,5

Время

безотказной

работы,

 

 

 

мин

 

3 8 ,5

3 7 ,0

4 4 ,0

4 5 ,5

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

9 .6 . Методом наименьших квадратов постройте наи­ лучшую прямую для зависимости 1 /S от 1 Т по данным, приведенным в первой части примера 9.4, при условии, что все случайные ошибки содержит только 1/5. График, изображенный на фиг. 9.10, построен на глаз. Сравните ваш результат с этой прямой, а также сопоставьте фор­ мулы для этих прямых. Считаете ли вы, что при построе­ нии прямой по этим данным целесообразно применять метод наименьших квадратов?

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

что ваш

выбор правилен.

 

 

 

 

 

9.7.

При испытании системы охлаждения получены

следующие данные:

 

 

 

 

 

 

Снижение

темпера­

18,9

16,9

14,9

12,9

10,9

8,9

туры ДТ

19,9

Время работы 6,

3,45

10,87

19,30

22,80

40,10

53,75

мин

0

9 .8 .

При исследовании

нелинейной

схемы

получена

следующая зависимость:

 

 

 

 

Напряжение, в

67,7

65,0

63,0

61,0

58,25

56,25

Ток, а

 

2,46

2,97

3,45

3,96

4,97

5,97

9.9.Получены следующие данные о падении давления

впатрубке:

Скорость, ж/се/е

0,37

0,44

0,49

0,52

0,59

0,76

Падение давления,

0,200

0,25

0,290

0,32

0,43

0,70

мм pm. cm.

9.10. Имеются следующие данные о производитель­ ности насоса:

Производительность,

0

7,57

15,14

22,71

30,28

37,85

45,43

л/мин

Напор, ж

29,9

32,3

33,1

34,7

35,7

33,1

32,0

9.11. Получена следующая зависимость удельной теп­ лоемкости сплава от содержания в нем компонента А:

Содержание компонента А

в

0,2

0,4

0,6

0,8

1,0

сплаве,

%

0

Удельная

теплоемкость

0,2

0,27

0,35

0,47

0,69

1,0

Указание: проверьте функцию (7) из табл. 9.1.

9 .1 2ч В какой-либо или во всех приведенных выше пяти задачах исследуйте функцию, которую вы нашли для крайних точек, и укажите любые другие значения, пред­ ставляющие теоретический интерес.

9.13. Покажите, что формулы (1)—(7) в табл. 9.1 позволяют построить линейный график для данных, при­ веденных в^задачах 9.7—9.11.

9.14. Найдите уравнение наилучшей прямой, представ­ ляющей зависимость доли скоростного напора, сохранив­ шегося после частичного перекрытия трубопровода гидрав­ лической системы, от процента перекрываемой площади поперечного сечения трубопровода с регулируемым пе­ рекрытием. Кроме того, решите, целесообразна ли экст­ раполяция этих данных.

Сохранившаяся

доля

 

ско­

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ростного напора

 

 

 

0 ,8 6

 

0 ,8 5

 

0 ,6 0

 

0 ,5 5

0 ,5 2

0 ,4 5

Процент перекрытия

 

 

 

10

 

20

 

30

 

40

 

50

 

60 1*

 

 

 

 

 

 

 

 

 

 

 

Л И Т Е Р А Т У Р А

 

 

 

 

 

 

 

 

1.

A c t o n

F .,

 

A n a ly sis

of S tr a ig h t-lin e

D a ta ,

W iley , N .Y .,

1959 .

2 .

B a c o n

R .,

«B est» S tra ig h t

L in e am ong

th e

P o in ts,

A m . J . Phys.,

 

21

(6),

428— 4 4 5

(S ep tem b er 1953).

 

 

 

 

 

 

 

 

 

3 .

6

о w

k e r

 

A .,

 

L i

e

b e

r

m a n

G .,

E n g in eerin g

S ta tistic s ,

 

C h . 9,

P ren tice -H a ll,

E n g lew o o d

C liffs,

 

N .J .,

1959.

 

 

 

4 .

G u e s t

P .,

 

F ittin g

a

S tra ig h t

L in e

by th e

M ethod

o f

G rouping,

 

A m . J . P hys.,

18

(5),

3 24 (M ay 1950).

 

 

 

 

 

 

 

 

 

5 . H o e l s c h e r

R. P. ,

A r n o l d

J . N. ,

 

P i e r c e

S . H .,

 

G rap h ic

A id s

 

in

E n g in eerin g

C o m p u ta tio n ,

C h . 2,

M cG raw -H ill,

 

N . Y . 1952 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6 . M

a c k e

у

С . О .,

G rap h ical

S o lu tio n s, C h .

5

and

6,

W iley ,

N .Y .,

 

1943 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7 .

W

i 1 s

о

n

E . B .,

A n

In tro d u ctio n to

S c ie n tific

R esearch,

C h . 8,

 

M cG raw -H ill,

 

N .Y ., 1952 .

 

 

 

 

 

 

 

 

 

 

 

 

8 .

W о r t h

i

n

g

A .

G .,

G

e

f f n e r

 

J .,

T reatm en t

of E x p eri­

 

m en tal

D a ta ,

C h s.

1— 3

and

11, W ile y ,

N .Y .,

1943.

 

 

 

МАТЕМАТИЧЕСКИЙ АНАЛИЗ ДАННЫХ

Невозможно провести четкую границу между чисто статистическим анализом и графическим анализом, ко­ торый в свою очередь отличается от математического анализа. Метод наименьших квадратов позволяет полу­ чить математическое выражение кривой, которое нахо­ дится чисто численными способами и имеет прочную ста­ тистическую основу. Почти все, что можно сказать в за­ щиту такого довольно произвольного деления, исполь­ зуемого в последних четырех главах данной книги, это то, что вывод некоторых аналитических методов основан главным образом на статистических методах, рассмотрен­ ных в гл. 8 . Другие методы связаны с графическими ,по­ строениями (они описаны в гл. 9). Вопросы, которые не были рассмотрены в этих двух главах,'собраны в дан­ ной главе и названы математическими методами, хотя в некоторых случаях лучше было бы назвать их численными методами. Кроме того, необходимо подчеркнуть, что в ли­ тературе накоплен обширный материал по математиче­ скому анализу экспериментальных данных; здесь мы мо­ жем очень кратко коснуться лишь некоторых наиболее важных и элементарных вопросов. Такие важные методы, как гармонический анализ, численное решение дифферен­ циальных уравнений, построение номограмм, и множест­ во других вопросов можно найти в соответствующих спра­ вочных пособиях. Значительную часть этого специального материала экспериментатор обычно никогда не усваивает, тем более что он всегда может обратиться за помощью к специалистам. Методы, изложенные в данной главе, яв­ ляются довольно простыми, носят общий характер и используются многими экспериментаторами.

1 0.1 . Значащие цифры

Вопрос о значащих цифрах, вероятно, следовало бы рассмотреть ранее, хотя именно при операциях с числами и вычислениях они приобретают наибольшую важность. Значащие цифры несут информацию о значении измеряе­ мой величины. Если например, измеренная емкость кон­ денсатора равна 0,213 мкф, т. е. известны три значащие цифры, то этот факт не изменится, если емкость записать как 0,21310_в ф или 213 000 пкф. При записи таких боль­ ших чисел, как последнее, возможна путаница, так как можно подумать, что емкость действительно известна с точностью до шестой цифры. Чтобы исключить такие недоразумения, лучше всего записывать эту величину как 213-103 пкф.

Как и следовало ожидать, число значащих цифр вы­ бирается исходя из результатов анализа неопределенности, рассмотренного в гл. 2 и 3. Если, например, неопределен­ ность измерения составляет 1 0 %, то в лучшем случае мы можем указать не более двух значащих цифр. Напри­ мер, такой отсчет, как 98, означает, что результат заклю­ чен между 97 и 99, но в то же время с вероятностью 10% он может лежать между 95 и 101 или даже в более широ­ ких пределах. При неопределенности, равной 1%, мы

можем указать

три цифры,

при 0 ,1 % — четыре и т. д.

В естественном

стремлении

несколько оптимистически

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

Как показано в гл. 3, при обработке данных ошибки определенным образом комбинируются. Если допустить небрежность, то можно потерять значащие цифры. Пред­ положим, что при испытаниях на коррозию определяется потеря веса. Перед испытаниями вес образца составлял 1,54 кГ, а после испытаний— 1,52 кГ. Потеря веса со­ ставляет 0,02 кГ, т. е. мы имеем только одну значащую цифру. Пусть требуется сравнить этот результат с резуль­ татом при другом испытании, который также оказался равным 0,02 кГ. Очевидно, что мы не можем прийти к ка­ кому-либо выводу. Значащие цифры практически исчез­ ли, поэтому в этом случае необходимо использовать весы,

которые позволяли бы получить четыре, а то и пять зна­ чащих цифр.

Если исследователю предстоит большой объем обработ­ ки данных с применением технических средств или вычис­ лительных машин, то обычно в запасе необходимо иметь хотя бы одну значащую цифру. Пусть необходимо воз­ вести в квадрат скорость 107 м/сек, точный результат равен 11 449 м2/сек%. Поскольку вначале мы имели три значащие цифры, то, округляя, получаем 11Г400. Допу­ стим теперь, что к 11 449 нужно прибавить 6540, точный результат равен 17 989, а прибавляя 6540 к округленному числу, получаем только 17 940. Таким образом, при округ­ лении мы потеряли одну единицу третьего разряда. Если

бы мы

округлили

11 449

до

четырех знаков, что

дает

11 450,

то,

прибавив к этому числу 6540, получили бы

17 990,

а

при округлении

окончательного результата

получили

бы 18 0

0 0 (это

число лучше записать

как

180*102). В цифровых вычислительных машинах числа мо­ гут иметь до 1 2 знаков, и такие ошибки округления экспериментатор может практически исключить. Как только последовательность вычислений будет закончена, необходимо округлить результат до намеченного числа цифр, с тем чтобы он не вызывал у читателя недоверия или подозрения.

Стандартная логарифмическая линейка длиной 25 см при выполнении операций умножения и деления дает максимальную неопределенность 0,3—0,5% (для 95% всех отсчетов). Для линейки длиной 50 см неопределен­ ность в два раза меньше, а для спиральной линейки, имеющей шкалу длиной 150 см, она меньше в шесть раз. На такой линейке можно выполнять вычисления с точ­ ностью до четырех знаков, а для малых чисел — с точ­ ностью до пяти знаков. Четырехзначные таблицы ло­ гарифмов имеют точность до четвертого знака, и один знак необходимо добавлять в случае интерполяции. Точ­ ность номограмм, диаграмм и графиков в каждом случае необходимо устанавливать отдельно. Если графики по­ строены с соблюдением рекомендаций, изложенных в разд. 9.4, то инженер может снимать результат с точ­ ностью до половины наименьшего деления шкалы. При проведении эксперимента нельзя применять приборы.

имеющие различную точность. Точно так же неразумно применять точные приборы, имеющие ошибку порядка 0,1%, если получаемые данные будут обрабатываться с помощью простой логарифмической линейки. Аналогично продол­ жительные и утомительные вычисления на настольной счетной машине, имеющей большое число разрядов, пред­ ставляют собой лишь напрасную потерю времени, если данные были получены с неопределенностью 5%.

10.2. Подбор многочленов по эмпирическим данным

В предыдущей главе было рассмотрено несколько спо­ собов нахождения эмпирических уравнений графиче­ ским путем. Все эти методы основаны на том, что строит­ ся прямая либо непосредственно по данным, либо по ре­ зультатам, полученным после алгебраического преобразо­ вания переменных. Во многих экспериментах преобразо­ вание данных не позволяет получить линейный график и необходимо подбирать более общие функции, чем те, которые рассматривались в гл. 9. Наиболее широко ис­ пользуется многочлен, имеющий п + 1 членов:

Y = a + bX + cX2+ dX3+ • • • + пХп.

(10.1)

Если используются только два первых члена, то имеем общее уравнение прямой. Три члена дают параболу, про­ верка которой может производиться, как показано в табл. 9.1. Увеличивая число членов, можно подобрать более сложные кривые.

Обычно при подборе многочлена типа (10.1) на кривой выбирается столько точек, сколько постоянных должно быть определено, а затем решается система (n + 1) урав­ нений. При п — 2 или больше такой способ становится довольно громоздким и уравнения часто проще решить с помощью определителей1. К сожалению, уже при п = 2 получаем определители четвертого порядка, при исполь­ зовании которых часто встречаются гораздо большие трудности, чем при непосредственном решении системы трех уравнений. Вообще говоря, определители удобны

1 Краткая сводка определителей, используемых при подборе кривых, приводится в приложении 1 книги [6].

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

Ф и г. 10.1. График функции с максимумом, построенный по четырем точкам (одна из них в начале координат). Посколь­ ку вправо от максимума име­ ется точка перегиба, две из четырех точек находятся в этой части и только одна — слева.

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

Независимо от того, решается система уравнений не­ посредственно или с помощью определителей, исследова­ тель должен выбрать две—три или большее число точек таким образом, чтобы подобрать наилучшую функцию. Допустим, что для кривой, изображенной на фиг. 10.1, необходимо подобрать уравнение при п = 3. Левая часть кривой с непрерывно убывающим наклоном представляет собой функцию простого вида, а правая часть имеет точ­ ку перегиба. Таким образом, можно предположить, что подобрать уравнение для правой части кривой труднее, поэтому постараемся взять побольше точек на этом участ­ ке кривой. После того как будет найдено уравнение для п = 3, естественно проверить остальные точки, чтобы убедиться в том, что получено удовлетворительное соот­ ветствие для всей кривой. Если получено недостаточно хорошее соответствие, то можно подобрать четыре другие точки или перейти к уравнению для п = 4. Заметим, что если кривая проходит через начало координат, то постоян

ная а в формуле (10.1) равна нулю и система будет иметь на одно уравнение меньше (либо порядок определителя будет Меньше на единицу). Такой удобный результат можно получить всегда. Для этого достаточно перенести систему координат, как показано на фиг. 10.2. В систе-

Ф и г. 10.2.

Перенос

начала

i

координат в одну из

экспери-

ментальных

точек. При этом

у

число уравнений уменьшается

 

на единицу.

 

 

ме координат (X V ) имеем X' = X АХ

и У —

— Y AY. Многочлен второй степени

 

У = а + Ь Х + сХ2

(10.2)

в этом случае принимает вид

 

Y AY= a' -j-b' (X—ДЛ) + с' (X—АХ)2.

(10.3)

В системе координат (X', Y') постоянная а = 0. После определения Ь' и с' соотношение (10.3) можно преобразо­ вать и вернуться к первоначальной системе координат

У = Д У — ЬАХ + сАХЧЬХ+ с(Х2—2ХАХ),

или

У= АУ—bАХ + сАХ2+ X (Ь—2АХ) + сХ2.

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

Пример 10.1. На фиг. 10.3 показана плавная кривая, изображающая зависимость анодного тока от анодного

напряжения в триоде при постоянном напряжении на сетке, равном —5 в. Требуется найти уравнение этой кривой.

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

Ф и г . 10.3. График зависимости анодного тока от анодного напря­ жения в триоде (данные из примера 10.1). Пронумерованные точки (/, 2 и Э) использовались для построения кривой, а остальные точ­ ки (контрольные) иллюстрируют соответствие подобранного урав­ нения экспериментальным данным.

многочлен второго порядка. Перенеся начало координат в точку (20;0), получаем V = V — 20 и Г = I. В этом случае

Г = Ь ' (F—20)+с' (К—20)8

(при этом а' = 0). На кривой можно выбрать две точки с одинаковыми приращениями координаты V, одна*о необходимо иметь в виду, что в нижней части кривая имеет большую кривизну и здесь труднее подобрать уравнение.