Рис. 2.4. Передаточная функция фильтра, построенного с помощью firwin График передаточной функции фильтра представлен на рис. 2.4. Заметим, что в данном случае мы легко подавили сигнал с помощью фильтра высоких частот. Это связано с тем, что частота стробирования в данном примере 22850 достаточно велика. В противном случае для подавления основной мощности сигнала может понадобиться использовать функцию firwin2, позволяющую построить фильтр, подавляющий сразу несколько интервалов частот. Перейдем к внедрению водяного знака. Он равен коэффициентам фильтра, записанным в обратном порядке. В данном случае, построенный фильтра является симметричным, но в остальных случаях это сделать необходимо. После этого проверим корреляцию различных участков файла с найденным шаблоном. Wtr = ВС: : - 1 ] # Запись элементов # массива в обратном порядке DatN = n p . f l o a t 3 2 (D a t ) LnWtr = len(Wtr) Coef = 21 DatN[Pos:Pos+LnWtr] += Coef * Wtr Corr = np . co r r e la t e (Da tN [Po s - 40 : Pos+LnWtr+40], W t r , ’v a l i d ’ ) p r i n t sum(Wtr*Wtr) p i t . p l o t ( r a n g e ( 8 0 ) ,C o r r [ : 8 0 ] , ’k ’ )
RkJQdWJsaXNoZXIy MTExODQxMg==