Анализ периодических и непериодических сигналов

контрольная работа: Коммуникации и связь

Документы: [1]   Word-104964.doc Страницы: Назад 1 Вперед

Контрольная работа №1

Спектральный анализ периодического и непериодического управляющих сигналов


Дано:

Шифр сигнала тФА 4 из табл. 1[1];;

Длительность периода тФА Т = 0,001 с = 1000 мкс ;

Соотношение между периодом и длительностью импульса тФА Т = 3ПД


Рис. 1 - Периодический сигнал


Задание:

1.Выполнить математическое описание заданного периодического сигнала, изобразить графически 2-3 периода сигнала, указав на рисунке параметры.

Математическое описание заданного периодического сигнала


Рис. 2

В результате подстановки данных варианта получаем униполярные прямоугольные периодические импульсы.

Период сигнала : Т = 0,001 с = 1000 мкс ;

Длительность импульса:


ПД* = 2ПД = 2В· Т/3 = ;


Временной интервал между импульсами:


ПД = Т/3 =;


Четная симметрия относительно моментов времени


t = nВ·T/2, где n = 0,В±1, В±2, В±3тАж;

;


Скважность импульсов:




Анализ временных свойств сигнала и формулировка обоснованных предположений о свойствах и особенностях спектрального состава сигнала.

Сигнал является четной функцией времени

Сигнал представляет собой знакопостоянную последовательность импульсов. Постоянная составляющая ряда Фурье равна:



В разложении сигнала в ряд Фурье будут присутствовать только косинусоидальные гармонические составляющие, т.е.:



Ряд Фурье можно преобразовать следующим образом:



Вычисление спектров амплитуд и фаз. Характер огибающей спектра амплитуд.

Производим расчет весовых коэффициентов аn:



Амплитуды гармоник



Фазы гармоник



Результаты оформляем в виде таблицы.


n

0

1

2

3

4

5

6

7

8

9

10

an

-0,66667

0,5513

0,2757

0

-0,1378

-0,1103

0

0,0788

0,0689

0

-0,0551

bn


0

0

0

0

0

0

0

0

0

0

An

0,66667

0,5513

0,2757

0

0,1378

0,1103

0

0,0788

0,0689

0

0,0551

ПЖn


0

0

-

-ПА

-ПА

-

0

0

-

-ПА

0,66667

0,27565

0,13785

0

0,0689

0,05515

0

0,0394

0,03445

0

0,02755






Построение оценки сигнала



Для N=4 :



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


t/T

0

-0,66667

0,5513

0,2757

-0,1378

0,022494

0,1

-0,66667

0,4460

0,0852

0,1115

-0,02394

0,2

-0,66667

0,1704

-0,2230

-0,0426

-0,76191

0,3

-0,66667

-0,1704

-0,2230

-0,0426

-1,10265

0,4

-0,66667

-0,4460

0,0852

0,1115

-0,91601

0,5

-0,66667

-0,5513

0,2757

-0,1378

-1,08016

0,6

-0,66667

-0,4460

0,0852

0,1115

-0,91601

0,7

-0,66667

-0,1704

-0,2230

-0,0426

-1,10265

0,8

-0,66667

0,1704

-0,2230

-0,0426

-0,76191

0,9

-0,66667

0,4460

0,0852

0,1115

-0,02394

1

-0,66667

0,5513

0,2757

-0,1378

0,022494



Расчет относительного значения среднеквадратической погрешности представления сигнала оценкой из гармонических колебаний.

Квадрат сигнала численно равен мгновенной мощности, рассеиваемой на сопротивлении нагрузки 1Ом. Средняя мощность сигнала прямо пропорциональна энергии, запасакмой за период и обратно пропорциональна периоду.



Мощность n-ного гармонического сигнала:



Уравнение погрешности представления периодического сигнала усеченным рядом Фурье:



n

0

1

2

3

4

5

6

7

8

9

10

An

0,66667

0,5513

0,2757

0

0,1378

0,1103

0

0,0788

0,0689

0

0,0551

Pn

0,4444

0,1520

0,0380

0,0000

0,0095

0,0061

0,0000

0,0031

0,0024

0,0000

0,0015

PN

0,4444

0,5964

0,6344

0,6344

0,6440

0,6500

0,6500

0,6531

0,6555

0,6555

0,6570

δ,%

33,3340

10,5367

4,8374

4,8374

3,4051

2,4932

2,4932

2,0279

1,6717

1,6717

1,4438


По полученным результатам строим график


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

Рассмотрим непериодический сигнал s1(t), изображенный на рисунке.



Его спектральная плотность



Сигнал s2(t) образован суммой сигналов s1(t), один из которых сдвинут вправо, а другой "ево на величину τ. Применяя теорему сдвига, получим:




Спектральная плотность - действительная функция частоты, т.к. мнимая составляющая равна нулю. Размерность спектральной плотности - ВтИЩс.

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

Учитывая, что , имеем :



Нули спектральной плотности находим, учитывая, что sin(πk)=0 и cos(π/2+πk)=0, k=0,±1,±2тАж



По горизонтальной оси откладываем номер гармоники, основной частоты:



Аргумент спектральной плотности будет равен:



0

0,25

0,5

1

1,5

1,75

2

2,25

2,5

3

3,5

3,75

4

4,25

4,5

5

5,5

5,75

6

6,25

6,5

0

0,375

0,75

1,5

2,25

2,625

3

3,375

3,75

4,5

5,25

5,625

6

6,375

6,75

7,5

8,25

8,625

9

9,375

9,75

-6,6667Е-4

-0,00046

-3,7E-20

0,000424

3,68E-20

-6,6E-05

-2,6E-20

5,1E-05

1,43E-19

-0,00014

-3,7E-20

3,06E-05

2,6E-20

-2,7E-05

-3,7E-20

8,49E-05

-6E-20

-2E-05

-2,6E-20

1,84E-05

1,19E-19

-0,6667

-0,4594

0,0000

0,4244

0,0000

-0,0656

0,0000

0,0510

0,0000

-0,1415

0,0000

0,0306

0,0000

-0,0270

0,0000

0,0849

0,0000

-0,0200

0,0000

0,0184

0,0000







Сопоставление спектров периодического и непериодического сигналов.

Сопоставление спектров произведем на основании соотношения





Сравнение спектров периодического и непериодического сигналов показывает, что гармоники, построенные на частотах, кратных ω1, и ограниченные спектральной плотностью непериодического сигнала со значениями Сn на спектральных диаграммах совпадают.

Определение энергии и средней мощности заданного сигнала на участке цепи с сопротивлением 1 Ом.

Определим энергию сигнала по временному представлению.


Контрольная работа №2

Расчет прохождения периодических и непериодических сигналов через линейные электрические цепи первого порядка


Дано:

Шифр периодического сигнала s1 тФА 4 из табл. 3[1];


Рис. 3


После подстановки значений параметров и масштабирования, получаем:


Рис. 4


Длительность периода тФА Т = 0,001 с = 1000 мкс ;

Соотношение между периодом и длительностью импульса тФА Т = 3ПД

Соотношение параметров цепи и сигнала:



Шифр цепи - 2 из табл. 4[1];


Рис. 5


Значения сопротивлений из табл. 1[1] - R1 = 2R; R2 = R

Задание:

Рассчитать и построить в масштабе АЧХ и ФЧХ интегрирующей и дифференцирующей цепей в диапазоне от нуля до 10 кГц, полагая (по шкале абiисс сделать градуировку частоты в кГц и в безразмерных величинах ωτц);

Рассчитать и построить в масштабе переходную и импульсную характеристики цепей от нуля до tmax = 3τц (по шкале абiисс сделать градуировку времени в мкс и в безразмерных величинах t/τц);

Проверить выполнение предельных соотношений между частотными и импульсными характеристиками.

Рассчитать спектр амплитуд и фаз на выходе заданной цепи при действии периодического сигнала s1(t).

Построить с учетом масштаба на общей спектрограмме спектры амплитуд и фаз входного и выходного сигналов при действии сигнала s2(t).

Дать представление входного сигнала с помощью функций Хевисайда.

Получить динамическое представление отклика заданной цепи на действие сигнала s2(t)(с помощью переходных характеристик).

Изобразить отклик цепи на интервале времени от нуля до tmax, в три раза превышающем длительность воздействия сигнала s2(t) .

Сделать выводы по результатам проведенного анализа.

Расчет частотных характеристик интегрирующей и дифференцирующей цепей.

Выполнение пунктов 1-3 задания оформляем в виде таблицы.


Табл. 1 - Анализ дифференцирующей и интегрирующей цепей.

Дифференцирующая цепь

Интегрирующая цепь

Рис. 6

Рис. 7

Вводим оператор дифференцирования р, такой, что

Передаточный коэффициент цепи:

Передаточный коэффициент цепи:

Находим комплексный передаточный коэффициент, заменяя р на jω


После анализа цепей находим частотные характеристики.


Табл. 2 - Частотные характеристики цепей

Дифференцирующая цепь

Интегрирующая цепь

Амплитудно-частотная (АЧХ)

Фазочастотная (ФЧХ)

f,Гц

0

100

700

2000

5000

10000

ω,рад/с

0

628,3

4398,2

12566,4

31415,9

62831,9

ωτц,рад

0,000

0,043

0,303

0,866

2,165

4,329

К(ω)

0,000

0,262

0,885

0,984

0,997

~ 1

φ(ω),рад

-1,57

-1,31

-0,48

-0,18

-0,07

~0

φ(ω),град

-90

-74,8

-27,7

-10,4

-4,2

~0

f,Гц

0

100

700

2000

5000

10000

ω,рад/с

0

628,3

4398,2

12566,4

31415,9

62831,9

ωτц,рад

0,000

0,043

0,303

0,866

2,165

4,329

К(ω)

1

0,96

0,46

0,18

0,07

~ 0

φ(ω),рад

0,00

-0,27

-1,09

-1,39

-1,50

~ -1,57

φ(ω),град

0

-15,2

-62,3

-79,6

-85,8

~ -90

Рис. 8. АЧХ ДЦ

Рис. 9. АЧХ ИЦ

Рис. 10. ФЧХ ДЦ

Рис. 11. ФЧХ ИЦ


Расчет временных характеристик дифференцирующей и интегрирующей цепей.

Находим временные характеристики операторным методом, пользуясь значением операторного коэффициента найденного в пункте 1.


Табл. 3 - Временные характеристики цепей

Дифференцирующая цепь

Интегрирующая цепь

Импульсные

Переходные

t,мкс

0

200

500

700

1000

1298,7

t/ ПДц

0

0,462

1,155

1,617

2,31

3

h(t)

-2310

-1455,3

-727,8

-458,52

-229,29

~ 0

g(t)

1

0,63

0,32

0,20

0,10

~ 0

t,мкс

0

200

500

700

1000

1298,7

t/ ПДц

0

0,462

1,155

1,617

2,31

3

h(t)

2310

1455,35

727,78

458,52

229,29

~ 0

g(t)

0

0,37

0,68

0,80

0,90

~ 1

Рис. 12. ИХ ДЦ

Рис. 13 ИХ ИЦ

Рис. 14. ПХ ДЦ

Рис. 15. ПХ ИЦ


Проверка соотношений между частотными и временными характеристиками дифференцирующей и интегрирующей цепей.

Предельные соотношения



Табл. 4 - Предельные соотношения

Дифференцирующая цепь

Интегрирующая цепь


Расчет спектра амплитуд и фаз на выходе заданной цепи при действии периодического сигнала s1(t).

По известной формуле из теории четырехполюсников находим передаточный коэффициент заданной цепи:



Находим комплексный передаточный коэффициент, заменяя р на jω



Найдем спектральную плотность непериодического сигнала на входе(s2(t) ) и выходе (s(t)) цепи, соответствующему периодическому сигналу s1(t) на протяжении одного периода.

Спектральная плотность непериодического сигнала s2(t)(см. к.р.№1):



Спектральную плотность выходного сигнала s(t) найдем по формуле:



Учитывая, что


и ,

а также


Учитывая, что: или, получаем:

Спектр амплитуд выходного периодического сигнала s(t) :



Спектр фаз:



Табл. 5 - Спектры входного и выходного периодического сигналов

n

0

A0=|ao/2|=0,667

1

A0=|ao/2|=0,667

1,04575

1

0,551

0

0,361

-0,162

2

0,276

0

0,130

-0,152

3

0

-

0

-

4

0,138

-1

0,054

-1,102

5

0,110

-1

0,042

-1,085

6

0

-

0

-

7

0,079

0

0,029

-0,063

8

0,069

0

0,026

-0,055

9

0

-

0

-

10

0,055

-1

0,021

-1,045


Построение спектров амплитуд и фаз входного(s2(t) ) и выходного(s(t)) непериодического сигналов

В соответствии с пунктом 2, имеем:

Амплитудные спектры



Учитываем, что при ω=0



Амплитуда - четная функция частоты

Фазные спектры



Где функция sign(ω)=1 при ω>0 и sign(ω)=-1 при ω<0

Фаза - нечетная функция частоты


Табл. 6 - Спектры входного и выходного непериодического сигналов

-10

0,0276

1

0,0094

1,04575

-9,75

0

-

0

-

-9,375

0,0184

0

0,0063

0,048648

-9

0

-

0

-

-8,625

0,0200

1

0,0068

1,052634

-8,25

0

-

0

-

-7,5

0,0849

0

0,0294

0,059951

-6,75

0

-

0

-

-6,375

0,0270

1

0,0095

1,069491

-6

0

-

0

-

-5,625

0,0306

0

0,0109

0,077593

-5,25

0

-

0

-

-4,5

0,1415

1

0,0518

1,093529

-3,75

0

-

0

-

-3,375

0,0510

0

0,0199

0,115996

-3

0

-

0

-

-2,625

0,0656

1

0,0275

1,135528

-2,25

0

-

0

-

-1,5

0,4244

0

0,2259

0,164845

-0,75

0

-

0

-

-0,375

0,4594

1

0,4151

1,096523

0

0,6667

-

0,6667

-

0,375

0,4594

-1

0,4151

-1,096523

0,75

0

-

0

-

1,5

0,4244

0

0,2259

-0,164845

2,25

0

-

0

-

2,625

0,0656

-1

0,0275

-1,135528

3

0

-

0

-

3,375

0,0510

0

0,0199

-0,115996

3,75

0

-

0

-

4,5

0,1415

-1

0,0518

-1,093529

5,25

0

-

0

-

5,625

0,0306

0

0,0109

-0,077593

6

0

-

0

-

6,375

0,0270

-1

0,0095

-1,069491

6,75

0

-

0

-

7,5

0,0849

0

0,0294

-0,059951

8,25

0

-

0

-

8,625

0,0200

-1

0,0068

-1,052634

9

0

-

0

-

9,375

0,0184

0

0,0063

-0,048648

9,75

0

-

0

-

10

0,0276

-1

0,0094

-1,04575


Рис. 16 - Амплитудный спектр входного и выходного периодического сигналов.

Рис. 17 - Фазовый спектр входного и выходного непериодического сигналов.


Представление входного непериодического s2(t) с помощью единичной функции σ(t) (функции Хевисайда).

Входной сигнал является суммой функций Хевисайда, сдвинутых по временной оси и умноженных на амплитуду сигнала Е.



Рис. 18 - Функции Хевисайда


Динамическое представление отклика заданной цепи на действие сигнала s2(t).

Заданная цепь является суммой интегрирующего и дифференцирующего звеньев.



Находим переходную характеристику заданной цепи



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


, где


Динамическое представление отклика:



Табл. 7 - Отклик заданной цепи на действие непериодического сигнала

t/T

-1/2-0

-1/2+0

-1/3

-1/6-0

-1/6+0

0

1/6-0

1/6+0

1/3

1/2-0

1/2+0

2/3

5/6

1

7/6

4/3

3/2

5/3

11/6

2

0

-0,33307

-0,6912

-0,84558

-0,85702

-0,9338

-0,9669

-0,96935

-0,98581

-0,9929

-0,99343

-0,99696

-0,99859

-0,99935

-0,9997

-0,99986

-0,99994

-0,99997

-0,99999

-1

0

0

0

0

0,33328

0,6913

0,845627

0,857068

0,933821

0,966905

0,969358

0,985812

0,993431

0,996958

0,998592

0,999348

0,999698

0,99986

0,999935

1

0

0

0

0

0

0

0

-0,33335

-0,69133

-0,84564

-0,85708

-0,93383

-0,96936

-0,98581

-0,99343

-0,99696

-0,99859

-0,99935

-0,9997

-1

0

0

0

0

0

0

0

0

0

0

0,333729

0,691508

0,857164

0,933865

0,969379

0,985822

0,993435

0,99696

0,998593

1

0

-0,33307

-0,6912

-0,84558

-0,52374

-0,2425

-0,12127

-0,44563

-0,74332

-0,87164

-0,54743

-0,25347

-0,11736

-0,05434

-0,02516

-0,01165

-0,00539

-0,0025

-0,00116

0



Рис. 19 - Изображение входного непериодического сигнала s2 и отклик цепи на него.


Выводы.

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

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

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

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

Операторный метод и теория обобщенных функций дает достаточно мощный аппарат для исследования цепей и сигналов.

Контрольная работа №3

Расчет прохождения непериодического сигнала сложной формы через линейную цепь второго порядка


Дано:

Шифр входного сигнала s(t) тФА N1= 4; N2= 3 из табл. 2[1];

N2= 3 - Соотношение параметров импульса : t2 = 1,5t1 из табл. 5[1];

N1= 4 - Номер рисунка из табл. 5[1] - 4;


Рис. 20


Шифр цепи - N3 N4= 44 из табл. 3[1];

Номер рисунка из табл. 6[1] N3 N4= 44;


Рис. 21


Соотношение параметров цепи и сигнала:



Задание:

Рассчитать и построить в масштабе АЧХ и ФЧХ цепи

Рассчитать и построить в масштабе переходную и импульсную характеристики цепи;

Проверить выполнение предельных соотношений между частотными и импульсными характеристиками.

Дать поинтервальное аналитическое представление сигнала по его графику;

Рассчитать операторным методом или методом временного интегрирования отклик на выходе линейной цепи и дать его поинтервальное описание

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

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

Расчет частотных характеристик линейной цепи второго порядка

Для того, чтобы найти операторную передаточную функцию, преобразуем схему


Рис. 22

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


Рис. 23


На холостом ходу:


Рис. 24


Находим операторный передаточный коэффициент:



Находим корни знаменателя К(р)



Получаем



Получаем комплексный передаточный коэффициент заменяя р на jω


АЧХ цепи:

ФЧХ цепи:


Табл. 8 - Расчет АЧХ и ФЧХ цепи

f,кГц

0

0,80

1,59

2,25

3,18

4,77

7,96

15,92

31,83

ω,рад/с

0

5000

10000

14142,14

20000

30000

50000

100000

200000

ωτц,рад

0

0,5

1

2

3

5

10

20

К(ω)

0

0,22

0,32

0,33

0,32

0,26

0,18

0,10

0,05

0

φ(ω)

900

49,400

18,430

0,000

-18,430

-37,870

-56,890

-72,980

-81,430

-900


Рис. 25 - АЧХ цепи второго порядка


Рис. 26 - ФЧХ цепи второго порядка


Расчет временных характеристик линейной цепи второго порядка

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



Импульсная характеристика интегрирующей цепи (см. к.р.№2)



Переходная характеристика интегрирующей цепи (см. к.р.№2)



Импульсная характеристика заданной цепи:



Переходная характеристика заданной цепи:



Табл. 9 - Расчет временных характеристик

t,мкс

0

10

20

50

60

69,315

80

100

300

500

1000

t/ ПДц

0

0,1

0,2

0,5

0,6

0,69315

0,8

1

3

5

10

h(t)

10000

7326,2

5219,1

1292,3

535,8

0,0

-455,4

-972,1

-448,3

-66,5

-0,5

g(t)

0,000

0,086

0,148

0,239

0,248

0,250

0,247

0,233

0,047

0,007

0,000


Рис. 27 - Импульсная характеристика цепи

Рис. 28 - Переходная характеристика цепи


Предельные соотношения между частотными и временными характеристиками.



Поинтервальное описание входного сигнала

Обозначим длительность сигнала - Т.

Тогда



Рис. 29 - Входной сигнал


На интервале t=[0; 0,4T]

На интервале t=[0,4T; 0,6T]

На интервале t=[0,6T; T]

Используя функцию включения можно записать представление входного сигнала в виде:



Расчет отклика на выходе цепи

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



На интервале t=[0; 0,4T]



На интервале t=[0,4T; 0,6T]



На интервале t=[0,6T; T]



На интервале


t=[T; ] , где:


Найдем отклик операторным методом. Входной сигнал равен:


.


Пользуясь преобразованием Лапласа


,


Обратное преобразование по Лапласу выражений вида:



Обратное преобразование по Лапласу выходного сигнала:



Окончательно:



На интервале t=[0; 0,4T]



На интервале t=[0,4T; 0,6T]



На интервале t=[0,6T; T]



На интервале t=[T; ]


,


Построение графика отклика сигнала на выходе цепи


Табл. 10 - Расчет отклика сигнала

t, мкс

0

50

100

150

200

250

300

350

400

450

500

550

600

650

700

750

t/T

0

0,2

0,4

0,6

0,8

1

1,2

1,4

1,6

1,8

2

2,2

2,4

2,6

2,8

3

Sв(t)

0

0,077

0,200

0,224

0,097

-0,080

-0,147

-0,125

-0,089

-0,059

-0,038

-0,023

-0,014

-0,009

-0,005

-0,003


Рис. 30 - Входной и выходной сигналы


Выводы

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

Сравнение временного и операторного методов для нахождения отклика заданной цепи говорит в пользу последнего, так как операторный метод позволяет более полно охватывать процесс вычислений и является более емким.

Отклик на входной сигнал данной цепи второго порядка является апериодическим, так как корни знаменателя операторного передаточного коэффициента действительные.

По графику отклика можно судить об интегрирующей природе цепи.

Страницы: Назад 1 Вперед