EduKG
← Назад

Оценка движения по видео и датчикам: от кадров к модели

Цели занятия

  • Итог

    Извлекать координаты объекта из последовательных кадров видео или отсчётов датчика и строить таблицу x(t)x(t) с известным шагом Δt\Delta t.

  • Итог

    Вычислять среднюю и мгновенную скорость по разностям координат и строить график v(t)v(t).

  • Итог

    Определять характер движения (равномерное или равноускоренное), проверяя постоянство первых и вторых разностей координат.

  • Итог

    Сопоставлять измеренное ускорение со свободной схемой сил и оценивать доминирующую силу через ma⃗=∑F⃗m\vec a = \sum \vec F.

Понять идею

Оценка движения по видео и датчикам: от кадров к модели

Текст

Видео и датчики дают не само движение, а набор чисел: кадры с интервалом Δt=1/f\Delta t = 1/f или отсчёты акселерометра с частотой fsf_s. Первый шаг — задать масштаб по известному объекту (линейке, мячу известного диаметра) и превратить пиксели в метры, получив таблицу координат x(t)x(t). Скорость оценивают разностями между соседними кадрами: vi+1/2=(xi+1−xi)/Δtv_{i+1/2} = (x_{i+1}-x_i)/\Delta t, поэтому точка графика v(t)v(t) лежит «между» кадрами. Ускорение — это наклон графика v(t)v(t) или вторая разность a≈(xi+1−2xi+xi−1)/Δt2a \approx (x_{i+1}-2x_i+x_{i-1})/\Delta t^{2}. Реальные данные шумны, и шаг Δt\Delta t конечен, так что полезно оценивать погрешность δv≈δx/Δt\delta v \approx \delta x/\Delta t и проверять, действительно ли вторые разности постоянны. Кинематика становится физической моделью только тогда, когда найденное ускорение согласуется со свободной схемой сил и вторым законом Ньютона ma⃗=∑F⃗m\vec a = \sum \vec F.

Разобранные примеры

Напрямую

Задача

Видео снято с частотой 25 кадров/с, масштаб задан линейкой: 1 м=250 px1\ \text{м} = 250\ \text{px}. Мяч катится по полу, и в кадрах 0, 1, 2, 3 его координаты равны 0, 25, 50, 75 px. Найдите скорость мяча и её погрешность, если координата измеряется с точностью δx=1 px\delta x = 1\ \text{px}.

Решение

Шаг 1. Интервал между кадрами: Δt=1/f=1/25=0.04 с\Delta t = 1/f = 1/25 = 0.04\ \text{с}.

Шаг 2. Переводим пиксели в метры, деля на 250:

кадр 0 1 2 3
xx, м 0 0.10 0.20 0.30

Шаг 3. Разности координат: Δx=0.10 м\Delta x = 0.10\ \text{м} на каждом шаге. Скорость:

v=ΔxΔt=0.100.04=2.5 м/с.v = \frac{\Delta x}{\Delta t} = \frac{0.10}{0.04} = 2.5\ \text{м/с}.

Шаг 4. Все разности одинаковы, значит движение равномерное и v(t)v(t) — горизонтальная прямая v=2.5 м/сv = 2.5\ \text{м/с}.

Шаг 5. Погрешность: δv≈δx/Δt=(1/250)/0.04=0.1 м/с\delta v \approx \delta x/\Delta t = (1/250)/0.04 = 0.1\ \text{м/с}. Ответ: v=2.5±0.1 м/сv = 2.5 \pm 0.1\ \text{м/с}.

Применение

Задача

Объект падает, видео снято с частотой 20 кадров/с, поэтому Δt=0.05 с\Delta t = 0.05\ \text{с}. Масштабированные координаты (м) по кадрам: 0, 0.012, 0.049, 0.110, 0.196. Определите характер движения, найдите ускорение и запишите модель x(t)x(t). Согласуется ли результат со свободной схемой сил, если сопротивлением воздуха можно пренебречь?

Решение

Шаг 1. Первые разности Δx=xi+1−xi\Delta x = x_{i+1}-x_i: 0.012; 0.037; 0.061; 0.086 м. Они растут, значит движение не равномерное.

Шаг 2. Вторые разности Δ2x=xi+1−2xi+xi−1\Delta^2 x = x_{i+1}-2x_i+x_{i-1}: 0.025; 0.024; 0.025 м. Они постоянны (в пределах шума ±0.001\pm 0.001 м), значит ускорение постоянно:

a≈Δ2xΔt2=0.0247(0.05)2≈9.9 м/с2.a \approx \frac{\Delta^2 x}{\Delta t^{2}} = \frac{0.0247}{(0.05)^{2}} \approx 9.9\ \text{м/с}^{2}.

Шаг 3. Средняя скорость за первый интервал равна v1/2≈0.012/0.05=0.24 м/сv_{1/2} \approx 0.012/0.05 = 0.24\ \text{м/с} и относится к середине интервала, t=0.025 сt = 0.025\ \text{с}. Поэтому начальная скорость при t=0t=0 равна v0=v1/2−aΔt/2=0.24−9.9⋅0.025≈−0.008 м/сv_0 = v_{1/2} - a\Delta t/2 = 0.24 - 9.9\cdot 0.025 \approx -0.008\ \text{м/с}, то есть в пределах погрешности v0≈0v_0 \approx 0.

Шаг 4. Модель с x0=0x_0 = 0 и v0≈0v_0 \approx 0: x(t)≈12(9.9)t2x(t) \approx \tfrac{1}{2}(9.9)t^{2} м. Она воспроизводит таблицу в пределах погрешности. Модель x(t)=0.24t+12(9.9)t2x(t)=0.24t+\tfrac{1}{2}(9.9)t^{2} не подходит: она дала бы x(0.05)=0.0244 мx(0.05)=0.0244\ \text{м}.

Шаг 5. Свободная схема сил: на тело действует только сила тяжести mg⃗m\vec g, поэтому a=g≈9.8 м/с2a = g \approx 9.8\ \text{м/с}^{2}. Измеренное 9.9±0.4 м/с29.9 \pm 0.4\ \text{м/с}^{2} совпадает с gg — модель «свободное падение» подтверждена.

Перенос

Задача

Датчик телефона, закреплённого на тележке, даёт вдоль наклонной плоскости ускорение a=2.6 м/с2a = 2.6\ \text{м/с}^{2} при угле наклона θ=20∘\theta = 20^{\circ}. Оцените коэффициент трения μ\mu между тележкой и плоскостью и укажите, какие ещё проверки нужны, чтобы модели можно было доверять.

Решение

Шаг 1. Свободная схема сил: сила тяжести mg⃗m\vec g, нормальная реакция N⃗\vec N, сила трения F⃗тр\vec F_{тр}. Вдоль плоскости (ось xx вниз по наклону):

ma=mgsin⁡θ−Fтр.ma = mg\sin\theta - F_{тр}.

Шаг 2. Перпендикулярно плоскости: N=mgcos⁡θN = mg\cos\theta, а для скольжения Fтр=μNF_{тр} = \mu N. Подставляем:

a=g(sin⁡θ−μcos⁡θ).a = g(\sin\theta - \mu\cos\theta).

Шаг 3. Выражаем μ\mu и считаем:

μ=gsin⁡θ−agcos⁡θ=9.8⋅0.342−2.69.8⋅0.940=0.7529.21≈0.08.\mu = \frac{g\sin\theta - a}{g\cos\theta} = \frac{9.8\cdot 0.342 - 2.6}{9.8\cdot 0.940} = \frac{0.752}{9.21} \approx 0.08.

Шаг 4. Проверки модели: (а) вторые разности скорости по данным датчика должны быть постоянны — тогда aa действительно постоянна; (б) нужно оценить погрешность δa\delta a и δθ\delta\theta, потому что при малых μ\mu разность gsin⁡θ−ag\sin\theta - a мала и результат чувствителен к ошибкам; (в) полезно повторить опыт при другом θ\theta: если μ\mu остаётся тем же, модель трения Fтр=μNF_{тр} = \mu N разумна, если нет — нужна более сложная модель (например, зависимость от скорости).

Практика

Текст

План самостоятельной работы (безопасно, без падений с высоты).

  1. Съёмка. Положите телефон на штатив или на стопку книг так, чтобы он не упал, и снимите замедленное видео (120–240 кадров/с) катящегося по полу мяча. Рядом положите линейку. Бросать предметы и снимать падение с высоты не нужно: свободное падение разберём в симуляции.
  2. Масштаб. В одном кадре измерьте длину линейки в пикселях и получите коэффициент kk (м/пиксель).
  3. Таблица и график. Отметьте положение мяча в 6–10 последовательных кадрах, переведите в метры и постройте графики x(t)x(t) и v(t)v(t) (скорость — по разностям).
  4. Симуляция. Откройте симуляцию движения с построением графиков (например, «Motion» / «Движение с постоянным ускорением»), задайте a=9.8 м/с2a = 9.8\ \text{м/с}^{2} и сравните форму графика с вашими данными.
  5. Вывод. Запишите, какая модель (x=x0+v0tx = x_0 + v_0 t или x=x0+v0t+12at2x = x_0 + v_0 t + \tfrac{1}{2}at^{2}) описывает ваши данные, и какая погрешность δv\delta v ожидается при вашем Δt\Delta t.

Визуальная подсказка

Назначение

Показать переход от дискретных кадров к непрерывной модели: как из отметок положения на кадрах получаются графики x(t)x(t) и v(t)v(t) и откуда берётся ускорение.

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

Схема + графики. Сверху — ряд кадров (прямоугольники) с отметками времени t0,t1,t2,…t_0, t_1, t_2, \dots и отмеченным положением мяча; под ними — ось xx с теми же точками и штриховкой Δx\Delta x. Слева внизу — график x(t)x(t): точки, соединённые линией, с касательной, наклон которой даёт скорость. Справа внизу — график v(t)v(t): точки, лежащие между кадрами, и прямая v=v0+atv = v_0 + at, наклон которой равен ускорению aa.

Примеры из жизни

Смартфон на 240 кадрах в секунду: скорость бегуна

Сценарий

На школьной дорожке бегун пробегает мимо отметки 60 м. Смартфон снимает его с частотой 240 кадров/с; интервал между кадрами Δt=1/240 s≈4.17 ms\Delta t = 1/240\,\text{s} \approx 4.17\,\text{ms}. На двух соседних кадрах бедро бегуна смещается на 0.034 m0.034\,\text{m} вдоль дорожки. Масштаб задан линией разметки длиной 1.00 m1.00\,\text{m}, которая занимает 100 пикселей.

Связь

Здесь работает переход от кадров к таблице «положение–время»: калибровка пикселей в метры, привязка к временным меткам и численная производная. Средняя скорость на интервале равна v≈Δx/Δt=0.034 m/0.00417 s≈8.2 m/sv \approx \Delta x/\Delta t = 0.034\,\text{m}/0.00417\,\text{s} \approx 8.2\,\text{m/s}, то есть наклон графика x(t)x(t) между двумя кадрами.

Акселерометр телефона при торможении автомобиля

Сценарий

Телефон закреплён на приборной панели автомобиля массой 1400 kg1400\,\text{kg}. При торможении с 50 km/h=13.9 m/s50\,\text{km/h}=13.9\,\text{m/s} до нуля за 3.2 s3.2\,\text{s} встроенный акселерометр записывает данные с частотой 50 Hz50\,\text{Hz}, то есть 160 отсчётов. Среднее ускорение вдоль движения a≈−13.9/3.2≈−4.3 m/s2a \approx -13.9/3.2 \approx -4.3\,\text{m/s}^2.

Связь

Оценка движения по датчику требует численного интегрирования ускорения: Δv=a Δt\Delta v = a\,\Delta t для каждого из 160 интервалов, затем сложение. Схема сил автомобиля (сила тяжести, нормальная реакция, сила трения торможения) показывает, что горизонтальная сила трения создаёт это ускорение; при m=1400 kgm=1400\,\text{kg} средняя тормозная сила около F=ma≈6.0×103 NF=ma \approx 6.0\times10^3\,\text{N}.

Датчик на пешеходном мосту: колебания после грузовика

Сценарий

На пешеходном мосту акселерометр измеряет вертикальное ускорение с частотой 100 Hz100\,\text{Hz} в течение 50 s50\,\text{s} после проезда грузовика. Пиковое ускорение 0.80 m/s20.80\,\text{m/s}^2, доминирующая частота 2.5 Hz2.5\,\text{Hz}. Для синусоидального режима амплитуда смещения A=amax/(2πf)2≈0.80/(2π⋅2.5)2 m≈3.2 mmA = a_{\text{max}}/(2\pi f)^2 \approx 0.80/(2\pi\cdot2.5)^2\,\text{m} \approx 3.2\,\text{mm}.

Связь

Это пример обработки данных датчика: удаление постоянного смещения и дрейфа, переход от ускорения к смещению двойным численным интегрированием. Полученный график y(t)y(t) показывает колебание моста; по нему можно прочитать амплитуду и период.