function run_model( name_model)
    Path = (@__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
    sleep(0.1)
    return model_output
end
run_model("PUCCH_model") 

using EngeeDSP
using Statistics, LinearAlgebra
using FFTW  
using StatsBase 

function unwrap(phase::Vector{Float64})
    unwrapped = copy(phase)
    for i in 2:length(phase)
        diff = phase[i] - phase[i-1]
        if diff > π
            unwrapped[i] -= 2π
        elseif diff < -π
            unwrapped[i] += 2π
        end
    end
    return unwrapped
end

function compute_papr(signal::AbstractArray{<:Complex})
    power = abs2.(signal)
    papr_db = 10 * log10(maximum(power) / mean(power))
    return round(papr_db, digits=2)
end

function compute_spectrum(signal::AbstractArray{<:Complex}, fs=1.0)
    N = length(signal)
    fft_result = fft(signal)
    freqs = fftfreq(N, fs)
    pwr = abs2.(fft_result) / N
    return freqs, pwr
end

function analyze_constellation(signal::AbstractArray{<:Complex}, title)
    p = scatter(real.(signal), imag.(signal),
                title=title, xlabel="In-phase (I)", ylabel="Quadrature (Q)",
                legend=false, markersize=4, markerstrokewidth=0, alpha=0.7,
                aspect_ratio=:equal)
    return p
end

function compute_signal_stats(signal::AbstractArray{<:Complex}, name)
    println("\n" * "═"^50)
    println("统计数字：$名称")
    println("═"^50)
    println("   数据类型:$(eltype(signal))")
    println("   尺寸:$(尺寸(信号))")
    println("   PAPR:$(compute_papr(signal))dB")
    println("   能量：△(round(sum(abs2.（信号）），数字=2））")
    println("   平均振幅：△(round(mean(abs.（信号）），数字=4））")
    println("   平均相位：△(round(mean(angle.（信号）），数字=4））高兴")
end

y = collect(simout["PUCCH_model/PUCCH1/# Parse and validate inputs.y"]).value[end]
oSeq1 = collect(simout["PUCCH_model/PUCCH1/Spreading Format-1.1"]).value[end]
final = collect(simout["PUCCH_model/PUCCH1/# Get the PUCCH format 1 sequence.Seq"]).value[end]

compute_signal_stats(y, "信号y(nrpucch1后)")
compute_signal_stats(oSeq1, "正交序列1")
compute_signal_stats(final, "最后的序列")

p1 = analyze_constellation(y, "PUCCH1：'y'信号的星座")
p2 = analyze_constellation(final, "Pucch1：最后的星座")
p3 = plot(abs.(final), title="Pucch1：振幅包络",
            xlabel="计数", ylabel="的振幅", legend=false, grid=true)
display(plot(p1, p2, p3, layout=(1,3), size=(1200, 400)))
savefig("pucch1_analysis.png")

prbs_bits = collect(simout["PUCCH_model/PUCCH2/PUCCH PRBS-1.1"]).value[end]
modulated = collect(simout["PUCCH_model/PUCCH2/Symbol Modulate-1.1"]).value[end]
println("加扰比特(PRB):")
println("   位数：$(length(prbs_bits))")
println("   余额（1/0）：$（sum（prbs_bits））/$（sum（。!prbs_bits))")

compute_signal_stats(modulated, "QPSK调制字符")
unique_symbols = unique(round.(modulated, digits=6))
println("\Qpsk星座的Punic点:")
for sym in unique_symbols
    println("   I=$(round(real(sym), digits=4)), Q=$(round(imag(sym), digits=4))")
end

phase = angle.(modulated)
instant_freq = diff(unwrap(phase)) 
p1 = histogram(instant_freq, bins=30, title="Pucch2:瞬时频率",
                xlabel="频率", ylabel="发生频率", legend=false)
p2 = analyze_constellation(modulated, "PUCCH2：QPSK星座")

display(plot(p1, p2, layout=(1,2), size=(1200, 400)))
savefig("pucch2_analysis.png")

converted = collect(simout["PUCCH_model/PUCCH3/Convert in Int.1"]).value[end]
final = collect(simout["PUCCH_model/PUCCH3/Transform Precoding-1.1"]).value[end]
println("\转换的整数数据:")
println("   范围:$(最小(转换))-$(最大(转换))")
println("   平均值：$（圆形（平均值（转换），数字=2））")
compute_signal_stats(final, "变换预编码后的信号")

N = length(final)
fs = 1.0 
fft_result = fft(final)
freqs = fftfreq(N, fs)
pwr = abs2.(fft_result) / N
p1 = plot(freqs, pwr, title="Pucch3:功率谱密度",
            xlabel="频率", ylabel="功率(dB)", legend=false, grid=true,
            xlim=(-0.5, 0.5))
p2 = analyze_constellation(final, "Pucch3：DFT后的星座")

max_pwr = maximum(pwr)
half_pwr = max_pwr / 2 
indices = findall(pwr .> half_pwr)
if !isempty(indices)
    min_freq = minimum(freqs[indices])
    max_freq = maximum(freqs[indices])
    bandwidth = max_freq - min_freq
    println("\车道宽度标记:")
    println("   带宽在3分贝水平:$(轮(带宽,数字=4))")
else
    println("我们还没有能够估计3dB级别的带宽")
end
display(plot(p1, p2, layout=(1,2), size=(900, 400)))
savefig("pucch3_analysis.png")

prbs = collect(simout["PUCCH_model/PUCCH4/PUCCH PRBS.1"]).value[end]
modulated = collect(simout["PUCCH_model/PUCCH4/Symbol Modulate.1"]).value[end]
after_spreading = collect(simout["PUCCH_model/PUCCH4/Spreading.Out1"]).value[end]
final = collect(simout["PUCCH_model/PUCCH4/Transform Precoding.1"]).value[end]
compute_signal_stats(modulated, "在按块传播之前")
compute_signal_stats(after_spreading, "传播后")
compute_signal_stats(final, "变换后预编码")
gain_factor = mean(abs2.(after_spreading)) / mean(abs2.(modulated))
println("按块传播的效率:")
println("   功率增益:$(round(gain_factor,digits=2))")

plots = []
push!(plots, analyze_constellation(modulated, "在传播之前"))
push!(plots, analyze_constellation(after_spreading, "传播后"))
push!(plots, analyze_constellation(final, "最终信号"))
p4 = plot(abs.(modulated), label="在传播之前", linewidth=1.5, alpha=0.8)
plot!(p4, abs.(after_spreading), label="传播后", linewidth=1, alpha=0.7)
plot!(p4, abs.(final), label="预编码后", linewidth=1, alpha=0.7,
        title="Pucch4：比较振幅", xlabel="计数", ylabel="的振幅", grid=true)
push!(plots, p4)
display(plot(plots..., layout=(2,2), size=(1000, 800)))
savefig("pucch4_analysis.png")