Необходимое условие гармонических колебаний. Колебания. Гармонические колебания. Уравнение гармонических колебаний. Уравнение затухающих колебаний
Простейшим видом колебаний являются гармонические колебания - колебания, при которых смещение колеблющейся точки от положения равновесия изменяется с течением времени по закону синуса или косинуса.
Так, при равномерном вращении шарика по окружности его проекция (тень в параллельных лучах света) совершает на вертикальном экране (рис. 1) гармоническое колебательное движение.
Смещение от положения равновесия при гармонических колебаниях описывается уравнением (его называют кинематическим законом гармонического движения) вида:
где х - смешение - величина, характеризующая положение колеблющейся точки в момент времени t относительно положения равновесия и измеряемая расстоянием от положения равновесия до положения точки в заданный момент времени; А - амплитуда колебаний - максимальное смещение тела из положения равновесия; Т - период колебаний - время совершения одного полного колебания; т.е. наименьший промежуток времени, по истечении которого повторяются значения физических величин, характеризующих колебание; - начальная фаза;
Фаза колебании в момент времени t. Фаза колебаний - это аргумент периодической функции, который при заданной амплитуде колебаний определяет состояние колебательной системы (смещение, скорость, ускорение) тела в любой момент времени.
Если в начальный момент времени колеблющаяся точка максимально смещена от положения равновесия, то , а смещение точки от положения равновесия изменяется по закону
Если колеблющаяся точка при находится в положении устойчивого равновесия, то смещение точки от положения равновесия изменяется по закону
Величину V, обратную периоду и равную числу полных колебаний, совершаемых за 1 с, называют частотой колебаний:
Если за время t тело совершает N полных колебаний, то
Величину 
, показывающую, сколько колебаний совершает тело за с, называют циклической (круговой) частотой
.
Кинематический закон гармонического движения можно записать в виде:
Графически зависимость смещения колеблющейся точки от времени изображается косинусоидой (или синусоидой).
На рисунке 2, а представлен график зависимости от времени смещения колеблющейся точки от положения равновесия для случая .
Выясним, как изменяется скорость колеблющейся точки со временем. Для этого найдем производную по времени от этого выражения:
где - амплитуда проекции скорости на ось х.
Эта формула показывает, что при гармонических колебаниях проекция скорости тела на ось х изменяется тоже по гармоническому закону с той же частотой, с другой амплитудой и опережает по фазе смешение на (рис. 2, б).
Для выяснения зависимости ускорения найдем производную по времени от проекции скорости:
где - амплитуда проекции ускорения на ось х.
При гармонических колебаниях проекция ускорения опережает смещение по фазе на к (рис. 2, в).
Аналогично можно построить графики зависимостей
Учитывая, что , формулу для ускорения можно записать
т.е. при гармонических колебаниях проекция ускорения прямо пропорциональна смещению и противоположна ему по знаку, т.е. ускорение направлено в сторону, противоположную смещению.
Так, проекция ускорения - это вторая производная от смещения , то полученное соотношение можно записать в виде:
Последнее равенство называют уравнением гармонических колебаний .
Физическую систему, в которой могут существовать гармонические колебания, называют гармоническим осциллятором , а уравнение гармонических колебаний - уравнением гармонического осциллятора .
Гармоническое колебание - явление периодического изменения какой-либо величины, при котором зависимость от аргумента имеет характер функции синуса или косинуса. Например, гармонически колеблется величина, изменяющаяся во времени следующим образом:
где х - значение изменяющейся величины, t - время, остальные параметры - постоянные: А - амплитуда колебаний, ω - циклическая частота колебаний, - полная фаза колебаний, - начальная фаза колебаний.
Обобщенное гармоническое колебание в дифференциальном виде
(Любое нетривиальное решение этого дифференциального уравнения - есть гармоническое колебание с циклической частотой )
Виды колебаний
Свободные колебания совершаются под действием внутренних сил системы после того, как система была выведена из положения равновесия. Чтобы свободные колебания были гармоническими, необходимо, чтобы колебательная система была линейной (описывалась линейными уравнениями движения), и в ней отсутствовала диссипация энергии (последняя вызвала бы затухание).
Вынужденные колебания совершаются под воздействием внешней периодической силы. Чтобы они были гармоническими, достаточно чтобы колебательная система была линейной (описывалась линейными уравнениями движения), а внешняя сила сама менялась со временем как гармоническое колебание (то есть чтобы зависимость от времени этой силы была синусоидальной).
Уравнение гармонических колебаний
|
Уравнение (1)
|
дает зависимость колеблющейся величины S от времени t; это и есть уравнение свободных гармонических колебаний в явном виде. Однако обычно под уравнением колебаний понимают иную запись этого уравнения, в дифференциальной форме. Возьмем для определенности уравнение (1) в виде
дважды продифференцируем его по времени:
Видно, что выполняется следующее соотношение:
которое и называется уравнением свободных гармонических колебаний (в дифференциальной форме). Уравнение (1) является решением дифференциального уравнения (2). Поскольку уравнение (2) - дифференциальное уравнение второго порядка, необходимы два начальных условия для получения полного решения (то есть определения входящих в уравнение (1) констант A и ); например, положение и скорость колебательной системы при t = 0.
Математи́ческий ма́ятник - осциллятор, представляющий собой механическую систему, состоящую изматериальной точки, находящейся на невесомой нерастяжимой нити или на невесомом стержне в однородном поле сил тяготения. Период малых собственных колебаний математического маятника длины l неподвижно подвешенного в однородном поле тяжести с ускорением свободного падения g равен
и не зависит от амплитуды и массы маятника.
Физический маятник - осциллятор, представляющий собой твёрдое тело, совершающее колебания в поле каких-либо сил относительно точки, не являющейся центром масс этого тела, или неподвижной оси, перпендикулярной направлению действия сил и не проходящей через центр масс этого тела.
Основы теории Максвелла для электромагнитного поля
Вихревое электрическое поле
Из закона Фарадея ξ=dФ/dt следует, что любое изменение сцепленного с контуром потока магнитной индукции приводит к возникновению электродвижущей силы индукции и вследствие этого появляется индукционный ток. Следовательно, возникновение э.д.с. электромагнитной индукции возможно и в неподвижном контуре, находящемся в переменном магнитном поле. Однако э.д.с. в любой цепи возникает только тогда, когда в ней на носители тока действуют сторонние силы - силы неэлектростатического происхождения (см. § 97). Поэтому возникает вопрос о природе сторонних сил в данном случае.
Опыт показывает, что эти сторонние силы не связаны ни с тепловыми, ни с химическими процессами в контуре; их возникновение также нельзя объяснить силами Лоренца, так как они на неподвижные заряды не действуют. Максвелл высказал гипотезу, что всякое переменное магнитное поле возбуждает в окружающем пространстве электрическое поле, которое
и является причиной возникновения индукционного тока в контуре. Согласно представлениям Максвелла, контур, в котором появляется э.д.с., играет второстепенную роль, являясь своего рода лишь «прибором», обнаруживающим это поле.
первое уравнение Максвелла утверждает, что изменения электрического поля порождают вихревое магнитное поле.
Второе уравнен ие
Максвелла выражает закон электромагнитной индукции Фарадея: ЭДС в любом замкнутом контуре равна скорости изменения (т. е. производной по времени) магнитного потока. Но ЭДС равна касательной составляющей вектора напряженности электрического поля Е, помноженной на длину контура. Чтобы перейти к ротору, как и в первом уравнении Максвелла, достаточно разделить ЭДС на площадь контура, а последнюю устремить к нулю, т. е. взять маленький контур, охватывающий рассматриваемую точку пространства (рис. 9,в). Тогда в правой части уравнения будет уже не поток, а магнитная индукция, поскольку поток равен индукции, помноженной на площадь контура.
Итак, получаем: rotE = - dB/dt.
Таким образом, вихревое электрическое поле порождается изменениями магнитного, что и подано на рис. 9,в и представлено только что приведенной формулой.
Третье и четвертое уравнения
Максвелла имеют дело с зарядами и порождаемыми ими полями. Они основаны на теореме Гаусса, утверждающей, что поток вектора электрической индукции через любую замкнутую поверхность равен заряду внутри этой поверхности.
На уравнениях Максвелла основана целая наука - электродинамика, позволяющая строгими математическими методами решить множество полезных практических задач. Можно рассчитать, например, поле излучения различных антенн как в свободном пространстве, так и вблизи поверхности Земли или около корпуса какого-либо летательного аппарата, например, самолета или ракеты. Электродинамика позволяет рассчитать конструкцию волноводов и объемных резонаторов - устройств, применяющихся на очень высоких частотах сантиметрового и миллиметрового диапазонов волн, где обычные линии передачи и колебательные контуры уже непригодны. Без электродинамики невозможно было бы развитие радиолокации, космической радиосвязи, антенной техники и многих других разделов современной радиотехники.
Ток смещения
ТОК СМЕЩЕ́НИЯ, величина, пропорциональная скорости изменения переменного электрического поля в диэлектрике или вакууме. Название «ток» связано с тем, что ток смещения, так же как и ток проводимости, порождает магнитное поле.
При построении теории электромагнитного поля Дж. К. Максвелл выдвинул гипотезу (впоследствии подтвержденную на опыте) о том, что магнитное поле создается не только движением зарядов (током проводимости, или просто током), но и любым изменением во времени электрического поля.
Понятие ток смещения введено Максвеллом для установления количественных соотношений между изменяющимся электрическим полем и вызываемым им магнитным полем.
В соответствии с теорией Максвелла, в цепи переменного тока, содержащей конденсатор, переменное электрическое поле в конденсаторе в каждый момент времени создает такое магнитное поле, какое создавал бы ток, (названный током смещения), если бы он протекал между обкладками конденсатора. Из этого определения следует, что J см = J (т. е., численные значения плотности тока проводимости и плотности тока смещения равны), и, следовательно, линии плотности тока проводимости внутри проводника непрерывно переходят в линии плотности тока смещения между обкладками конденсатора. Плотность тока смещения j см характеризует скорость изменения электрической индукции D во времени:
J см = + ?D/?t.
Ток смещения не выделяет джоулевой теплоты, его основное физическое свойство - способность создавать в окружающем пространстве магнитное поле.
Вихревое магнитное поле создается полным током, плотность которого j , равна сумме плотности тока проводимости и тока смещения?D/?t. Именно поэтому для величины?D/?t и было введено название ток.
Гармоническим осциллятором называется система, которая совершает колебания, описываемые выражением вида d 2 s/dt 2 + ω 0 2 s = 0 или
где две точки сверху означают двукратное дифференцирование по времени. Колебания гармонического осциллятора есть важный пример периодического движения и служат точной или приближенной моделью во многих задачах классической и квантовой физики. В качестве примеров гармонического осциллятора могут быть пружинный, физический и математический маятники, колебательный контур (для токов и напряжений настолько малых, что можно было бы элементы контура считать линейными).
Гармонические колебания
Наряду с поступательными и вращательными движениями тел в механике значительный интерес представляют и колебательные движения. Механическими колебанияминазывают движения тел, повторяющиеся точно (или приблизительно) через одинаковые промежутки времени. Закон движения тела, совершающего колебания, задается с помощью некоторой периодической функции времени x = f (t ). Графическое изображение этой функции дает наглядное представление о протекании колебательного процесса во времени.
Примерами простых колебательных систем могут служить груз на пружине или математический маятник (рис. 2.1.1).
Механические колебания, как и колебательные процессы любой другой физической природы, могут быть свободными и вынужденными . Свободные колебания совершаются под действием внутренних сил системы, после того, как система была выведена из состояния равновесия. Колебания груза на пружине или колебания маятника являются свободными колебаниями. Колебания, происходящие под действием внешних периодически изменяющихся сил, называются вынужденными Простейшим видом колебательного процесса являются простые гармонические колебания , которые описываются уравнением
Частота колебаний f показывает, сколько колебаний совершается за 1 с. Единица частоты – герц (Гц). Частота колебаний f связана с циклической частотой ω и периодом колебаний T соотношениями:
дает зависимость колеблющейся величины S от времени t ; это и есть уравнение свободных гармонических колебаний в явном виде. Однако обычно под уравнением колебаний понимают иную запись этого уравнения, в дифференциальной форме. Возьмем для определенности уравнение (1) в виде
дважды продифференцируем его по времени:
Видно, что выполняется следующее соотношение:
которое и называется уравнением свободных гармонических колебаний (в дифференциальной форме). Уравнение (1) является решением дифференциального уравнения (2). Поскольку уравнение (2) - дифференциальное уравнение второго порядка, необходимы два начальных условия для получения полного решения (то есть определения входящих в уравнение (1) констант A и j 0); например, положение и скорость колебательной системы при t = 0.
Сложение гармонических колебаний одного направления и одинаковой частоты. Биения
Пусть совершаются два гармонических колебания одного направления и одинаковой частоты
Уравнение результирующего колебания будет иметь вид
Убедимся в этом, сложив уравнения системы (4.1)
Применив теорему косинусов суммы и сделав алгебраические преобразования:
Можно найти такие величины А и φ0 , чтобы удовлетворялись уравнения
Рассматривая (4.3) как два уравнения с двумя неизвестными А и φ0, найдем, возведя их в квадрат и сложив, а затем разделив второе на первое:
Подставляя (4.3) в (4.2), получим:
Или окончательно, используя теорему косинусов суммы, имеем:
Тело, участвуя в двух гармонических колебаниях одного направления и одинаковой частоты, совершает также гармоническое колебание в том же направлении и с той же частотой, что и складываемые колебания. Амплитуда результирующего колебания зависит от разности фаз (φ2-φ1) сгладываемых колебаний.
В зависимости от разности фаз (φ2-φ1):
1) (φ2-φ1) = ±2mπ (m=0, 1, 2, …), тогда A= А1+А2, т. е. амплитуда результирующего колебания А равна сумме амплитуд складываемых колебаний;
2) (φ2-φ1) = ±(2m+1)π (m=0, 1, 2, …), тогда A= |А1-А2|, т. е. амплитуда результирующего колебания равна разности амплитуд складываемых колебаний
Периодические изменения амплитуды колебания, возникающие при сложении двух гармонических колебаний с близкими частотами, называются биением.
Пусть два колебания мало отличаются по частоте. Тогда амплитуды складываемых колебаний равны А, а частоты равны ω и ω+Δω, причем Δω намного меньше ω. Начало отсчета выберем так, чтобы начальные фазы обоих колебаний были равны нулю:
Решим систему
Решение системы:
Результирующее колебание можно рассматривать как гармоническое с частотой ω, амплитуда А, которого изменяется по следующему периодическому закону:
Частота изменения А в два раза больше частоты изменения косинуса. Частота биений равна разности частот складываемых колебаний: ωб = Δω
Период биений:
Определение частоты тона (звука определенной высоты биений эталонным и измеряемым колебаниями - наиболее широко применяемый на метод сравнения измеряемой величины с эталонной. Метод биений используется для настройки музыкальных инструментов, анализа слуха и т. д.
Похожая информация.
На хабре было несколько статей по преобразованию Фурье и о всяких красивостях типа Цифровой Обработки Сигналов (ЦОС), но неискушённому пользователю совершенно не понятно, зачем всё это нужно и где, а главное как это применить.
АЧХ шума.
Лично мне после прочтения этих статей (например, этой) не стало понятно, что это и зачем оно нужно в реальной жизни, хотя было интересно и красиво.
Хочется не просто поглядеть красивые картинки, а так сказать, ощутить нутром, что и как работает. И я приведу конкретный пример с генерацией и обработкой звуковых файлов. Можно будет и послушать звук, и поглядеть его спектр, и понять, почему это так.
Статья не будет интересна тем, кто владеет теорией функций комплексной переменной, ЦОС и прочими страшными темами. Она скорее для любопытствующих, школьников, студентов и им сочувствующих:).
Сразу оговорюсь, я не математик, и многие вещи могу даже сказать неправильно (поправляйте личным сообщением), и данную статью пишу, опираясь на собственный опыт и собственное понимание текущих процессов. Если вы готовы, то поехали.
Пару слов о матчасти
Если мы вспомним школьный курс математики, то для построения графика синуса мы использовали круг. В общем-то так и получается, что вращательное движение можно превратить в синусоиду (как и любое гармоническое колебание). Самое лучшая иллюстрация этого процесса приведена в википедии
Гармонические колебания
Т.е. фактически график синуса получается из вращения вектора, который описывается формулой:
F(x) = A sin (ωt + φ),
Где A - длина вектора (амплитуда колебаний), φ - начальный угол (фаза) вектора в нулевой момент времени, ω - угловая скорость вращения, которая равна:
ω=2 πf, где f - частота в Герцах.
Как мы видим, что зная частоту сигнала, амплитуду и угол, мы можем построить гармонический сигнал.
Магия начинается тогда, когда оказывается, что представление абсолютно любого сигнала можно представить в виде суммы (зачастую бесконечной) различных синусоид. Иначе говоря, в виде ряда Фурье.
Я приведу пример из английской википедии . Для примера возьмём пилообразный сигнал.
Пилообразный сигнал
Его сумма будет представлена следующей формулой:
Если мы будем по очерёдно суммировать, брать сначала n=1, затем n=2 и т.д., то увидим, как у нас гармонический синусоидальный сигнал постепенно превращается в пилу:
Наверное красивее всего это иллюстрирует одна программа, найденная мной на просторах сети. Выше уже говорилось, что график синуса является проекцией вращающегося вектора, а как же быть в случае более сложных сигналов? Это, как ни странно, проекция множества вращающихся векторов, а точнее их суммы, и выглядит это всё так:
Вектора рисуют пилу.
Вообще рекомендую сходить самим по ссылке и попробовать самим поиграться с параметрами, и посмотреть как меняется сигнал. ИМХО более наглядной игрушки для понимания я ещё не встречал.
Ещё следует заметить, что есть обратная процедура, позволяющая получить из данного сигнала частоту, амплитуду и начальную фазу (угол), которое называется Преобразование Фурье.
Разложение в ряд Фурье некоторых известных периодических функций (отсюда)
Я детально на нём останавливаться не буду, но покажу, как это можно применить по жизни. В списке литературы порекомендую то, где можно почитать подробнее о матчасти.
Переходим к практическим упражнениям!
Мне кажется, что каждый студент задаётся вопросом, сидя на лекции, например по матану: зачем мне весь этот бред? И как правило, не найдя ответа в обозримом будущем, к сожалению, теряет интерес к предмету. Поэтому я сразу покажу практическое применение данных знаний, а вы эти знания уже будете осваивать сами:).
Всё дальнейшее я буду реализовывать на сях. Делал всё, конечно, под Linux, но никакой специфики не использовал, по идее программа будет компилироваться и работать под другими платформами.
Для начала напишем программу для формирования звукового файла. Был взят wav-файл, как самый простой. Прочитать про его структуру можно .
Если кратко, то структура wav-файла описывается так: заголовок, который описывает формат файла, и далее идёт (в нашем случае) массив 16-ти битных данных (остроконечник) длиной: частота_дискретизации*t секунд или 44100*t штук.
Для формирования звукового файла был взят пример . Я его немного модифицировал, исправил ошибки, и окончательная версия с моими правками теперь лежит на гитхабе тут
Сгенерируем двухсекундный звуковой файл с чистым синусом частотой 100 Гц. Для этого модифицируем программу таким образом:
#define S_RATE (44100) //частота дискретизации
#define BUF_SIZE (S_RATE*10) /* 2 second buffer */
….
int main(int argc, char * argv)
{
...
float amplitude = 32000; //берём максимальную возможную амплитуду
float freq_Hz = 100; //частота сигнала
/* fill buffer with a sine wave */
for (i=0; i Обращаю внимание, что формула чистого синуса соответствует той, о которой мы говорили выше. Амплитуда 32000 (можно было взять 32767) соответствует значению, которое может принимать 16-ти битное число (от минус 32767 до плюс 32767). В результате получаем следующий файл (можно его даже послушать любой звуковоспроизводящей программой). Откроем данный файл audacity и увидим, что график сигнала в действительности соответствует чистому синусу: Поглядим спектр этого синуса (Анализ->Построить график спектра) Виден чистый пик на 100 Гц (логарифмический масштаб). Что такое спектр? Это амплитудно-частотная характеристика. Существует ещё фазочастотная характеристика. Если помните, выше я говорил, что для построения сигнала надо знать его частоту, амплитуду и фазу? Так вот, можно из сигнала получить эти параметры. В данном случае у нас график соответствий частот амплитуде, при чём амплитуда у нас не в реальных единицах, а в Децибелах. Я понимаю, что чтобы объяснить, как работает программа, надо объяснить, что такое быстрое преобразование Фурье, а это как минимум ещё на одну некислую статью. Для начала алокируем массивы: C = calloc(size_array*2, sizeof(float)); // массив поворотных множителей
in = calloc(size_array*2, sizeof(float)); //входный массив
out = calloc(size_array*2, sizeof(float)); //выходной массив
Скажу лишь, что в программе мы читаем данные в массив длиной size_array (которое берём из заголовка wav-файла). While(fread(&value,sizeof(value),1,wav)) {
in[j]=(float)value;
j+=2;
if (j > 2*size_array) break;
}
Массив для быстрого преобразования Фурье должен представлять собой последовательность {re, im, re, im,… re, im}, где fft_size=1<< p - число точек БПФ. Объясняю нормальным языком: Int p2=(int)(log2(header.bytes_in_data/header.bytes_by_capture));
Логарифм от количество байт в данных, делённых на количество байт в одной точке. После этого считаем поворотные множители: Fft_make(p2,c);// функция расчёта поворотных множителей для БПФ (первый параметр степень двойки, второй алокированный массив поворотных множителей).
И скармливаем наш считанный массив в преобразователь Фурье: Fft_calc(p2, c, in, out, 1); //(единица означает, что мы получаем нормализованный массив).
На выходе мы получаем комплексные числа вида {re, im, re, im,… re, im}. Для тех, кто не знает, что такое комплексное число, поясню. Я не зря начал эту статью с кучи вращающихся векторов и кучи гифок. Так вот, вектор на комплесной плоскости определяется действительной координатой a1 и мнимой координатой a2. Или длиной (это у нас амплитуда Am) и углом Пси (фаза). Обратите внимание, что size_array=2^p2. Первая точка массива соответствует частоте 0 Гц (постоянная), последняя точка соответствует частоте дискретизации, а именно 44100 Гц. В результате мы должны рассчитать частоту, соответствующей каждой точке, которые будут отличаться на частоту дельта: Double delta=((float)header.frequency)/(float)size_array; //частота дискретизации на размер массива.
Алокируем массив амплитуд: Double * ampl;
ampl = calloc(size_array*2, sizeof(double));
И смотрим на картинку: амплитуда - это длина вектора. А у нас есть его проекции на действительную и мнимую ось. В результате у нас будет прямоугольный треугольник, и тут мы вспоминаем теорему Пифагора, и считаем длину каждого вектора, и сразу пишем её в текстовый файл: For(i=0;i<(size_array);i+=2) {
fprintf(logfile,"%.6f %f\n",cur_freq, (sqrt(out[i]*out[i]+out*out)));
cur_freq+=delta;
}
…
11.439514 10.943008
11.607742 56.649738
11.775970 15.652428
11.944199 21.872342
12.112427 30.635371
12.280655 30.329171
12.448883 11.932371
12.617111 20.777617
...
Теперь скармливаем получившейся программе тот звуковой файл синуса ./fft_an ../generate_wav/sin\ 100\ Hz.wav
format: 16 bits, PCM uncompressed, channel 1, freq 44100, 88200 bytes per sec, 2 bytes by capture, 2 bits per sample, 882000 bytes in data
chunk=441000
log2=18
size array=262144
wav format
Max Freq = 99.928 , amp =7216.136
И получаем текстовый файл АЧХ. Строим его график с помощью гнуплота Скрипт для построения: #! /usr/bin/gnuplot -persist
set terminal postscript eps enhanced color solid
set output "result.ps"
#set terminal png size 800, 600
#set output "result.png"
set grid xtics ytics
set log xy
set xlabel "Freq, Hz"
set ylabel "Amp, dB"
set xrange
#set yrange
plot "test.txt" using 1:2 title "AFC" with lines linestyle 1
Обратите внимание на ограничение в скрипте на количество точек по X: set xrange . Частота дискретизации у нас 44100, а если вспомнить теорему Котельникова, то частота сигнала не может быть выше половины частоты дискретизации, следовательно сигнал выше 22050 Гц нас не интересует. Почему так, советую прочитать в специальной литературе. Обратите внимание на резкий пик на частоте 100 Гц. Не забывайте, что по осям - логарифмический масштаб! Шерсть справа, как я думаю, ошибки преобразования Фурье (тут на память приходят окна). А давайте! Давайте поглядим спектры других сигналов! Double d_random(double min, double max)
{
return min + (max - min) / RAND_MAX * rand();
}
Она будет генерировать случайное число в заданном диапазоне. В результате main будет выглядеть так: Int main(int argc, char * argv)
{
int i;
float amplitude = 32000;
srand((unsigned int)time(0)); //инициализируем генератор случайных чисел
for (i=0; i Сгенерируем файл , (рекомендую к прослушиванию). Поглядим его в audacity. Поглядим спектр в программе audacity. И поглядим спектр с помощью нашей программы: Хочу обратить внимание на очень интересный факт и особенность шума - он содержит в себе спектры всех гармоник. Как видно из графика, спектр вполне себе ровный. Как правило, белый шум используется для частотного анализа пропускной способности, например, аудиоаппаратуры. Существуют и другие виды шумов: розовый, синий и другие
. Домашнее задание - узнать, чем они отличаются. А теперь давайте посмотрим другой интереснейший сигнал - меандр. Я там выше приводил табличку разложений различных сигналов в ряды Фурье, вы поглядите как раскладывается меандр, выпишите на бумажку, и мы продолжим. Для генерации меандра с частотой 25 Гц мы модифицируем в очередной раз наш генератор wav-файла: Int main(int argc, char * argv)
{
int i;
short int meandr_value=32767;
/* fill buffer with a sine wave */
for (i=0; i В результате получим звуковой файл (опять же, советую послушать), который сразу надо посмотреть в audacity Не будем томиться и поглядим его спектр: Пока не очень что-то понятно, что такое… А давайте поглядим несколько первых гармоник: Совсем другое дело! Ну-ка поглядим табличку. Смотрите-ка, у нас есть только 1, 3, 5 и т.д., т.е. нечётные гармоники. Мы так и видим, что у нас первая гармоника 25 Гц, следующая (третья) 75 Гц, затем 125 Гц и т.д., при этом у нас амплитуда постепенно уменьшается. Теория сошлась с практикой! Эта картинка прям как картинка из википедии , где для примера меандра берутся не все частоты, а только первые несколько. Сумма первых гармоник, и как меняется сигнал
Меандр так же активно используется в радиотехнике (надо сказать, что - это основа всей цифровой техники), и стоит понимать что при длинных цепях его может отфильтровать так, что, родная мама не узнает. Его так же используют для проверки АЧХ различных приборов. Ещё интересный факт, что глушилки телевизоров работали именно по принципу высших гармоник, когда сама микросхема генерировала меандр десятки МГц, а его высшие гармоники могли иметь частоты сотни МГц, как раз на частоте работы телевизора, и высшие гармоники успешно глушили сигнал вещания телевизора. Вообще тема подобных экспериментов бесконечная, и вы можете теперь сами её продолжить. Для тех, кто нифига не понял, что мы тут делаем, или наоборот, для тех, кто понял, но хочет разобраться ещё лучше, а так же для студентам, изучающим ЦОС, крайне рекомендую эту книгу. Это ЦОС для чайников, которым является автор данного поста. Там доступным даже для ребёнка языком рассказываются сложнейшие понятия. В заключении хочу сказать, что математика - царица наук, но без реального применения многие люди теряют к ней интерес. Надеюсь, данный пост подстегнёт вас к изучению такого замечательного предмета, как обработка сигналов, и вообще аналоговой схемотехнике (затыкайте уши, чтобы не вытекали мозги!). :) «Физика - 11 класс»
Ускорение - вторая производная координаты по времени.
Мгновенная скорость точки - это производная координаты точки по времени. где х" - вторая производная координаты по времени. При свободных колебаниях координата х
изменяется со временем так, что вторая производная координаты по времени прямо пропорциональна самой координате и противоположна ей по знаку. Гармонические колебания
Из математики: вторые производные синуса и косинуса по их аргументу пропорциональны самим функциям, взятым с противоположным знаком, и никакие другие функции таким свойством не обладают. Периодические изменения физической величины в зависимости от времени, происходящие по закону синуса или косинуса, называются гармоническими колебаниями
. Амплитуда колебаний
Амплитудой
гармонических колебаний называется модуль наибольшего смещения тела от положения равновесия. Амплитуда определяется начальными условиями, а точнее энергией, сообщаемой телу. График зависимости координаты тела от времени представляет собой косинусоиду. х = x m cos ω 0 t
Тогда уравнение движения, описывающее свободные колебания маятника: Период и частота гармонических колебаний.
При колебаниях движения тела периодически повторяются. Частота колебаний - это число колебаний в единицу времени. В Международной системе единиц (СИ) единица частоты называется герцем
(Гц) в честь немецкого физика Г. Герца. Число колебаний за 2π с равно: Величина ω 0 - это циклическая (или круговая) частота колебаний. Частоту свободных колебаний называют собственной частотой
колебательной системы. Зависимость частоты и периода свободных колебаний от свойств системы.
1. для пружинного маятника
Собственная частота колебаний пружинного маятника равна: Она тем больше, чем больше жесткость пружины k, и тем меньше, чем больше масса тела m. Период колебаний равен: Период колебаний пружинного маятника не зависит от амплитуды колебаний. 2. для нитяного маятника
Собственная частота колебаний математического маятника при малых углах отклонения нити от вертикали зависит от длины маятника и ускорения свободного падения: Период же этих колебаний равен Период колебаний нитяного маятника при малых углах отклонения не зависит от амплитуды колебаний. Период колебаний возрастает с увеличением длины маятника. От массы маятника он не зависит. Чем меньше g, тем больше период колебаний маятника и, следовательно, тем медленнее идут часы с маятником. Так, часы с маятником в виде груза на стержне отстанут за сутки почти на 3 с, если их поднять из подвала на верхний этаж Московского университета (высота 200 м). И это только за счет уменьшения ускорения свободного падения с высотой.
Чистый ламповый синус
График спектра
это массив комплексных чисел. Я даже боюсь представить, где используется комплексное преобразование Фурье, но в нашем случае мнимая часть у нас равна нулю, а действительная равна значению каждой точке масива.
Ещё одна особенность именно быстрого преобразования Фурье, что оно обсчитывает массивы, кратные только степени двойки. В результате мы должны вычислить минимальную степень двойки:
Вектор на комплексной плоскости
В результате получаем файл примерно такого вида:Пробуем!
Итак (барабанная дробь), запускаем скрипт и лицезреем:
Спектр нашего сигнала
А давайте побалуем?
Вокруг шум…
Для начала построим спектр шума. Тема про шумы, случайные сигналы и т.п. достойна отдельного курса. Но мы её коснёмся слегка. Модифицируем нашу программу генерации wav-файла, добавим одну процедуру:
Сигнал в audacity
Спектр
Наш спектр
А компот?
Его величество - меандр или меандр здорового человека
Спектр меандра
Первые гармоники
А теперь внимание! В реальной жизни сигнал меандра у нас имеет бесконечную сумму гармоник всё более и более высокой частоты, но как правило, реальные электрические цепи не могут пропускать частоты выше какой-то частоты (в силу индуктивности и ёмкости дорожек). В результате на экране осциллографа можно часто увидеть вот такой сигнал:
Меандр курильщика
Книга
Заключение
Удачи!
Ускорение точки - это производная ее скорости по времени, или вторая производная координаты по времени.
Поэтому уравнение движения маятника можно записать так:
Поэтому:
Координата тела, совершающего свободные колебания, меняется с течением времени по закону синуса или косинуса.
Промежуток времени Т, за который система совершает один полный цикл колебаний, называется периодом колебаний
.
Если одно колебание совершается за время Т то число колебаний за секунду
Через промежуток времени, равный одному периоду, колебания повторяются.
Часто для краткости циклическую частоту называют просто частотой.
Жесткая пружина сообщает телу большее ускорение, быстрее меняет скорость тела, а чем тело массивнее, тем медленнее оно изменяет скорость под влиянием силы.
