книги / Теория инженерного эксперимента
..pdfсложным, он все же довольно прост, и с помощью линейки и циркуля (для деления отрезков в отношении 1 : 2 ) за несколько минут можно обработать 20—30 точек. Для проверки построение прямой можно произвести с друго го конца, откладывая отрезки, соответствующие V3 S.
Фи г . 9.4. Пример построения прямой графическим методом наи меньших квадратов (методом Асковица).
О экспериментальные точки; X точки, соответствующие 2/sS; ф точки на пря мой, построенной методом наименьших квадратов.
Третья точка также будет лежать на прямой. Это построе ние не показано на фиг. 9.4, но читатели могут попытать ся выполнить его самостоятельно.
Метод Асковица не является приближенным, если по оси X откладываются совершенно одинаковые отрезки, но он имеет те же недостатки, что и классический метод, а именно : только переменная Y может иметь ошибку, гра фик; должен быть линейным и все данные должны иметь одинаковую точность. Этот важный метод построения прямой игнорируется в большинстве учебников по мате-
магической статистике, очевидно, на том основании, что он позволяет лишь найти прямую и не дает какого-либо количественного показателя точности. Применяя простой метод, описанный в предыдущем примере, легко найти величину отклонения точек от прямой и вычислить сред нюю квадратическую ошибку. Данный метод вполне за служивает того, чтобы инженеры освоили его.
Пример] 9.3. Используя данные из примера 9.1, с по мощью графического метода Асковица найдите прямую наименьших квадратов.
Решение. Применение этого способа для построения прямой возможно благодаря тому, что каждый раз уве личение переменной L происходит на одинаковую вели чину. Построение прямой показано на фиг. 9.5. Полу чаем в точности такую же прямую, что и при использо вании численного метода, в пределах ошибок графиче ского построения.
В настоящее время известен ряд других методов пост роения прямой, не связанных с нахождением наимень ших квадратов [1J. Наиболее простым из них является так называемый метод группировки. Для применения это го метода необходимо знать хотя бы одну точку прямой (это может быть начало координат), так как определяется лишь угловой коэффициент, и не обязательно, чтобы слу чайную ошибку имела только переменная Y. Построение прямой производится очень быстро.
Если задан ряд экспериментальных точек с коорди натами X и Y, то угловой коэффициент прямой, соеди няющей эти точки, можно вычислить по формуле
(9 1 ) 2
где Y и X — координаты т точек, сгруппированных в одной части графика, a Y' и X ' — координаты т точек, находящихся в другой части графика. Нас интересует, сколько точек следует брать из общего числа полученных данных, чтобы составить две группы по т точек.
С помощью статистического анализа [4] можно пока зать, что если интервалы между значениями X примерно одинаковы, то имеющиеся данные необходимо разделить
Ф и г . 9 |
.5 . П остроение прям ой методом А сковица п о |
данны м из |
|||
прим ера |
9 .3 . П рим енение |
м етода |
в озм ож но, |
п оскол ьк у |
эти данны е |
имеют равны е приращ ения |
п о оси |
X . Чтобы |
не усл ож н я ть чертеж , |
||
здесь п ок азан о н ахож ден и е только одной из д в у х |
точек . |
О эксперим ентальны е точки; X точки, соответствую щ ие */«£; ф точки н а п р я мой, построенной методом наим еньш их квадратов.
на три равные группы и использовать группу точек с координатами (X, У), находящихся в верхней части графика, и группу точек с координатами (X ', У ), находящихся в нижней части графика, исключив среднюю группу то чек. В этом случае нарушается общее практическое пра вило статистического анализа: «Статистический критерий оказывается наиболее эффективным при использовании всех данных». Можно сказать, что если до проведения эксперимента известно, что будет использован метод груп пировки, то не следует тратить усилия на получение дан ных в средней части интервала значений X , а необходимо сосредоточить внимание на первой^и последней'третях этого интервала. Однако на практике могут потребовать ся данные для всего интервала значений, с тем чтоб#
Ф и г . 9 .6 . П остроение прям ой |
методом группировки для данны х |
||
об изм енении тем пературы из |
прим ера 9 .1 . П ри таком |
небольш ом |
|
числе данны х предпочтительно строить прям ую |
на |
гл аз. |
|
А — п рям ая, построенная методом |
наим еньш их квадратов; |
Б — п р ям ая, по |
|
строен н ая методом гр уп п и ровки . |
|
|
исследователь мог убедиться, что зависимость действи тельно линейная. Затем, если желательно применить ме тод группировки, средние точки можно отбросить. По* скольку данный метод не является точным, то нецелесо образно полностью полагаться только на него.
Пример 9.3а. Используя метод группировки, опреде* лите угловой коэффициент наилучшей прямой, соответст вующей данным из примера 9.1.
Решение. В формулу (9.12) можно подставить данные для двух первых и двух последних точек, опуская единст венную среднюю точку, так как общее число точек не делитсяна три. Затем находим
(20 + 22) - (5 + 7)
В данном случае Y = АТ и X — L.
На фиг. 9,.6 изображена прямая наименьших квадра тов, полученная в примере 9.1 или 9.2, и для сравнения проведена прямая через ту же точку L = 0, угловой коэф фициент которой найден методом группировки. При таком небольшом числе точек, очевидно, лучше всего провести прямую на глаз, чем применять этот приближенный ме тод. Однако при наличии 20—30 точек метод группиров ки обеспечивает гораздо меньшее отклонение от прямой наименьших квадратов. При большом числе точек он наи более удобен для быстрого, построения прямой. И нако нец, необходимо сделать одно краткое, но важное замеча ние. Ни один инженер не станет применять какой-либо статистический метод в ущерб наглядности или здраво му смыслу. Непринятие прямой, построенной методом группировки, означало бы нарушение всех канонов ин женерной целесообразности.
9.3. Исследование функций графическими методами
В предыдущих главах было показано, что различные функциональные соотношения можно преобразовать та ким образом, чтобы получить графики в виде прямой, что удобно при выборе интервалов между значениями пере менной (разд. 6.1), при экстраполяции (разд. 7.3) и при менении метода наименьших квадратов (разд. 9.1 и 9.2). В экспериментальной работе часто может возникать обрат ная ситуация. Имея данные, которые в определенной си стеме координат могут быть представлены в виде прямой, мы должны найти ее уравнение. Такое уравнение обычно называется эмпирическим. Это означает, что его находят на основе эмпирических данных, а не из теоретических соображений. Эмпирические соотношения получают при изучении физических, химических и электрических свойств в теплотехнике, гидромеханике и во многих но вых областях, еще не имеющих прочной теоретической базы.
Независимо от характера эксперимента основной зада чей является выбор и преобразование системы координат таким образом, чтобы полный набор данных давал по воз
можности прямую линию. Если полученные данные не образуют прямой на линейной графической бумаге, то можно попытаться построить график в логарифмических координатах (или наносить логарифмы значений X и Y на линейную графическую бумагу). В логарифмических координатах график простой, но важной функции
Y = k X a |
(9.13) |
имеет вид прямой. Переходя к логарифмам, получаем
logK =log£+alogX , |
(9.14) |
где k и а — согласующие постоянные. После линейного соотношения эта функция, по-видимому, наиболее часто встречается в экспериментальной работе и широко приме няется при анализе размерностей (см. разд. 4.6).
Имеется также третий тип графической бумаги — полулогарифмическая, когда одна шкала является лога рифмической, а другая — линейной. В этом случае по лучается прямая, если данные подчиняются закону
Y = k(W y * .
Данное выражение идентично следующему:
У =£г2-3026аХ.
После преобразования этой функции имеем
lo g r= lo g £ + o X . |
(9.15) |
Чтобы получилась прямая, шкала по оси Y должна быть логарифмической, а по оси X — линейной.
Иногда встречается бумага специального вида (напри мер, с тремя осями координат, гиперболическая), однако фактически в ее применении нет необходимости. Напри мер, гиперболическую функцию
Y = |
X |
|
(9.16) |
|
й -f~ |
ЬХ. |
|||
|
|
можно представить в виде прямой, построив в линейных координатах зависимость Х/Y от X или MY от MX.
В некоторых случаях получают U-образную или коло колообразную кривую в линейных координатах. Это сви детельствует о том, что имеет место параболическая функция или более общая полиномиальная функция
Y = a + b X + cX2. |
(9.17) |
При наличии у кривой нескольких экстремумов эта функция может иметь дополнительные члены. В случае уравнения параболы, состоящего из трех членов, воз можны несколько способов получения линейного графика. Возьмем на гладкой кривой, построенной по полученным данным, произвольную точку (Xlt Ух). Тогда
УЬ Х ^ -|-
Вычитая это уравнение из (9.17), получаем
У—Yi=b (Х - Х у) + с (Х2-* ? ).
Разделив обе части данного уравнения на ( X —Х х) , имеем
|
| = £ = |
6 + с ( Х |
+ Х |
х), |
где (b + |
сХi) — постоянная. Если это уравнение удовле |
|||
творяет |
полученным |
данным, |
то |
график зависимости |
(У — У х) /( Х — Х х) от X будет иметь вид прямой. Другой способ состоит в том, что дифференцируется
уравнение (9.17):
ж = ь + ь х .
Если значения X берутся через равные интервалы, то график зависимости AY/AX от X будет иметь вид прямой.
Эти и некоторые другие способы построения линейных графиков приведены в табл. 9.1. Доказательство этих правил оставим для упражнений.
Следует заметить, что в таких экспериментах нет не обходимости пользоваться только фактическими данными. Наилучшая методика состоит в том, что вначале данные наносятся на график в линейных координатах, а затем через точки проводится плавная кривая. После этого вы-
Таблица 9 J
Ф ункция |
Что необходимо сделать, |
|
чтобы получить на графике прямую |
||
|
1.У = аХ + b
2.Y = kXa
3 у __ £(Ю)аХ или Y = keaX
4. |
Y = |
Xj(a + ЬХ) или |
|
5. |
1/У = a/X + |
b |
|
Y = |
a + b X + |
cX2 |
|
6. |
Y = |
X/(a + ЬХ) + c |
7.Y = k(lO)bx+cX* или y = kebx+ cx*
В линейных координатах строится за висимость Y от X
На логарифмической бумаге строится зависимость Y от X
На полулогарифмической бумаге стро ится зависимость Y от X (ось X — линейная, ось Y —логарифмическая)
В линейных координатах строится за висимость 1/Y от \/Х или Х /Y от X
Влинейных координатах строится за висимость (У—Yi)/(X—Xj) от X
Влинейных координатах строится за
висимость (X—X x)/(Y—Ух) от X
В линейных координатах строится за висимость (log У—log Ух) от X
бирается наиболее подходящая функция и на кривой берутся произвольные точки для проверки соответствия принятой функции. Заметим, что такая алгебраическая проверка должна проводиться для всей области данных. В некоторых экспериментах для одной части данных гра фик имеет вид прямой, а в остальной части он является криволинейным. Однако даже в этих случаях целесооб
разно пытаться получить линейный график. Даже если |
|
для части данных можно построить |
линейный график, |
то это значительно удобнее, чем в том |
случае, когда для |
всей области данных график является криволинейным. Реальный пример такого анализа показан на фиг. 9.7 и 9.8. Студенты старшего курса провели эксперимент,
в процессе которого непрерывно регистрировалась тем пература в различных точках трубы отопительной систе мы. На фиг. 9.7 изображен обычный график зависимости температуры от времени. Температура измерялась с по мощью термопары, устанавливаемой в средней части тру бы. Между 260 и 316 °С наблюдается небольшой участок, где температура остается без изменения или даже несколь ко понижается. Чтобы получить линейный график для данных о неустановившемся тепловом потоке, часто
Ф и г . .9.7. График, построенный по данным об изменении темпера туры в дымоходе. Изображенная кривая построена на глаз с по мощью лекала. Стрелкой показан момент перехода на другую шкалу потенциометра. Читатель может убедиться, возможно ли обнаружить на таком графике отклонение результата.
Ф и г.*9.8. График, построенный по тем же данным, что и |
на |
фиг. 9.7, но время отложено на логарифмической оси. Стрелкой |
по |
казан момент перехода на другую шкалу потенциометра. В этом случае можно легко обнаружить резкое изменение наклона прямой.
строится график зависимости температуры от логариф ма времени (фиг. 9.8). На этом полулогарифмическом гра фике отчетливо видно резкое изменение наклона кривой при 288 °С. При наличии таких конкретных результатов не представляет труда проверить исходные данные и об наружить, что при этой температуре потенциометр, даю щий показания от нескольких термопар, переключается на другую шкалу. Таким образом, с достаточной уве ренностью можно утверждать, что данный прибор дает ошибочные показания, во всяком случае в одном из диа пазонов, и должен пройти калибровку.
Кривые, характеризующие физические и химические свойства, часто представляют собой результат взаимодей ствия двух или большего числа факторов; изменения наклона графика наблюдаются в тех случаях, когда один эффект преобладает над другим. Таким образом, наиболь ший теоретический интерес представляет точка, где одна кривая переходит в другую. Наглядным примером яв ляется график прочности при растяжении для мягкой стали, который до предела текучести имеет вид прямой, а затем становится криволинейным. Уортинг и Геффнер [8 ] описали случай из оптики, когда изучение точки соеди нения двух кривых привело к открытию двух различных оптических рецепторов, существующих в глазу челове ка, — палочек и колбочек.
Пример] 9.4. В примере 6.2 был описан эксперимент, спланированный по методу латинского квадрата, когда изучалась работа двигателя автомобиля «Додж» при раз личных значениях нагрузки, скорости и температуры. Одним из полученных графиков была зависимость удель ного расхода горючего от температуры двигателя. С плав ной кривой были сняты следующие данные:
Температура |
|
|
|
|
|
|
|
|
двигателя, |
49 |
54 |
66 |
66 |
71 |
77 |
82 |
88 |
°С |
||||||||
Удельный |
0,545 |
0,532 |
0,523 |
0,518 |
0,513 |
0,510 |
0,507 |
0,505 |
расход |
Найдите уравнение кривой.