function run_model(name_model)
    Path = string(@__DIR__) * "/" * name_model * ".engee"
    
    if name_model in [m.name for m in engee.get_all_models()] # Проверка условия загрузки модели в ядро
        model = engee.open( name_model ) # Открыть модель
        model_output = engee.run( model, verbose=true ); # Запустить модель
    else
        model = engee.load( Path, force=true ) # Загрузить модель
        model_output = engee.run( model, verbose=true ); # Запустить модель
        engee.close( name_model, force=true ); # Закрыть модель
    end
    return model_output
end

Eb_No = 13;
Select_index = 1;

run_model("QPSK_timing_sync"); # Запуск модели.
println(simout)
println()
println("BER: $(collect(BER_QPSK).value[end][1])")

# Подключение библиотек
neededLibs = ["FFTW", "DSP", "Statistics"]
for lib in neededLibs
    try
        eval(Meta.parse("using $lib"))
    catch ex
        Pkg.add(lib)
        eval(Meta.parse("using $lib"))
    end
end

Out_Transmit_Filter_value = vcat(simout["QPSK_timing_sync/RX_out"].value...)
fs = 800  
power_spectrum_raw = abs.(fft(Out_Transmit_Filter_value)).^2 / (length(Out_Transmit_Filter_value) * fs)
power_spectrum_db_raw = 10*log10.(power_spectrum_raw)
window_size = 20  
kernel = ones(window_size) / window_size 
power_spectrum_smoothed = conv(power_spectrum_db_raw, kernel)[window_size÷2+1:end-window_size÷2]
freqs = fftfreq(length(power_spectrum_smoothed), fs)

plot(fftshift(freqs), fftshift(power_spectrum_smoothed),
     title="Энергетический спектр (скользящее среднее, окно $window_size)",
     xlabel="Частота, Гц",
     ylabel="Спектральная плотность мощности, дБ/Гц",
     legend=false,
     linewidth=2,
     grid=true)

ber_qpsk = zeros(8);
RX_out = Vector{Vector{ComplexF64}}(undef, 8)
Selection_data = Vector{Vector{ComplexF64}}(undef, 8) 
fs = 800 
global Select_index

@time for i in 1:8
    println("Запуск модели для Select_index = $i")
    Select_index = Int32(i)
    run_model("QPSK_timing_sync")
    ber_qpsk[i] = collect(BER_QPSK).value[end][1]
    RX_out[i] = vcat(simout["QPSK_timing_sync/RX_out"].value...)
    Selection_data[i] = simout["QPSK_timing_sync/Выбор.1"].value
end

println("\nРезультаты BER:")
for i in 1:8
    println("Select_index = $i: BER = $(ber_qpsk[i])")
end
plot(1:8, ber_qpsk, 
     title="Зависимость BER от Select_index",
     xlabel="Select_index", 
     ylabel="BER",
     marker=:circle,
     linewidth=2,
     grid=true,
     legend=false,
     yscale=:log10)

best_idx = argmin(ber_qpsk)
worst_idx = argmax(ber_qpsk)
best_selection = Selection_data[best_idx]
worst_selection = Selection_data[worst_idx]
downsample_factor = 20
best_constellation = best_selection[1:downsample_factor:end]
worst_constellation = worst_selection[1:downsample_factor:end]
function compute_spectrum_q(signal, fs)
    q_signal = imag(signal)
    N = length(q_signal)
    power_spectrum = abs.(fft(q_signal)).^2 / (N * fs)
    power_spectrum_db = 10*log10.(power_spectrum)
    freqs = fftfreq(N, fs)
    return fftshift(freqs), fftshift(power_spectrum_db)
end
freqs_best, spectrum_best = compute_spectrum_q(RX_out[best_idx], fs)
freqs_worst, spectrum_worst = compute_spectrum_q(RX_out[worst_idx], fs)
p_combined = plot(layout=(3,1), size=(800, 900))

t = 1:200
plot!(p_combined[1], t, imag(best_selection[t]), 
      label="Лучший (Index $best_idx)", linewidth=2, color=:blue)
plot!(p_combined[1], t, imag(worst_selection[t]), 
      label="Худший (Index $worst_idx)", linewidth=2, color=:red, linestyle=:dash)
title!(p_combined[1], "Временные данные Q компонента")

plot!(p_combined[2], freqs_best, spectrum_best, 
      label="Лучший", linewidth=2, color=:blue)
plot!(p_combined[2], freqs_worst, spectrum_worst, 
      label="Худший", linewidth=2, color=:red, linestyle=:dash)
title!(p_combined[2], "Спектры Q компонентов")

scatter!(p_combined[3], 1:length(best_constellation), imag(best_constellation), 
         label="Лучший", markersize=4, alpha=0.7, color=:blue,
         markerstrokewidth=0, marker=:circle)
scatter!(p_combined[3], 1:length(worst_constellation), imag(worst_constellation), 
         label="Худший", markersize=4, alpha=0.7, color=:red,
         markerstrokewidth=0, marker=:x)
title!(p_combined[3], "Q компоненты созвездий")
plot!(p_combined, titlefontsize=12, legendfontsize=10)