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

Методическое пособие 370

.pdf
Скачиваний:
4
Добавлен:
30.04.2022
Размер:
1.04 Mб
Скачать

Рис. 6.3. График импульса по условиям (6.3)

Условия (6.3) получаются из совместного рассмотрения уравнения прямых линий и определения необходимых коэффициентов, т.е.:

y1 kt,

y2 kt b,

0 kt b, b kT,

y2 k t T .

Аналогично можно получить выражение для функции, имеющей график, показанный на рис.6.4:

y1 kt, y2 kt b, 0 kT b, (6.4) b kT, y2 kt kT k T t .

Рис. 6.4. График пилообразного импульса

150

Симметричный импульс в соответствии с рис. 6.4 показан на рис. 6.5.

Рис. 6.5. Симметричный пилообразный импульс

Представление импульса (рис. 6.5) в записи для расчета на ЭВМ, с учетом (6.3) и (6.4), будет иметь вид

 

 

T

 

 

 

T

 

 

 

 

u t : if t

 

 

, y

, if t 3

 

, y

2

t

, y

t ,

 

 

 

4

1

 

 

4

 

 

3

 

 

 

 

 

 

 

 

 

где y1 kt, y2

k T t ,

y3 k t T .

 

 

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

Рис. 6.6. Импульс в форме трапеции

Описание импульса (рис. 6.6) следует производить в программном варианте, т.е.

151

 

 

 

u1t if

t

 

T

 

 

 

 

 

 

 

 

 

 

 

 

C N

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u2 t if t C 1

 

 

 

 

 

 

u t :

N C

 

 

 

 

 

(6.5)

 

 

 

T

 

 

 

 

 

 

T

 

 

 

 

 

u3 t if

 

 

 

t C 1

 

 

 

 

 

 

 

 

 

 

N C

 

 

 

 

 

 

 

 

N C

 

 

 

 

 

 

 

 

 

0 if t

T

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

где u t p t,

 

t p

 

T

 

 

t p

T

 

 

u

t

 

 

,

u

 

.

(6.6)

 

 

C N

1

2

 

 

 

N

3

 

 

 

 

 

Значения величин в (6.6) могут задаваться в различных ва-

риантах. При этом следует иметь в виду, что N T - скваж- tu

ность импульса с длительностью tu и периодом T , а величина

C - коэффициент наклона (фронт) боковых граней импульса.

6.2. Определение спектра сигнала через ряды Фурье

Порядок определения спектра периодического сигнала на основе рядов Фурье приведен в п. 5.2.

Будем считать, что для наглядного представления спектральных составляющих сигнала заданы:

количество гармоник n, период колебаний Т, частота сигнала f 1T ,

анализируемая функция u t , время анализа t.

С учетом этих обозначений определение коэффициентов ряда Фурье должно осуществляться по следующим формулам:

 

2

T

 

a n :

u t cos 6,28 f t n dt,

(6.7)

T

 

0

 

152

 

2

 

T

 

 

 

 

 

 

 

b n :

 

u t sin 6,28 f t n dt,

T

 

 

0

 

 

 

 

 

 

 

M n :

 

 

 

 

 

,

 

 

a2 n b2 n

 

 

 

b n

 

180

 

Y n : atan

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a n

 

 

D n : if a n 0,Y n ,Y n 180 .

(6.8)

(6.9)

(6.10)

(6.11)

График спектра амплитуд строится по формуле (6.9), а график спектра фаз по выражению (6.11).

Приведем примеры построения спектров с помощью про-

граммы Mathcad.

Пример 6.1. Построить спектр амплитуд и спектр фаз для прямоугольного импульса, определяемого формулой

 

T

 

u t : if t

 

, 0.9, if t T, 0.9, 0 .

2

 

 

Порядок набора формул в программе Mathcad для определения спектра следующий:

n: 0,1 10

T : 0.005

t : 0, 0.00001 0.05

1

 

 

 

 

 

 

T

 

 

 

 

 

 

f :

 

 

 

 

 

u t : if t

 

, 0.9, if t T, 0.9,

0

 

T

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

T

 

 

 

 

 

 

 

 

 

 

a n :

 

 

u t cos 6,28 f t n dt

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

T

 

 

 

 

 

 

 

 

 

 

b n :

 

 

u t sin 6,28 f t n dt

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

b n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

180

 

M n :

 

a2 n b2 n

 

 

 

 

Y n : atan

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a n

 

 

D n : if a n 0,Y n ,Y n 180 .

153

Затем следует, используя правила программы Mathcad, построить графики сигнала и спектров амплитуд и фаз

(рис. 6.7-6.9).

1

u(t) 0

1

0 0.002 0.004 0.006

t

Рис. 6.7. График входного сигнала

1.146

1.5

 

 

 

1

 

 

M(n)

 

 

 

 

0.5

 

 

0

0

5

10

 

0

 

0

n

10

 

Рис. 6.8. График спектра амплитуд

154

 

200

 

 

H(k)

0

 

 

 

200

5

10

 

0

 

 

k

 

Рис. 6.9. График спектра фаз

Восстановить сигнал по его спектру (рис.6.10) можно, если сложить, вычисленные ранее, гармоники с учетом их фаз и амплитуд, т.е.

3

 

F t M n cos 6,28 f n t D n .

(6.12)

n 1

n

2

F(t)

0 u(t)

2

0 0.002 0.004 0.006

t

Рис. 6.10. График восстановленного сигнала по значениям трех его гармоник

155

Пример 6.2. Используя представление сигнала в форме (6.5), проследить изменения в его спектре при различных параметрах импульса.

Зададимся следующими значениями в (6.5):

С : 100

N : 4 T : 0.008

t 0, 0.00001..0.008

f :

1

 

T

p : 6,28 f C

 

 

 

 

 

u1t

p t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

T

 

 

 

u2 t : p t

 

 

 

n: 0,1 20

u3 t : p

 

 

 

 

 

 

 

N

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

u1t if

t

 

T

 

 

 

 

 

 

 

 

 

 

C N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u2 t if t C 1

T

 

 

 

 

 

 

 

 

 

N C

 

 

 

 

 

(6.13)

 

u t :

 

T

 

 

 

T

 

 

 

 

u3 t if

 

t C 1

 

 

 

 

 

 

 

 

 

N C

 

 

 

 

 

 

 

 

 

N C

 

 

 

 

 

 

 

0 if t T . N

График сигнала для этого случая показан на рис. 6.11.

2

ut( ) 1

0

0

0.005

0.01

t

Рис. 6.11. График прямоугольного импульса со скважностью N 4

156

Вычисление спектра амплитуд (рис. 6.12) осуществляется по формулам:

 

2

 

T

 

 

 

a n :

 

u t cos 6,28 f n t dt

 

 

 

T

 

 

 

 

0

 

 

 

 

2

 

T

 

 

 

b n :

u t sin 6,28 f n t dt

M n :

a2 n b2 n

 

T

 

0

 

 

 

0.777

M (n)

8.195 10 3

1

 

 

0.5

 

 

0

10

20

0

0

n

20

trace 1

Рис. 6.12. График спектра амплитуд прямоугольного импульса при N 4

Если в (6.13) принять N 2, то график спектра амплитуд для сигнала (рис. 6.13) примет иной вид (рис. 6.14).

4

u(t) 2

0

0

0.005

0.01

t

Рис. 6.13. График сигнала при N 2

157

3.109

4

 

 

 

 

 

M ( n )

2

 

 

0.032

0

10

20

 

0

 

0

n

20

trace 1

Рис. 6.14. График спектра амплитуд прямоугольного импульса при N 2

Если в (6.13) принять С : 4 при N : 2, то график сигнала приобретет форму, показанную на рис.(6.15), где склоны импульса станут пологими.

4

u(t)

2

0

0

0.005

0.01

t

Рис. 6.15. График импульса в виде трапеции

График спектра амплитуд (рис. 6.16) у импульса с пологими склонами (рис. 6.15) существенно отличается от случая спектра импульса с крутыми склонами (рис. 6.14).

158

3

 

 

2

 

 

M(n)

 

 

1

 

 

0

10

20

0

 

n

 

Рис. 6.16. График спектра импульса с пологими склонами

Отличия в графиках по рис. 6.14 и рис. 6.16, прежде всего в том, что амплитуды высокочастотных гармоник в спектре импульса с пологими склонами значительно ниже, чем у импульса с крутыми склонами, т.е. спектр содержит меньшее число высокочастотных составляющих. Число гармоник под первым лепестком при N 4 (рис. 6.12) равно трем, в то время как при N 2 имеется лишь одна гармоника (рис. 6.14).

6.2.1. Определение спектра сигнала через быстрое преобразование Фурье

Вычисление спектра сигнала связано с вычислением интегралов (5.8 – 5.10), подынтегральные функции в которых быстро осциллируют, что существенно затрудняет вычисление таких интегралов с заданной точностью и ведет к значительным затратам времени. В связи с этим были разработаны специальные методы быстрого (или дискретного) преобразования Фурье (БПФ или ДПФ) [10], а в дальнейшем, применительно к ним, созданы специальные программные средства в системе Mathcad [4], использование которых в сотни раз ускоряет процедуру определения спектра сигнала и особенно это заметно при воссоздании приближения функции рядом Фурье, то есть

159