Анализ линейной электрической цепи во временной и частотной областях

Тип работы:
Курсовая
Предмет:
Физика


Узнать стоимость

Детальная информация о работе

Выдержка из работы

СОДЕРЖАНИЕ

Введение

Техническое задание

1. Спектральный анализ

1.1. Вычисление переходной характеристики цепи

1.2. Определение реакции цепи на импульс заданной формы с помощью интеграла Дюамеля

1.3. Вычисление спектра сигнала на входе и на выходе цепи

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

2. Временной анализ

2.1. Вычисление дискретного сигнала на выходе цепи U2(n)

2.2. Спектральные характеристики дискретизированного сигнала

2.3. Синтез схемы дискретной цепи

2.4. Корректор

Заключение

Список литературы

ВВЕДЕНИЕ

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

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

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

ТЕХНИЧЕСКОЕ ЗАДАНИЕ

L

R R

R

Рисунок 1. Схема анализируемой цепи

Рисунок 2. Входной импульс

Задание на курсовую содержит схему анализируемой цепи и входной сигнал в виде одиночного импульса, параметры которого указаны на рисунке. Все резисторы схемы имеют сопротивление R=1КОм, индуктивность катушек L=1Гн, ёмкость конденсаторов С= 1 мкФ.

В процессе выполнения работы необходимо получить следующие результаты:

1. Рассчитать и построить график напряжения на выходе цепи u2(t);

2. Вычислить передаточную функцию цепи H (jщ), спектральную плотность сигнала на входе и выходе цепи, построить соответствующие графики;

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

4. провести дискретизацию входного сигнала и импульсной характеристики цепи, вычислить отсчёты дискретного сигнала на входе и построить графики спектра дискретизированного сигнала и АЧХ дискретизированной цепи;

5. по отсчётам входного сигнала вычислить его спектральную плотность;

6. составить схему дискретной цепи, выполнив Z-преобразование дискретизированной импульсной характеристики;

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

1. СПЕКТРАЛЬНЫЙ АНАЛИЗ

1.1 Вычисление переходной характеристики цепи

Вычисляется переходная характеристика цепи, как реакция на входное воздействие в виде единичной функции 1(t).

(1. 1)

(1. 2) при t> 0-

рисунок 1.1. Схема в момент времени t = 0-

IL (0-)= IL (0+)=0

В момент времени t=0- индуктивность заменяем на провод.

при t=0+

рисунок 1.2. Схема в момент времени t = 0+

В момент времени t=0+ индуктивность заменяем на разрыв. Тогда напряжение на выходе:

IR3= = =0.5 мА

UR2=0

UR3=0.5 В

При t=?

рисунок 1.3. Схема в момент времени t = ?

Rэкв= + R3= +1 = 1.5 кОм

IR3= = = 0. 667 мА

IR1=IR2=0. 5*IR3=0. 3335 Ма

UR2=0. 3335 В

UR3=0. 667 В

p= = -1500

U2= g (t)=Uпр+Uсв = Uпр +(U (0+)-Uпр)=1 — 0. 5

Далее, используя переходную характеристику, находим реакцию цепи с помощью интеграла Дюамеля.

1. 2 Определение реакции цепи на импульс заданной формы с помощью интеграла Дюамеля

Весь отрезок времени 0? t<? разбивается на два интервала. Границы интервалов приходятся на моменты времени t0 = 0, t1 = 4мс.

U11, 0? t< t1

U (t) =

U12, t? t2

-2500t+5, 0? t< t1

U (t) =

0, t? t2

U11`= -2500; U12`= 0.

U21(t) = U11(0)*g (t) + =5*(1 — 0.5 + =5 — 2.5 —

— 2500 = 5 — 2.5 — 2500t + (-) = 5 — 2500t — 2.5 + 0. 833 — 0. 833 =

5. 833 — 2500t — 3. 333

U22(t) = U11(0)*g (t) + + (U12(t2) — U11(t-2)) * g (t — t2) + += 5(

Таблица 1.1. Результаты вычислений методом интеграла Дюамеля

t, мс

0

0,2

0,4

0,6

0,8

1

1,2

1,4

1,6

1,8

U2, В

2. 5

2. 864

3. 004

2. 978

2. 83

2. 59

2. 282

1. 925

1. 531

1. 109

t, мс

2

2,2

2,4

2,6

3

3,2

3,4

3,6

3,8

4

U2, В

0. 667

0. 21

-0. 258

-0. 735

-1. 217

-1. 704

-2. 194

-2. 687

-3. 182

-3. 678

Рисунок 1.4. График зависимости выходного сигнала от времени

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

Выполним проверку:

t=1 мс

U21=

t=4 мс

U22=

1.3 ВЫЧИСЛЕНИЕ СПЕКТРА СИГНАЛА НА ВХОДЕ И НА ВЫХОДЕ ЦЕПИ

Для нахождения спектральной плотности входного сигнала функция U1(t) представляется в виде суммы четырех «простейших» функций:

Рисунок 1.5. Входной импульс

Рисунок 1.6. Разложение входного сигнала на простейшие

f1(t) = 5 > F1(p) =;

f2(t) =; f2(t) > F2(p) = -;

f3(t) =; f3(t) > F3(p) =;

f4(t) =; f4(t) > F4(p) =

Изображение входного сигнала записывается как сумма изображений «простейших» функций:

F (p) = F1(p) + F2(p) + F3(p) + F4(p) = + + =

U1(jщ) = = =

= (1. 3)

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

U1(щ) = (1. 4)

Фазовая характеристика спектральной плотности входного сигнала:

ц1(щ) =arctg ((1. 5)

Передаточная функция по напряжению цепи

Рисунок 1.7. Схема анализируемой цепи

H (jщ)= (1. 6)

H (jщ) = = (1. 7)

Амплитудно-частотная характеристика:

H (щ) = (1. 8)

Фазо-частотная характеристика:

Ф (щ) = (1. 9)

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

U2(щ) = U1(щ)*H (щ) (1. 10)

Фазовая характеристика спектральной плотности сигнала на выходе цепи записывается в виде:

ц2(щ) = ц1(щ)*Ф (щ) (1. 11)

Вычисление модулей и аргументов спектральных плотностей на входе и выходе цепи, а так же АЧХ и ФЧХ ее производится с помощью программы «FREAN».

Таблица 1.2. Результаты вычислений с помощью программы «FREAN»

F, кГц

U1(щ), мВ*с

ц1(щ), градус

H (щ)

Ф (щ)

U2(щ), мВ*с

ц2(щ), градус

0

0

0

1

0

0

0

0,1

7,125

342

0,942

-10,9

6,714

331,1

0,2

8,302

53,91

0,831

-17,23

6,9

36,68

0,3

3,474

125,9

0,736

-19,36

2,555

106,6

0,4

1,975

17,86

0,669

-19,23

1,32

-1,363

0,5

3,186

89,82

0,624

-18,17

1,988

71,65

0,6

0,495

161,8

0,594

-16,83

0,294

144,9

0,7

1,987

53,73

0,573

-15,48

1,139

38,25

0,8

1,501

125,7

0,558

-14,23

0,837

111,5

0,9

0,686

17,68

0,547

-13,1

0,375

4,578

1

1,593

89,64

0,539

-12,11

0,859

77,53

1,2

1,12

53,54

0,528

-10,46

0,591

43,09

1,4

0,404

17,5

0,521

-9,163

0,21

8,335

1,6

0,27

161,4

0,516

-8,137

0,139

153,3

1,8

0,699

125,4

0,513

-7,308

0,358

118,1

2

0,797

89,27

0,51

-6,627

0,407

82,64

2,5

0,637

89,09

0,507

-5,363

0,323

83,72

3

0,531

88,9

0,505

-4,498

0,268

84,41

3,5

0,455

88,72

0,503

-3,87

0,229

84,85

Рисунок 1.8. Амплитудная характеристика на входе цепи

Рисунок 1.9. Фазовая характеристика на входе цепи

Рисунок 1. 10. Амплитудно-частотная характеристика цепи

Рисунок 1. 11. Фазо-частотная характеристика цепи

импульс цепь дискретный сигнал

Рисунок 1. 12. Амплитудная характеристика на выходе цепи

Рисунок 1. 13. Фазовая характеристика на выходе цепи

Выполним проверку:

F=1 кГц ==> щ=6280 рад/с

U1(щ) =

H (щ) =

U2(щ) = U1(щ)*H (щ) = 1,603*0,539 = 0,864 мВ*с

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

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

Вычислим импульсную характеристику цепи:

h (t) = g (o)*д (t) + g`(t) (1. 12)

h (t) = 0. 5*д (t) + 750 (1. 13)

H (jщ) = dt (1. 14)

H (jщ) =)dt = 0.5 =

= (1. 15)

Полученный результат совпадает с результатом, полученным в пункте 1.3.

2. ВРЕМЕННОЙ АНАЛИЗ

2.1 Вычисление дискретного сигнала на выходе цепи U2(n)

Максимум модуля спектральной плотности среди значений U1(n):

U1max = 8,302 мВ? с.

Мы выбрали верхнюю границу спектра входного сигнала fв = 2,5 кГц. Тогда частота дискретизации берется равной fд = 5кГц. Соответственно период дискретизации Тд = = 0,2 мс.

Составляем аналитическое выражение для U1(t)

0, t < 0

5 — 2500t, 0? t < t1

U1(t) = - 2500(t — t1), t1? t < t2

0, t > t2

Подставляя вместо t последовательность моментов дискретизации, вычисляем значения дискретных отсчетов входного сигнала U1(n). Аналогичным образом вычисляются значения дискретных отсчетов импульсной характеристики цепи H (n) на интервале времени 0? t < t2.

Дискретные значения импульсной характеристики вычисляются по формуле:

H (n) = 0. 5*д (n) + Т*750 (2. 1)

Таблица 2.1. Дискретные значения функции входного сигнала и импульсной характеристики

t, мс

0

0,2

0,4

0,6

0,8

1

1,2

1,4

1,6

1,8

n

0

1

2

3

4

5

6

7

8

9

U1(n)

5

4,5

4

3,5

3

2,5

2

1,5

1

0,5

H (n)

0,5

0,11 112

0,8 232

0,6 099

0,4 518

0,3 347

0,2 479

0,1 837

0,1 361

0,0101

t, мс

2

2,2

2,4

2,6

2,8

3

3,2

3,4

3,6

3,8

4

n

10

11

12

13

14

15

16

17

18

19

20

U1(n)

0

-0,5

-1

-1,5

-2

-2,5

-3

-3,5

-4

-4,5

-5

H (n)

0,747

0,0055

0,004

0,003

0,0023

0,002

0,001

0,0009

0,68

0,0005

0,37

U2(m) = (2. 2)

Таблица 2.2. Дискретные значения функции выходного сигнала

n

0

1

2

3

4

5

6

7

8

9

t, мс

0

0,2

0,4

0,6

0,8

1

1,2

1,4

1,6

1,8

U2(n)

2,5

2,80 561

2,91 166

2,86 987

2,71 855

2,48 609

2,19 352

1,85 643

1,48 635

1,9 183

n

10

11

12

13

14

15

16

17

18

19

t, мс

2

2,2

2,4

2,6

2,8

3

3,2

3,4

3,6

3,8

U2(n)

0,6792

0,25 317

-0,1828

-0,6261

-1,0749

-1,528

-1,984

-2,44 161

-2,90 128

-3,36 218

Рисунок 2.1. Значение дискретного сигнала на входе цепи

Рисунок 2.2. Значение дискретного сигнала на выходе цепи

Рисунок 2.3. Значение дискретных отсчетов импульсной характеристики цепи

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

2. 2 Спектральные характеристики дискретизированного сигнала

Спектральные характеристики дискретизированного сигнала U1(n) рассчитываются по формуле:

U1(jщ) = T* (2. 3)

На частотах:

щ = 0 => f = 0

щ = р/4Т = 3925 рад/с; => f = 0. 625 кГц

щ = р/2Т = 7850 рад/с; => f = 1,250 кГц

щ = р/Т = 15 700 рад/с => f = 2,5 кГц

Таблица 4.3. Спектральная плотность дискретизированного сигнала

F, кГц

0

0. 625

1,25

2,5

U1(щ), мВ*с

0

1,332

1,01

0,0173

Рисунок 2.4 График спектральной плотности дискретизированного сигнала

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

2.3 Синтез схемы дискретной цепи

Z — преобразование импульсной характеристики цепи записывается в виде:

H (z) = 0.5 + = (2. 4)

Y (z) = X (z)*H (z) = X (z) *

Y (z) *() = X (z) * (

Y (z) = 0. 65X (z) — 0. 37X (z) + 0. 74Y (z)

a0 = 0. 65; a1 = - 0. 37; b1 = 0. 74.

Рисунок 2.5. Схема дискретной цепи

X (n) Y (n)

Рисунок 2.6. Канонический вид схемы дискретной цепи

Коэффициенты передачи масштабных усилителей a0, а1, b1 те же, что и в предыдущей схеме. T — элемент памяти с задержкой на один период дискретизации.

Передаточная функция дискретной цепи

H (jщ) = (2. 5)

АЧХ дискретной цепи имеет вид:

H (щ) = (2. 6)

Таблица 2.4. Значения АЧХ дискретной цепи

f, кГц

0

1

2

3

4

5

H (щ)

1,7 692

0,89 148

0,65 578

0,65 517

0,89 105

1,0769

Рисунок 2.5. Амплитудно-частотная характеристика дискретной цепи

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

2.4 Корректор

Z — преобразование передаточной функции корректора H'(Z) находится как величина, обратная H (Z) исходной цепи, вычисленной в пункте 4.3.

H`(z) = (2. 7)

H`(z) = = (2. 8)

Y (z) = 1. 538X (z) — 1. 138X (z) + 0. 569Y (z)

a0 = 1. 538; a1 = - 1. 138; b1 = 0. 569.

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

U2`= (2. 9)

Отсчеты импульсной характеристики корректора находятся путем деления полинома числителя H'(Z) на его знаменатель.

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

n

0

1

2

3

4

5

6

7

8

9

t, мс

0

0,2

0,4

0,6

0,8

1

1,2

1,4

1,6

1,8

U2(n)

2,5

2,806

2,912

2,87

2,719

2,486

2,194

1,856

1,486

1,092

H'(n)

1,538

-0,263

-0,15

-0,085

-0,048

-0,027

-0,015

-0,009

-0,005

-0,003

U2'(n)

3,845

3,658

3,365

3,015

2,631

2,228

1,815

1,393

0,968

0,54

n

10

11

12

13

14

15

16

17

18

19

t, мс

2

2,2

2,4

2,6

2,8

3

3,2

3,4

3,6

3,8

U2(n)

0,6792

0,25 317

-0,1828

-0,6261

-1,0749

-1,528

-1,984

-2,442

-2,9

-3,36

H'(n)

-0,002

-0,0011

-0,0006

-0,0003

-0,0002

0

0

0

0

0

U2'(n)

0,10 972

-0,3209

-0,7518

-1,1828

-1,6139

-2,045

-2,476

-2,907

-3,34

-3,77

Схема дискретной цепи, реализующая функцию корректора в соответствии с формулой H`(z) = в канонической форме имеет вид:

X (n) Y (n)

Рисунок 2.6. Канонический вид схемы дискретного корректора

Рисунок 2.7. Дискретные значения импульсной характеристики корректора

Рисунок 2.8. Дискретные значения сигнала на выходе корректора

Амплитудно-частотная характеристика корректирующей цепи

H`(jщ) = (2. 10)

H`(щ) = (2. 11)

Таблица 2.6. Значения АЧХ корректирующей цепи

f, кГц

0

1

2

3

4

5

H (щ)

0,92 807

1,20 804

1,60 551

1,60 642

1,20 887

0,92 809

Рисунок 2.9. Амплитудно-частотная характеристика корректирующей цепи

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

ЗАКЛЮЧЕНИЕ

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

Этот анализ заключался в следующем:

Мы вычислили переходную характеристику и рассчитали с ее помощью интеграл Дюамеля.

Получили аналитические выражения для АЧХ и ФЧХ цепи, спектральных плотностей входного и выходного сигналов, а так же вычислили их значения и получили их графическое изображение при помощи программы «FREAN».

Применили формулу преобразования Фурье, чтобы проверить связь между временными и частотными характеристиками, при этом результат вычисления совпал с формулой Н (jщ) полученной при нахождении передаточной функции цепи.

Определил частоту и период дискретизации входного сигнала. Затем при помощи теоремы свертки определил дискретный сигнал на выходе цепи. Далее используя аналогию между Z — преобразованием и преобразованием Фурье вычислил значения спектральной плотности входного дискретного сигнала. При этом результат расчета практически совпал с результатами спектральной плотности аналогового сигнала.

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

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

СПИСОК ЛИТЕРАТУРЫ

1. Тихобаев В. Г. Расчет электрических цепей при импульсном воздействии (методические указания к курсовой работе). Новосибирск: СибГУТИ, 2001

2. Шебес М. Р. Каблукова М.В. Задачник по теории линейных электрических цепей. М.: Высшая школа, 1990

3. Бакалов В. П. Игнатов А.Н. Крук Б. И. «Основы теории электрических цепей и электроники». М.: Радио и связь, 1989

ПоказатьСвернуть
Заполнить форму текущей работой