In [ ]:
Pkg.add(["DSP"])
In [ ]:
using Plots, DSP; # Подключение библиотек PLots и DSP
现在,让我们使用 DSP 库为滤波器生成系数。
In [ ]:
st = 0.0001;
fc = [500 1200]'; # Частоты сигнала
fs = 6000; # Частота дискретизации сигнала
t = [0:1/fs:0.001;]; # Диапазон времени сигнала
x = cos.(2*pi*fc[1]*t) + cos.(2*pi*fc[2]*t); # Пример сигнала для фильтрации
responsetype = Lowpass(2000; fs); # Определение полосы пропускания
print("Размер окна фильтрации:")
N = parse(Int,readline())
designmethod = FIRWindow(hanning(N)); # Определение размеров окна
x_filt =filt(digitalfilter(responsetype, designmethod), x); # Фильтрация сигнала
c=digitalfilter(responsetype, designmethod); # FIR-фильтр
coef=length(c); # Количество коэффициентов фильтра
Out[0]:
让我们绘制原始信号和滤波后信号的功率谱密度图。
In [ ]:
p1 = DSP.periodogram(x);
plot(freq(p1), power(p1), xlabel="частота, Гц", ylabel="спектральная плотность мощности", label="исходный сигнал")
p2 = DSP.periodogram(x_filt);
plot!(freq(p2), power(p2), label="отфильтрованный сигнал")
Out[0]:
现在,让我们创建一个模型,在其中生成滤波器本身。
In [ ]:
print("Напишите название вашей модели:")
name_model = readline()
Path = (@__DIR__) * "/" * name_model * ".engee"
if isdir(Path)
rm(Path;force = true, recursive = true)
end
engee.create(name_model) # Создать модель
Out[0]:
首先,让我们设置模型的输入和输出端口。
In [ ]:
engee.add_block("/Basic/Ports & Subsystems/In1", name_model*"/"); # Создание входного порта для подсистемы
engee.add_block("/Basic/Ports & Subsystems/Out1", name_model*"/"); # Создание выходного порта для подсистемы
现在让我们声明一个创建 FIR 滤波器模型的循环。
In [ ]:
for n in 1:coef-1
name_gain="Gain-"*string(n); # Задаем имя блока для Gain
engee.add_block("/Basic/Math Operations/Gain", name_model*"/"*name_gain); # Добавим в модель Gain
engee.set_param!(name_model*"/"*name_gain, "Gain" => c[n]); # Установим значения коэффициентов фильтра
name_delay="Delay-"*string(n); # Задаем имя блока для Delay
engee.add_block("/Basic/Discrete/Delay", name_model*"/"*name_delay); # Добавим в модель Delay
engee.set_param!(name_model*"/"*name_delay, "DelayLength" => 1); # Установим длину задержки равной 1
engee.set_param!(name_model*"/"*name_delay, "SampleTime" => st); # SampleTime для Delay
name_add="Add-"*string(n); # Задаем имя блока для Add
engee.add_block("/Basic/Math Operations/Add", name_model*"/"*name_add); # Добавим в модель Add
if n==1
engee.add_line(name_gain*"/1", name_add*"/1"); # Соединим Gain и Add 1 вход
end
if n>1
name_delay_1="Delay-"*string(n-1); # Задаем имя блока для предыдущего Delay
engee.add_line(name_delay_1*"/1", name_delay*"/1"); # Соединим Delay n-1 и Delay n
engee.add_line(name_delay_1*"/1", name_gain*"/1"); # Соединим Delay n-1 и Gain n
name_add_1="Add-"*string(n-1); # Задаем имя блока для предыдущего Add
engee.add_line(name_add_1*"/1", name_add*"/1"); # Соединим Add n-1 и Add 1 вход
engee.add_line(name_gain*"/1", name_add_1*"/2"); # Соединим Gain n-1 и Add 2 вход
end
if n==coef-1
name_gain="Gain-"*string(n+1); # Задаем имя блока для Gain
engee.add_block("/Basic/Math Operations/Gain", name_model*"/"*name_gain); # Добавим в модель Gain
engee.set_param!(name_model*"/"*name_gain, "Gain" => c[n+1]); # Установим значения коэффициентов фильтра
engee.add_line(name_delay*"/1", name_gain*"/1"); # Соединим Delay и Gain
engee.add_line(name_gain*"/1", name_add*"/2"); # Соединим Gain и Add 2 вход
engee.add_line(name_add*"/1", "Out1/1"); # Соединим Add и Out1
end
end
engee.add_line("In1/1", "Gain-1/1"); # Соединим In1 и Gain-1
engee.add_line("In1/1", "Delay-1/1"); # Соединим In1 и Delay-1
将结果保存到模型中并更改建模参数。
In [ ]:
engee.save(Path)
model = engee.load(Path, force=true ) # Загрузить модель
engee.set_param!(model, "StopTime" => 5, "FixedStep" => 0.1) # меняем фиксированный размер шага и время окончания симуляции
param = engee.get_param(model) # Получение параметров текущей модели
结论¶
在本例中,我们展示了一个使用命令控制自动生成模型的例子。