Моделирование узкополосного низкочастотного фильтра

Описанная процедура может быть легко реализована в виде модели Simulink, для чего нам потребуются следующие блоки:

генератор синусоидальных колебаний Sine Wave (папка Simulink\ Sources);

генератор шума Random Number (Simulink\Sources);

сумматор Sum (Simulink\Math);

полифазный фильтр-дециматор FIR Decimation (DSP Blockset\Filtering\ Multirate Filters);

полифазный фильтр-интерполятор FIR Interpolation (DSP Blockset\ Filtering\Multirate Filters);

анализатор спектра Buffered FFT Frame Scope (DSP Blockset\DSP Sinks);

осциллограф Scope (Simulink\Sinks).

Правила построения моделей подробно изложены в справочной системе MATLAB, с ними также можно ознакомиться, обратившись к соответствующим источникам [6-8].
Построенная модель, где блоки (F+D)1 и (F+D)2 - это фильтры-дециматоры, а блоки (I+F)1 и (I+F)2 - фильтры-интерполяторы. Перед тем, как её запустить, необходимо задать параметры блоков и режима работы модели. В нашем примере параметры блоков имеют следующие значения:

Блок (F+D)1. FIR filter coefficients = b1, Decimation factor = 20;

Блок (F+D)2. FIR filter coefficients = b2, Decimation factor = 2;

Блок (I+F)2. FIR filter coefficients = b2, Interpolation factor = 2;

Блок (I+F)1. FIR filter coefficients = b1, Interpolation factor = 20;

Анализаторы спектра 1 и 5. Параметры блоков одинаковы и указаны на рис. 13;

Анализаторы спектра 2 и 4. FFT length и Buffer size = 256, Sample time of original time series = 20/8000, остальные - как на рис. 13;

Анализатор спектра 3. FFT length = = 256 и Buffer size = 128, Sample time of original time series = 20·2/8000, остальные - как на рис. 13;

Генератор 1. Amplitude = 1, Frequency = 90, Sample Mode = Discrete, Sample Time = 1/8000;

Генератор 2. Amplitude = 1, Frequency = 105, Sample Mode = Discrete, Sample Time = 1/8000;

Генератор шума. Mean = 0, Variance = 3e-1, Sample Time = 1/8000.

Параметры режима работы модели устанавливаются в окне Simulation Parameters.