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

В физике и инженерии графический анализ колебательных процессов является одним из фундаментальных методов определения параметров системы. В контексте математического маятника под «поиском подклада на графике» чаще всего понимают определение положения статического равновесия (нулевой линии, базовой линии или «подклада» сигнала), а также вычисление таких характеристик, как период, амплитуда, логарифмический декремент затухания и длина маятника. Данная статья подробно рассматривает теоретические основы, методы компьютерного моделирования и алгоритмы обработки экспериментальных данных для точного извлечения этих параметров из временных диаграмм φ(t) или ω(t).

Теоретическая модель математического маятника

Рассмотрим классическую модель: материальная точка (шарик) массы m, подвешенная на невесомой нерастяжимой нити длины l в однородном поле тяжести с ускорением g. Сопротивление среды в первой приближении пренебрегаем.

Уравнение движения и интеграл энергии

Движение маятника описывается вторым законом Ньютона для вращательного движения:

I ε = M_гр → m l² (d²φ/dt²) = ー m g l sin φ

Отсюда получаем дифференциальное уравнение:

d²φ/dt² + (g/l) sin φ = 0

Используя закон сохранения механической энергии (кинетическая + потенциальная = const), можно получить точную формулу для угловой скорости ω = dφ/dt в зависимости от угла отклонения φ при заданной амплитуде φ₀ (начальная скорость равна нулю):

ω(φ) = ± √[ (2g/l) (cos φ ー cos φ₀) ]

Используя тригонометрические тождества (cos φ = 1 ー 2 sin²(φ/2)), формулу можно привести к виду, удобному для интегрирования:

ω(φ) = ± 2 √(g/l) √[ sin²(φ₀/2) ౼ sin²(φ/2) ]

Точная формула периода колебаний

Период T вычисляется как четвертую часть времени полного цикла (от 0 до φ₀ и обратно), умноженное на 4. Интеграл приводит к полному эллиптическому интегралу первого рода K(k):

T(φ₀) = 4 √(l/g) K(k), где k = sin(φ₀/2)

Разложение K(k) в ряд дает известную формулу зависимости периода от амплитуды:

T ≈ T₀ [ 1 + (1/4) sin²(φ₀/2) + (9/64) sin⁴(φ₀/2) + … ]

где T₀ = 2π √(l/g), период малых колебаний (φ₀ → 0).

Важно: Формула T₀ = 2π √(l/g) применима только при малых углах (обычно φ₀ < 5°...10°). Задача компьютерного эксперимента — определить границу применимости этой формулы, анализируя график T(φ₀).

Компьютерный эксперимент: Постановка и алгоритм

Схема проведения числового эксперимента мало отличается от натурного. Основные этапы:

1. Построение модели: Численное интегрирование уравнения движения (метод Рунге-Кутты 4-го порядка, Верле или Эйлера-Кромера) с заданными l, g, φ₀.
2. Генерация данных: Получение массивов значений времени tᵢ, угла φᵢ и угловой скорости ωᵢ.
3. Построение графиков: φ(t), ω(t), фазовая траектория ω(φ).
4. Обработка графика (Поиск подклада и параметров): Анализ временного ряда для извлечения физических констант.

Что такое «Подклад» на графике φ(t)?

В экспериментальных данных (реальных или смоделированных с шумом) график угла φ(t) часто смещен относительно нуля из-за:

• Ошибки калибровки датчика угла (нулевой сдвиг).
• Невыверенности установки (вертикаль не совпадает с нулем датчика).
• Наличия постоянной внешней силы (например, электрического поля или ветра).

Поиск подклада — это процедура определения константы C (смещения базовой линии), такой что φ_ист(t) = φ_изм(t) ー C симметричен относительно нуля (среднее значение за период равно нулю).

Алгоритм поиска подклада (нулевой линии):

1. Метод среднего значения: Для затухающих колебаний подклад равен среднему значению огибающей (среднее между максимумами и минимумами). C = (φ_max + φ_min) / 2.
2. Метод наименьших квадратов (МНК): Аппроксимация данных функцией φ(t) = A e^-βt cos(ωt + φ₀) + C. Параметр C и есть искомый подклад.
3. Фазовый портрет: На графике ω(φ) траектория должна быть симметрична относительно оси ω (φ=0). Смещение центра эллипса/спирали по оси φ дает значение подклада.

Определение параметров маятника по графику

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

Период T и частота ω₀

• По переходам через ноль: Измеряют время между соседними переходами графика φ(t) через нуль в одну сторону (или между пиками одного знака). T = t_i+n ౼ t_i (усредняют по n периодам для точности).
• Спектральный анализ (БПФ/FFT): Быстрое преобразование Фурье временного ряда φ(t). Пик спектра дает точную частоту f₀ = 1/T. Этот метод устойчив к шуму и позволяет выделить гармоники (важно для больших амплитуд, когда движение несинусоидально).

Длина маятника l и ускорение g

Зная точный период T (для малых амплитуд) или используя точную формулу с эллиптическими интеграми (для больших амплитуд), вычисляют:

l = g (T / 2π)² (для малых углов)

Если l известно точно, определяют g (гравиметрия).

Логарифмический декремент затухания λ и коэффициент сопротивления

Если в модели учтено сопротивление (сила F_сопр = -b v), амплитуда затухает по экспоненте: A(t) = A₀ e^-βt.

• Строят огибающую пиков амплитуды (линии, соединяющие максимумы).
• В полулогарифмических координатах (ln A от t) график становится прямой линией.
• Наклон прямой дает β, откуда находят коэффициент затухания γ = 2β и качество Q = ω₀ / γ.

Анализ зависимости периода от амплитуды T(φ₀)

Это центральная задача предложенной лабораторной работы. Результаты компьютерного эксперимента оформляют в виде таблицы и графика.

Порядок проведения серии измерений

График T(φ₀) и критерий «малого угла»

На графике (ось X, φ₀ в градусах или радианах, ось Y — T или T/T₀) видна нелинейность.

• Строят теоретическую кривую T_теор(φ₀) = 4√(l/g) K(sin(φ₀/2)).
• Наносят экспериментальные точки с погрешностями.
• Определяют предельную амплитуду φ_гр, при которой отклонение от формулы малых колебаний превышает допустимую погрешность прибора (например, 0.5% или 1%).

Фазовый портрет: ω(φ) как инструмент диагностики

График угловой скорости от угла (фазовая траектория) дает исчерпывающую информацию о системе без привязки ко времени.

• Консервативная система (без трения): Замкнутая кривая, симметричная относительно осей. Площадь, ограниченная кривой, пропорциональна полной энергии E = (1/2) m l² ω² + m g l (1 ౼ cos φ).
• Диссипативная система (с трением): Спираль, сворачивающаяся в точку равновесия (0, 0).
• Поиск подклада на фазовом портрете: Если центр спирали смещен по оси φ на значение C, то C и есть искомый подклад (смещение равновесия).
• Нелинейность: Форма кривой отклоняется от эллипса (который соответствует линейному осциллятору) при больших амплитудах, становясь «свернутым» по вертикали (угловая скорость بالقرب экстремумов угла ниже, чем предсказывает линейная теория).

Практическое задание: Обработка данных в Python (пример)

Ниже приведен фрагмент кода для численного моделирования и поиска параметров (включая подклад) на графике.

import numpy as np
from scipy.integrate import solve_ivp
from scipy.signal import find_peaks
from scipy.optimize import curve_fit
import matplotlib.pyplot as plt

g = 9.81
l = 1.0
phi0 = np.deg2rad(30) # 30 градусов
omega0 = 0.0
beta = 0.05 # Коэффициент затухания (для реализма)
C_offset = 0.02 # Искусственный "подклад" (смещение нуля датчика)

def pendulum_eq(t, y):
 phi, omega = y
 dphi = omega
 domega = -2betaomega ー (g/l)*np.sin(phi)
 return [dphi, domega]

sol = solve_ivp(pendulum_eq, [0, 40], [phi0, omega0], max_step=0.01, dense_output=True)
t = np.linspace(0, 40, 4000)
phi_true = sol.sol(t)[0]
phi_measured = phi_true + C_offset # Данные с "подкладом"

peaks, _ = find_peaks(np.abs(phi_measured), distance=300)
t_peaks = t[peaks]
phi_peaks = phi_measured[peaks]

def envelope(t, A, beta, C):
 return A * np.exp(-beta * t) + C

popt, _ = curve_fit(envelope, t_peaks, phi_peaks, p0=[phi0, 0.05, 0])
C_found = popt[2]
print(f"Истинный подклад: {C_offset:.4f}, Найденный подклад: {C_found:.4f}")

phi_corrected = phi_measured ౼ C_found

zero_crossings = np.where(np.diff(np.sign(phi_corrected)))[0]
if len(zero_crossings) > 1:
 T_measured = 2 * (t[zero_crossings[2]] ー t[zero_crossings[0]]) # 2 полупериода
 print(f"Измеренный период: {T_measured:.3f} с")
 print(f"Теория (малые): {2np.pinp.sqrt(l/g):.3f} с")

fig, axes = plt.subplots(2, 2, figsize=(10, 8))
axes[0,0].plot(t, phi_measured, label='Сырые данные (с подкладом)')
axes[0,0].axhline(C_found, color='r', linestyle='--', label=f'Найденный подклад={C_found:.3f}')
axes[0,0].legend; axes[0,0].set_ylabel('φ, рад'); axes[0,0].grid

axes[0,1].plot(t, phi_corrected, label='Скорректированные данные', color='green')
axes[0,1].legend; axes[0,1].set_ylabel('φ, рад'); axes[0,1].grid

axes[1,0].plot(phi_corrected, sol.sol(t)[1], color='purple')
axes[1,0].set_xlabel('φ, рад'); axes[1,0].set_ylabel('ω, рад/с')
axes[1,0].set_title('Фазовый портрет (ω vs φ)'); axes[1,0].grid

axes[1,1].plot(t_peaks, phi_peaks, 'o', label='Пики амплитуды')
axes[1,1].plot(t, envelope(t, *popt), 'r-', label='Аппроксимация огибающей')
axes[1,1].legend; axes[1,1].set_xlabel('t, с'); axes[1,1].set_ylabel('Амплитуда'); axes[1,1].grid

plt.tight_layout
plt.show

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

Ключевые итоги:

• Поиск подклада (смещения нуля) является обязательным предварительным этапом обработки данных. Наиболее надежно он находится через аппроксимацию огибающей амплитуд (МНК) или анализ фазового портрета.
• Зависимость периода от амплитуды T(φ₀) отлично описывается полным эллиптическим интегралом первого рода. Отклонение от формулы T₀ = 2π√(l/g) становится заметным уже при φ₀ > 5°–10°.
• Компьютерный эксперимент (численное интегрирование уравнений движения) позволяет смоделировать реальные условия (затухание, шумы, подклад) и отработать алгоритмы обработки данных перед натурным опытом.
• График в фазовых координатах ω(φ) служит наглядным индикатором сохранения энергии (замкнутые кривые) или диссипации (спирали), а также мгновенно показывает наличие систематического смещения (подклада) по оси абсцисс.

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