|
|
F07~=0.443*Fdiskr/N. И спад АЧХ (боковиков) 6 дб на октаву.
Аналогичный спад АЧХ (уже не боковиков) имеет IIR первого порядка с разностным уравнением y(n)=G*x(n)+R*y(n-1). R от 0 до 1.
Означенное звено имеет перед. ф-ию H(Z)=G/(1-R*Z^(-1)).
Находите квадрат модуля АЧХ этого звена, т.е. кв. модуля указанной перед. ф-ии на единичной окружности в Z-плоскости, т.е.
при Z=exp(j*2pi*f/Fdiskr). Припавниваете 3-х дб полосу звена к F07=0.443*Fdiskr/N и находите отсюда R. G -из нормировки.
Обычно берут немного измненное звено, не как выше приведено. Вводят нуль на половине частоты дискретизации.
Перед. ф-ия такого звена H1(Z)=(G1*(1+Z^(-1)))/(1-R1*Z^(-1)).
Разностное уравнение такого звена y(n)=G1(x(n)+x(n-1))+R1*y(n-1)
Аналогично - находите 3-х дб.полосу и прравниваеие ее к F07, откуда находите R1. Его значение близко к R будет.
Звено, аналогичное H1, автоматом получается из обычной интегрирующей цепочки? с передаточной ф-ей HA(Sn)=1/(1+Sn), через билинейное z-преобр.
Sn - это нормированная лапласовкая переменная (n не индекс:-)).
Sn=S*tau. Из равенства 3-х дб полос: tau=N/(2pi*0.443*Fdiskr).
Билинейное z-преобразование: Sn=B*(1-Z^(-1))/(1+Z^(-1)).
К-т B=ctg(pi*F07/Fdiskr)=ctg(pi*0.443/N).
Бил. Z-пр. даст H2(Z)=G2*(1+Z^(-1)))/(1-R2*Z^(-1)).
Ну и получается, что R2=(B-1)/(B+1), а G2=1/(B+1).
R2 близко к R1 и к R
Ну и на G2 можно дополнит. условия нормировки наложить.
|
|
E-mail: 
info@telesys.ru
