Дискретный КИХ-фильтр
В данном примере реализованы модель фильтрации синусоиды и сравнение работы КИХ-фильтра в MATLAB и Engee. Блок Discrete FIR Filter независимо фильтрует каждый канал входного сигнала с заданным цифровым КИХ-фильтром.
Для реализации фильтрации используем стандартные настройки блока Discrete FIR Filter.
Модель Engee
Модель Simulink
using Plots
using MATLAB
using CSV
using DataFrames
using Statistics
plotlyjs();
mat"start_simulink"
mat"p = genpath('/user/start/examples'); addpath(p);"
Запуск модели Engee
Посмотрим на результаты моделирования в Engee:
# Запуск модели
modelName = "FIR"
#engee.close( modelName, force=true ) # Для повторного запуска модели
PID_model = modelName in [m.name for m in engee.get_all_models()] ? engee.open( modelName ) : engee.load( "$(@__DIR__)/$(modelName).engee");
results = engee.run( PID_model )
Dict{String, DataFrame} with 2 entries:
"FIR_Out" => 101×2 DataFrame…
"FIR_In" => 101×2 DataFrame…
# Результаты моделирования
FIR_In_x = results["FIR_In"].time;
FIR_In_y = results["FIR_In"].value;
FIR_Out_x = results["FIR_Out"].time;
FIR_Out_y = results["FIR_Out"].value;
# Построение графиков
plot(FIR_Out_x , FIR_Out_y, legend = false)
plot!(FIR_In_x , FIR_In_y, legend = false)
plot!(title = "Результаты моделирования в Engee", ylabel = "Данные", xlabel="Время, c")
Запуск модели Simulink
Сравним результаты моделирования в Engee и Simulink. Для этого сначала запустим симуляцию модели.
mat"run_test_model('FIR')";
Посмотрим на результаты в Simulink:
# Чтение выходов модели
x1 = mat"In.Data";
x2 = mat"Out.Data";
y = mat"Out.Time";
# Построение графиков
plot(y, x2, legend = false)
plot!(y, x1, legend = false)
plot!(title = "Результаты моделирования в Simulink", ylabel = "Отклик", xlabel="Время, c")
Теперь построим график наложения выходных данных Simulink и Engee.
plot(y, x2, legend = false) # Simulink
plot!(FIR_Out_x , FIR_Out_y, legend = false) # Engee
Как мы видим из графиков представленных выше, модели работают одинаково.