Формулу (8.75) называют интегралом Дюамеля.

С помощью интеграла Дюамеля можно найти не только ток, но и любую другую физическую величину, например напряжение. В этом случае в формулу вместо переходной проводимости g(/) будет входить переходная функция h(t), если на входе цепи действует источник ЭДС

(напряжения), и переходное сопротивление /?(/), если на входе цепи дей ствует источник тока.

§ 8.54 Последовательность расчета с помощью интеграла Дюа­меля. Расчет с помощью интеграла Дюамеля проводят в четыре этапа:

1) определение переходной проводимости g(/) (переходной функции А(г)) для исследуемой цепи;

2) нахождение g(/ - т) (Л(/ - т)). С этой целью в формуле для g(t) заменяют / на (/ - т);

3) определение Для этого находят производную от заданного напряжения u(t) по времени /ив полученном выражении заменяют / на т;

4) подстановка найденных на этапах 1-3 функций в формулу (8.75), интегрирование по переменной т и подстановка пределов.

Пример 101. Найти ~ f(i) и и2 = f(t) при замыкании ключа на схеме рис. 8.40. Напряжение источника ЭДС «(/) = 100(1-е-а/) В; а~0,25с~}; Я» 0,5 Ом; /,|=1Гн;

М = 0.5 Гн.

Рис. 8.40

Решение. Переходная проводимость цепи, состоящей из последовательно включенных R и L.

g<0 = 1 (j _ Q-h<),

 


Л
где

рмуле (8.75) выпадает, так как ц(0) = 0. При этом

»'(/) = —100(1-е'я')= 100 о с'; dt

ц'(т) = 100ое"ят;

'](/)= р'(т) g{t - т) dx - |е'о: (l-e"A('”:))</t.

• к *

о о

При интегрировании учитываем, что о г г не зависит:

/,(/) = 200 (1 + е’°-5' - 2 е"025') А.

Напряжение на зажимах вторичной обмотки

м2(Г) = м = 50 (е“°-25' - е’°-5') В.

§ 8.55 Применение интеграла Дюамеля при сложной форме на­пряжения. Пусть напряжение u(t) изменяется во времени по сложному закону, например в соответствии с рис. 8.41. Начальное напряжение равно и(0). В интервале от / = 0 до /=/, напряжение плавно растет, и закон

Рис. 8.41

его изменения u^t). В момент t = t} оно ме­няется скачком от иа до uh, а затем снова изменяется, но уже по другому закону и2 (/) во времени. При t-t2 напряжение скачком уменьшается от ис до нуля.

Требуется найти ток в каждом из трех ин­тервалов времени. Под первым интервалом будем понимать интервал времени от t = 0 до t - z, (не включая скачка напряжения от иа до иь)\ под вторым — от до 12, включая скачок от иа до uhi но не включая скачок от ис до 0; под третьим — при t > t2, включая скачок от ис до 0.

Интегрирование по-прежнему проводим по т, понимая под t фикси­рованный момент времени, в который требуется найти ток. На основании принципа наложения ток в любой момент времени t определится как сум­ма токов от всех напряжений, воздействовавших на цепь до момента /.

В первый интервал времени

/(/) = и(0) g(t) + JWj'(t) g(r - г) dx.
о

Во второй интервал времени

h

i(t) = w(0) g(0 + Jw[(t) g(t- t) dx + (uh о

где слагаемое (ub -ua) g(l~/t) обусловлено скачком напряжения от иа и иь в момент времени

В третий интервал времени

/(Г) = w(0) g(t) + Ji/,'(т) g(t - т) dx + (uh - иа) g(t -/j) +
о

'j

+ /«гСО g(/-x)^T + (0-wj g(/-/2)‘

'i

Пример 102. В электрической цепи (см. рис. 8.40) в момент времени / = 0 замы кает - ся ключ и напряжение «(/) изменяется в соответствии с рис. 8.41; и(0) = 50 В. В первый интервал времени от / = 0 до / = /(=4с напряжение w((/) = 150-100е-<,/1, где а = 0,25 с"1. Во второй интервал времени от / = /1=4с до / = /2=6с и2(/) = 50 + 100е"с('_/|), где с = 0.4с'1. Параметры схемы (см. рис. 8.40) /? = 0,5Ом; L] = I Гн (вторичная цепь разомкнута).

Найти закон изменения тока во времени для обоих интервалов времени, а также значения тока it при г, равном 2 и 5 с.

Решение. В соответствии с § 8.54 переходная проводимость

g(f) = 1 (1 “ е’*'); b = у = 0.5 с"1; g(t - т) = ~ (1 - е-А(/т))-

К L к

В первый интервал времени и'(т) = 100 се'"1. Поэтому

/| (О = «(0) g(/) + jw'(T) g(t - х) dx =
о

= + [е'"т(1-е'А(,-т))<Уг =

R R *

о

= I00(!-e-0>s') + 200(l + e’°-s' -e'025').

Приг = 2 c /, = l00(l-e”l) + 200(l + e"1-2e-0,5) = 94,9 A.

Во второй интервал времени (включая скачок иь - иа = 36,9 В )

Ъ /

'1(0 = и(0) g(/) + Jw[(r) g(t - т) dx + (uh - u„) g(i J«2(T) g(t ~x)dx\
0 'I

^(т) = -Ю0се‘‘гсс'1;

/, U) = 100 (I - e '°-5') + 200 (0,632 - 1.718 e -A5') + (1 - e'03 >) -

0 j 5

—e’c' +Az£c-‘'i +e"c'i c-c('-'i))e-<'».

(b-c)R с c

При I = 5 c /) = 204.32 A.

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

Однако классический метод физически более прозрачен, чем опера­торный, в котором решение уравнений во многом формализовано.

Если при сравнении методов исходить из объема вычислительной работы, то решение уравнений первого, второго, а иногда и третьего порядков для источников постоянной (синусоидальной) ЭДС или тока целесообразно проводить классическим методом, а решение уравнений более высоких порядков — операторным. Объясняется это тем, что чем выше порядок характеристического уравнения, тем более громоздкой и трудоемкой оказывается операция нахождения постоянных интегрирова­ния в классическом методе. Операторный метод имеет перед классичес­ким явное преимущество при решении задач, в которых определение принужденной компоненты искомой величины оказывается затруднитель­ным вследствие сложного характера вынуждающей силы, а также при решении уравнений в частных производных (см. § 12.13-12.15). Если воздействующее напряжение изменяется во времени, например линейно или в виде всплеска одной или нескольких экспонент, рекомендуется применять операторный метод или интеграл Дюамеля. Но основной об­ластью применения интеграла Дюамеля являются случаи, когда напря­жение изменяется по сложному закону во времени, например при нали­чии скачков напряжения (см. § 8.55), или когда переходная проводимость g(r) и (или) воздействующее на схему напряжение заданы графически

(в последнем случае интеграл Дюамеля берется путем численного интег­рирования).

Рассматриваемый в § 8.66 метод расчета переходных процессов, получивший название метода пространства состояний, используется глав­ным образом, когда расчет осуществляется с применением ЭВМ. Для ручного счета этот метод громоздок.

Классический и операторный методы, а также метод пространства состояний в аналитической форме и интеграл Дюамеля имеют общий недостаток: необходимость определения всех корней характеристического уравнения, что для уравнений высоких степеней (например, 5-, 6-, 7-й, ...) требует много времени. В этих случаях может быть рекомендо­вано числовое решение на ЭВМ уравнений, составленных по методу про­странства состояний; может быть применен и спектральный метод в том виде, в каком он рассмотрен, например, в гл. 9. Кроме того, в этих слу­чаях используют моделирующие установки.

§ 8.57 Дифференцирование электрическим путем. Для четырехпо­люсников (рис. 8.42) при определенных условиях выходное напряжение и2(/) пропорционально производной от входного напряжения ^(Г),

 


 
а б Рис. 8.42

 

т. е. м2(/)« dux(f)ldt. Схему на рис. 8.42, а применяют чаще схемы на рис. 8.42, б, так как при практическом осуществлении она обладает мень­шими габаритами, массой и более удобна при регулировке.

Если u}(t)=U}(p), то du\(t)/dt = р Отсюда следует, что четы­рехполюсник осуществляет дифференцирование, если для него

Для схемы на рис. 8.42, а 6/2(р) " Ц(р)

Чтобы схема осуществила дифференцирование, необходимо выполнить условие | Л С р |« 1, тогда t/2(p)« RC р Ц(р). Для синусоидального процесса заменим р на j со и тогда схема на рис. 8.42, а будет выпол­нять свои функции, если со RC « I.

Аналогично доказывается, что для схемы на рис. 8.42, б необходимо выполнить условие (со L/ R) « 1. Если ux(t) — несинусоидальная пери­одическая функция, то эти условия должны выполняться для наивысшей частоты функции Wj(z).

При дифференцировании импульсных воздействий длительностью /н параметры схем на рис. 8.42 должны удовлетворять условиям RC « tH и L!R « гн. Эти условия получим из двух предыдущих, если в первом приближении будем считать, что поступление на вход четырехполюсни­
ка импульса длительностью гн соответствует воздействию на вход од­ной полуволны синусоиды частотой (о = 2л/(2/н) = л:/гн.

§ 8,58. Интегрирование электрическим путем. Для четырехполюс­

ников (рис. 8.43) при определенных условиях выходное напряжение

 


     
а б
 
Рис. 8.43

 

 

Схема на рис. 8.43, а предпочтительнее схемы на рис. 8.43, б по при­чинам, упомянутым в § 8.57.

Если то р](г) р' Отсюда следует, что схема

выполняет свои функции, если соотношение между ее параметрами обес­печивает выполнение соотношения (У2(/?) = U\{pVР-

JXnn схемы на рис. 8.43, a U2(p') = Ux(p')/(R С р +1), т. е. для нее дол­жно быть | R С р | » 1. Заменив р на j со, найдем условие w R С » 1, при котором схема на рис. 8.43, а будет выполнять функции интегриру­ющего звена при синусоидальном процессе. Для схемы на рис. 8.43, б со Ы R » 1.

При интегрировании импульсных воздействий длительностью zH дол­жны быть выполнены следующие условия: /?С»гн для схемы на рис. 8.43, а и L/ R » tH для схемы на рис. 8.43, б.

Напряжение с выхода интегрирующего (дифференцирующего) устрой­ства подается для наблюдения (записи) на электронный осциллограф.

§ 8.59 Передаточная функция четырехполюсника на комплексной частоте. Под передаточной функцией четырехполюсника К(р) на ком­плексной частоте р понимают отношение выходного напряжения 4/2(р) ко входному U\(p} (рис. 8.44, а):

         
 

 

К(р) зависит от схемы четырехполюсника, числового значения элемен­тов схемы и от частоты р. Для четырехполюсника (рис. 8.43, б)

К(р) =------------- . Из уравнения (8.76) следует, что

R + р L

U2(p) = K(p)U}(p).

Под передаточной функцией четырехполюсника для синусоидально­го процесса на частоте со понимают отношение

         
 
   


 

X'(yco) получают из K(p) заменой р на у со, |/С(усо)|— модуль, а

ср(со)— аргумент K(j со). Для схемы на рис. 8.43, б K(j со) ------------------------------------------- ------------------------------------------------------------- -

ср(со) = arctg

Зависимости | K(J со) и ф(со) изображены на рис. 8.44, б, в. Если не-

сколько четырехполюсников, например три, соединены каскадно (рис. 8.44, г) и известны передаточные функции каждого четырехполюс­ника, то передаточная функция каскада в соответствии с формулой (8.77) равна произведению передаточных функций этих четырехполюсников:

     


 

 

Пример 103. На рис. 8.45 изображена замкнутая система (система с обратной связью). Она состоит из основного четырехполюсника с передаточной функцией К(р) и четырехполюсника обратной связи с Кос(р) Функцию последнего часто выполняет уси­литель, работающий в режиме пропорционального усиления. Вывести формулу передаточ­ной функции всей системы Кзс(р).

Решение. На вход основного четырехпо­люсника поступает основной сигнал G't(p) и сигнал с выхода четырехполюсника обратной свя­зи, поэтому

U2(p) = (U](p)±Uoc(p))K(p). (8.80)

Кроме того,

Uoc(P>= Koc(p)U2(p). (8 81)

Подставим (8.81) в (8.80). Получим

Л. (Р) —--------- — —-------------------------

и^р) \±К{р)Кос(р)

Если 1 - К(р) КО(:(р) = 0, то в системе возникнут автоколебания, амплитуда их будет ограничиваться нелинейностью системы. Плюс в формуле (8.80) и минус в формуле (8.82) соответствуют положительной обратной связи. Минус в формуле (8.80) и плюс в (8.82) — отрицательной.

§ 8.60 Переходные процессы при воздействии импульсов напря­жения. Ток в любой схеме при действии на нее импульса напряжения (ри- c. 8.46, а) можно найти, например, тремя способами:

1) применяя интеграл Дюамеля;

2) определяя ток при t <t} так же, как от действия постоянного на­пряжения £/; при t < Zj действующее на систему напряжение равно нулю. Следовательно, система освобождается от вынуждающих ЭДС и по ней протекают свободные токи, обусловленные запасом энергии в индуктив­ных и емкостных элементах системы;

3) представляя импульс в виде двух постоянных напряжений. Поло­жительное напряжение U действует начиная с t = 0, отрицательное — начиная с / = При t < токи в цепи определяются одним напряжени­ем U\ при i > Zj — обоими напряжениями с учетом сдвига второго на­пряжения на время

и к

Рис. 8.46

 

Рассмотрим третий способ. Положим, что требуется найти ток в цепи при подключении ее к источнику напряжения, имеющего форму равно­бедренного треугольника (рис. 8.46, б). Задача решается в три приема.

Сначала определяем ток в интервале времени от / = 0 до Z</| от действия напряжения щ - к / (рис. 8.46, в). Затем для интервала времени l2 находим ток в цепи от действия двух напряжений (рис. 8.46, в, г): от продолжающего действовать напряжения = к t и от вступающего в действие при t -tx дополнительного напряжения и2 = —2 к (t —).

Для интервала времени / > t2 ток определяется действием трех напря­жений: продолжающих действовать напряжений щ и и2 и вновь всту­пающего в действие при t =t2 напряжения и2 - к (I -t2) (при t >t2 сум­ма напряжений wb и2 и иу (рис. 8.46, д) даст нуль).

Из трех перечисленных способов наиболее экономным является пер­вый.

При воздействии серий импульсов переходный процесс рассчитыва­ют часто операторным методом.

Пример 104. На последовательно соединенные R и L поступает серия прямоугольных импульсов напряжения единичной амплитуды; длительность импульса х и длительность паузы также т (рис. 8.46, е). Используя третий способ в сочетании с теоремой запаздыва­ния (см. § 8.40), определить ток в цепи.

Решение. Найдем изображение напряжения:

t/(p) = l(l-e-^ +е-2^-е’3рт + ...).

Р

Выражение в скобках представляет собой бесконечную геометрическую прогрессию со знаменателем — е"7”. Сумма членов ее равна ——- -—. Изображение тока

     
 


 

     
 


Применим формулу разложения. Корни знаменателя:

р' = 0; p‘ = -R!L\ х Pt ={ak + J Ь^)х = j п(2 к+ 1) (-a> <£<<»).

Группируя член к=0 с к--\. член А = ] с членом к = -2 и т. д., получим

х $in(n(2fc + 1)--Ф2* + |)

+-Z--------------------------- —

Л *=0 *2[§] [**]1

 


(ад

-Ат / 2


dt равна единице. Единицей измерения 8-функции


является се­


 

кунда в минус первой степени. Если импульс действует при некотором

 


 
а б в Рис. 8.47

времени t = /1} то он обозначается как 50* - б), т. е. импульс действует, когда аргумент 5-функции равен нулю.

Единичной функцией 1(/) (рис. 8.47, б) называют функцию, равную единице при />0 и равную 0 при /<0. Единичной функцией ](-/) (рис. 8.47, в) называют функцию, равную 1 при t <0 и равную 0 при Г>0. Функции 1(г) и 1(-г) имеют нулевую размерность. Рассмотрим свойства перечисленных функций:

— производная по времени функция !(/) равна S-функции, т. е.

®

at

5-функция обладает фильтрующим действием:

f^^t-t.^f^b^t-t^

— изображение 5-функции по Лапласу равно единице: а изображение функции 5(f - ) на основании теоремы смещения рав­но е“р/‘;

 


   
ОО
 
0
 


 
dt = \>

— единичные функции 1(/) и 1(-/) также обладают фильтрующим действием.

Умножение произвольной функции /(/) на ](/) обращает произведе­ние f(t) 1(г) в f(t) при />0 и в нуль при / < 0. Аналогично произве­дение /(/)!(-/) обращается в нуль при />0 ив /(г) при 1<0. Им­пульсное (игольчатое) напряжение или ток в виде 5-функций единичной площади записывают так: 5(/) • 1. Здесь единица имеет размерность вольт- секунда или ампер-секунда соответственно.

Если на электрическую цепь, входная проводимость которой равна g(/), при нулевых начальных условиях воздействует единичный импульс напряжения 5(/)-1 Вс (рис. 8.47, а), то для определения вызванного им тока в цепи представим импульс в виде двух прямоугольных напряже­ний (положительного и отрицательного) с одинаковыми числовыми значе­ниями ] / Дт, сдвинутых во времени на величину Дт -> 0 (рис. 8.47, г).

Ток в цепи после окончания действия импульса

>(О = 4- (g(r + Дт / 2) - g(Z - Дт / 2)).

Дт

Разложим функции £(/ + Дт/2) и £(г-Дт/2) в ряд Тейлора по сте­пеням Дт/2, учтем малость Дт и получим

'■(О=4" («(')+г'(') - 8(0+-V = I ■«’(/).

Дт у 2 2 J

Здесь и далее 1 имеет размерность Вс.

Величину g'(/) = называют импульсной переходной проводимо- d t

стью. При />Дт (Дт-»0) g'(0» умноженная на 1 Вс, численно равна току в цепи при воздействии на нее единичного импульса напряжения

$(/)•! Вс. Аналогично величину /?*(/) =

переходной функцией четырехполюсника. При t > Дт (Дт неженная на 1 Вс, численно равна напряжению на выходе четырехпо­люсника при воздействии на его вход единичного импульса напряжения 5(/) • 1 В с. В интервале времени от / = -Дт / 2 до / = Дт / 2 (во время дей­ствия импульса 5(z)l Вс) напряжение на выходе четырехполюсника

и2(О = Л(0)б(/)-1 + Л'(О-1,

а ток на входе двухполюсника

/(/) = g(0) 6(0 • 1 + g'(0 • 1 •

Наряду с понятиями «переходная проводимость» g(0 и «импульсная переходная проводимость» g'(f) применяют дуальные им понятия: пе­реходное сопротивление г(Г) и импульсное переходное сопротивление Переходное сопротивление rafr(z) численно равно напряжению на входе цепи uQh(t) при воздействии на ее вход единичного тока:

"«л(О = Ц^)^(О.

Импульсное переходное сопротивление rj5(z) численно равно напря­жению на входе цепи uaf)(t), после того как на ее вход воздействовал импульс тока в виде 5-функции единичной площади:

uuh(0 = 5(0- 1(А с)-г'аЬ(1).

Величины r(z) и r'(t) могут быть входными и взаимными, однако g(t) и R(t) не являются взаимно-обратными величинами; g(z) опреде­ляется при питании схемы от источника ЭДС, а Я(/) — при питании схемы от источника тока.

Подчеркнем, что в литературе по переходным процессам в зависимо­сти от рассматриваемого вопроса под одним и тем же названием — им­пульсная переходная функция — понимают функцию либо либо A5(z). Между этими функциями имеется зависимость

Л5(/) = Л(0,)6(0 + Л'0);

Л'(/) характеризует реакцию четырехполюсника (его выходное напряже­ние) после окончания воздействия на его вход единичным импульсом напряжения l-5(z)B-c, а /?5(z) —напряжение на выходе четырехполюс­ника и во время действия импульса, и после окончания.

Аналогичные соотношения существуют между двумя импульсными переходными проводимостями:


s5W=g(ot)S(o+g'(/)


и между двумя импульсными переходными сопротивлениями:


Л5(/) = Я(0ф)8(г) + Ш


при воздействии на вход схемы единичным импульсом тока. С помощью h5(t) интеграл Дюамеля запишется так:


«г(0 = |“(^)Л8(/-т) с/т.
о


Здесь Л 5 (/-?) = Л(0) 5(г - т) + h'(J - т).


Формулу интеграла Дюамеля в математических работах называют формулой свертки двух функций, в данном случае функций г/(/) и As(z).


§ 8.62. Определение h{t) и А6(/) через К{р). Как упоминалось, при воздействии на вход четырехполюсника единичного напряжения «((/) = ](/) напряжение на выходе его u2(t) = h(t). Если это положение записать относительно изображений, учитывая, что !(/) = —- и обозначив р изображение Л(/) через Н(р\ то Н(р') = К(р)/р.

Отсюда


К(р) = р Н{р).


(8.83)


Определим теперь Л(/) через К(р). Поскольку h(t)^H(p\ а Н(р') определено предыдущей строкой, то


(8.84)


При воздействии на вход четырехполюсника единичным импульсом напряжения ut(f)- 1 ♦§(/) = ] =щ(р) напряжение на выходе его


«2(о=/>8(о=<адад=1-ад,

таким образом,


*5(0=ад-


(8.85)


Пример 105. Запишем h(t\ hb(t) для схемы на рис. 8.42, а.


hb(t)=K(p) =


ПС


 

 

.63. Метод пространства состояний. Метод пространства состо­яний (метод переменных состояния) представляет собой упорядоченный способ нахождения состояния системы в функции времени, использу­ющий матричный метод решения системы дифференциальных уравнений первого порядка, записанных в форме Коши (в нормальной форме). При­менительно к электрическим цепям под переменными состояния пони­мают величины, определяющие энергетическое состояние цепи, т. е. токи через индуктивные элементы и напряжения на конденсаторах. Значения этих величин полагаем известными к началу процесса. Переменные со­стояния в обобщенном смысле назовем х. Так как это некоторые функ­ции времени, то их можно обозначить х(Г).

Пусть в системе п переменных состояния. Матрицу-столбец перемен­ных состояния в n-мерном пространстве состояний обозначим [х] = т выходных величин (токи, напряжения) обозначим у, матрицу-столбец

У]

выходных величин [у] =

L-M

Источники воздействий (источники ЭДС и тока) будем именовать z.

Матрица-столбец источников воздействий [z] =

Для электрических цепей можно составить матричные уравнения вида

     
 


 

[у] = И И+[0] И»

где [Л/], (/V], [Р], [£>] — некоторые матрицы, определяемые структу­рой цепи и значениями ее параметров.

На основании принципа наложения решение (8.86)

[x(r)] = eIM,f [х(0)]+ |е(ЛУ]ит) [W][z(t)J dx,
о

где (х(0)]— матрица начальных значений х.

Первое слагаемое в формуле (8.88) описывает свободные процессы в системе, второе — принужденные и свободные при нулевом исходном состоянии (вывод формулы (8.88) см. в конце параграфа).

Из (8.87) и (8.88) находим

МО) = И е'"1' МО)] + j[P] e|M l('-”[W] [z(t)] dt + [Q] [z(/)J. (8.89)

0

Рис. 8.48

Систему уравнений 8.86 составляют либо на основании законов Кирхгофа, либо путем при­ведения схемы к резистивной (без элементов L и С). Как это можно сделать, пояснено в при­мере 100.

         
 


[М] = -у; [W) = |; И = Е;

Поясним формулу (8.88) на простом приме­ре. Ток в схеме на рис. 8.48 до коммутации был /(0_) = Е/(2 Я). Уравнение состояния для этой схемы dildi = -(Я/L) i + (Е/Е), т. е.

 

 

Матричную функцию e^' в формуле (8.88) вычисляют по формуле (теореме) Сильвестра:

=e^^j1] + e^/ [^l + .^ + e’1"' [Л„],

где

п in

[Л]=П(["]-Х>[Ч) /1Ж-Ъ);

7 = 1 / 7-1

f г

Здесь Хг— собственные значения (характеристические числа) квад­ратной матрицы [Л/], т. е. корни уравнения

det([A/]-X[i]) = 0. (8.92)

Из уравнения (8.92) следует, что уравнение относительно X состав­ляют, приравнивая к нулю определитель матрицы [М], в котором все элементы этой матрицы атт = I,..., и), расположенные по главной диагонали, заменяют на элементы атт -X.

Характеристические числа X— это не что иное, как корни характе­ристического уравнения послекоммутационной схемы. Запись решения в виде ряда (8.90) предполагает, что все характеристические числа раз­личны (нет кратных корней).

Если же среди корней уравнения det([A/]-X[l]) = 0 будет кратный корень X, кратности 5, то составляющая е^', обусловленная этим кор­нем, имеет вид

         
 


где A dj (Х[1]-[Л/])— присоединенная матрица к матрице Х[1}-[Л/]. В ней все элементы заменены на алгебраические дополнения, а за­тем проведено транспонирование. Составляющие решения по формуле (8.93) соответствуют части решения по формуле разложения (см. § 8.50), учитывающей кратные корни.

При машинном счете функцию е^', подсчитывают разложением в ряд:

е1"1' =[1] + [М]; + ^ЦД- + ....

3.

Пример 106. Методом пространства состояний исследовать переходный процесс в схеме на рис. 8.49, а. До коммутации был установившийся режим:Е = 4В; J = 1 А; К = 2 Ом; Л = 1Гн; С = 1 Ф.

Решение. Обозначим токи и напряжения в соответствии с рис. 8.49, а. До комму­тации

     


 
Е 2 R
 


     
uc(0_) = R
 
2R)

 

В качестве переменных состояний выбираем ток /( и напряжение на конденсаторе и(-.

Известно несколько способов составления уравнений состояния. Рассмотрим иногда применяемый, основанный на сведении послекоммутационной схемы к резистивной с ис­точниками ЭДС и тока. С этой целью индуктивные элементы в послекоммутационной схеме заменяем на источники тока, которые доставляют ток в том же направлении, что и в ис­ходной схеме (в рассматриваемом примере L заменяем на источник тока ij с напряжени­ем на нем L di} Idt), а конденсатор С — на источник ЭДС, причем в соответствии с тео­ремой компенсации ЭДС этого источника должна быть направлена встречно току в ветви с конденсатором, т. е. встречно напряжению на конденсаторе (в рассматриваемом при­мере конденсатор С с напряжением на нем ис заменен на источник ЭДС Е) =и(-).

В результате схема окажется без индуктивных и емкостных элементов (чисто резис­тивной). но с дополнительными источниками тока и ЭДС (рис. 8.49. б).

 


 
а а

 

а б

Рис. 8.49

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

     


 

По известным потенциалам узлов рассчитаем напряжения на источниках тока Lkdikldt эквивалентируюших индуктивные элементы Lk, и токи im=C„duCmldt че­рез источники ЭДС, эквивалентирующие емкостные элементы емкостью Ст.

Для первой ветви схемы на рис. 8.49, б

Фа = (/! + J) Я + = £ -> /( R~ L ■

dt

Отсюда

di, 2 R . и/- Е R

—L ---------- f _ _s_ +-----

dt L ' L L L

Ток второй ветви i2 можно определить по первому закону Кирхгофа или по закону Ома для участка цели с источником ЭДС:

^а~ис <i}+J)R + uc-uc

R R

Следовательно,

du(-/ dt = (ц/С) + (J/С).

Таким образом, уравнения переменных состояния для послекоммутаиионной схемы на рис. 8.49, а таковы:

 


         
 


     
 


 
[Л/] =
                 
 


Составим уравнение для определения характеристических чисел X:

 


 

Таким образом.

2?+4Х + 1 = 0; Л|=-0,27; Х2 = -3,73с“‘.

По формуле (8.91),

 


 
-0,078 -0,289 0,289 1,077
         
 


 


 
По формуле (8 88),

 

Выполнив подсчеты, получим

/j =-1+0,75 е'027' +0.75 е-3'73' Л:
ие = 6-2.8е“°-27/-0,2е-3,73' В.

Если за выходную величину^ принять напряжение между точками d и/, то

 


 
[«</] = [-*
                           
 
 


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

ч(/)] = [<р(/-т)] [«(/)] (х(т)].

Общее решение

[х(П] = [Ф(/ - т)] [x(t)J + [Ф(/ - г)] [«(/)] (х(/)] =
= [ф(/ - т)] (1] + [»(/)] [х(т)] = [<р(/ - т)] [Л(0],

где Л(г) требуется определить.

Подставим

(х(г)] = [ф(/ - т)1 {/?(/)] (8.96)

в уравнение (8.86):

((<i>0 - t)J - ] [Ф(г - т)]) (Я(Г)] + [Ф(/ - т)] (Я(/)] = [.л/] [z]. (897)

Поскольку (<р(/-т)| есть матрица, столбцы которой являются решением уравнения (8.94), то первый член выражения (8.97) — нулевая матрица. Следовательно.

(ад = [<р(/-т)Г1 bVH-J (898)

Проинтегрируем (8.98) от т до/:

(/?(/)]-[ВД]= /[ф(Х-т)Г? [.VjUJrfk. (8.99)

т

Из уравнений (8.96) и (8.99) следует

[Ф(/ - т)Г' КО] = (Ф(О)Г' (х(т)] + |(<р(Х - т)]'1 [;V][z(X)J d\, (8.100)

X

но [ф(0)} = [1]. Умножая (8.100) слева на [<р(/ - т)] и учитывая, что

[ф(/-т)Нф('-т)Г’ =е[Л/1('-,) e'(W](Vtl1Л/](,

получим

[х(/)) = [ф(/-т)][х(т)1+ |[ф(Х-т)ЦУ)(с(Х)]агХ. (8.101)

X

Полагая в (8.101) т = 0 и заменяя затем переменную X на т, получим формулу (8.88).

§ 8.64 Дополняющие двухполюсники. Два двухполюсника, содер­жащие элементы R, L, С, называют дополняющими, если входное сопро­тивление при их последовательном (параллельном) соединении оказы­вается чисто активным, не зависящим от комплексной частоты р. Так, двухполюсник из параллельно соединенных L и R2 и двухполюсник из параллельно соединенных С и R} (рис. 8.50, а) являются дополняющи­ми при их последовательном соединении и выполнении условия

r]=r2=r = /l i с.

 


 
a б в Рис. 8.50

Двухполюсники Я2, Си Л,, L при их параллельном соединении (рис. 8.50, б) являются дополняющими при том же условии.

Элементы двух дополняющих двухполюсников взаимно-дуальны. Эле­ментам С,, /?! одного соответствуют такие дуальные элементы С2, L2, Я2 дополняющего, что произведение сопротивлений двух взаимно дуальных элементов должно быть равно R2, где R — произвольное активное сопротивление.

Последовательное соединение и С} в исходном двухполюснике заменяют на параллельное соединение С2 = Ц/ R2 и L2 = Q R2 0 допол­няющем.

Параллельное соединение С] и Д, в исходном двухполюснике заме­няют на последовательное соединение L2 = С( R2 и С2 = L^l R2 в допол­няющем.

§ 8.65 Системные функции и понятие о видах чувствительнос­ти. Системные функции Н(р) — это обобщенное название функций, характеризующих четырехполюсник. Ими могут быть, например, пере­даточная функция напряжения U2(p)IU\(p), передаточная функция тока /2(р)/ АС/7) и т. п. Если какой-либо параметр (Я, L, С) в схеме четырех­полюсника изменяется, то изменяются модуль и аргумент системной функции и можно говорить о чувствительности системной функции к изменению этого параметра.

Под классической чувствительностью понимают отношение относи­тельного изменения функции &Н(р)/Н(р} к относительному измене­нию параметра Дх/х:

             


Применительно к установившемуся синусоидальному режиму рассматривают чувствительность модуля и чувствительность аргумента //(/ со.)

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

Понятие чувствительности используют главным образом в задачах синтеза. Электрические цепи стремятся сформировать так, чтобы они
были по возможности малочувствительны к изменению параметра. Если Н(р) зависит от многих параметров и все они могут изменяться, то верх­ней границей возможной ошибки считают сумму модулей чувствитель- ностей по всем параметрам. При определении классической чувствитель­ности можно воспользоваться теоремой вариаций (см. § 2.19) и теоремой Теллегена (см. § 3.43).

§ 8.66 Обобщенные функции и их применение к расчету переход­ных процессов. Обобщенными функциями (ОФ) называют функции времени fit), которые терпят разрыв, например при г = 0. Значение функции при t <0 обозначим А (О, при />0 А (О (рис. 8.50, в).

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

f W = f- (О К-0 + Ш ко.

В общем случае /(/) может содержать также 8-функцию и ее про­изводные, Производная от /(/)

А/ м/ ч r'/м/ч г / ч

= А (О к-0 + А (О КО + А (О А (0 ~~ =

dt d(-t) dt dt

= а'(О к-0 + л’ (О КО + 5(r) (A (0) - A (0)L

Используя ОФ, можно решать задачи на переходные процессы, о ко­торых говорилось в § 8.28, а также задачи на импульсные воздействия. В этом случае необходимо составить уравнения для послекоммутацион- ной схемы, выразить токи, напряжения и их производные через ОФ и, воспользовавшись фильтрующим свойством l(-r), l(z), и 8(/), прирав­нять в этих уравнениях коэффициенты, содержащие только !(-/)> толь­ко !(/) и только 8(/), затем решить их совместно.

Пример 107. Путем использования обобщенных функций решить задачу примера 86 (см. рис. 8.26).

Решение. В уравнении для послекоммутационной схемы подставим

     
 


“л = *л-(О Ц“0 + «л+0) 1(0;

и ('2 = ис2-{{У U-0 + исг*(О КО;

u ci - К_О + т,{I) КО + 5(О (0,) - U(~।(0_)];

«С2 - «С2(О Ц-0 + «гзЛО1(0 + S(O 1«С2(°+) - «С2(°-)]•
£ = £!(-/) + £![/].

Коэффициенты яри 1(~г), 1(0, и 6(/) дают три уравнения:

Л (С| W^)+(0 + Cj «Г2-(0)+ Uci-V) =

Я(С( »cu(0 + C2 Uc2+(0) + «Cb(/)= Е>
u Ci (0* ) (Q + С'2);1 uCl(Q-) + С-2 и('2^-У

Из (8.103) u(?1_ (/) = £, из (8.105) Wc)(0*) = C] £/(C)+C2); далее решаем (8.104) клас­сическим или операторным методом, имея в виду, что (/) = иГ2 </). В результате получаем тот же ответ, что и в примере 86.

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

иДО = I sin о i - lme7<D/.

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

и2 (0 = 1т((Л(0) + j*’(T) е‘,<3 ' dx) е'w/) = 1т(а(ш,г) е'“'). (8.106)

Здесь о

<

о(<о,0 - Л(0)+ J/f(r) е"7"1 dx = /«(©,/) + j л(о.О = . (8.107)

о

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

1/|(г)=1л1((/я1(г)е')

и определим

и2(/) = 1гп((/?(0) (/„,(/) + р’(т)U„,(t) (/ - т)e-/ttdx)е'"').
О

Учтем, что —*-■ -Л'(т)е'1шг ~а'(и).т) и Л(0) = о(е,0). Тогда

dt

«>(/) = 1т(Л(и.г) eywZ). (8.108)

где

>4(aj) = a(0,O)t/w(r)+ |о'(<о,т) Um(J~ *)dx\ (8.109)

о

/l(w,z) — огибающая выходного напряжения. Формулу (8.105) называют интегралом Дю­амеля для огибающей, она позволяет, не вдаваясь в мелкие детали, выявить макрострук­туру переходного процесса.

Пример 108. Определим огибающую тока в цепи, когда на вход последовательно соединенных R и L воздействует напряжение и,(/) = к t sin <о t.

1

Вместо h(t) используем g(t) = — (l-e ). R

В соответствии с формулой (8.107)

о(<й,г) = 5(0)+ fg'(x) е-?от dx~-—!---------------------------------------------------------------- (1-е-*');---------------------------------------------------------------- = —+ >ш.

J R + j L L

Учтем, что g(0) = 0, = -J-e ч\ Um (t-т) = к (/ - т).

Огибающая тока в пели по формуле (8.109):

.. . к t. . -q х j к? к Lt —at к I*

= — (г - т) с * dx - ----------- --------- - +-------------------- (е ? -1) = —-------------------------------

L ' R + jaL (/? + ;<о£)2 Л2+(соО

I, ~ *' ч2 z К>

Il R t **------------- I J----------------------- —

У — + е 1- cos(<o t + 2 <p) - cos 2 <p + |wz + e L sin(co t + 2 (p) - sin 2 ф

x A J \

Rf

tof + e L sin(w / + 2 ф) - sin 2 q>

Rt

R t -b~~

— + e f- cos(co t + 2 ф) - cos 2 ф

Вопросы для самопроверки

1. Дайте определение переходном}' процессу. 2. Что понимают под принужденными и свободными токами и напряжениями? 3. Сформулируйте законы (правила) коммутации. 4. Дайте определение зависимым и независимым начальным условиям. 5. Какие вы знае­те способы составления характеристического уравнения? 6. Объясните, почему при состав­лении характеристического уравнения путем приравнивания к нулю входного сопротив­ления Z(p) = ЛЧр)/Л/(Г) в общем случае нельзя сокращать числитель и знаменатель дроби на общий множитель. 7. Чем определяется число корней характеристического урав­нения? 8. Изложите сущность классического метода расчета и принцип составления урав­нений для определения постоянных интегрирования. 9. Переходный процесс в некоторой цепи сопровождается биениями. О чем это может свидетельствовать? 10. Дайте обосно­вание обобщенным законам коммутации. 11. Запишите известные вам соотношения меж­ду /(О и Г(р), а также теоремы операторного метода и предельные соотношения. 12. По­чему/? называют комплексной частотой? 13. Охарактеризуйте этапы расчета операторным методом. 14. В чем особенности расчета переходных процессов операторным методом при синусоидальном источнике и ненулевых начальных условиях? 15. Охарактеризуйте свой­ства единичной функции !(/) и свойства дельта-функции 5(f). 16. Определите переход­ную и импульсную переходную проводимости (сопротивления). Укажите, с какой целью они используются. 17. Охарактеризуйте идею расчета с помощью интеграла Дюамеля. 18. Прокомментируйте известные вам формы записи интеграла Дюамеля. 19. Какими спо-

     
 


     
в г д
 
Рис. 8.51


 


     
 
1\ Ri Е
Рис. 8.52 Рис. 8.53

собами можно определить отзвук системы, когда на нее воздействует импульс напряже­ния или тока? 20. Поясните принцип работы интегрирующих и дифференцирующих це­пей. Запишите условия, при которых эти цепи выполняют свои функции. 21. Чем следует руководствоваться при формировании дополняющих двухполюсников? 22. Поясните идею расчета переходных процессов с помощью обобщенных функций. 23. Перечислите основ­ные этапы расчета методом переменных состояния. 24. Как составляют уравнения пере­менных состояния путем сведения послекоммутационной схемы к чисто резистивной?

25. Охарактеризуйте сильные и относительно слабые стороны известных вам методов рас­чета переходных процессов. 26. Что понимают под системными функциями? Какие виды чувствительности системных функций вы знаете? 27. В схеме на рис. 8.51, а с источни­ком тока Jq в момент / = 0 одновременно размыкается ключ Кг и замыкается По­казать, что заряды, протекающие через сопротивление и Я2 за время от 0 до не

зависят от емкостей Q и С2. Определить величины этих зарядов. (Ответ:

.) 28. В схеме рис. 8,6. а при размыкании ключа происходит переходный

процесс. Определить законы изменения во времени напряжений и( ) и «сг на конденса­торах. Задано j(t) = I sin(wt+ 90°) А, /? = 1/юС = 1Ом. со = 100рад/с. (Ответ: uCi = 0,447 sin(<o I + 63°27') -0,253 - 0.15 е-'200' В; «Г2 = О-447 $in«»' + Ы°2?') + 0,253 - -О,15е‘у200' В.) 29. Покажите, что в симметричной мостовой схеме (рис. 8.51, в), в

I

которой выполняется условие £,/<’=/?2, переходная функция Л(/)»-у+е 4 . 30. В схе­ме рис. 8.51, б R = L - С = 1. Покажите, что входная переходная проводимость равна t е"'.

31. Покажите, что энергия, запасаемая в L схемы рис. 8.51, е (начальные условия нулевые), равна тепловым потерям в R. 32. Первичная обмотка трансформатора рис. 8.51,0 при ну­левых начальных условиях подключается к источнику постоянной ЭДС Е, R} - R2 = R: L}=L2=M. Определите q(0^), /2(0„). (Ответ: /((0+) = -/2(0+) = £/ (2 R).) 33. Опре­делите степень характеристического уравнения для схемы рис. 8.52. (Ответ: пятая.)

34. Как определить К(р) через h(t) и через Л&(/)?35. По h(r) = — (1 + 2е [6]- ) четы­

                         
 
 


/? = 0,2Ом. С = 5Ф, £ = 1Гн. (Ответ: Л(/) = 1,62е-0,724'-0.62 е’0,276/.) 37. На вход че-

J fc)

тырехполюсника с K(j со) =------------------ воздействует единичный импульс напряжения в виде

5-функции. Определите напряжение на выходе четырехполюсника после окончания дей­ствия импульса. (Ответ: 0,25 е"^5'.) 38. В схеме на рис. 8.53, а до коммутации был уста­новившийся режим; £ = 20В, Jim р/](р) = 2А и lim р /((р) = 5 А. Определите Л, />->0

и R2. (Ответ: Я, = 4 Ом, Я2 = 6 Ом.) 39. В схеме на рис. 8.53, б до коммутации был ус­тановившийся режим; £ = 10 В, limp/2(p) = 2A и Нт р 12(р) = 1,428 А. Определите />-»» р~»О

Я, и Я3, еслиЯ3 = Я2. (Ответ: Я, =4 Ом, Я2 = 3 Ом.) 40. Решите задачи 11.4; 11.12; 11.15; 11.26; 11.29; 11.32; 11.38; 11.40; 11.47; 11.50; 11.55; 11.57.


Глава девятая

ИНТЕГРАЛ ФУРЬЕ.

СПЕКТРАЛЬНЫЙ МЕТОД. СИГНАЛЫ


§ 9.1. Ряд Фурье в комплексной форме записи. Как известно из предыдущего (см. § 7.2), в ряд Фурье можно разложить любую периоди­ческую функцию /(/), удовлетворяющую условиям Дирихле.

Обозначим период функции Г, а основную частоту — соо = 2 л !Т. Ряд Фурье можно записать двояко.

Первая форма записи:


00

Л0=4 + Х4 sin<*


вторая форма записи:


/(/) - Ло + (Л£ sin к o)0 t + Ак cos к со0 /),

*=|


(9-D

(9-2)


где — постоянная составляющая ряда; Ак амплитуда ^-гармони­ки ряда; — начальная фаза А-гармоники;


А'к - Ак cos ук; Ак - Ак sin ;


i 772

4) = 7 f ЛО

1 -тн


Т!2

]/(/) sin к w0 t dt\

-TH


(9-3)

(9-4)


772

]"/(/) COS А (00 1 <*1-

-TH


(9-5)


~&~JX

Из курса математики известно, что sin х =------------------- . Следовательно,

2 У


e J(k ®0+ч/4) _e-(*®o+Vi) sin(A соо t + у J -------------------------------------- —---------------


(9.6)


Подставив правую часть формулы (9.6) в выражение (9.1), получим


(9.7)


 

 

 

Обозначим

А.к =-лк e,4t.

Тогда ряд (9.7) можно записать так:

(9-10)

Формула (9.10) представляет собой комплексную форму записи ряда Фурье. Текущий индекс к может принимать все целые числовые значе­ния от -со до +<», но не может равняться нулю, так как постоянная составляющая ряда выделена в виде отдельного слагаемого.

Пример 109. Представить функцию f(t) = 2 + 3 sin(o0 t + 30°) + 2 sin(2 co0 /- 45°) в комплексной форме записи.

Решение. Ло=2; Д =3eyJ0’; 4, = -Зе*'34 Л2 = 2 е’ 42=-2еу45‘;

у (у) _ 2 4. 1 (3 1 + _ 3 g”? (иу “ +30*) 2 f-45*) _ 2 /-45”)

2 У

Составим выражение для комплексной амплитуды Ак. По определе­нию (см. формулу (9.8)),

4 = 4 e>Vx = Ак C0SVa + J 4 sin = 4 + j Л”к, (9.11)

где 4 определяется формулой (9.4), 4 — формулой (9.5).

Подставим правые части формул (9.4) и (9.5) в формулу (9.11):

2 772

4=~ J7(0 (sin к ©о t + j cos к ©0 /) d{ = ' -77 2

772

Jf (t) (cos к o)0 t - J sin к ©0 t) dt,

-7/2

или

Ak dl. (9.12)

' -772

Подставим правую часть формулы (9.12) в формулу (9.10):

7/2

7(') = Л+ Z ^*“"'1 dt. (9.13)

к » -ОО / -7’/2

§ 9.2. Спектр функции и интеграл Фурье. Ряд Фурье — это триго­нометрический ряд, представляющий собой изображение периодической функции суммой синусоид, амплитуды которых конечны, а аргументы кратны основной частоте ©0.

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

Формулу для интеграла Фурье получают из формулы для ряда Фурье (из формулы (9.13)) предельным переходом при стремлении периода Т к бесконечности.

На функцию /(/) при представлении ее интегралом Фурье наклады- -нх>

вают ограничение, а именно, полагают, что $f(f)dt есть величина

конечная (не бесконечно большая). Это серьезное ограничение. Ряд фун­кций этому условию не удовлетворяет'’. । т/2

Так как по определению (см. формулу (9.3)), Л = ~ f/(0 а при +со -7/2

Т —> оо [/(f) dt есть величина конечная, то Ао = 0.

] 7/2

Преобразуем выражение — j/(f) е" ' с/f, стоящее под знаком -772

суммы в формуле (9.13). С этой целью произведение ко заменим на со (под со будем понимать изменяющуюся (текущую) частоту). В ряде Фурье разность двух смежных частот Дсо=©о =2 л/Т. Следовательно, МТ = До)/(2 л).

При Т -> оо, заменив Део дифференциалом с/со, получим

1 '(/О)Л = Л.

1 -Т/2 71 -77 2

Обозначим

5(/о) =

-со

Формула (9.14) дает возможность преобразовать функцию времени /(f) в функцию частоты S(J со); преобразование (9.14) называют пря­мым преобразованием Фурье, a S(j о) — спектром функции f(t). Это комплексная величина, зависящая от вида функции /(f). В соответствии

772 j

f/(f) dt на — S(y to) и учтем, что

•Лл 2 л

*’ Средн функций /(/}, для которых интеграл ^f(t)di расходится, наиболее важной для практики является функция /(/) = А, где А — постоянное число. Для того чтобы эту функцию представить интегралом Фурье, пользуются следующим приемом. Находят интеграл Фурье для функции /(/) = Ле’₽/. где р>0 и /(/) = 0 при t<Q. Для этой функ­ции J/(z) dt сходится, поэтому она может быть представлена интегралом Фурье.

Далее в полученном выражении устремляют р к нулю.

при изменении к от -оо до +<ю со = к со0 также изменяется от -оо до -ню. Следовательно,

i (0- 4-00

ЛО = у- X S(jo)е‘J‘' А».

(О ш -ОС

Заменив сумму интегралом, найдем

1 -rw

Л0 = — Js(> а>) rfw.

Формула (9.15) представляет собой запись интеграла Фурье (форму­лу обратного преобразования Фурье). Она выражает непериодическую функцию f(t) в виде бесконечно большого числа синусоидальных колебаний с бесконечно близкими частотами и бесконечно малыми амп­литудами S(j(&)d(j3 (S(j со) конечно, но произведение S(j d® бес­конечно мало, так как бесконечно мало значение с/со).

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

Преобразования (9.14) и (9.15) являются взаимно обратными.

Отметим, что представление функции f(t) в комплексной форме в виде интеграла Фурье (формулы (9.15)) привело к необходимости фор­мально ввести отрицательную угловую частоту. При этом сумма слага­емых подынтегральной функции (9.15) при ±со дает синусоидальные колебания частоты со.

Сопоставим формулу (9.14) с формулой преобразования по Лапласу:

со

F(p) = f/(r) е'" Л, (9.16)

о

когда /(г) = 0 при t <0.

Если учесть, что f(t) = O при /<0, и заменить р на у со, то (9.16) переходит в (9.14). Следовательно, формулы для спектра функции SV о) могут быть получены из соответствующих выражений для изображений по Лапласу, если в последнихр заменить на у со.

Пользуясь соотношениями § 8.39, найдем спектр функции f(t) = е~а/, полагая, что /(/) = 0 при t <0.

Изображение ее по Лапласу 1 / (а + р). Заменим р на j со и получим спектр S(j со) = I /(а + j со); S(j со) есть комплексная величина, равная S(co)ey<‘. Модуль ее равен 1/Vet2 ч- со2 , аргумент фл. = arctg(-co/a). Графики для экспоненциального импульса изображены на рис. 9.1, а, б.

Пример 110. Найти S(o) и ф<со) для прямоугольного импульса (рис. 9.1, в) ампли­тудой Л и длительностью

 


Рис. 9 1
                           
 
   


S(co) = —sin = A /
co /„ 2

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

cos со -1 со /..

tg Ф., = —:—-— = - tg-r21- sin со /„ 2

График Ф, дан на рис. 9.1, д. При значениях со/,, = л. Зя.... фл возрастает скачком на л.

Обратим внимание на то, что при определении S(J со) путем замены р на j со в формуле для F'(p') следует соблюдать некоторую осторож­ность, если функция /(/) имеет импульсный характер, иначе можно потерять импульсную компоненту в S(J со) в виде дельта-функции. На­пример, изображение функции 1(0 по Лапласу равно \/ р, тогда как

I

спектр £(/со) функции 1(f) равен не l/yco, а л5(со) +------------------------------------------------------------- . Чтобы

показать это, определим спектр функции 1(/) е”р' (Р > 0), а затем устре­мим р->0:

е-и/ е-у«>/ dt 1-------------- р у m

P + J<0 Р~ +(0“ Р2+СО2

Первое слагаемое правой части при 0 -> 0 и при со -> 0 стремится к бесконечности, т. е. имеет вид дельта-функции а 8(со), второе слагаемое правой части при (3—>0 равно 1//<о. Чтобы вычислить коэффициент а, проинтегрируем р/(02 + со2) = а 5(со) по со от -оо до +а>:

 


 
d(a = a JS(co) da). — СО
                                     
 
 
   


 

л5(со) выпадают.

§ 9.3 Спектр функции, смешенной во времени. Спектр суммы

функций времени. Если функции времени /(/) соответствует спектр

9(7 со), то функции соответствует спектр е”уыт S(j со), что сле­дует из теоремы смещения в области оригиналов (см. § 8.40), если заме­нить р на j со.

Так как модуль функции е"ушт равен единице, то модуль спектра функции /0-т) равен модулю спектра функции /(/), т. е. равен 5(<о), однако аргумент спектра функции f{t - т) отличается от аргумента спек­тра функции /(г) на -сот.

Если /(/) представляет собой сумму нескольких функций времени, например /(/) = /(/) + /2(0» а каждая из них имеет спектр соответствен­но (у со) и S2(j то спектр S(J со) функции /(г) равен сумме спек­тров этих функций, т. е. S(J со) - S}(j со) + S2(y со). Это следует из линей­ности преобразования (9.14). Однако модуль S(co) * ^(со) + 52((о) и аргумент срА.(со) * срд| (со) + <рд2(со).

§ 9.4 Теорема Рейли. Теорему Рейли (Релея) записывают следующим образом:

(9.17)

Функция /(/) = 0 при г<0; S(co) представляет собой модуль спект­ра 5(7 со) функции /(/):

5(7 а) = J/We-'"' dt.

-СО

Если принять, что f(f) есть напряжение, приложенное к активному сопротивлению 1 Ом, то левая часть в (9.17) представляет собой энер­гию, выделяющуюся в этом сопротивлении.

Таким образом, площадь квадрата модуля спектра 5(со), разделенная на я, является энергией, рассеиваемой в активном сопротивлении, на которое воздействует /(/).

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

4

7(0 = — р(/о>) е'“'Ло.

2 It *

-со

Умножим обе части последнего равенства на /(/) и проинтегрируем по t от -со до 4-ос:

 


     
+оО
dm dt.
             
 


 


 
dt = S(-j со),

-со

следовательно,

+со | +О0

\S(J со) S(-j w) dсо - — [S2 (со) dm.

J 2 л <

-00 -00

Для перехода к формуле (9.17) учтем, что при / < 0 функция /(/) = 0. Это дает возможность заменить в левой части нижний предел с -оо на 0. Приняв во внимание, что квадрат модуля 3’2(со) есть четная функция +00 +00

частоты, заменим J в правой части последнего уравнения на 2 В

-сс О

результате получим формулу (9.17). Величину 52(<о) называют спектральной плотностью энергии сиг­нала, а функцию S2(co):= /(со)—энергетическим спектрам.

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

Положим, что на вход некоторого четырехполюсника с передаточной функцией К(у со) - К"(со) е7Ф<<0) при нулевых начальных условиях воздей­ствует сигнал /j(/), имеющий спектр 5вх(у со). На выходе четырехпо­люсника появится сигнал /2(/), спектР которого

^ВЫХ(jg)) = /C(;co)Sbx (у св),

4-00

где S,x0 <o)= (ЛОТd’-

-СО

Так как сигнал /2(/) может отличаться от сигнала по значению (по амплитуде), положим в а раз, и запаздывать на некоторое время /0, но по форме должен быть таким же, как и /(f), то можно записать, что /2<О = я/1('-'о)-

Если к функции /2(/) применить преобразование Фурье, то окажет­ся, что спектр функции /2(г) равен

aS„xO<o)e^“'“.

Рис. 9.2

Сравнивая (9.19) и (9.20), замечаем, что

K(j со) = К(со) е7 Ф(С1>) = а е"7°'('.

Следовательно, для прохождения импуль­са или модулированного колебания через четырехполюсник без искажения формы необходимо, чтобы модуль передаточной функции был постоянен (не зависел от час­тоты), а аргумент <р(<в) = -со/о линейно из­менялся в функции частоты (рис. 9.2, а).

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

/С(со) к минимальному равно J1 (рис. 9.2, б). Такой характеристикой обладает, например, схема рис. 3.42, а. Для этой полосы приближенно полагают, что = const; <р(со) = -со /0.

Для того чтобы сигнал при прохождении через четырехполюсник не изменил своей формы, необходимо, чтобы важнейшие гармонические составляющие частотного спектра сигнала находились внутри полосы пропускания четырехполюсника. Для импульсных сигналов треугольной, трапецеидальной, прямоугольной, колоколообразной и некоторых других форм принимают, что они занимают полосу частот от со = 0 до со = 2 л//и, где — длительность импульса.

Если же необходимо передать через четырехполюсник основную часть энергии сигнала (например, 90 % энергии сигнала), то полосу частот можно сузить примерно до 0-г1//и.

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

Первый способ состоит в непосредственном применении прямого и обратного преобразований Фурье.

Основные этапы этого способа таковы:

1) нахождение спектра Ц(усо) входного сигнала w,(/);

2) определение передаточной функции четырехполюсника K(j со);

3) получение спектра выходного сигнала U2(j со) = K(j со) Ц(/ со);

4) вычисление w2(r) по (ЛД/со).

Последнюю операцию можно осуществить с помощью форму­лы (9.15), но практически ее удобнее выполнить, используя таблицу изоб­ражения по Лапласу, заменив на р в U2(ja)).

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

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

Для того чтобы выяснить влияние отдельных элементов схемы на искажение формы импульса, прежде всего составляют полную схему замещения четырехполюсника, учитывая в ней все факторы, влияющие на частотные свойства (паразитные емкости ламп, импульсных трансфор­маторов, индуктивности рассеяния трансформаторов, емкостные свойства р—«-переходов транзисторов, зависимость коэффициентов усиления транзисторов от скорости процесса (от частоты <л)).

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

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

В каждой из этих расчетных схем с учетом упрощений, рассмотрен­ных в § 8.16, число оставшихся индуктивных элементов и конденсаторов оказывается значительно меньше, чем в полной схеме замещения.

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

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

Проиллюстрируем сказанное примером. На рис. 9.3, а изображена схема лампового усилителя, где RH нагрузочное сопротивление; Ср — относительно большая разделительная емкость (через нее прохо­дит только переменная составляющая выходной величины); С2 —отно­сительно малая емкость нагрузки и (или) емкость второго каскада уси­ления. Штриховой линией показаны источник анодного напряжения £н и малые по сравнению с Ср (по нескольку пикофарад) межэлектродные емкости анод—сетка Са^с, сетка—катод Сс_к и С| (емкость анод—-ка­тод и емкость монтажа). В дальнейшем емкости Сс_а и Сс_к не учиты­ваем, как оказывающие малое влияние на работу схемы.

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

           
 


 
Рис. 9.3
         


РЛС £,] ---------- +------- +------- .

*/ *к

• <1 п

Для схемы рис. 9.3, г

8,2 =

р С, ию (р)
рС
<?э2

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