сдвиги на этой частоте различались как можно больше. Как и следовало ожидать, это имеет место для параметров, отличающихся больше всего. Пусть это <ц = 0.9 и 02 = —0.9. Следующий шаг нахождение фазовых сдвигов в терминах отсчетов. Пусть исходный контейнер записан с частотой стробирования Fr = 22050. Как и в скрипте 18, вычисляем значения функции f (n) = sin(2 * л * 800 * Т * п), Т = 1 / Ғ г и находим смещения отфильтрованной последовательности по отношению к оригиналу. Оказалось, что эти значения равны 5 и 15. Переходим к процедуре внедрения ЦВЗ. Скрипт 19. Внедрение ЦВЗ в контейнер de f i n s B i t ( B i t ,C o e l ,C o e 2 ,O r i g , F i t , I n t e r v ) : У У У Для согласования на стыках надо знать последние значения в оригинальном интервале (Orig) и отфильтрованном ( F i t ) для функции l f i l t i c I ! 5 i f B i t == ’ 1 ’ : Сое = Coel e l s e : Сое = Coe2 A = [ 1 . , -Coe] В = [Coe, - 1 . ] Zg = s g n . l f i l t i c ( B , A , [ F i t ] , [Orig]) 0ut ,Z = s g n . l f i l t e r ( B , A , I n t e r v , z i = Zg) re turn Out ####################################### def in s e r tW t r ( I n , B i t s , Coe1, Coe2, LnFragm): LnBits = l e n ( B i t s ) LnFi le = l e n ( I n ) i f LnFi le / LnFragm < LnBits: p r i n t ’Big Lnfragm’ return 0 Beg = 0
RkJQdWJsaXNoZXIy MTExODQxMg==