Лекция 12. Из принципа суперпозиции следует, что реакция ДЛС на произвольное воздействие может быть вычислена по
формуле свертки.
N 1
y nT h vT X nT vT
v 0
Нас интересует гармоническое входное воздействие:
|
|
|
|
|
|
|
Aвх |
cos nT DLS Aвых cos nT 0 |
|
||
X nT |
|
|
|
||
|
Aвх |
sin nT DLS Aвых |
sin nT 0 |
|
|
|
|
|
|
|
|
Составим из пар входов и выходов комплексные сигналы:
|
nT |
|
A |
cos |
nT |
|
|
|
|
isin |
|
nT |
|
|
Aвых |
ФЧХ |
|||||||||
|
|
0 |
0 |
|
|||||||||||||||||||||
y |
|
вых |
|
|
|
|
|
|
|
|
|
|
|
|
ei 0 H |
||||||||||
X |
|
nT |
|
|
|
|
вх |
|
nT |
|
isin |
|
nT |
|
|
|
|
||||||||
|
|
|
|
|
|
|
вх |
||||||||||||||||||
|
|
|
|
|
A cos |
|
|
|
|
|
|
|
|
A |
|||||||||||
АЧХ
Докажем, что h vT F 1 H
В формулу свертки подставим комплексный гармонический входной сигнал:
|
N 1 |
|
|
i nT vT |
|
|
|
i nT |
N 1 |
|
i vT |
|
|
|
||
|
|
|
e |
e |
h vT e |
|
|
|
||||||||
y nT h vT |
|
|
|
|
|
|
|
|||||||||
|
v 0 |
|
|
|
|
|
|
|
|
v 0 |
|
|
|
|
|
|
|
|
|
|
|
|
N 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
nT |
|
|
ei nT h vT |
e i vT |
N 1 |
|
|
|
|||||
|
|
|
|
v 0 |
|
|
|
|
|
|
|
|
|
|||
H |
|
y |
|
|
|
|
|
h vT e i vT F h vT |
, ч.т.д. |
|||||||
X |
|
|
e |
i nT |
|
|||||||||||
|
|
nT |
|
|
|
|
|
|
v 0 |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
*
Соотношение (*) нужно для доказательства. На практике снять комплексную характеристику с реального устройства не представляется возможным.
0 2 з f 2 з
T
Пример расчета КИХ ФНЧ.
На входе заданная частотная характеристика фильтра. Поскольку фильтр цифровой – периодический.
АЧХ |
|
|
40дб |
fн |
fS |
60дб |
|
|
|
ФЧХ |
|
Если мы хотим сделать действительный фильтр, надо сделать исходную АЧХ периодической и симметричной относительно нулевой частоты – относительно fS .
Задание: |
fП |
fПT 0.25 |
N 8 |
|
|||
|
fS |
|
|
Im0
1
2 3 4 5 6 7
Re 
© Дорошенко Е., Подкопаев И. |
1 |
Если исходная функция действительная и симметричная – спектр реальный и симметричный. Убедимся в том, что фильтр готов:
Если N увеличить до 1024, то получим:
Попробуем вырезать центральную часть и сделать БПФ для N 8 , получим
Если опять взять N 1024, то картина сожмется к краям:
Появились колебания Гиббса.
Поскольку мы обнулили серединные точки, остальные значения образуют симметричную функцию относительно n 0, следовательно ДПФ опять будет реальной, т.е. Im 0, 0, т.е. остальные значения должны
быть симметричными. |
|
|
|
|
|
|
||||
Сместим импульсную характеристику на |
NФ 1 |
дискретов, |
где |
NФ |
- число оставшихся необнуленными |
|||||
|
||||||||||
2 |
|
|
|
|
|
|
|
|||
значений фильтра. В нашем случае |
5 1 |
2. |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
||||
2 |
|
|
|
|
|
|
|
|
|
|
При этом: |
|
|
|
|
|
|
||||
- АЧХ не изменится |
|
|
|
|
|
|
||||
- ФЧХ получит дополнительное смещение e i з |
i |
N 1 |
T |
i |
|
|
||||
e |
2 |
|
e |
2 |
|
|||||


2
Получили линейную ФЧХ, линейность не изменяет входной сигнал.
X |
nT |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
D |
|
|
D |
|
D |
|
|
D |
|
|
|
|
|
||||||
|
|
|
|
1 |
|
|
2 |
|
|
|
3 |
|
|
4 |
|
|
|
|
|
|
1 |
|
2.4 |
|
|
|
5 |
|
|
|
2.4 |
|
|
|
1 |
||||
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
y nT |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
© Дорошенко Е., Подкопаев И. |
2 |
Снижение колебаний Гиббса с помощью взвешенной усеченной импульсной характеристики.
|
|
n |
1 |
|
|
|
|
|
Как и в случае ЛЧМ воспользуемся окном Хемминга g v 0.54 0.46cos2 |
|
|
, |
n 0..N |
Ф |
1 |
||
|
|
|||||||
|
|
Nф 1 |
2 |
|
|
|
||
|
|
|
|
|
|
|
||
h vT 1 2.4 |
5 2.4 1 |
|
|
|
|
|
|
|
g v 0.08 0.5410.54 0.08
hH vT 0.081.3 5 1.3 0.08
Если этими коэффициентами зарядить фильтр АЧХ примет вид
Кривая завалилась, но колебания Гиббса уменьшились.
8.5. Расчет КИФ ФВЧ, ФП, ФР.
ФВЧ
0 1 2 3 4 5 6 7
ФП
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
FД1
0 1 2 3 4 5 6 7
FД1
0 1 2 3 4 5 6 7
ФР
FД1
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
Задача АЧХ.
1.АЧХ в дискретном виде.
2.Вычисление ОБПФ – получаем h vT .
3.Выбираем длину фильтра по заданной переходной характеристике фильтра
4. Смещаем импульсную характеристику по кругу на Nф 1 шагов. 2
5.Домножаем h vT g v
6.Вычисляем результат через БПФ.
© Дорошенко Е., Подкопаев И. |
3 |