а дискретизації вихідного сигналу, і спектр отриманого сигналу масштабований по осі абсцис щодо спектра вихідного сигналу. p align="justify"> Якщо вихідний сигнал містить частоти, що перевищують частоту Найквіста децімірованного сигналу, то при децимації буде мати місце алиасинг (накладення спектрів).
Таким чином, для збереження спектру необхідно до децимації видалити з вихідного сигналу частоти, що перевищують частоту Найквіста децімірованного сигналу. Ця операція проводиться цифровими фільтрами. p align="justify"> Інтерполяція це процес зворотний децимації (збільшення частоти дискретизації)
Дано
№ за списком в журналі 17.
Fd = 8400 Гц, f1 = 2000 Гц, f2 = 3800, l = 4, m = 5
F 1 , F 2 < span align = "justify"> - несучі частоти сігнла,
l - коеффіціентінтероляціі,
m-коефіцієнт децимації,
F d - частота діскрітізаціі.
В
Малюнок 4.1 - Структурна схема проектованого фільтру.
.1 Задамо сигнал у вигляді суміші гаромніческіх і шуму:
clear all = 8400 = 2000 = 3800 = 1 = 1/Fs = 0: td: t1 = sin (2 * pi * f1 * t) + cos (2 * pi * f2 * t) +0.1 * randn (size (t))
В
Малюнок 4.2 - Тимчасова діаграма нашої суміші
В
Малюнок 4.3 - Частотний спектр суміші, за допомогою програми sptool
4.2 Проводимо інтерполяцію: збільшуємо частоту дискретизації в 3 рази командою
= upfirdn (x, 1,4,1)
В
Малюнок 4.4 - Частотний спектр, після інтерполяції
.3 Проектуємо фільтр і готуємо сигнал для децимації
Вибираємо частоту зрізу виходячи зі співвідношення: fd * 3/16 => 2662 Гц
На основі цього проектуємо фільтр із заданою частотою дискретизації і частотою зрізу за допомогою команди fdatool
В
Малюнок 4.5 - Амплітудно-частотна характеристика, необхідного ФНЧ
В
Малюнок 4.6 - Імпульсна характеристика нашого ФНЧ
.4 Деціміруем сигнал і використовуємо розроблений фільтр, для виділення першої гаромонікі
Зменшення частоти дискретизації у 8 разів), за допомогою команди: w1 = upfirdn (w, Num, 1,5), де Num - коефіцієнти іспульсной характеристики нашого ФНЧ.
В
Малюнок 4.7 - Тимчасова діаграма сигналу і його спектр, після децимації та фільтрації
На спектрі видно, що тепер наша частота f1 = 2000 Гц, передається вж...
