Курсовая ТОЭ 3 курс
.pdfПрименим метод контурных токов.
|
|
|
|
|
|
|
|
iк |
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
1 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
к |
|
|
|
к |
|
к |
|
к |
|
|
uC |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
R1i1 |
R1i2 0i3 |
0i4 |
|
|
(4.1) |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
iк i |
L1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
0iк |
0i |
к |
0i |
к R iк |
|
|
u |
C2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
1 |
|
|
2 |
3 |
|
|
|
|
|
2 4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
i к |
i |
|
(4.2) |
|
|
|
|
iк |
R1i1 uC1 |
|
(4.3) |
|
iк i |
|
(4.4) |
i к |
uC2 |
(4.5) |
||||||||||||||||||||
|
|
|
|
|
|
|
|
L |
|
|||||||||||||||||||||||||||||
1 |
|
|
1 |
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
R1 |
|
|
|
|
|
3 |
|
|
|
4 |
R2 |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
к |
|
|
к |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
iC1 |
i2 |
i3 |
i1 |
|
|
|
|
|
uC1 |
|
|
iL1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
R1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
к |
|
|
к |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
iC |
i3 i4 iL |
|
|
|
|
|
|
|
uC |
|
|
(4.6) |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
R2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
2 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
uL |
|
uC |
uC |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
1 |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Используя соотношения i |
|
C u , |
i |
C u |
, |
u |
L1 |
L i получим уравнения состояния: |
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
C1 |
|
1 C1 |
C2 |
2 C2 |
|
|
1 L1 |
|
|
|
||||
|
|
u u |
0u |
C2 |
i |
L1 |
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
C1 |
|
|
|
C1 |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
0uC uC |
|
iL |
|
|
0i1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
(4.7) |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
uC |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
i |
|
|
|
|
|
1 |
|
|
|
2 |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
0.5u |
C1 |
0.5u |
C2 |
0i |
L1 |
0i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
L1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Уравнения состояния в матричной форме: |
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
1 |
|
0 |
|
|
|
|
1 |
|
uC |
|
1 |
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
uC |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
0 |
1 |
|
|
1 |
uC2 |
|
0 i1 |
(4.8) |
|
|
|
|
|
|
|
|
|||||||||||||||||
uC2 |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
i |
|
|
|
0.5 |
0.5 |
0 |
|
|
i |
L1 |
|
0 |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
L1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Характеристический полином: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
1 p |
|
|
|
|
|
0 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
det A p E |
|
|
|
0 |
|
|
|
|
|
1 p |
|
|
1 |
p3 2 p2 |
2 p 1 0 |
(4.9) |
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0.5 |
|
|
|
|
|
|
0.5 |
|
|
p |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Корни характеристического полинома – частоты собственных колебаний цепи:
p1 1 |
p 0.5 j0.5 |
3 0.5 j0.866 |
|
2,3 |
|
Так как корни характеристического полинома равны полюсам передаточной функции – уравнения состояния цепи составлены правильно.
11
3. Проектирование ЛЦФ методом соответствия переходных характеристик
3.1. Определение частоты дискретизации
Частота дискретизации д 2 определяется с учётом четырёх критериев для максимальной частоты из учитываемых частот:
1)при использовании однопроцентного критерия ширины спектра импульсной характеристики > ∆ с = 4.64151;
2)по критерию удовлетворительного описания минимального временного интервала
процессов цепи, то есть |
|
> 2 ⁄ |
|
= 31,41593; |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
||
3) при использовании условия устойчивости ЛЦФ в случае комплексного полюса |
в ПФ |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
фильтра-прототипа (для |
явной и |
|
смешанной |
форм алгоритма |
Эйлера) |
|
|
> |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
2 | |2⁄(−2 ( )) = 6,28321. |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
Таким образом, |
= 40. |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
Выбираем |
«удобную» |
частоту |
дискретизации |
= 1256 2 |
, то есть |
период |
||||||
|
|
|
|
|
|
|
д |
|
|
|
|
|
дискретизации |
= 2 ⁄ |
= 0,005 = |
|
|
⁄40. |
|
|
|
|
|
|
|
д |
д |
|
|
|
|
|
|
|
|
|
||
|
3.2. Расчёт дискретной передаточной функции |
|
|
|
|
|
||||||
Определяем |
ПФ ЛЦФ |
методом |
|
совпадения аналоговой 1а( ) |
и |
дискретной |
1д( ) |
переходных характеристик в дискретные моменты времени = , где = д = 0,005. Для ≥ 0 имеем
1а( ) = 0,5 − 0,5 − + 0,57734 −0,5 cos(0,86603 + 1,5708) =
=0,5 − 0,5(1 − ) + 0,57734(1 − 0,5 ) cos(0,86603 + 1,5708)
=0,5 − 0,5 ∙ 0,995 + 0,57734 ∙ 0,9975 cos(0,00433 + 1,5708) = 1д( )
Далее находим
0,5 0,51( ) = − 1 − − 0,995 +
+0,57734 |
( − 0,9975 ∙ cos(0,00433)) cos(1,5708) − 0,9975 sin(0,00433) sin(1,5708) |
||||||||
|
|
|
|
||||||
|
|
|
|
|
2 − 2 ∙ 0.99749 + 0,99752 |
|
|
||
|
1( ) = |
0,5 |
− |
0,5 |
+ |
−0,00249 |
|
= |
|
|
|
|
2 − 1.99498 + 0.995006 |
||||||
|
|
|
− 1 |
− 0,995 |
|
|
|||
|
= |
|
0,00001 3 − 0,00001 2 + 0,00001 |
= |
|
||||
|
( 3 − 2,98998 2 + 2,98001 − 0,99003)( − 1) |
|
0,00001 3 − 0,00001 2 + 0,00001 = ( − 1)( − 0,995)( 2 − 1.99498 + 0.995006).
После чего определяем ПФ ЛЦФ
12
( ) = |
− 1 |
1( ) = |
0,00001 2 − 0,00001 + 0,00001 |
= |
||
|
( − 0,995)( 2 |
− 1.99498 + 0.995006) |
||||
|
|
|
0,00001 2 − 0,00001 + 0,00001 = 3 − 2,98998 2 + 2,98001 − 0,99003.
3.3. Построение схемы ЛЦФ
Построение схемы ЛЦФ осуществляется на основании выражения для передаточной функции дискретной системы. Схема ЛЦФ изображена на рисунке 10.
0,00001 2 − 0,00001 + 0,00001( ) = 3 − 2,98998 2 + 2,98001 − 0,99003 =
0,00001 −1 − 0,00001 −2 + 0,00001 −3 = 1 − 2,98998 −1 + 2,98001 −2 − 0,99003 −3.
Рис.10. Схема ЛЦФ.
3.4. Численный контроль переходной характеристики ЛЦФ
Осуществим контроль переходной характеристики ЛЦФ как по выражению, так и по передаточной функции.
При контроле по выражению, используем следующую формулу
1д( ) = 0,5 − 0,5 ∙ 0,995 + 0,57734 ∙ 0,9975 cos(0,00433 + 1,5708).
13
При контроле по передаточной функции используем разностное уравнение, получаемое из передаточной функции:
( ) = |
( ) |
0,00001 −1 − 0,00001 −2 + 0,00001 −3 |
||
|
= |
|
→ |
|
( ) |
1 − 2,98998 −1 + 2,98001 −2 − 0,99003 −3 |
→ ( )(1 − 2,98998 −1 + 2,98001 −2 − 0,99003 −3) = = ( )(0,00001 −1 − 0,00001 −2 + 0,00001 −3).
Применяя обратное z-преобразование, получаем
( ) = 0,00001 ( − 1) − 0,00001 ( − 2) + 0,00001 ( − 3) +
+2,98998 ( − 1) − 2,99151 ( − 2) + 0,99003 ( − 3)
В таблице три приведены значения переходной характеристики для 30 первых точек.
Таблица 3
|
1( ) |
1д( ) |
∆1 |
Разностное |
∆2 |
|
уравнение |
||||||
|
|
|
|
|
||
|
|
|
|
|
|
|
1 |
0 |
0 |
0 |
0 |
0 |
|
2 |
0 |
0.00001 |
-0.00001 |
0.00001 |
-0.00001 |
|
3 |
0 |
0.00001 |
-0.00001 |
0.00001 |
-0.00001 |
|
4 |
0 |
0.00002 |
-0.00002 |
0.00002 |
-0.00002 |
|
5 |
0 |
0.00003 |
-0.00002 |
0.00003 |
-0.00002 |
|
6 |
0 |
0.00003 |
-0.00003 |
0.00003 |
-0.00003 |
|
7 |
0 |
0.00004 |
-0.00004 |
0.00004 |
-0.00004 |
|
8 |
0 |
0.00005 |
-0.00004 |
0.00005 |
-0.00004 |
|
9 |
0.00001 |
0.00005 |
-0.00005 |
0.00005 |
-0.00005 |
|
10 |
0.00001 |
0.00006 |
-0.00005 |
0.00006 |
-0.00005 |
|
11 |
0.00001 |
0.00007 |
-0.00006 |
0.00007 |
-0.00006 |
|
12 |
0.00001 |
0.00008 |
-0.00007 |
0.00008 |
-0.00007 |
|
13 |
0.00002 |
0.00009 |
-0.00007 |
0.00009 |
-0.00007 |
|
14 |
0.00002 |
0.0001 |
-0.00008 |
0.0001 |
-0.00008 |
|
15 |
0.00003 |
0.00011 |
-0.00008 |
0.00011 |
-0.00008 |
|
16 |
0.00003 |
0.00012 |
-0.00009 |
0.00012 |
-0.00009 |
|
17 |
0.00004 |
0.00014 |
-0.00009 |
0.00014 |
-0.00009 |
|
18 |
0.00005 |
0.00015 |
-0.0001 |
0.00015 |
-0.0001 |
|
19 |
0.00006 |
0.00016 |
-0.00011 |
0.00016 |
-0.00011 |
|
20 |
0.00007 |
0.00018 |
-0.00011 |
0.00018 |
-0.00011 |
|
21 |
0.00008 |
0.0002 |
-0.00012 |
0.0002 |
-0.00012 |
|
22 |
0.00009 |
0.00021 |
-0.00012 |
0.00021 |
-0.00012 |
|
23 |
0.00011 |
0.00023 |
-0.00013 |
0.00023 |
-0.00013 |
|
24 |
0.00012 |
0.00025 |
-0.00013 |
0.00025 |
-0.00013 |
|
25 |
0.00014 |
0.00027 |
-0.00014 |
0.00027 |
-0.00014 |
|
26 |
0.00015 |
0.0003 |
-0.00014 |
0.0003 |
-0.00014 |
|
27 |
0.00017 |
0.00032 |
-0.00015 |
0.00032 |
-0.00015 |
|
28 |
0.00019 |
0.00035 |
-0.00015 |
0.00035 |
-0.00015 |
|
29 |
0.00021 |
0.00037 |
-0.00016 |
0.00037 |
-0.00016 |
|
30 |
0.00024 |
0.0004 |
-0.00016 |
0.0004 |
-0.00016 |
|
|
|
|
14 |
|
|
Как видно из таблицы 3, полученные значения практически совпадают, а, значит, найденные
нами выражения верны.
4. Проектирование ЛЦФ методом использования уравнений
состояния
4.1. Определение дискретной передаточной функции ЛЦФ, синтезированного явным
алгоритмом Эйлера
Расчёт ЛЦФ произведём явным алгоритмом Эйлера.
( ) = ( )| =( −1)⁄ =200( −1) =
0,5 = (200( − 1))3 + 2(200( − 1))2 + 2(200( − 1)) + 1 =
0.0000000625 = 3 − 2,99 2 + 2,98005 − 0,99005.
Как видно из формулы, передаточная функция ЛЦФ, синтезированного методом Эйлера,
отличается от передаточной функции, полученной из переходной характеристики.
Переходная характеристика ЛЦФ
|
1( ) = ( ) |
|
|
= |
|
0,0000000625 |
|
|
|
|
. |
|
||||||||||
|
− 1 |
( − 1)( 3 − 2,99 2 + 2,98005 − 0,99005) |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
Определим полюса звена третьего порядка |
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
3 − 2,99 2 + 2,98005 − 0,99005 = 0 |
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
1 = 0,995; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2,3 = 0.9975 ± 0.00433; |
|
|
|
|
|
|
|
|
|||
= + |
1 |
+ |
2 |
|
+ |
3 |
|
+ |
4 |
|
|
|
. |
|||||||||
|
|
|
− 0,995 |
− 0.9975 − 0.00433 |
− 0.9975 + 0.00433 |
|||||||||||||||||
0 |
|
− 1 |
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
0 = 1(0) = 0. |
|
|
|
|
|
|
|
|
|
1 = lim |
|
|
|
|
|
|
|
|
|
|
0,0000000625 |
|
|
|
|
|
|
= 0,5. |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
→1 ( − 0,995)( − 0.9975 − 0.00433)( − 0.9975 + 0.00433) |
|
|
|
||||||||||||||||||
2 = |
lim |
|
|
|
|
|
|
|
|
|
0,0000000625 |
|
|
|
|
|
= −0,5. |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
→0,995 ( − 1)( − 0.9975 − 0.00433)( − 0.9975 + 0.00433) |
|
|
|
|
|||||||||||||||||
3 = |
|
|
lim |
|
|
|
|
|
|
|
0,0000000625 |
|
|
|
|
= 0.28868. |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
→0.9975+ 0.00433 ( − 1)( − 0,995)( − 0.9975 + 0.00433) |
|
|
|
|
|
||||||||||||||||
4 = |
|
|
lim |
|
|
|
|
|
|
|
0,0000000625 |
|
|
|
= − 0.28868. |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
→0.9975− 0.00433 ( − 1)( − 0,995)( − 0.9975 − 0.00433) |
|
|
|
|
|
|
Обратное z-преобразование суммы при наличии комплексно-сопряжённых полюсов может
быть найден следующим образом:
15
−1 |
[ |
3 |
|
|
+ |
4 |
|
] = |
|
− (− + ) |
− (− − ) |
||||||||
|
|
|
|
=2| 3||− + | cos(arg(− + ) + arg 3) 1( ).
Внашем случае | 3| = | 0.28868| = 0.28868, |0.9975 + 0.00433| = 0.99751, arg(0.9975 +
0.00433) = 0.00434, arg ̇3 = 1.5708.
Таким образом, переходная характеристика имеет
1( ) = (0,5 − 0,5 ∙ 0,995 + 0.57736 ∙ 0.99751 cos(0,00434 + 1.5708)) 1( ).
4.2. Построение схемы ЛЦФ, полученной на основе явного алгоритма Эйлера
Построение схемы ЛЦФ осуществляется на основании выражения для передаточной функции дискретной системы. Схема ЛЦФ изображена на рисунке 11.
0.0000000625( ) = 3 − 2,99 2 + 2,98005 − 0,99005 =
0,0000000619 −3 = 1 − 2,99 −1 + 2,98005 −2 − 0,99005 −3.
Рис.11. Схема ЛЦФ, синтезированного явным алгоритмом Эйлера.
16
4.1. Определение дискретной передаточной функции ЛЦФ, синтезированного смешанным
алгоритмом Эйлера
Расчёт ЛЦФ произведём смешанным алгоритмом Эйлера.
( ) = ( )| =( −1)⁄ =200( −1) =
0,5 = (200( − 1))3 + 2(200( − 1))2 + 2(200( − 1)) + 1 =
0.0000000625 = 3 − 2,99 2 + 2,98005 − 0,99005.
Как видно из формулы, передаточная функция ЛЦФ, синтезированного методом Эйлера,
отличается от передаточной функции, полученной из переходной характеристики.
Переходная характеристика ЛЦФ
1( ) = ( ) |
|
= |
|
0,0000000625 2 |
. |
||
− 1 |
( − 1)( 3 |
− 2,99 2 |
+ 2,98005 − 0,99005) |
||||
|
|
|
Определим полюса звена третьего порядка
|
|
|
3 − 2,99 2 + 2,98005 − 0,99005 = 0 |
|
|||||||
|
|
|
|
|
1 = 0,995; |
|
|
|
|||
|
|
|
|
2,3 = 0.9975 ± 0.00433; |
|
|
|
||||
= + |
1 |
+ |
2 |
+ |
3 |
|
|
+ |
4 |
. |
|
− 1 |
− 0,995 |
− 0.9975 − |
0.00433 |
− 0.9975 + 0.00433 |
|||||||
0 |
|
|
|
|
0 = 1(0) = 0.
0,0000000625
1 = lim = 0,5.→1 ( − 0,995)( − 0.9975 − 0.00433)( − 0.9975 + 0.00433)
2 = |
lim |
|
0,0000000625 |
|
|
= −0,5. |
||
|
|
|
|
|
||||
|
|
|
|
|
||||
|
→0,995 ( − 1)( − 0.9975 − 0.00433)( − 0.9975 + 0.00433) |
|
||||||
3 = |
lim |
0,0000000625 |
|
= 0.28868. |
||||
|
|
|
||||||
|
|
|
||||||
|
→0.9975+ 0.00433 ( − 1)( − 0,995)( − 0.9975 + 0.00433) |
|
|
|||||
4 = |
lim |
0,0000000625 |
= − 0.28868. |
|||||
|
|
|||||||
( − 1)( − 0,995)( − 0.9975 − 0.00433) |
||||||||
→0.9975− 0.00433 |
|
|
|
Обратное z-преобразование суммы при наличии комплексно-сопряжённых полюсов может
быть найден следующим образом:
−1 |
[ |
3 |
|
|
+ |
4 |
|
] = |
|
− (− + ) |
− (− − ) |
||||||||
|
|
|
|
=2| 3||− + | cos(arg(− + ) + arg 3) 1( ).
Внашем случае | 3| = | 0.28868| = 0.28868, |0.9975 + 0.00433| = 0.99751, arg(0.9975 +
0.00433) = 0.00434, arg ̇3 = 1.5708.
Таким образом, переходная характеристика имеет
17
1( ) = (0,5 − 0,5 ∙ 0,995 + 0.57736 ∙ 0.99751 cos(0,00434 + 1.5708)) 1( ).
4.2. Построение схемы ЛЦФ, полученной на основе смешанного алгоритма Эйлера
Построение схемы ЛЦФ осуществляется на основании выражения для передаточной функции дискретной системы. Схема ЛЦФ изображена на рисунке 10.
0.0000000625( ) = 3 − 2,99 2 + 2,98005 − 0,99005 =
0,0000000625 −2 = 1 − 2,99 −1 + 2,98005 −2 − 0,99005 −3.
Рис.12. Схема ЛЦФ, синтезированного явным алгоритмом Эйлера.
18
4.5. Сравнение данных расчёта переходных характеристик ЛЦФ и фильтра-прототипа
Произведём сравнение переходных характеристик, полученных непосредственно из переходной характеристики аналогового прототипа, явного метода Эйлера и смешанного метода Эйлера. В таблице 4 показаны значения для первых 30 отсчётов переходной характеристики, а на рисунке 13 графики переходных характеристик.
Таблица 4
|
1д( ) |
Явный |
∆ |
Смешанный |
∆ |
|
|
|
|
|
|
0 |
0 |
0 |
0 |
0 |
0 |
4 |
0.00003 |
0 |
0.00003 |
0 |
0.00003 |
8 |
0.00005 |
0 |
0.00005 |
0 |
0.00005 |
12 |
0.00009 |
0.00001 |
0.00008 |
0.00001 |
0.00008 |
16 |
0.00014 |
0.00003 |
0.0001 |
0.00003 |
0.0001 |
20 |
0.0002 |
0.00007 |
0.00013 |
0.00007 |
0.00013 |
24 |
0.00027 |
0.00012 |
0.00015 |
0.00012 |
0.00015 |
28 |
0.00037 |
0.00019 |
0.00018 |
0.00019 |
0.00018 |
32 |
0.00049 |
0.00029 |
0.00021 |
0.00029 |
0.00021 |
36 |
0.00064 |
0.00041 |
0.00023 |
0.00041 |
0.00023 |
40 |
0.00082 |
0.00056 |
0.00026 |
0.00056 |
0.00026 |
44 |
0.00103 |
0.00074 |
0.00029 |
0.00074 |
0.00029 |
48 |
0.00127 |
0.00096 |
0.00031 |
0.00096 |
0.00031 |
52 |
0.00156 |
0.00122 |
0.00034 |
0.00122 |
0.00034 |
56 |
0.00188 |
0.00151 |
0.00036 |
0.00151 |
0.00036 |
60 |
0.00224 |
0.00185 |
0.00039 |
0.00185 |
0.00039 |
64 |
0.00264 |
0.00223 |
0.00041 |
0.00223 |
0.00041 |
68 |
0.00309 |
0.00265 |
0.00044 |
0.00265 |
0.00044 |
72 |
0.00359 |
0.00312 |
0.00046 |
0.00312 |
0.00046 |
76 |
0.00413 |
0.00365 |
0.00049 |
0.00365 |
0.00049 |
80 |
0.00473 |
0.00422 |
0.00051 |
0.00422 |
0.00051 |
84 |
0.00537 |
0.00484 |
0.00053 |
0.00484 |
0.00053 |
88 |
0.00607 |
0.00552 |
0.00056 |
0.00552 |
0.00056 |
92 |
0.00683 |
0.00625 |
0.00058 |
0.00625 |
0.00058 |
96 |
0.00764 |
0.00703 |
0.0006 |
0.00703 |
0.0006 |
100 |
0.0085 |
0.00788 |
0.00063 |
0.00788 |
0.00063 |
104 |
0.00943 |
0.00878 |
0.00065 |
0.00878 |
0.00065 |
108 |
0.01041 |
0.00974 |
0.00067 |
0.00974 |
0.00067 |
112 |
0.01145 |
0.01076 |
0.00069 |
0.01076 |
0.00069 |
116 |
0.01255 |
0.01184 |
0.00071 |
0.01184 |
0.00071 |
19
Рис. 13. Графики переходных характеристик ЛЦФ Из графиков для переходных характеристик мы видим, что переходные характеристики не
отличаются. Это значит, что методом Эйлера и методом билинейного преобразования достаточно точно можно синтезировать по прототипу ЛЦФ.
20