Формула для периода малых колебаний математического маятника имеет вид $
T=2 \pi \sqrt{\cfrac{l}{g}}.
$
Из уравнения закона сохранения механической энергии для тела маятника (нулевой уровень потенциальной энергии принят на уровне точки подвеса)
\[
\frac{m l^{2} \omega^{2}}{2}=m g l\left(\cos \varphi-\cos \varphi_{0}\right)
\]
следует формула для угловой скорости движения маятника
\[
\omega=\sqrt{\frac{2g}{l}\left(\cos \varphi-\cos \varphi_{0}\right)}.
\]
Разобьем весь участок движения от \( \varphi_{0} \) до нуля на бесконечно малые интервалы \( d \varphi . \) Время \( d t, \) за которое маятник пройдет этот интервал, равно
\[
d t=\frac{d \varphi}{\omega}=\frac{d \varphi}{\sqrt{\frac{2g}{l}\left(\cos \varphi-\cos \varphi_{0}\right)}}
\]
Тогда время движения \( t_{1} \) определяется как сумма малых интервалов, что сводится к простому интегрированию
\[
t_{1}=\int_{0}^{\varphi_{0}} d t=\int_{0}^{\varphi_{0}} \frac{d \varphi}{\sqrt{\frac{2g}{l}\left(\cos \varphi-\cos \varphi_{0}\right)}}
\]
Период колебаний в 4 раза больше найденного времени \( t_{1}, \) то есть $T=4 t_{1}$.
Угловая скорость в безразмерных единицах выражается следующим образом
\[
\widetilde{\omega}=\frac{d \varphi}{d \tau}=\frac{d \varphi}{\sqrt{\frac{g}{l}} d t}=\sqrt{\frac{l}{g}}\omega
\]
Период малых колебаний в безразмерных единицах равен
\[
\widetilde{T}=\sqrt{\frac{g}{l}} T=2 \pi
\]
Зависимость угловой скорости \( \widetilde{\omega} \) от угла отклонения \( \varphi \) имеет вид
\[
\widetilde{\omega}(\varphi)=\sqrt{2(\cos \varphi-\cos \varphi_{0})}
\]
Интервал разбиения равен $\Delta \varphi=\frac{\varphi_{0}}{N}$.
Сначала следует задать «нулевой» угол отклонения \( \varphi_{0} \), а координаты остальных точек разбиения задаются формулой
\[
\varphi_{k}=\varphi_{k-1}-\Delta \varphi
\]
Угловая скорость \( \omega_{k} \) в точке \( \varphi_{k} \) рассчитывается по формуле
\[
\omega_{k}=\sqrt{2(\cos \varphi_{k}-\cos \varphi_{0})}
\]
В данном случае (при \( \varphi_{0}=\cfrac{\pi}{2} \) ) эта формула еще больше упрощается
$$
\omega_{k}=\sqrt{2\cos \varphi_{k}}
$$
В рекомендуемом приближении равноускоренного движения средняя скорость на выбранном интервале равна среднему арифметическому скоростей на концах интервала
\[
\langle\omega\rangle=\frac{1}{2}\left(\omega_{k-1}+\omega_{k}\right)
\]
Поэтому время \( \Delta t_{k} \) прохождения \( k \) -го интервала от \( \varphi_{k-1} \) до \( \varphi_{k} \) рассчитывается по формуле
\[
\Delta t_{k}=\frac{2 \Delta \varphi}{\omega_{k-1}+\omega_{k}}
\]
Время \( t_{k} \) прохождения точки \( \varphi_{k} \) находится по формуле
\[
t_{k}=t_{k-1}+\Delta t_{k}
\]
при начальном условии \( t_{0}=0 \).
Период колебаний \( T_{N} \) при разбиении на \( N \) интервалов рассчитывается по формуле $
T_{N}=4 t_{N}
$.
Результаты расчетов угловых скоростей, времен и периодов колебаний при указанных значениях числа интервалов разбиения приведены в Таблице 1
Для расчета графика закона движения в приближении малых колебаний необходимо
воспользоваться формулой
\[
\varphi(t)=\frac{\pi}{2} \cos t
\]
Результаты расчетов по этому закону представлены в таблице 2 и на графике. Там же приведен график рассчитанного закона
движения. Интересно отметить, что в первом случае задаются значения времен, а рассчитываются соответствующие углы отклонения; а во втором, наоборот: задаются углы отклонения и проводится расчет времен.
См. график 3.2
Результаты расчетов погрешностей $\varepsilon_{N}$ периодов колебаний при различном числе интервалов разбиения \( N \) приведены в Таблице 3.
Для определения параметров зависимости \( \varepsilon_{N}=\cfrac{C}{N^{\gamma}}, \) ее необходимо представить в двойном логарифмическом масштабе \( \ln \varepsilon_{N}=\ln C-\gamma \ln N \).
На рисунке показан график этой зависимости, построенный по данным Таблицы 3. Линейность полученной зависимости показывает применимость формулы для зависимости погрешности от числа \( N \). Параметры этой линейной зависимости, рассчитанные по методу наименьших квадратов равны: коэффициент наклона \( a \approx-2,3 \) и сдвиг \( b \approx-1,65 \). Следовательно, искомые параметры зависимости равны
$$
\gamma \approx-a=2.3 \\
C=\exp (b) \approx 0.19
$$
Из формулы в п. 3.4 несложно найти число интервалов разбиений, необходимых для достижения
погрешности \( \varepsilon=0.002 \)
\[
N_{\min }=\left(\frac{C}{\varepsilon}\right)^{1 / \gamma} \approx 8 .
\]
Итак, все дальнейшие расчеты можно проводить при $ N=N_\min=8 $.
Результаты расчетов периодов колебаний для различных амплитуд колебаний показаны в Таблице 4.
При малых колебаниях формула для периода колебаний
\[
T\left(\varphi_{0}\right)=T_{0}\left(a+\frac{\varphi_{0}^{2}}{b}\right)
\]
должна совпадать с формулой для периода колебаний математического маятника, откуда следует, что параметр \( a=1 . \) Для проверки применимости формулы выше к описанию результатов расчетов необходимо построить график зависимости величины \( \left(\frac{T}{T_{0}}-1\right) \) от квадрата амплитуды \( \varphi_{0}^{2}, \) который показан на рисунке.
Коэффициент наклона полученной зависимости равен \( 0,0706, \) поэтому параметр \( b, \) входящий в формулу п. 4.2, примерно равен \( b \approx 14 \).
При допустимой погрешности реального эксперимента в \( 5 \% \) отклонение величины периода от периода малых колебаний не будет заметно, если выполняется неравенство
\[
\frac{\varphi_{0}^{2}}{b}<0.05
\]
откуда следует, что «малыми» углами могут считаться углы, меньшие, чем \( \varphi_{0}<45^{\circ} \).