2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




Начать новую тему Ответить на тему На страницу 1, 2, 3  След.
 
 Помогите разобраться с выводом формул по колебаниям
Сообщение10.08.2015, 17:54 
Аватара пользователя


13/02/13
777
♍ — ☉ — ⊕
Здравствуйте! Заинтересовался я колебаниями и (как обычно) решил что-нибудь вывести по данной теме. Выбор пал на амплитуду колебаний груза, подвешенного на пружине.
1) Итак, я начинаю с уже готовой формулы (потому что её вывод очевиден): $x''=-{\omega_0}^2x$
2) Решим ДУ, чтобы получить формулу координаты груза.

$x''+{\omega_0}^2x=0$

${\lambda}^2+{\omega_0}^2=0$

$\lambda=\pm\omega_0i$

$x=e^{\alpha t}(C_1\cos\beta t+C_2\sin\beta t)$

$x= C_1 \cos\omega_0 t+C_2\sin\omega_0 t$

3) Теперь надо определить значения констант. Для этого решим задачу Коши для данного ДУ со следующими условиями - $x(0)=x_0$ и $v(0)=v_0\Rightarrow x'(0)=v_0$

$$\left\{
\begin{array}{lcl}
 x=C_1\cos\omega_0 t+C_2\sin\omega_0 t \\
 x'=-\omega_0 C_1 \sin\omega_0 t +\omega_0 C_2 \cos \omega_0 t \\
\end{array}
\right.
\Rightarrow
\left\{
\begin{array}{lcl}
 x_0=C_1 \cos 0+C_2 \sin 0=C_1\\
 v_0=- \omega_0 C_1\sin 0 + \omega_0 C_2 \cos 0=\omega_0 C_2\\
\end{array}
\right.
\Rightarrow
\left\{
\begin{array}{lcl}
 C_1=x_0\\
 C_2=\frac{v_0}{\omega_0}\\
\end{array}
\right.
$$

4) Тогда уравнение координаты запишется так: $x=x_0 \cos \omega_0 t + \frac{v_0}{\omega_0}\sin \omega_0 t$

5) По определению амплитуда - это максимальное смещение от положения равновесия. То есть $A=\max(|x-x_0|)$
Но так как $x_0=\operatorname{const}$, то данная разность окажется максимальной тогда, когда $x$ примет максимальное значение.

6) Теперь требуется найти максимальное значение функции $x=x_0 \cos \omega_0 t + \frac{v_0}{\omega_0}\sin \omega_0 t$
Чтобы его найти, нужно выяснить при каком значении $t$ функция $x(t)$ принимает наибольшее значение. Для этого найдём нули производной функции.
$$t=\frac{1}{\omega_0}\arctg\frac{v_0}{x_0\omega_0}+\frac{\pi}{\omega_0}k \text{ , где } k\in\mathbb{Z}$$
Но легко видеть, что $t(k)$ - функция монотонная. Поэтому не получается найти такое $t$, при котором $x(t)$ примет максимальное значение. Выходит, что колебаний вообще быть не должно.

Скажите, пожалуйста, где и в чём я ошибаюсь.

 Профиль  
                  
 
 Re: Помогите разобраться с выводом формул по колебаниям
Сообщение10.08.2015, 18:22 


19/07/15
74
Atom001 в сообщении #1043981 писал(а):
6) Теперь требуется найти максимальное значение функции $x=x_0 \cos \omega_0 t + \frac{v_0}{\omega_0}\sin \omega_0 t$
Чтобы его найти, нужно выяснить при каком значении $t$ функция $x(t)$ принимает наибольшее значение. Для этого найдём нули производной функции.


Отвечу только по этому пункту.

Можно воспользоваться преобразованием выражения вида $Acos(y) + Bsin(y)$ в выражение вида $Csin(y + \phi)$. Константы $C$ и $\phi$ находятся по формулам $C = \sqrt{A^2 + B^2}, \phi = atan2(A, B)$. Имея в виду такое преобразование, сразу ясно, что
$$|x|_{max} = \sqrt{x_0^2 + \frac{v_0^2}{\omega_0^2}$$

 Профиль  
                  
 
 Re: Помогите разобраться с выводом формул по колебаниям
Сообщение10.08.2015, 18:25 
Аватара пользователя


13/02/13
777
♍ — ☉ — ⊕
Kephe в сообщении #1043992 писал(а):
Отвечу только по этому пункту.

Спасибо! Это достаточно полезное замечание.


Но всё-каки хочется добраться до такого же результата моим путём.

 Профиль  
                  
 
 Re: Помогите разобраться с выводом формул по колебаниям
Сообщение10.08.2015, 18:30 
Заслуженный участник


20/07/09
4026
МФТИ ФУПМ
Atom001 в сообщении #1043981 писал(а):
Но легко видеть, что $t(k)$ - функция монотонная. Поэтому не получается найти такое $t$, при котором $x(t)$ примет максимальное значение.
Я не понимаю это "поэтому".

 Профиль  
                  
 
 Re: Помогите разобраться с выводом формул по колебаниям
Сообщение10.08.2015, 18:35 
Аватара пользователя


13/02/13
777
♍ — ☉ — ⊕
Nemiroff в сообщении #1044000 писал(а):
Я не понимаю это "поэтому".

Мы ведь должны найти знаки постоянства $x'(t)$. Но на всей области определения $x'(t)\geqslant 0$. А это значит, что $x_{\text{max}}=x(\infty)$.

 Профиль  
                  
 
 Re: Помогите разобраться с выводом формул по колебаниям
Сообщение10.08.2015, 18:47 
Заслуженный участник


20/07/09
4026
МФТИ ФУПМ
Atom001 в сообщении #1044002 писал(а):
Мы ведь должны найти знаки постоянства $x'(t)$
Какое постоянство? Вы же сами нашли корень $$t=\frac{1}{\omega_0}\arctg\frac{v_0}{x_0\omega_0}+\frac{\pi}{\omega_0}k \text{ , где } k\in\mathbb{Z}$$

Atom001 в сообщении #1044002 писал(а):
Но на всей области определения $x'(t)\geqslant 0$.
Это ещё почему?

 Профиль  
                  
 
 Re: Помогите разобраться с выводом формул по колебаниям
Сообщение10.08.2015, 18:58 
Аватара пользователя


13/02/13
777
♍ — ☉ — ⊕
Nemiroff в сообщении #1044006 писал(а):
Какое постоянство?

Прошу прощения. "Промежутки знакопостоянства".

Nemiroff в сообщении #1044006 писал(а):
Это ещё почему?

Это я погорячился.

Nemiroff в сообщении #1044006 писал(а):
Вы же сами нашли корень

Всё, я понял. Спасибо!
Просто запутался во всех этих функция и их производных.

 Профиль  
                  
 
 Re: Помогите разобраться с выводом формул по колебаниям
Сообщение10.08.2015, 19:11 
Заслуженный участник


20/07/09
4026
МФТИ ФУПМ
Ну ведь если бы вы полином третьей степени на экстремумы проверяли --- вы бы нашли производную, а это полином второй степени. Приравняли нулю --- два корня. И эти два корня в общем случае разные, и могут "монотонно зависеть от номера корня". И это ничего не значит. Тут то же самое, только корней счётное число.
Atom001 в сообщении #1044013 писал(а):
Всё, я понял. Спасибо!
Отлично.

 Профиль  
                  
 
 Re: Помогите разобраться с выводом формул по колебаниям
Сообщение10.08.2015, 22:55 
Заслуженный участник
Аватара пользователя


30/01/06
72407
Atom001 в сообщении #1043981 писал(а):
Но легко видеть, что $t(k)$ - функция монотонная. Поэтому не получается найти такое $t$, при котором $x(t)$ примет максимальное значение. Выходит, что колебаний вообще быть не должно.

Всё вы пишете правильно, кроме этого пункта. На самом деле, $x(t_\mathrm{zero}(k))$ все одинаковы для всех $k,$ так что наибольшее значение достигается в каждой из этих точек. (Вспомните, в определении наибольшего значения говорится, что значение функции в этой точке должно быть больше или равно, чем в любой другой точке - то есть, например, у константной функции в любой точке наибольшее значение.)

На будущее, если у вас есть периодическая функция (а ваша $x(t),$ как легко заметить, периодическая), то для анализа на экстремумы достаточно рассмотреть только один период этой функции.

 Профиль  
                  
 
 Re: Помогите разобраться с выводом формул по колебаниям
Сообщение11.08.2015, 11:31 
Аватара пользователя


13/02/13
777
♍ — ☉ — ⊕
Munin в сообщении #1044120 писал(а):
На будущее, если у вас есть периодическая функция (а ваша $x(t),$ как легко заметить, периодическая), то для анализа на экстремумы достаточно рассмотреть только один период этой функции.

Спасибо! Буду иметь ввиду.
А вообще, в школе исследованию периодических функций на экстремумы уделялось ну очень мало внимания, да и сам я очень редко с подобными задачами сталкивался. А тут вот столкнулся и сразу же запутался. Я-то решил заняться колебаниями, потому что начал изучать ДУ. Мне сразу же захотелось поприменять ДУ на практике, а в колебаниях это сделать как раз и можно. И я думал, что именно в ДУ где-нибудь ошибусь, но оказалось, что ошибся в поиске наибольшего значения.


Так. Амплитуду нашли. Теперь нужно найти начальную фазу (как известно, она подобно амплитуде зависит от начальных условий).
Но, если с амплитудой было всё просто (амплитудой мы назвали наибольший $x(t)$), то с начальной фазой не так всё очевидно. Какое условие нужно поставить, чтобы из уравнения координаты $x(t)$ найти и выразить через начальные условия начальную фазу колебаний?
Возможно, $\phi$ в формуле Kephe как раз и является этой самой $\varphi_0$. Но нельзя ли получить начальную фазу другим способом?

 Профиль  
                  
 
 Re: Помогите разобраться с выводом формул по колебаниям
Сообщение11.08.2015, 13:25 
Заслуженный участник


25/02/08
2961
Atom001
Начальная фаза и появляется именно в выражении вида $\[A\sin (\omega t + \varphi )\]$. Вообще здесь 2 типа начальных условий. Вы можете задавать начальную координату и начальную скорость (как это обычно делается в механике) и получить $\[{x_0}\cos \omega t + \frac{{{v_0}}}{\omega }\sin \omega t\]$, либо задавать в качестве начальных условий амплитуду и фазу (как это более удобно в колебаниях), и тогда получите $\[A\sin (\omega t + \varphi )\]$. Переход от одних к другим уже описал Kephe
$\[a\sin x + b\cos x = \sqrt {{a^2} + {b^2}} \sin (x + {\mathop{\rm arctg}\nolimits} \frac{b}{a})\]$
Арктангенс не всегда даёт верную "чётверть", где находится фаза, поэтому можно использовать
$\[\left\{ \begin{array}{l}
\sin \varphi  = \frac{b}{{\sqrt {{a^2} + {b^2}} }}\\
\cos \varphi  = \frac{a}{{\sqrt {{a^2} + {b^2}} }}
\end{array} \right.\]$
Можете сразу искать амплитуду и фазу через $ \[{x_0}\]$ и $\[{v_0}\]$ записав $\[x\]$ в виде $\[A\sin (\omega t + \varphi )\]$. Тут найдите скорость и получите два уравнения.

 Профиль  
                  
 
 Re: Помогите разобраться с выводом формул по колебаниям
Сообщение11.08.2015, 13:28 
Аватара пользователя


13/02/13
777
♍ — ☉ — ⊕
Ms-dos4
Ясно. Спасибо!

 Профиль  
                  
 
 Re: Помогите разобраться с выводом формул по колебаниям
Сообщение11.08.2015, 13:37 
Заслуженный участник
Аватара пользователя


30/01/06
72407
На самом деле, у слов "амплитуда" и "начальная фаза" нет такого общего универсального смысла, как вы думаете. Их нельзя использовать по отношению к какой-то произвольной функции колебаний, и всегда находить аналитическими способами. Скорее, смысл этих слов более условный. Когда колебание имеет вид синусоиды, то просто условились называть амплитудой и начальной фазой подчёркнутые константы в формуле $\underline{A}\sin(\omega t+\underline{\varphi})$ (чаще пишут формулу с косинусом $\underline{A}\cos(\omega t+\underline{\varphi}),$ причём вы сразу видите, что здесь начальной фазой называют другую величину).

Если есть какое-то другое периодическое колебание, то понятия "амплитуды" и "начальной фазы" применяются реже. Но иногда можно их понимать так:
- амплитуда - это размах от наибольшего до наименьшего значения за период;
- начальная фаза - это задержка функции по времени от какого-то условленного, отдельно оговорённого "начального положения", которая происходит без изменения, или почти без изменения, формы самой функции.

Ещё нередко встречается формула затухающих колебаний: $\underline{A}\,e^{-t/\tau}\cos(\omega t+\underline{\varphi})$ - тогда тоже можно называть эти константы $A$ и $\varphi$ амплитудой и фазой, а $\tau$ называется характе́рным временем затухания.

-- 11.08.2015 13:41:47 --

(Оффтоп)

Ms-dos4 в сообщении #1044410 писал(а):
Арктангенс не всегда даёт верную "чётверть", где находится фаза

Для этого как раз используется условная функция $\operatorname{atan2}(x,y)$ от двух переменных, которая не только даёт арктангенс, но и размещает его в нужной четверти, смотря по знакам аргументов. В программировании, такую функцию математические пакеты предоставляют почти всегда (и она даже реализована как команда микропроцессора).

 Профиль  
                  
 
 Re: Помогите разобраться с выводом формул по колебаниям
Сообщение11.08.2015, 13:47 


19/07/15
74
На всякий случай уточню, что под $\varphi = atan2(A, B)$ имел в виду стандартное в компьютерной сфере (но не уверен, что всем понятное) обозначение для арктангенса от двух параметров: https://en.wikipedia.org/wiki/Atan2
Такая функция позволяет получить правильную фазу для всех возможных сочетаний знаков множителей при $\cos$ и $\sin$.

Об этой тонкости с выбором правильной четверти чуть выше уже упомянул Ms-dos4.

edit: Munin опередил.

 Профиль  
                  
 
 Re: Помогите разобраться с выводом формул по колебаниям
Сообщение11.08.2015, 14:11 
Аватара пользователя


13/02/13
777
♍ — ☉ — ⊕
Munin
Понятно. Спасибо!

Munin в сообщении #1044419 писал(а):
Для этого как раз используется условная функция $\operatorname{atan2}(x,y)$ от двух переменных

А что это такое? Я вообще-то думал, что Kephe ошибся и вместо $\arctan\frac{B}{A}$ написал $\operatorname{atan2}(A,B)$. Эта функция используется только в программировании или в математике тоже? Я просто в математике никогда такого не встречал.

Kephe в сообщении #1044423 писал(а):
https://en.wikipedia.org/wiki/Atan2

Наверное всё-таки в математике данная функция не используется, потому что куда проще пользоваться системой из синуса и косинуса.

-- 11.08.2015, 19:11 --

Я ещё нашёл период колебаний. Всё сошлось. Может быть подкинете ещё пару задачек по этой теме (в идеале, чтобы было связано с ДУ).

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 35 ]  На страницу 1, 2, 3  След.

Модераторы: photon, whiterussian, profrotter, Jnrty, Aer, Парджеттер, Eule_A, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group