EduKG
← Назад

Проверка упрощённой модели движения: от гипотезы к данным (необязательное занятие)

Цели занятия

  • Итог

    Записать упрощённую модель движения в виде xmodel(t)=x0+v0t+12at2x_{\text{model}}(t)=x_0+v_0t+\frac{1}{2}at^2 и перечислить не менее двух её допущений.

  • Итог

    Вычислить остатки ri=xi−xmodel(ti)r_i=x_i-x_{\text{model}}(t_i) по таблице данных и сравнить их с заявленной погрешностью измерения σx\sigma_x.

  • Итог

    По графику остатков определить, случайно расхождение или систематично: например, растёт как tt или как t3t^3.

  • Итог

    Принять, уточнить или отклонить модель, назвав нарушенное допущение и предложив проверяемое уточнение.

Понять идею

Проверка упрощённой модели движения: от гипотезы к данным

Текст

Упрощённая модель движения заменяет тело материальной точкой и предсказывает координату по формуле

xmodel(t)=x0+v0t+12at2,x_{\text{model}}(t)=x_0+v_0t+\frac{1}{2}at^2,

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

Данные видеотрекинга или датчика дают пары (ti,xi)(t_i,x_i) с погрешностью σx\sigma_x, поэтому проверять модель «на глаз» нельзя — нужно вычислять остатки ri=xi−xmodel(ti)r_i=x_i-x_{\text{model}}(t_i).

Если ∣ri∣|r_i| сравнимы с σx\sigma_x и не имеют тренда, модель адекватна в пределах разрешения измерений.

Систематический тренд указывает на нарушенное допущение: остаток, растущий как tt, — это сдвиг начала отсчёта времени, а растущий как t3t^3 — слабое линейное сопротивление воздуха.

Проверка заканчивается решением: принять модель, уточнить её или отклонить, и формулировкой следующей измеримой гипотезы.

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

Напрямую

Задача

Модель постоянной скорости задана как xmodel(t)=0.50+1.20 tx_{\text{model}}(t)=0.50+1.20\,t (единицы СИ). Видеотрекинг дал координаты xx = 0.52, 1.68, 2.82, 3.94 м в моменты tt = 0, 1, 2, 3 с. Погрешность измерения σx=0.03\sigma_x=0.03 м. Вычислите остатки и решите, согласуется ли модель с данными.

Решение

Прогноз модели: xmodelx_{\text{model}} = 0.50, 1.70, 2.90, 4.10 м при tt = 0, 1, 2, 3 с.

Остатки ri=xi−xmodel(ti)r_i=x_i-x_{\text{model}}(t_i): +0.02+0.02, −0.02-0.02, −0.08-0.08, −0.16-0.16 м.

Сравнение с погрешностью: ∣r1∣,∣r2∣≤σx|r_1|,|r_2|\le \sigma_x, но ∣r3∣≈2.7 σx|r_3|\approx 2.7\,\sigma_x и ∣r4∣≈5.3 σx|r_4|\approx 5.3\,\sigma_x.

Остатки не случайны: их модуль монотонно растёт, и они отрицательны на больших временах, значит реальная скорость падает (например, из-за трения).

Вывод: модель постоянной скорости отклоняется; данные требуют модели с замедлением.

Применение

Задача

Тележка скатывается по наклонному жёлобу. Данные: tt = 0, 0.2, 0.4, 0.6, 0.8 с; xx = 0.000, 0.038, 0.150, 0.336, 0.596 м; σx=0.005\sigma_x=0.005 м. Оцените ускорение по вторым разностям, постройте модель x=12at2x=\frac{1}{2}at^2 (при x0=0x_0=0, v0≈0v_0\approx 0), вычислите остатки и проверьте модель.

Решение

Первые разности Δx\Delta x: 0.038, 0.112, 0.186, 0.260 м; вторые разности: 0.074, 0.074, 0.074 м — постоянны, что и требуется при постоянном ускорении.

При Δt=0.2\Delta t=0.2 с: a≈0.074(0.2)2=1.85a\approx\frac{0.074}{(0.2)^2}=1.85 м/с², откуда модель xmodel(t)=0.925 t2x_{\text{model}}(t)=0.925\,t^2.

Прогноз: 0.037, 0.148, 0.333, 0.592 м.

Остатки: +0.001+0.001, +0.002+0.002, +0.003+0.003, +0.004+0.004 м — все меньше σx=0.005\sigma_x=0.005 м.

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

Перенос

Задача

Шар отпускают без начальной скорости, модель xmodel(t)=12gt2x_{\text{model}}(t)=\frac{1}{2}gt^2 с g=9.8g=9.8 м/с². Остатки: rr = −0.001-0.001, −0.008-0.008, −0.027-0.027, −0.064-0.064 м при tt = 0.2, 0.4, 0.6, 0.8 с; σx=0.001\sigma_x=0.001 м. Одногруппник утверждает: «это просто сдвиг начала отсчёта времени». Проверьте это утверждение, выяснив, как остаток масштабируется с tt, t2t^2 и t3t^3; решите, какой механизм подтверждается, и сформулируйте уточнённую гипотезу вместе с одним измерением для её проверки.

Решение

Сдвиг времени δt\delta t даёт r≈−v(t) δt∝tr\approx -v(t)\,\delta t\propto t; постоянная ошибка датчика даёт r=constr=\text{const}; линейное сопротивление воздуха при малых tt даёт r≈−g6τt3r\approx -\frac{g}{6\tau}t^3.

Проверяем масштабирование: r/tr/t = −0.005-0.005, −0.020-0.020, −0.045-0.045, −0.080-0.080 (не постоянно); r/t2r/t^2 = −0.025-0.025, −0.050-0.050, −0.075-0.075, −0.100-0.100 (не постоянно); r/t3r/t^3 = −0.125-0.125, −0.125-0.125, −0.125-0.125, −0.125-0.125 м/с³ (постоянно).

Следовательно, r∝t3r\propto t^3, и версия о сдвиге времени не подтверждается.

Уточнённая гипотеза: сила сопротивления пропорциональна скорости, a=g−v/τa=g-v/\tau; тогда коэффициент −g6τ=−0.125-\frac{g}{6\tau}=-0.125 м/с³ даёт оценку τ≈13\tau\approx 13 с.

Проверка: повторить опыт с шаром той же формы, но большей массы — при линейном сопротивлении τ∝m\tau\propto m, и наклон графика rr от t3t^3 должен уменьшиться; либо удлинить дорожку и убедиться, что отношение r/t3r/t^3 остаётся постоянным.

Практика

Текст

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

  1. Выберите модель: постоянная скорость для игрушечной машинки на гладком полу или равноускоренное движение для шарика, скатывающегося по низкому жёлобу (высота не более 10 см), либо используйте симуляцию (PhET, готовые примеры Tracker).
  2. Снимите видео на телефон под прямым углом к плоскости движения, положите в кадр линейку для масштаба; извлеките пары (ti,xi)(t_i,x_i) с шагом 0.2 с.
  3. В таблице вычислите xmodel(ti)x_{\text{model}}(t_i), остатки ri=xi−xmodel(ti)r_i=x_i-x_{\text{model}}(t_i) и постройте график rr от tt.
  4. Сравните ∣ri∣|r_i| с σx\sigma_x и определите тип тренда: случайный, ∝t\propto t или ∝t3\propto t^3.
  5. Запишите решение: принять модель, уточнить допущение (например, добавить сопротивление) или отклонить, и предложите следующий измеримый тест.

Не бросайте тяжёлые предметы, не используйте лестницы и не проводите измерения рядом с движущимся транспортом.

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

Назначение

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

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

Схема из двух панелей. Сверху — график xx от tt: точки данных и кривая модели xmodel(t)=x0+v0t+12at2x_{\text{model}}(t)=x_0+v_0t+\frac{1}{2}at^2. Снизу — график остатков ri=xi−xmodel(ti)r_i=x_i-x_{\text{model}}(t_i) от tt с полосой ±σx\pm\sigma_x, на котором слева направо показаны три схематических случая: случайное облако (модель адекватна), линейный рост r∝tr\propto t (сдвиг времени) и рост r∝t3r\propto t^3 (сопротивление воздуха). Каждый случай подписан соответствующим механизмом.

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

Спринт на 100 м: проверка модели равномерного движения

Сценарий

Приложение в телефоне фиксирует время бегуна на отметках 20, 40, 60, 80 и 100 м: 2.9 с, 5.0 с, 7.0 с, 9.0 с, 10.5 с. Простая модель считает скорость постоянной, v=10.0 м/сv = 10.0\ \text{м/с}, то есть s=vts = vt, и предсказывает 2.0, 4.0, 6.0, 8.0 и 10.0 с.

Связь

Разности «измерение минус модель» (невязки) равны +0.9+0.9, +1.0+1.0, +1.0+1.0, +1.0+1.0 и +0.5+0.5 с: знак у них одинаковый, а не разбросанный, значит модель систематически ошибается на старте. Это указывает на разгонный участок s=12at2s = \frac12 a t^2: модель уточняют по данным, а не отбрасывают.

Поездка на автомобиле: модель постоянной скорости против данных

Сценарий

Автомобилю нужно проехать 150 км150\ \text{км} по трассе; модель задаёт постоянную скорость 100 км/ч100\ \text{км/ч}, откуда t=s/v=1.5 ч=90 минt = s/v = 1.5\ \text{ч} = 90\ \text{мин}. Часы и одометр дают 118 мин118\ \text{мин}, из них 12 мин12\ \text{мин} — стоянка на заправке; в движении спидометр колеблется от 9090 до 110 км/ч110\ \text{км/ч}.

Связь

Невязка +28 мин+28\ \text{мин} распадается на слагаемое, которого в модели нет (12 мин12\ \text{мин} стоянки), и остаток 16 мин16\ \text{мин}, объясняемый разбросом скорости вокруг 100 км/ч100\ \text{км/ч}. Проверка модели отделяет пропущенный член от случайного разброса.

Мяч и телефон: проверка модели свободного падения замедленной съёмкой

Сценарий

Мяч роняют с высоты h=1.80 мh = 1.80\ \text{м}, телефон снимает падение с частотой 240 кадр/с240\ \text{кадр/с}, поэтому соседние кадры разделены 1/240 с≈4.2 мс1/240\ \text{с} \approx 4.2\ \text{мс}. Модель свободного падения без сопротивления воздуха h=12gt2h = \frac12 g t^2 даёт t=2h/g=2⋅1.80/9.81=0.606 сt = \sqrt{2h/g} = \sqrt{2\cdot 1.80/9.81} = 0.606\ \text{с}. По кадрам получается 0.63±0.01 с0.63 \pm 0.01\ \text{с}.

Связь

Невязка 0.024 с0.024\ \text{с} (около 4%4\%) больше интервала между кадрами и погрешности момента отпускания, значит модель без сопротивления воздуха проверку не проходит и требует члена сопротивления (∼v2\sim v^2). Именно здесь работают оценка движения по видео и её временное разрешение.