4 Функции Simulink для выполнения анализа спектра сигналов
В программе Simulink можно производить анализ спектра методом ДПФ (дискретного преобразования Фурье). При этом возможно анализировать спектры только непериодических сигналов по конечным наборам отсчетов.
Для анализа спектра периодических сигналов более целесообразно использовать функции классического языка MATLAB. Однако в этом случае моделирование спектра напоминает работу в программе Mathcad.
Работа в программе MATLAB осуществляется с помощью симулятора работы виртуального прибора Simulink. Запуск пакета Simulink можно произвести из командного окна MATLAB, нажав пиктограмму
При запуске Simulink открываются два окна: пустое окно untitled (окно для создания блок–диаграммы модели) и окно Library Simulink (библиотека) с перечнем основных разделов библиотеки.
В открывшееся окно untitled необходимо добавить блоки, моделирующие работу источников сигналов, измерительных приборов и аналоговых систем.
Для изменения параметров блоков необходимо выполнить двойной щелчок на пиктограмму блока. При этом должно открыться окно настройки параметров Block Parameters.
Создание модели фильтра
Для проведения спектрального анализа аналогового фильтра требуется составить следующую структурную схему (модель), рис. 5.
Рис. 5. Структурная схема для проведения спектрального анализа фильтра
Модель аналогового фильтра (в данной работе исследуется фильтр Баттерворта) выглядит так, как показано на рис. 6, и создается с помощью следующих блоков:
аналоговый фильтр с единичным усилением Analog Filter Design
(DSP blockset/Filtering/Filter Design/Analog Filter Design), рис. 7;
усилитель Gain (Simulink/Math/Gain), рис. 8.
Рис. 6. Модель фильтра Баттерворта в программе MATLAB
Рис. 7. Расположение блока Analog Filter Design
Рис. 8. Расположение блока усилителя Gain
Исходные данные необходимые для синтеза аналогового фильтра:
порядок фильтра N;
пульсации в полосе пропускания / заграждения (не требуется для фильтров Баттерворта);
коэффициент усиления K.
Исходные данные задаются в блоке параметров фильтра Block Parameters: Analog Filter Design (рис. 9), где:
Design method – вид аппроксимации;
Filter order – порядок фильтра.
Коэффициент усиления задается отдельно в блоке параметров усилителя Block Parameters: Gain (рис. 10).
Рис. 9. Блок параметров Analog Filter Design
Рис. 10. Окно настройки параметров блока Gain
Блок Power Spectral Density (анализатор спектра)
Для просмотра спектральной плотности сигнала используют анализаторы спектра Power Spectral Density (Simulink Extras/Additional Sinks/Power Spectral Density), рис.11.
Рис. 11. Расположение блока Power Spectral Density (а) и результаты анализа спектра (б)
В окне настройки Block Parameters: Power Spectral Density, рис. 12, задаются следующие параметры анализатора спектра:
Length of buffer – длина буфера (по умолчанию 128);
Number of points for fft – число анализируемых точек (по умолчанию 512);
Plot after how many points– количество точек, после которого производится построение графика (по умолчанию 64);
Sample time – период дискретизации.
Все параметры, кроме периода дискретизации, должны быть кратны 2 N , где N – целое число.
Рис. 12. Окно настройки параметров блока Power Spectral Density
Результаты анализа спектра, выполняемого блоком Power Spectral Density, представляются в трех окнах (рис. 11.б):
исследуемый сигнал (Time history);
амплитудный спектр (Power Spectral Density);
фазовый спектр (Power Spectral Density (phase)).
Чтобы просмотреть спектр сигнала нужно запустить работу модели. Это можно сделать, вызвав меню Simulation и нажав команду Start, либо нажав сочетание клавиш Ctrl + T или же щелкнув пиктограмму
Блок Signal Generator (генератор сигналов)
Для задания гармонического сигнала, используется блок Signal Generator
(Simulink/ Sources/Signal Generator), рис. 13.
Рис. 13. Расположение блока Signal Generator
В окне настройки блока Signal Generator задаются следующие параметры:
Wave form – форма сигнала:
sine– синусоидальный сигнал;
square – прямоугольный сигнал;
saw tooth – пилообразный сигнал;
random – случайный сигнал (шум);
Amplitude и Frequency – амплитуда и частота сигнала;
Units – единица измерения частоты (Hertz – Герцы или rad/sec – рад/сек).
Рис. 14. Окно настройки параметров блока Signal Generator
Для создания более сложных сигналов используется сумматор (Simulink/ Math Operations/Sum), рис. 15. Блок Sum выполняет суммирование входных сигналов. Пример использования сумматора для сложения двух входных сигналов показан на рис. 16.
Рис. 15. Расположение блока Sum
Рис. 16. Пример использования сумматора для сложения двух сигналов
Сумматор может использоваться в двух режимах:
Сложение входных сигналов (в том числе с разными знаками);
Суммирование элементов вектора, поступающих на вход блока.
Окно настройки параметров блока Sum изображено на рис. 18. Вводя значения в поле List of sings (список знаков), можно управлять режимами работы блока Sum.
Значения могут задаваться одним из трех способов:
в виде последовательности знаков «+» и «–», причем число знаков определяет число входов блока, а сами знаки – полярности соответствующих входных сигналов. При большом количестве слагаемых их целесообразно разбить на несколько групп, отделяя одну группу от другой символом │ (например: + +│– +);
в виде целой положительной константы (больше 1), значение которой равно числу входов блока, а все входы считаются положительными (например, ввод константы 4 аналогичен вводу «списка знаков» в форме + + + +).
ввод значении 1 означает вычисление суммы элементов входного вектора (в этом случае внутри блока выводится символ ?).
Рис. 17. Окно настройки параметров блока Sum
Два других параметра настройки блока имеют следующий смысл:
раскрывающийся список Icon shape (форма значка) позволяет выбрать форму блока: round (окружность) или rectangular (прямоугольник);
флажок Show additional parameters (показать дополнительные параметры) дает возможность задать необходимость «усечения» результата сложения, если он превышает диапазон, установленный для целочисленных значений.
Режимы моделирования систем в программе Simulink (меню Simulation)
Посредством этих команд входящих в меню Simulation (моделирование), рис. 18, разработчик получает возможность не только динамически управлять сеансом моделирования, но и изменять многие важнейшие параметры модели, такие, например, как способ изменения модельного времени, алгоритм расчета и формат представления результатов моделирования.
Рис. 18. Меню Simulation
Элементы управления на вкладке Solver окна Simulation Parameters (параметры моделирования), рис. 19, собраны в три группы.
Simulation time (интервал моделирования) – выбор интервала моделирования посредством указания начального (Start time) и конечного (Stop time) значений модельного времени.
Рис. 19. Окно установки параметров моделирования
Solver options (параметры расчета) – выбор метода реализации (расчета) модели.
Output options (параметры вывода) – параметры вывода выходных параметров моделируемой системы (при моделировании с переменным шагом).
Под выбором метода реализации модели имеется в виду следующее. Имея структуру исследуемой системы в виде блок–диаграммы, разработчик может выбрать метод отображения хода моделирования. С помощью двух раскладывающихся списков Type (тип) система может быть реализована в следующих формах:
с дискретными состояниями и дискретным временем перехода из одного состояния в другое;
с дискретными состояниями и непрерывным временем перехода;
с непрерывными состояниями и дискретным временем переходов;
с непрерывными состояниями и непрерывным временем переходов.
Первый список (слева) позволяет выбрать способ изменения модельного времени:
Variable – step (переменный шаг) – моделирование с переменным шагом;
Fixed – step (фиксированный шаг) – моделирование с фиксированным шагом.
Второй список (справа) позволяет выбрать метод расчета нового состояния системы. Первый вариант (discrete) обеспечивает расчет дискретных состояний системы. Остальные пункты списка обеспечивают выбор метода расчета нового состояния для непрерывных систем. Эти методы различаются для переменного (Variable – step) и для фиксированного (Fixed – step) шага времени, но основаны на единой методике – решение обыкновенных дифференциальных уравнений(ode).
Ниже двух раскрывающихся списков Type находится поле, название которого изменяется в зависимости от выбранного способа изменения модельного времени.
Параметры других вкладок используется по умолчанию.
Источник
Документация
Оцените спектр мощности в Simulink
Спектр мощности (PS) сигнала временной области является распределением степени, содержавшей в сигнале по частоте, на основе конечного множества данных. Представление частотного диапазона сигнала часто легче анализировать, чем представление временного интервала. Много приложений обработки сигналов, таких как подавление помех и система идентификации, основаны на специфичных для частоты модификациях сигналов. Цель степени спектральная оценка состоит в том, чтобы оценить спектр мощности сигнала от последовательности выборок времени. В зависимости от того, что известно о сигнале, методы оценки могут включить параметрические или непараметрические подходы и могут быть основаны на временном интервале или анализе частотного диапазона. Например, общий параметрический метод включает подбор кривой наблюдениям к авторегрессивной модели. Общий непараметрический метод является периодограммой. Спектр мощности оценивается с помощью методов преобразования Фурье, таких как валлийский метод и метод набора фильтров. Для сигналов с относительно маленькой длиной подход набора фильтров производит спектральную оценку с более высоким разрешением, более точным уровнем шума и peaks, более точным, чем валлийский метод, с низким или никакой спектральной утечкой. Эти преимущества прибывают за счет увеличенного расчета и более медленного отслеживания. Для получения дополнительной информации об этих методах смотрите Спектральный анализ. Можно также использовать другие методы, такие как максимальный энтропийный метод.
В Simulink ® , можно выполнить спектральный анализ в реальном времени динамического сигнала с помощью блока Spectrum Analyzer . Можно просмотреть спектральные данные в спектре анализатор. Чтобы получить последние спектральные данные для последующей обработки, создайте SpectrumAnalyzerConfiguration возразите и запуститесь getSpectrumData функция на этом объекте. Альтернативно, можно использовать блок Spectrum Estimator из dspspect3 библиотека, чтобы вычислить спектр мощности и блок Array Plot , чтобы просмотреть спектр.
Оцените спектр мощности Используя спектр Анализатор
Можно просмотреть спектр мощности (PS) сигнала с помощью блока Spectrum Analyzer . PS вычисляется в режиме реального времени и меняется в зависимости от входного сигнала, и с изменениями в свойствах блока Spectrum Analyzer . Можно изменить динамику входного сигнала и видеть то, что вызывает те изменения, имеют на спектре сигнала в режиме реального времени.
Модель ex_psd_sa кормит шумным синусоидальным сигналом блок Spectrum Analyzer . Синусоидальный сигнал является суммой двух синусоид: один на частоте 5 000 Гц и другом на частоте 10 000 Гц. Шум во входе является Гауссовым с нулевым средним значением и стандартным отклонением 0,01.
Откройте и смотрите модель
Чтобы открыть модель, введите ex_psd_sa в MATLAB ® командная строка.
Источник