ж) Пример определения кинематических характеристик по стробоскопическим фотографиям

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

 

Задание 1. Найти кинематический закон движения точки.

Спроецируем точки на координатные оси с учётом масштаба и выпишем в таблицу 1 значения координат точки, считая, что фотографирование началось при t = 0 и движение происходит по часовой стрелке. Измерения координат x и y прямые, поэтому оценим их погрешности по методике, изложенной в пункте б). Поскольку в данном случае нет особого смысла много раз измерять координаты, ибо мы будем получать всё время один и тот же результат, то следует предположить, что ∆хразбр. = ∆yразбр. = 0. Это не значит, конечно, что случайных ошибок нет – просто они меньше точности используемых инструментов. Приборная погрешность при измерении стандартной линейкой длиной 200 мм составляет ∆хпр. = 0,2∙ =0,13 мм. Погрешность отсчёта и округления при округлении координат до 1 мм составит 0,5 мм. Следует учесть неидеальность процедуры проектирования, которая также приводит к погрешности отсчёта и округления и составляет примерно 0,5 мм (подумайте, почему!) результирующая погрешность будет равна по формуле (10):

 

 

(

Рисунок 6 – Стробоскопическая фотография движения материальной точки

 

Таблица 1

t, cек 0,10 0,20 0,30 0,40 0,50 0,60 0,70 0,80 0,90 1,0
x, см 1,0 2,0 3,0 4,0 5,0 6,0 7,0 8,0 9,0 10,0 11,0
y, см 8,0 8,9 9,6 10,1 10,5 10,6 10,3 10,1 9,5 8,8 8,0

 

Для установления вида функциональной зависимости x = x(t) изобразим данные таблицы 1 на рисунке 7, откладывая время по горизонтали, координату x – по вертикали (в том же масштабе, что и на рисунке 6, руководствуясь при этом правилами, изложенными в пункте д).

При этом учитываем, что погрешность ∆t задана неявно (она равна 0,005 с). Из рисунка 7 сразу видно, что искомая функциональная зависимость x = x(t) линейная. Задача, следовательно, состоит в том, чтобы провести по точкам на рисунке 7 прямую, наилучшим в некотором смысле образом соответствующую этим точкам. Можно, конечно, это сделать графически, однако это не даёт полной уверенности, что прямая – наилучшая.

 

 

Рисунок 7 – Зависимость x = x(t)

 

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

Идея метода состоит в следующем. Пусть уравнение искомой прямой имеет вид x = at + b, где a и b – постоянные, подлежащие определению. При каждом значении времени ti (0,1,…10) найдём величину(ati + b – xi)2, представляющую квадрат разности между экспериментальным значением величины xi и значением (ati + b), вычисленным по формуле, выражающей ожидаемую линейную зависимость. Образуем далее сумму . Прямая x = at + b будет соответствовать экспериментальным точкам наилучшим образом, если мы найдём такие значения a и b, при которых достигается минимум суммы S. Условия минимума имеют вид , , что даёт систему уравнений:

 

и

 

Система, может быть переписана в виде:

 

 

Подставляя численные значения и решая систему, получим a = 10,0 см/с,
b = 1,0 см, так что искомая зависимость x(t) имеет вид

 

(37)

 

Существуют формулы, позволяющие определить погрешности значений a и b, однако мы их приводить не будем, а воспользуемся приближённым методом: относительная погрешность измерения координаты больше всего для наименьшего значения xmin = 1,0 см и составляет 10% и меньше всего для максимального значения xmax = 11,0 см (примерно 1%). Поэтому значения , вычисляемые по формуле (36) должны получаться с такой же относительной погрешностью. Именно поэтому мы имеем a = (10,0 0,5) см/с, b = (1,00 0,05) см. Приводим на рисунке 7 наилучшую прямую.

Для нахождения вида функциональной зависимости y = y(t) поступим аналогично, изобразив данные из таблицы 1 на координатной плоскости (y,t) (рис. 8).

 

Рисунок 8 – Зависимость y = y(t)

Из рисунка 8 не вытекает, однако, с определённостью предположение о виде зависимости y = y(t). В таких случаях обычно выдвигаются гипотезы о том, какому классу функций (полиномов, показательных, тригонометрических и т. д.) принадлежит искомая зависимость, а затем эти гипотезы принимаются или отвергаются. Чаще всего выдвигается гипотеза о принадлежности неизвестной функции y(t) к классу полиномов некоторой степени n: . Степень полинома обычно берется вначале минимальной, совместимой с характером расположения экспериментальных точек. Из рисунка 8 сразу видно, что зависимость y(t) нелинейная, то есть n ≠ 1. Таким образом, мы берём функцию y(t) = aоt2 + a1t + a2 и ищем значения параметров aо, a1, a2, при которых эта функция наилучшим образом соответствует экспериментальным точкам рисунка 8. Задача решается на основе метода наименьших квадратов. Условия минимума суммы:

 

дают:

 

 

 

Подставляя численные значения и решая систему уравнений, находим после округления aо = 10,0; a1 = 10,0; a2 = 8,0 (количество значащих цифр, в значениях a выбрано исходя из того, что относительная погрешность в определении координаты составляет примерно 1%). Таким образом, зависимость координаты
y = y(t) имеет вид

y = -10t2 + 10t + 8,0 (см) (38)

 

На рисунке 8 построена кривая (парабола), соответсвующая уравнению (38). Как видно, кривая достаточно хорошо проходит через экспериментальные точки. Следует, однако, помнить, что предположение о полиноминальной зависимости y = y(t) является лишь гипотезой. Ведь вполне возможно, что функция вида y = abt + c, где постоянные подобраны с помощью метода наименьших квадратов, или полином степени большей 2 значительно лучше соответствуют экспериментальным точкам рисунка 8. Иными словами, возникает вопрос, насколько оправдана гипотеза о полиноминальной зависимости степени 2, то есть насколько функция (38) соответствует экспериментальным точкам.

На первый взгляд, естественным представляется следующий путь. С помощью метода наименьших квадратов определим значения a, b, c для функции вида y = abt + c, при которых она наилучшим образом соответствует экспериментальным точкам, затем для этих значений a, b, вычислим сумму квадратов разностей, фигурирующих в методе наименьших квадратов, и сравним её с суммой для полиноминальной зависимости (38). Естетственнно, что та зависимость, для которой эта сумма меньше, лучше отвечает экспериментальным точкам. Ясно, однако, что этот путь, хотя и возможен, но трудоёмок и малоперспективен, поскольку существует множество функций времени, которые могли бы, в принципе, соответствовать экспериментальным точкам рисунка 6. Например, зависимость y = Asin(Bt + C) + D с надлежаще подобранными константами A, B, C, D. Поэтому вопросы совместимости гипотезы о той или иной зависимости (в нашем случае зависимости 38) с экспериментальными данными решаются с помощью так называемых критериев согласия (другое название – критерии значимости). Одним из наиболее удобных критериев является так называемый «критерий χ2» (читается хи-квадрат) или критерий Пирсона. В методе наименьших квадратов вычисляется величина χ2:

 

, (38)

 

то есть сумма квадратов отклонений экспериментальных значений yi от вычисленных по формуле (38), деленная на квадрат погрешности измерения величины y. В нашем случае χ2 = 3. Найденное значение χ2 должно быть сопоставлено с теорией. Это делается с помощью таблицы распределения χ2, фрагмент которой приведён в таблице 2. В данной таблице n – это число степеней свободы распределения χ2, равной числу измерений минус увеличенное на единицу число параметров, определяемых из эксперимента. В нашем случае число измерений равно 11 и с помощью метода наименьших квадратов было определено 3 параметра, так что n = 11 – (3 + 1) = 7. Число Р в таблице – вероятность, выражаемая в процентах. По найденному значению χ2 = 1,3 и числу степеней свободы n = 7 находим, что P ≈ 98%. Это означает, что если гипотеза о зависимости (38) справедлива, то найденное или большее значение χ2 должно встречаться примерно в 98% случаев. Следовательно, на уровне доверительной вероятности 98% мы подтвердили зависимость (38). Если, например, при тех же условиях
χ2 = 14,1, то это означало бы, что при справедливости гипотезы (38) такие большие отклонения встречались бы лишь в 5% случаев, так что наше найденное значение χ2 = 14,1 свидетельствовало бы о ненадёжности гипотезы, и это заставило бы искать другую зависимость y(t), например, в виде полинома третьей степени и т. д.

Выпишем окончательно найденный кинематический закон движения:

 

x = 10,0t + 1,0 (см) (39)

y = -10t2 + 10t + 8,0 (см)

 

На рис. 8 построена кривая (парабола), соответствующая уравнению (38). Как видно, кривая достаточно хорошо проходит через экспериментальные точки.

 

Задание 2. Найти среднюю скорость частицы в интервале времениt (0,4;0,8),её модуль и углы с осями координат.

 

Имеем: x(0,4) = 5,0 см; x(0,8) = 9,0 см y(0,4) = 10,4 см; y(0,8) = 9,6 см

Тогда ∆x = 4,0 см, ∆y = – 0,8 см, ∆t = 0,4 с

;

; ;

 

Обратите внимание, что для того, чтобы получить правильное значение для , нужно вычислить значение модуля средней скорости более точно, чем записано выше:

 

Задание 3. Найти модуль мгновенной скорости в момент t1 = 0,7 сек и углы с осями координат.

;

; ;

 

Задание 4. Найти ускорение частиц в тот же момент и углы, составляемые вектором ускорения с осями координат.

;

,

;

 

Вектор скорости и ускорения изображены на рисунке 5.

 

Задание 5. Найти тангенциальное и нормальное ускорения в тот же момент времени.

Направлен вектор так же, как и . Изображаем его на рисунке. Вектор может быть найден геометрически: . Этот вектор также показан на рисунке 6.

Задание 6. Найти радиус кривизны траектории в точке, соответствующей тому же моменту времени.

Используя формулу (34), находим:

 

Во многих случаях оказывается полезным приближенный графический способ нахождения радиуса кривизны. Для этого точку на траектории, соответствующего моменту времени t1 = 0,70 c,соединим прямолинейными отрезками с соседними точками, соответсвующими моментам t2 = 0,60 c и t3 = 0,80 c. Из середины этих отрезков восстанавливаем перпендикуляр до их перечисления в точке 0.

Точка 0 примерно совпадает с центром соприкасающейся окружности, соответствующей участку траектории вблизи точки, для которой велось построение. Радиус окружности примерно равен R1.

 

Задание 7. Найти зависимость пройденного пути S от времени t, то есть функцию S = S(t).

Имеем:

Проведите вычисление самостоятельно, используя табличный интеграл

 

Получим:

Задание 8. Написать уравнение траектории точки.

Исключая время t из уравнений (39), имеем:

 

Или y = – 0,1x2 + 1,2x + 6,9 (см). Это и есть уравнение траектории.

 

Задание 9. Найти скорость удаления частицы от начала координат в момент времени t1 = 0,7 c.

Расстояние от частицы до начала координат согласно формуле (39) равно:

 

 

Поэтому:

(См. формулу 36 стр.15)

Теперь достаточно подставить значениеt1 = 0,7 c.

Задание 10. Найти относительную скорость двух частиц в момент времени

t1 = 0,7 c.

Допустим, что спустя время tо = 0,2 c из той же начальной точки с той же начальной скоростью по той же траектории движется другая частица. Требуется найти относительную скорость частиц в момент t1 = 0,7 c. Скорость второй частицы в момент t1 = 0,7 c будет такая же, как скорость первой в момент 0,5 с. Поэтому, так же, как и в задании 3, можно определить скорость частицы в момент 0,5 с.

 

,

 

Тогда Vотнх = 0,0 см/с, Vотну = – 4,0 см/с и скорость первой частицы относительно второй частицы в момент времени 0,7 секунд равна 4 см/с и направлена вниз.

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

 

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

1. Какие ошибки (пункт а)) имели место при выполнении работы и как они учитывались?

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

3. Как можно проверить отсутствие промахов в серии наблюдений?

4. Изложите методику расчёта погрешностей при измерении объёма цилиндра штангенциркулем.

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

6. Запишите выражение для векторов скорости и нормального ускорения в указанный преподавателем момент времени и проверьте выполнение условий .

7. Нанесите экспериментальные точки и постройте теоретическую кривую зависимости от времени той координаты, для которой она нелинейная, откладывая вдоль оси абсцисс значения , а вдоль оси ординат – значения этой координаты. Сделайте выводы.

 

Таблица 2 χ2 распеделения

 

N/P%
0,8 0.7 1,1 1,6 2,2 3,4 6,0 9,5
0,6 1,1 1,6 2,3 3,0 4,4 7,3 11,1
0,9 1,6 2,2 3,1 3,8 5,3 8,6 12,6
1,3 2,2 2,8 3,8 4,7 6,3 9,8 14,1
1,6 2,7 3,5 4,6 5,5 7,3 11,0 15,5
2,1 3,3 4,2 5,4 6,4 8,3 12,2 16,9
2,6 3,9 4,9 6,2 7,3 9,3 13,4 18,3
3,1 4,6 5,6 7,0 8,3 10,3 14,6 19,7
3,6 5,2 6,3 7,8 9,0 11.3 15,8 21,0
4,1 5,9 7,0 8,6 9,9 12,3 17,0 22,4