close

Вход

Забыли?

вход по аккаунту

?

Resnichenko Analiz elektr

код для вставкиСкачать
Санкт-Петербургский государственный
архитектурно-строительный университет
Факультет инженерной экологии и городского хозяйства
Кафедра электроэнергетики и электротехники
АНАЛИЗ ЭЛЕКТРИЧЕСКИХ СХЕМ
Методические указания
Санкт-Петербург
2017
Министерство образования и науки
Российской Федерации
Санкт-Петербургский государственный
архитектурно-строительный университет
Факультет инженерной экологии и городского хозяйства
Кафедра электроэнергетики и электротехники
АНАЛИЗ ЭЛЕКТРИЧЕСКИХ СХЕМ
Методические указания
Санкт-Петербург
2017
0
1
УДК 658.382.3:621.3
ЦЕЛЬ И ЗАДАЧИ КУРСОВОЙ РАБОТЫ
Рецензент канд. техн. наук, доцент В. П. Демидов (СПбГАСУ)
Анализ электрических схем: метод. указания / сост. В. В. Резниченко, В. Я. Соколов; СПбГАСУ. – СПб., 2017. – 28 с.
Содержат краткое изложение вопросов теории и пояснение к оформлению работы.
Предназначены для выполнения курсовой работы по теоретическим
основам электротехники для студентов специальности 13.03.02 «Электроэнергетика и электротехника», целью которой является закрепление знаний
об основных методах расчета линейных цепей во временной и частотной
областях.
Табл. 1. Ил. 8. Библиогр.: 8 назв.
Целью данной курсовой работы являются закрепление
и углубление знаний, полученных на лекциях, практических и лабораторных занятиях по курсу «Общая электротехника», а также
приобретение необходимых навыков и опыта при самостоятельном освоении различных методов количественного и качественного анализа процессов передачи сигналов через линейные электрические цепи.
Расчет начинается с нахождения входной и передаточной
функций заданной цепи. Затем определяются переходная и импульсная характеристики для выходного сигнала. Реакция на выходе цепи находится:
а) при воздействии сигнала в виде одиночного импульса;
б) при периодическом несинусоидальном режиме.
Для использования вычислительных средств в курсовой работе введен раздел о формировании описания заданной электрической цепи по методу пространства состояний.
В целом с учетом содержания курса «Общая электротехника
и электроника» и специфики последующих дисциплин, изучаемых
студентами специальности 13.03.02, материал курсовой работы
охватывает следующие разделы теории цепей: анализ линейных
электрических цепей во временной области; переходные и импульсные характеристики схемы; анализ цепей при установившемся гармоническом режиме методом комплексных амплитуд; теория
четырехполюсников; методы расчета сложных цепей; использование операторного и частотных методов анализа.
© Санкт-Петербургский государственный
архитектурно-строительный университет, 2017
2
3
ЗАДАНИЕ К КУРСОВОЙ РАБОТЕ
Z4
Таблица
Вариант
Нечетный
Z2
Рис. 3
Четный
В курсовой работе исследуется цепь второго порядка
с наиболее распространенной лестничной структурой. Для заданной цепи находятся входная и передаточная функции, затем определяются переходная и импульсная характеристики для выходного
сигнала. Формируются уравнения цепи через переменные состояния. Определяются спектр входного воздействия, его ширина и частотная характеристика цепи. После этого определяется спектр реакции цепи на единичное воздействие. Реакция цепи находится:
а) при воздействии сигнала в виде одиночного импульса;
б) при периодическом несинусоидальном режиме.
Задание. Резистивная нагрузка R = 1 Ом соединена с источником при помощи пассивного линейного четырехполюсника,
представляющего собой цепь типа RL, RC или RLC. Возможные
схемы цепей показаны на рис. 1 и 2. В таблице указаны параметры
схем, форма и длительность одиночного импульса для одного из
вариантов расчета. Формы импульсов представлены на рис. 3.
Характеристики
сигналов
Характеристики элементов
1
2
3
4
Форма
τ, с
(рис. 3)
АмT
плиn=
туда
τ
А/В
Схема (рис. 1)
Рис. 1
Рис. 2
4
1
3
5
7
9
11
13
15
17
19
21
23
25
L=4
R=1
L=2
R=2
R=2
C=1
R=1
C=3
R=4
R = 0.5
L = 0.5
R = 0.5
L = 0.5
R=1
L=4
R=2
L=1
R=1
R=1
C=4
R=1
C = 0.5
R = 0.5
R=1
L=2
R=2
L=2
R=2
R=1
L=4
L=1
C=3
R=2
R=2
C=1
C=2
C=2
R=2
R=2
R=2
L=1
L=1
R= 2
L=2
R=2
C=1
C=1
R=1
C = 0.5
R=1
R=2
C=1
5
а
е
в
г
д
а
б
в
е
д
а
б
в
0.1
5
2
2
3
20
1
5
2
2
1
1
0.5
5
1
10
5
1
10
5
1
10
10
1
10
10
2
2
2
2
2
2
2
2
2
2
2
2
2
Окончание таблицы
Четный
Нечетный
Вариант
2
4
6
8
10
12
14
16
18
20
22
24
26
Характеристики
сигналов
Характеристики элементов
1
C=4
R=1
C=2
R = 0.5
R = 0.5
L=I
R=I
L=3
R = 0.25
R = 0.25
C=2
R=1
C=2
2
R=1
L=4
R = 0.5
C=1
R=1
R=1
L=4
R=1
L=2
R=1
R=1
C=4
R=1
3
4
C=1
R = 0.5
R=1
C=4
C=2
L=2
R = 0.5
R = 0.5
L=1
L=1
L=2
R = 0.2
R = 0.2
R = 0.5
C=1
C=1
R = 0.5
C=1
R = 0.5
L=1
L=1
R=1
L=2
R = 0.25
L=2
L=4
Форма
τ, с
(рис. 3)
д
а
г
в
а
е
в
а
е
в
д
в
г
0.1
10
3
2
2
30
0.5
10
2
1
2
1
0.5
АмT
плиn=
туда
τ
А/В
20
10
2
20
10
2
20
10
2
20
20
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
Единицы измерения величин: L – Гн; С – Ф; R – Oм.
1. ОПРЕДЕЛЕНИЕ ВХОДНЫХ И ПЕРЕДАТОЧНЫХ
ФУНКЦИЙ ЦЕПИ, ИХ НУЛЕЙ И ПОЛЮСОВ
Рассматриваемые функции входного сопротивления или
входной проводимости, а также передач по напряжению и току
могут быть определены любыми изученными методами анализа
цепей [5, 6, 8]. Для цепей лестничной структуры наиболее простой
результат дает метод пропорциональных величин [8]. При этом за
исходные величины принимаются напряжение Uн(S) и ток Iн(S)
в нагрузочном сопротивлении Rн и перемещающиеся справа налево по цепи.
6
Операторные сопротивление и проводимость цепи будут
1
U (S )
,
Z11 ( S ) = 11 ;
G11 ( S ) =
I11 ( S )
Z11 ( S )
соответствующие передаточные функции по напряжению и току:
U (S )
I (S )
.
HU ( S ) = н
H I (S ) = н
;
U1 ( S )
I1 ( S )
Для выполнения расчетов вначале получают выражения
в общем виде, а затем осуществляют численные подстановки заданных параметров. Выражения общего вида понадобятся в последующих расчетах. При использовании метода контурных токов
(МКТ) для определения передаточных функций по току и напряжению можно руководствоваться следующими выражениями:
U ( S ) Z н ( S ) ⋅ I н ( S ) Rн ⋅ Δ Z 12 ( S )
.
HU ( S ) = н
=
=
U1 ( S )
U1 ( S )
Δ Z (S )
Здесь принято, что единственный источник напряжения
U1(S), действующий на входе цепи, вызывает в первом, входном,
контуре ток I1(S), а в выходном – ток Iн(S); ΔZ(S) – определитель
системы по МКТ; ΔZ12(S) – eгo алгебраическое дополнение.
Для передаточных функций по току и для входных сопротивлений получим:
U ( S ) Δ Z 11 ( S )
I ( S ) Δ Z 12
.
=
=
;
H1 (S ) = н
Z1 (S ) = 1
I1 ( S ) Δ Z 11
I1 (S )
Δ Z (S )
Расчеты по методу узловых напряжений (МУН) приведут
к следующим соотношениям:
Δ (S )
U ( S ) Δ G12 ( S )
I (S )
HU ( S ) = н
=
= G12
, H I (S ) = н
,
U1 ( S ) Δ G11 ( S )
I1 ( S ) Δ G ( S ) ⋅ R2
Δ (S )
G1 ( S ) = G12
.
Δ G (S )
Здесь ΔG(S) – определитель системы по МУН; ΔG12 – алгебраическое дополнение этой системы. В этом случае предполагается,
что на входе цепи действует источник тока.
При выполнении расчетов искомые функции цепи должны
быть проверены с помощью метода, не используемого при первоначальном расчете. Например, если расчет произведен по методу
7
контурных токов, то проверку можно произвести по методу узловых напряжений.
Корни полиномов числителя передаточных функций цепи
называются нулями, а корни знаменателя – полюсами. Все полюса
и нули для входных функций цепи располагаются в левой полуплоскости комплексной плоскости. В случае входных функций
они должны чередоваться. Собственные частоты наносятся в виде
крестиков (полюсы) и кружочков (нули) комплексной плоскости
с обозначением масштабов.
2. ОПРЕДЕЛЕНИЕ ПАРАМЕТРОВ ЧЕТЫРЕХПОЛЮСНИКА
И ИХ СВЯЗЬ С ФУНКЦИЯМИ ЦЕПИ
Рассмотрим четырехполюсную
цепь, представленную на рис. 4, подключенную к нагрузочному сопротивлению Rн.
Описание четырехполюсника
(ЧП) с помощью параметров проводимостей [G] и сопротивлений
Рис. 4
[Z] может быть получено различными методами [8]. Так, например,
используя МУН для четырехполюсника с двумя внешними источниками тока Iн(S), I1(S), получим матрицу G-параметров следующего
вида:
G
G =  11
G 21
G12 
.
G 22 
Из (1) определяется матрица параметров передачи [A].
G
1 

− 22

a
a

G 21
G 21 
12 
.
A =  11
=

a 21 a 22  G − G 22 ⋅ G11 G11 
 21
G 21
G 21 
8
(1)
Полезно проверить правильность расчетов матрицы A путем каскадного соединения простейших четырехполюсников с известными А-параметрами.
Матрица А-параметров в свою очередь позволяет определить
найденные ранее входные и передаточные функции цепи.
a ⋅R +a
−1
G11 ( S ) = 21 11 22 = Z11
(S ) ,
a11 ⋅ Rн + a12
Rн
HU ( S ) =
,
a11 ⋅ Rн + a22
1
H I (S ) =
.
a21 ⋅ Rн + a22
3. ПЕРЕХОДНАЯ И ИМПУЛЬСНАЯ ХАРАКТЕРИСТИКИ
ЦЕПИ
Переходная характеристика h1(t) представляет собой реакцию
цепи при нулевых начальных условиях на воздействие сигнала
в виде единичной ступенчатой функции δ1(t). Так, например, в цепи с источником напряжения на входе переходная характеристика
для выходного напряжения имеет изображение HU(S)/S, откуда
 H (S ) 
h1 (t ) = L−1  U
,
 S 
(2)
где L–1 – символ обратного преобразования Лапласа.
Импульсная характеристика h(t) представляет собой реакцию
цепи при нулевых начальных условиях на сигнал в виде единичной импульсной функции δ(t), изображение которой L{δ(t)} = 1.
Таким образом, импульсная характеристика для выходного напряжения u2(t) в цепи с источником напряжения будет иметь вид
h(t ) = L−1{HU ( S )}.
(3)
В цепи с источником тока в выражения (2), (3) следует подставить HI(S) вместо HU(S).
Для проверки полученных соотношений можно использовать
выражение обобщенной производной, связывающее h(t) и h1(t),
9
h(t ) =
d
⋅ [h1 (t )] .
dt
(4)
В выражении (4) степени полинома числителя и знаменателя
могут оказаться одинаковыми. В этом случае h(t) может содержать
в качестве слагаемого импульсную функцию.
diL u L (t )
=
.
L
dt
Рассмотрим пример формирования уравнений по методу пространства состояний. Обратимся к цепи, показанной на рис. 5.
R1
4. СИСТЕМА УРАВНЕНИЙ ПО МЕТОДУ
ПРОСТРАНСТВА СОСТОЯНИЙ
dx
= F ( x, u ),
dt
(5)
y = Cx + Du ,
(6)
где x(t) – вектор переменных состояний; y(t) – вектор реакций цепи; u(t) – вектор входных воздействий, называют системой уравнений по методу пространства состояний.
Уравнение (5) называют уравнением состояния, а уравнение (6) – выходным уравнением.
В линейном случае уравнение (5) можно записать в виде
(7)
В качестве переменных состояний принимают непрерывные
функции (напряжение на емкостях или ток на индуктивностях).
Одним из простейших методов формирования уравнения (7) является следующий прием [1]: на основании теоремы замещения индуктивности заменяют на источники тока iL(t), а емкости – на источники напряжения uС(t). В результате получаем резистивную
цепь с источниками и внешними воздействиями. Затем проводится
анализ цепи, при этом определяются напряжения на индуктивностях uL(t) и токи на емкостях iС(t). Производя перегруппировку
членов, находят из уравнения (7) следующие ключевые выражения
duС iС (t )
=
,
dt
C
10
IL
UL
Описание цепи в виде системы дифференциальных уравнений
dx(t )
= Ax(t ) + Bu (t ) .
dt
(9)
(8)
C1
I C1
C2
IC 2
R2
Рис. 5
Заменим ветви с реактивными элементами соответствующими источниками напряжения и тока (рис. 6).
R1
I0
I L UL
UC1
I C1
UC2
R2
IC2
Рис. 6
На основании законов Кирхгофа получим следующую систему уравнений:
u − uC1
iC1 = iR1 − iL =
− iL ,
R1
u
iC 2 = iL − iR 2 = − C1 + iL ,
R2
u L = u C1 − u C 2 .
11
реакции. Исходя из определения импульсной характеристики,
то есть при x(0), u(t) = δ(t), получим из (10) и (5):
Или, используя (8), (9), получим:
duC1
u
i
u
= − C1 − L +
,
dt
R1C1 C1 C1 + R1
du C 2
u
i
= − C2 + L ,
dt
R2 C 2 C 2
diL uC1 iC 2
=
+
,
dt
L
L
x(t ) = e At ⋅ B ;
Таким образом, в приведенном примере импульсная характеристика для тока будет следующей:

h(t ) = 0

откуда матрицы А и В из выражения (7) для цепи на рис. 5 имеют вид
 1
− R C
 1 1
A= 0

 1

 L
0
1
−
R2C1
1
−
L
1
C1 

1 
,
C2 

0 

−
h(t ) = C ⋅ e At ⋅ B + D ⋅ δ(t ) .
 1 
R C 
 1 1
B =  0 .
 0 


1
R2

 1
0 e At 

 R1C1
T

0 0 .

Преобразуем по Лапласу (6), (7) при нулевых начальных
условиях:
SX ( S ) = AX ( S ) + BU ( S ) ,
Y ( S ) = CX ( S ) + DU ( S ) .
Функция цепи связана с описанием по методу переменных
состояний следующим образом:
H ( S ) = D + C ( SI − A) −1 ⋅ B = L{h(t )}.
Если принять, что выходной реакцией является, например,
ток в резисторе R2, то система (6) сводима к уравнению
uC 2 

1
= 0
0 ⋅ [u C1 u C 2 i L ]T + 0 ⋅ u .
R2 
R2

Таким образом, матрицы C и D имеют вид
iR2 =
5. РАСЧЕТ РЕАКЦИИ ЦЕПИ ПРИ ОДИНОЧНЫХ
ВХОДНЫХ СИГНАЛАХ
В данном разделе рекомендуется использовать операторный
метод анализа электрических цепей [5]. Изображение искомой реакции на выходе F2(S) определяем по выражению


1
C = 0
0 ,
D = 0.
R2


Решение уравнения (6) во временной области имеет вид
F2 ( S ) = F1 ( S ) ⋅ H ( S ) .
Последнему соотношению во временной области соответствует интеграл наложения
t
t
X (t ) = e At ⋅ x(0) +  e A(t − τ) ⋅ Bu (t )dt ,
(10)
0
где eAt – матричная экспонента; x(0) – вектор начальных условий.
Первое слагаемое (10) отвечает реакции при нулевом входе,
а второе – при нулевых начальных условиях. Таким образом, выражение (10) есть сумма свободной и вынужденных составляющих
12
f 2 =  f1 (τ) ⋅ h(t − τ)dt .
0
Для получения изображения одиночных импульсов может
быть использована теорема запаздывания.
Рассмотрим, например, прямоугольный импульс напряжения
(рис. 7).
13
 F при 0 ≤ t ≤ τ ,
f1 =  M
 0 при t > 0 .
H ( jω) =
M ( jω) M (ω) ⋅ e jϕ1 ( ω)
=
=
N ( jω) N (ω) ⋅ e jϕ2 ( ω)
M (ω) j
=
⋅ e (ϕ1 (ω) − φ ⋅ ϕ2 (ω)).
N (ω)
(11)
M (ω)
, назыN (ω)
вается амплитудно-частотной характеристикой цепи, а ϕ(ω) =
= ϕ M (ω) − ϕ N (ω) – фазочастотной характеристикой цепи.
Зависимость модуля H ( jω) от ω, т. е. H ( S ) =
7. ОПРЕДЕЛЕНИЕ СПЕКТРОВ ВХОДНОГО И ВЫХОДНОГО
СИГНАЛОВ В ВИДЕ ОДИНОЧНЫХ ИМПУЛЬСОВ
Рис. 7
f1 (t ) = FM ⋅ δ1 (t ) − F1 ⋅ δ1 (t − τ) .
Преобразование по Лапласу этой суммы с учетом теоремы запаздывания имеет вид
F1 ( S ) =
FM
(1 − e −tS ) .
S
Реакция на выходе цепи находится с помощью теоремы разложения для конкретного случая полюсов F2(S).
6. ОПРЕДЕЛЕНИЕ АМПЛИТУДНО-ЧАСТОТНЫХ
И ФАЗОЧАСТОТНЫХ ХАРАКТЕРИСТИК ЦЕПИ
Входные и передаточные функции цепи представляют собой
рациональные функции переменной S
H (S ) =
M (S )
,
N (S )
где M(S), N(S) – полиномы. При подстановке S = jω получим
14
Спектр импульсного сигнала находят по изображению сигнала, используя связь преобразований Лапласа и Фурье. Для одиночного импульса напряжения u1(t), имеющего изображение U1(s),
комплексный спектр
U1 ( jω) = U1 ( s ) |s = jω .
Модуль функции U1(jω) определяет амплитудный спектр
A1(ω), а аргумент U1(jω) – фазовый спектр φ1(ω).
Так, например, для сигнала в виде прямоугольного импульса
получим комплексный спектр:
π
F
2 F − jω
 τ
F1 ( jω) = M (1 − e − jωt ) = M e 2 ⋅ sin  ω  .
ω
jω
 2
Амплитудный спектр
2F
 π
А1 (ω) = M ⋅ sin  ω 
ω
 2
2πk
при ω =
, где k = ±1, ±2, …
τ
Фазовый спектр
τ
 τ
ϕ1 = arg F1 ( jω) = −ω + arg sin  ω  .
2
 2
Графики амплитудного и фазового спектра показаны на рис. 8.
15
8. РАСЧЕТ ВЫНУЖДЕННОГО РЕЖИМА
ПРИ НЕСИНУСОИДАЛЬНОМ ПЕРИОДИЧЕСКОМ
ВОЗДЕЙСТВИИ
Рис. 8
Спектр выходного сигнала определяется с учетом частотных
характеристик цепи.
F2 ( jω) = F1 ( jω) ⋅ H ( jω) ,
или
F2 ( jω) = A2 (ω) ⋅ e jϕ2 (ω) = A1 (ω) ⋅ A(ω) ⋅ e (ϕ( ω) + ϕ1 (ω)) ,
где A1(ω), ϕ1(ω) – амплитудный спектр и фазовый спектр выходного сигнала соответственно.
Отсюда следует, что
A2 (ω) = A1 (ω) ⋅ A(ω) ,
ϕ 2 (ω) = ϕ1 (ω) + ϕ(ω) .
При построении графиков спектров их целесообразно изображать на одном рисунке.
16
При рассмотрении этого режима сигнал на входе цепи представляют в виде периодической последовательности заданных
в условии задачи импульсов. Период Т принимают равным nτ, где
τ – длительность импульса, а величина n и форма импульса указаны в таблице вариантов.
Для нахождения комплексных амплитуд Аk периодического
сигнала f1(t) можно использовать соотношение, связывающее величины Аk с изображением одиночных импульсов F1(S), совпадающее с периодической функцией f1(t) в пределах периода [0;T].
Известно, что комплексные амплитуды определяются соотношением
2
Ak = F1 ( S ) |S = jkω1 ,
T
2π
где ω1 =
.
T
Очевидно, что связь Аk с комплексным спектром будет аналогичной:
2
Ak = ⋅ F1 ( jω) |S = kω1 .
T
Определим комплексные амплитуды для периодической последовательности прямоугольных импульсов длительностью
T
t= .
2
При k = 0 A = τF ,
π
2τFM − j 2
⋅e ,
π
k = 2 A0 = 0 ,
k = 1 A0 =
π
2τFM − j 2
⋅e
и т. д.
3π
Переходя во временную область, получим ряд Фурье в виде
2τ
τ
π  2τ
π


f1 = FM + FM cos  ω1t −  + FM cos  ω1 3t −  .
2
2  3π
2
π


k = 3 A0 =
17
По полученным выражениям комплексных амплитуд могут
быть построены амплитудный и фазовый спектры периодического
сигнала.
Для нахождения действующего значения напряжения используют формулу [8]
A12 + A2 2 + A3 2 + ...
.
2
Для получения комплексного спектра выходного сигнала используют соотношение
F2 ( jkω1 ) = F1 ( jkω1 ) ⋅ H ( jkω1 ) .
F1 = A2 +
9. АНАЛИЗ ЦЕПИ НА ЭВМ
Анализ цепи необходимо провести, используя различные
прикладные программы: MathCad [2], Electronic WorkBench [3],
MatLab [4, 7]. Целью этого раздела является ознакомление с методами приближенного определения реакций цепи и сопоставление
полученных результатов с аналитическими расчетами, полученными в предыдущих разделах.
Для проведения анализа в среде MathCad необходимо определить матрицы уравнений состояний. Для этого разработана программа MPS; алгоритм программы изложен в [1]. Работа с программой осуществляется в диалоговом режиме. Рассмотрим ее на
примере схемы рис. 2. Зададимся параметрами этой цепи:
R1 = 2 Ом, R2 = 1 Ом, С1 = 1 Ф, С2 = 2 Ф, L = 1 Гн, U = 1 В.
Перед работой с ЭВМ необходимо пронумеровать узлы схемы. Каждый элемент задан параметром, узлом входа и выхода.
Пример диалога
Количество элементов R?
2<>
Параметры и узлы входа и выхода первого элемента R?
2 1 2<>
Параметры и узлы входа и выхода второго элемента R?
1 3 4<>
18
Количество элементов С?
2<>
Параметры и узлы входа и выхода первого элемента С?
1 2 4<>
Параметры и узлы входа и выхода второго элемента C?
2 3 4 <>
Количество элементов L?
1<>
Параметры и узлы входа и выхода первого элемента L?
1 2 3 <>
Количество элементов I ?
0<>
Количество элементов U ?
1<>
Параметры и узлы входа и выхода первого элемента U ?
1 4 1 <>
Результаты расчета
Матрица А
− 0,50 0,00 1,00 
 0,00 − 0,50 0,50


 1,00 − 1,00 1,00 
Матрица В
0,50
0,00


0,00
Выходные уравнения для токов через резисторы R1 и R2
Матрица С
− 1,00 0,00 0,00
 0,00 1,00 0,00


Для реакции y = iR2 в матрицах C и D выбирают последние строки.
19
Имея описание цепи, с помощью уравнений состояний легко
осуществить численный расчет реакции схемы на различные воздействия, так что численное решение уравнений состояния сведется к их численному интегрированию.
Рассмотрим уравнение (7). Чтобы перейти к численной схеме
его решения, в простейшем случае можно воспользоваться методом Эйлера.
xk − xk −1
= Axk −1 + Bu k
h
или итерационная процедура – вид
0
TOL = 0,01,
j = −1,
ω
(12)
Очевидно, что выражение (12) необходимо дополнить начальными условиями x0 = 0; t0 = 0.
Выражение (12) можно принять как базовое для численного
решения уравнений, описывающих цепь. Наиболее просто составить программу для решения дифференциальных уравнений, используя систему MathCad.
Для задания входного воздействия необходимо взять функцию
системы MathCad Хевисайда δ1(t). Для получения функции δ(t) используют приближенное выражение
δ(t ) = 100(δ1 (t ) + δ1 (t − 0,01) ,
что соответствует узкому прямоугольному импульсу единичной
площади.
Отметим, что при D, не равном нулю, невозможно получить
точное значение импульсной характеристики из-за наличия в ней
слагаемого Dδ(t).
При D = 0 рекомендуется учитывать эту матрицу в конце расчета.
Образец программы с пояснениями представлен в прил. 1.
Расчет спектральных характеристик и сигналов в системе
MathCad осуществляется с помощью интегрирования функций
комплексной переменной.
20
inf
F ( f , ω) =  exp(− jωt ) f (t )dt ,
где inf = τ – длительность сигнала; TOL – точность интегрирования, задаваемая от площади сигнала; f(t) – представление сигнала
в программе.
Так, для прямоугольного сигнала f (t ) = A(δ1 (t ) − δ1 (t − τ) .
Пример выполнения программы представлен в прил. 2.
В программе реализован контроль правильности полученных
результатов с помощью обратного преобразования Фурье
dx Δx xk − xk −1
=
=
,
dt Δt
h
где h – шаг по времени, принимаемый постоянным.
Тогда уравнение (7) примет вид
xk = xk −1 + h( Axk −1 + Bu k ) .
Зададим интеграл Фурье в виде
2 гр
f (t ) =  A(ω) cos(ω t )dt ,
π 0
где ωгр – верхняя граница амплитудного спектра.
Расчет частотных характеристик цепи можно произвести
в системе MathCad, используя выражение
H ( jω) = H ( S ) S = jω ,
где H ( S ) = C ( IS − A) −1 ⋅ B + D .
Амплитудно-частотная и фазочастотная характеристики будут иметь вид A(ω) = |H(jω)|; Ф(ω) = arg H ( jω) .
Целесообразно сопоставить полученные результаты с результатами моделирования в системе Electronic WorkBench [3] при
применении Боде-плоттера.
10. РЕКОМЕНДАЦИИ К ВЫПОЛНЕНИЮ
И ОФОРМЛЕНИЮ КУРСОВОГО РАСЧЕТА
Курсовую работу оформляют в виде пояснительной записки
с титульным листом, на котором указывают министерство, институт, кафедру, наименование курсовой работы, исполнителя и фамилию преподавателя, а также место и год выполнения. На первом
21
листе пояснительной записки должно быть оглавление, на втором –
техническое задание; порядок расчета соответствует разделам
данного методического руководства; в конце пояснительной записки следует привести заключение (краткие выводы по курсовой
работе с анализом полученных результатов) и список использованной литературы.
Расчеты следует сопровождать ясными и краткими комментариями. Пояснительная записка должна иметь рубрикацию с выделением соответствующих разделов работы. К каждому разделу
прилагаются листинги программ на ЭВМ с соответствующими результатами расчетов.
Работа выполняется на листах формата А4 в печатном виде.
22
Рекомендуемая литература
1. Чуа Л. О. Машинный анализ электронных схем: алгоритмы и вычислительные методы / Л. О. Чуа, Лин Пен-Мин. – М. : Энергия, 1980. – 640 с.
2. Дьяконов В. П. Система MathCad / В. П. Дьяконов. – М. : Радио
и связь, 2001. – 127 с.
3. Электротехника и электроника в экспериментах и упражнениях.
Практикум на Electronic WorkBench / под общ. ред. Д. И. Панфилова. – М. :
Додэка, 1999. – Т. 1. – 304 с.
4. Новгородцев А. Б. Расчет электрических цепей в MatLab: учеб. курс /
А. Б. Новгородцев. – СПб. : Питер, 2004. – 250 с.
5. Новгородцев А. Б. 30 лекций по теории электрических цепей /
А. Б. Новгородцев. – СПб. : Политехника, 1995. – 520 с.
6. Демирчян К. С. Теоретические основы электротехники. – 4-е изд. /
К. С. Демирчян, Л. Р. Нейман, Н. Б. Коровкин, В. Л. Чечурин. – СПб. : Питер, 2003. – Т. 1. – 424 с.
7. Дьяконов В. П. MatLab 6. – СПб. : Питер, 2001. – 592 с.
8. Бычков Ю. А. Основы теории электрических цепей / Ю. А. Бычков,
В. М. Золотницкий, Э. П. Чернышев. – СПб. : Лань, 2004. – 464 с.
23
ПРИЛОЖЕНИЯ
Приложение 2
Приложение 1
Программа преобразования Фурье
Программа расчета переходного процесса
− 0.19 − 0.06
0.19
, B=
A=

 , C = [− 0.75
− 0.25 − 0.75
0.25
определяющие цепь в пространстве состояний
h : = 0.01
n : = 2000
T : = 0.02
i : = 0…n – 1
– шаг интегрирования
– количество шагов
– параметр воздействия
– начало цикла интегрирования
функция воздействия:

 T 
u (t ) = (t ) − t −   ⋅ 100
 2 

F1 ( x1, x 2, t ) = A0,0 ⋅ x1 + A0,1 ⋅ x 2 + B0 ⋅ u (t )
F2 ( x1, x 2, t ) = A1,0 ⋅ x1 + A1,1 ⋅ x 2 + B0 ⋅ u (t )
задание начальных условий:
 t 0  0
 x1  : = 0 ,
 0  
 x 20  0
− 0.20] = 0.25 – матрицы,
inf : = 1
– верхняя граница интегрирования
TOL : = 0.05 – точность интегрирования
– задание комплексной единицы
j = −1
inf
Fur ( f , ω) =  exp(− jωt ) f (t )dt
0
n : = 300
dw : = 0.1
k : = 0…n
Wk : = k dw
T:=1
Форма сигнала:


 T  T
f (t ) = 10δ(t )t − 2δ t −  ⋅ t −  + δ(t − T )(t − T ) ,
 2  2


yk = Fur[ f , wk ] ,
uk = if [ [ yk
численное интегрирование по методу Эйлера:
t+n
 t i +1  

 x1  : =  x1 + h ⋅ F [x1 , x 2 , t ]  ,
1
i
i
 i +1  

 x 2i +1   x 2 + h ⋅ F2 [x1i , x 2i , t ]
]] > 0,
arg[ yk ], 0 ,
bk = Re [ yk ] .
Амплитудный спектр
определение реакции:
 x1 
y (t ): = C  i  + D ⋅ u[t i ] .
 x2 i 
24
25
Фазовый спектр
Оглавление
3.2
u
k
-3.2
0
w
30
k
26
Цель и задачи курсовой работы…………………………………………..
Задание к курсовой работе………………………………………………...
1. Определение входных и передаточных функций цепи,
их нулей и полюсов…………………………………………………….
2. Определение параметров четырехполюсника и их связь
с функциями цепи……………………………………………………..
3. Переходная и импульсная характеристики цепи…………………
4. Система уравнений по методу пространства состояний…………
5. Расчет реакции цепи при одиночных входных сигналах…………
6. Определение амплитудно-частотных и фазочастотных
характеристик цепи……………………………………………………
7. Определение спектров входного и выходного сигналов
в виде одиночных импульсов…………………………………………
8. Расчет вынужденного режима при несинусоидальном
периодическом воздействии…………………………………………..
9. Анализ цепи на ЭВМ………………………………………………..
10. Рекомендации к выполнению и оформлению курсового
расчета………………………………………………………………….
Рекомендуемая литература…….………………………………………….
Приложения ………………………………………………………………..
27
3
4
6
8
9
10
13
14
15
17
18
21
23
24
Учебное издание
АНАЛИЗ ЭЛЕКТРИЧЕСКИХ СХЕМ
Методические указания
Составители: Резниченко Виктор Васильевич,
Соколов Владимир Яковлевич
Редактор А. В. Афанасьева
Корректор М. А. Молчанова
Компьютерная верстка И. А. Яблоковой
Подписано к печати 28.12.2017. Формат 60×84 1/16. Бум. офсетная.
Усл. печ. л. 1,6. Тираж 100 экз. Заказ 175. «С» 134.
Санкт-Петербургский государственный архитектурно-строительный университет.
190005, Санкт-Петербург, 2-я Красноармейская ул., д. 4.
Отпечатано на ризографе. 190005, Санкт-Петербург, ул. Егорова, д. 5/8, лит. А.
28
Документ
Категория
Без категории
Просмотров
0
Размер файла
391 Кб
Теги
analiz, resnichenko, elektro
1/--страниц
Пожаловаться на содержимое документа