Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Metodichka_SMIFV

.pdf
Скачиваний:
34
Добавлен:
23.05.2015
Размер:
5.45 Mб
Скачать

 

 

 

 

 

11

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

dx(t)/dt 0

 

 

 

 

 

 

 

 

 

 

-1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

-2

 

 

 

 

 

 

 

 

 

 

-3

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-4

 

 

 

 

 

 

 

 

 

 

-1

-0.8

-0.6

-0.4

-0.2

0

0.2

0.4

0.6

0.8

1

 

 

 

 

 

x(t)

 

 

 

 

 

 

 

1. Фазовая траектория для НУ x(0)=0,5;

 

 

 

 

 

dx(0)/dt=0;5

 

 

 

 

2. Фазовая траектория для НУ x(0)=0,5;

 

 

 

 

 

 

dx(0)/dt=0;5

 

 

 

Рис. 6. Фазовые траектории для решений уравнения (3) при различных начальных условиях.

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

Пусть теперь у маятника появилось трение в точке подвеса. Тогда уравнение (3) надо переписать в виде:

(6) d 2 x 2 dx 2 x 0 dt2 dt

От уравнения маятника без трения оно отличается членом 2 dxdt ,

содержащим первую производную от координаты (скорость). Его решение

имеет вид (рис. 7):

(7) x xme t cos( t )

 

 

 

 

 

 

12

 

 

 

 

 

 

 

 

 

 

 

0.4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x (син.), 0

 

 

 

 

 

 

 

 

 

 

 

 

 

dx/dt (красн.)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-0.2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-0.4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-0.6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-0.8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.2

0.4

0.6

0.8

1

1.2

1.4

1.6

1.8

2

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

x 104

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис.

7. Колебания маятника с трением в точке

 

 

 

подвеса. Начальные условия: x=0,3; dx/dt=0,5.

 

 

 

Из-за трения положение маятника стало гармонической функцией,

затухающей по экспоненте. – коэффициент затухания.

 

 

 

 

 

Фазовая траектория для такого осциллятора есть спираль (рис. 8). Из

любого пункта фазового пространства спираль придет в начало координат, в

положение равновесия, в нуль. Эта точка «притягивает» фазовую

траекторию, которая закончится в ней откуда бы она не исходила. Поэтому

 

 

 

 

 

 

 

 

 

 

 

эта

 

точка

для

0.6

 

 

 

 

 

 

 

 

 

 

траектории является

 

 

 

 

 

 

 

 

 

 

 

аттрактором.

 

 

 

 

 

 

 

 

 

 

 

 

0.4

 

 

 

 

 

 

 

 

 

 

(attraction (англ.)-

 

 

 

 

 

 

 

 

 

 

 

притягивать).

 

0.2

 

 

 

 

 

 

 

 

 

 

 

Теперь можно

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

скомпенсировать

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

потери

 

на

трение

-0.2

 

 

 

 

 

 

 

 

 

 

(как это происходит

-0.4

 

 

 

 

 

 

 

 

 

 

в

 

механических

 

 

 

 

 

 

 

 

 

 

часах,

 

основой

 

 

 

 

 

 

 

 

 

 

 

 

-0.6

 

 

 

 

 

 

 

 

 

 

устройства которых

 

 

 

 

 

 

 

 

 

 

 

также

 

является

-0.8

 

 

 

 

 

 

 

 

 

 

маятник.

В

часах

 

 

 

 

 

 

 

 

 

 

 

-1

 

 

 

 

 

 

 

 

 

 

 

 

компенсацию

 

 

 

 

 

 

 

 

 

 

осуществляет

гиря

-0.3

-0.2

-0.1

0

 

0.1

0.2

0.3

 

0.4

 

Рис. 8. Фазовая траектория маятника с затуханием

или пружина).

 

 

 

 

 

Чтобы

x 10 4

13

скомпенсировать потери на трение, в осциллятор нужно периодически вносить энергию, или как говорят физики, его возбуждать. Уравнение маятника с трением и периодическим возбуждением имеет вид:

(8)d 2 x 2 dx 2 x Fm cos t dt 2 dt m0

где Fm- возбуждающая сила. - частота собственных колебаний осциллятора, - частота вынуждаемых силой колебаний. Такой осциллятор имеет принципиальное отличие от идеального маятника: при наличии возбуждения и начальных уловиях, лежащих в определенной области, колебания стремятся выйти на вполне определенный уровень (рис. 9).

10

5

x(спл),

dx/dt(пункт) 0

-5

-

0

0.

1

1.

2

2.

3

 

 

 

 

 

 

x 10 4

x(спл), dx/dt(пункт)

5

0 -5

0

0.

1

1.

2

2.

3

t

Рис. 9. Колебания маятника с трением в точке подвеса и с периодическим возбуждением. Верхний график: начальные условия: x=2; dx/dt=10. Нижний: начальные

условия: x=1; dx/dt=1. Для обоих случаев колебания с течением времени выходят на один и тот же режим.

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

Особенно наглядно это свойство реального осциллятора демонстрирует следующий рисунок (рис. 10), на котором приведены фазовые траектории реального маятника. Из разных начальных условий траектории приходят на один эллипс. Этот эллипс является аттрактором.

14

10

5

dx/dt 0

-5

-10

 

 

 

 

 

 

 

-2

-1

0

1

2

3

 

 

x

 

 

 

 

Рис. 10. Фазовые траектории маятника с трением в точке подвеса и с периодическим возбуждением

Системы, которые существуют за счет рассеивания или диссипации энергии, как в рассмотренном примере, носят название диссипативных систем. Диссипативная система в процессе своей эволюции всегда стремится попасть на аттрактор. Для системы может существовать несколько аттракторов. Например, для рассмотренного примера, аттрактором является также точка с координатами (0,0).

Следующие рисунки показывают, что происходит с осциллятором, когда он еще сильнее приближен к реальности, когда он становится нелинейным. Для этого достаточно чтобы в последнем члене левой части уравнения (8) коэффициент 2 умножался не на переменную x, а на sin(x).

 

 

 

 

 

0

 

 

 

d 2 x

 

dx

 

F

(9)

 

 

2

 

2 sin x

m

cos t

 

2

 

 

 

dt

 

dt

0

m

Поведение маятника меняется в корне (рис. 11 – зависимость положения и скорости от времени, рис. 12 – фазовая траектория).

15

60

50

40

x (син) 30 dx/dt (кр) 20

10

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

-10

 

 

 

 

 

 

 

 

-20 0

1

2

3

4

5

6

7

4

 

 

 

t

 

 

 

x 10

 

 

 

 

 

 

 

 

Рис. 11. Зависимость положения и скорости от времени

30

 

 

 

 

 

 

 

25

 

 

 

 

 

 

 

20

 

 

 

 

 

 

 

15

 

 

 

 

 

 

 

10

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

dx/dt 0

 

 

 

 

 

 

 

-5

 

 

 

 

 

 

 

-10

 

 

 

 

 

 

 

-15

 

 

 

 

 

 

 

-20

25

30

35

40

45

50

55

20

x

Рис. 12. Фазовая траектория «нелинейного» маятника

16

Генератор (осциллятор) Ван-дер-Поля2

Представляет собой математическую модель электронного генератора:

d 2 x

0 (1

x2

)

dx

2 x

 

 

 

2

2

 

, где: x- координата; t- время;

0 -

dt

 

 

dt

 

 

x0

 

 

 

коэффициент затухания; – частота собственных колебаний; x0- начальная амплитуда.

Параметр затухания 0 (1 x2 ) отрицателен при малых амплитудах

x02

x, положителен при больших. Тогда при <0 ограничивается рост энергии колебаний и при >0 компенсируются потери на трение (затухание).

У такого осциллятора колебания с малой амплитудой усиливаются, а с большой – затухают (рис. 13,14).

Van der Pol

 

50

 

 

 

 

 

 

 

40

 

 

 

 

 

 

 

30

 

 

 

 

 

 

 

20

 

 

 

 

 

 

 

10

 

 

 

 

 

 

d(THETA)/dt

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-10

 

 

 

 

 

 

 

-20

 

 

 

 

 

 

 

-30

 

 

 

 

 

 

 

-40

 

 

 

 

 

 

 

-50

 

 

 

 

 

 

 

-5

0

5

10

15

20

25

 

 

 

 

THETA

 

 

 

Рис. 13 Колебания с малой амплитудой (синяя) усиливаются, а с большой

(красная) – затухают. Фазовое пространство. Выход из разных начальных

условий происходит на аттрактор. 0=5; =4.176.

 

 

 

2 Van der Pol, B., "On relaxation-oscillations", The London, Edinburgh and Dublin Phil. Mag. & J. of Sci., 2(7), 978-992 (1927).

 

 

 

 

 

 

17

 

 

 

 

 

25

 

 

 

 

 

 

 

 

 

 

20

 

 

 

 

 

 

 

 

 

 

15

 

 

 

 

 

 

 

 

 

THETA

10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

-50

0.5

1

1.5

2

2.5

3

3.5

4

4.5

 

 

 

 

 

 

time

 

 

 

4

 

 

 

 

 

 

 

 

 

 

x 10

 

50

 

 

 

 

 

 

 

 

 

d(THETA)/dt

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-500

0.5

1

1.5

2

2.5

3

3.5

4

4.5

 

 

 

 

 

 

time

 

 

 

4

 

 

 

 

 

 

 

 

 

 

x 10

Рис. 14 Колебания с малой амплитудой усиливаются, а с большой – затухают. Вверху амплитуда, внизускорость

Ниже (рис. 15-20) приводятся решения в фазовом пространстве, полученные для разных значений и 0.

18

Van der Pol

 

30

 

 

 

 

 

 

 

 

20

 

 

 

 

 

 

 

 

10

 

 

 

 

 

 

 

d(THETA)/dt

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-10

 

 

 

 

 

 

 

 

-20

 

 

 

 

 

 

 

 

-30

 

 

 

 

 

 

 

 

-10

-5

0

5

10

15

20

25

 

 

 

 

 

THETA

 

 

 

Рис. 15 0=0,5; =4.176. Фазовая траектория

 

 

 

 

30

 

 

 

 

 

 

 

20

 

 

 

 

 

 

THETA

10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

-100

0.5

1

1.5

2

2.5

3

 

 

 

 

time

 

 

4

 

 

 

 

 

 

 

x 10

 

30

 

 

 

 

 

 

 

20

 

 

 

 

 

 

d(THETA)/dt

10

 

 

 

 

 

 

0

 

 

 

 

 

 

-10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-20

 

 

 

 

 

 

 

-300

0.5

1

1.5

2

2.5

3

 

 

 

 

time

 

 

4

 

 

 

 

 

 

 

x 10

Рис. 16 0=0,5; =4.176. Координата (сверху), скорость (снизу)

19

Van der Pol

 

25

 

 

 

 

 

 

 

20

 

 

 

 

 

 

 

15

 

 

 

 

 

 

 

10

 

 

 

 

 

 

 

5

 

 

 

 

 

 

THETA

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-5

 

 

 

 

 

 

 

-10

 

 

 

 

 

 

 

-15

 

 

 

 

 

 

 

-20

 

 

 

 

 

 

 

-25

 

 

 

 

 

 

 

-5

0

5

10

15

20

25

 

 

 

 

time

 

 

 

Рис. 17 0=5; =4.176 Фазовая траектория

 

 

 

 

30

 

 

 

 

 

20

 

 

 

 

THETA

10

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

-100

0.5

1

1.5

2

 

 

 

 

time

4

 

 

 

 

 

x 10

 

20

 

 

 

 

d(THETA)/dt

10

 

 

 

 

0

 

 

 

 

-10

 

 

 

 

 

 

 

 

 

 

-20

 

 

 

 

 

0

0.5

1

1.5

2

 

 

 

 

time

4

 

 

 

 

 

x 10

Рис. 18 0=5; =4.176 Координата (сверху), скорость (снизу)

20

Van der Pol

 

60

 

 

 

 

 

 

 

 

 

40

 

 

 

 

 

 

 

 

 

20

 

 

 

 

 

 

 

 

d(THETA)/dt

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-20

 

 

 

 

 

 

 

 

 

-40

 

 

 

 

 

 

 

 

 

-60

 

 

 

 

 

 

 

 

 

-15

-10

-5

0

5

10

15

20

25

 

 

 

 

 

THETA

 

 

 

 

Рис. 19 =4.176; 0=.1; Фазовая траектория

 

 

 

 

30

 

 

 

 

 

 

 

20

 

 

 

 

 

 

THETA

10

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-10

 

 

 

 

 

 

 

-200

1

2

3

4

5

6

 

 

 

 

time

 

 

4

 

 

 

 

 

 

 

x 10

 

60

 

 

 

 

 

 

 

40

 

 

 

 

 

 

d(THETA)/dt

20

 

 

 

 

 

 

0

 

 

 

 

 

 

-20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-40

 

 

 

 

 

 

 

-600

1

2

3

4

5

6

 

 

 

 

time

 

 

4

 

 

 

 

 

 

 

x 10

Рис. 20 =4.176; 0=.1; Координата (сверху), скорость (снизу)

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]