function ssp_motor_and_load()
    # 护照：P2轴动力。
    P2 = 15_000.0; Uline = 380.0; f = 50.0; nN = 2930.0
    eta = 0.894; cosphi = 0.89; IN = 28.6
    ki = 7.5; mP = 2.2; mM = 2.3; polepairs = 1
    # 相电压、滑差和力矩。
    Uphase = Uline / sqrt(3)
    n0 = 60f / polepairs
    sN = (n0 - nN) / n0
    sK = sN * (mM + sqrt(mM^2 - 1))
    omega0 = 2pi*f / polepairs
    omegaN = pi*nN / 30
    MN = P2 / omegaN
    MP = mP*MN
    MM = mM*MN
    Pmechloss = 0.05P2
    # 转子的阻力，考虑到机械损耗。
    Rr = (P2 + Pmechloss)*sN / (3IN^2*(1-sN))
    C = 1.02
    Rs = Lsigma = Lm = 0.0
    # 替换方案中c系数的五个改进。
    for _ in 1:5
        Rs = Uphase*cosphi*(1-eta)/IN - C^2*Rr - Pmechloss/(3IN^2)
        Lsigma = Uphase/(4pi*f*(1+C^2)*ki*IN)
        L1 = Uphase/(2pi*f*IN*sqrt(1-cosphi^2) - (2/3)*(2pi*f*MM*sN)/(polepairs*Uphase*sK))
        Lm = L1 - Lsigma
        C = 1 + Lsigma/Lm
    end
    # 负载的全电功率和惯量。
    S = P2/(eta*cosphi)
    Jrotor = 0.05
    Jfan_assumed = 0.15
    J = Jrotor + Jfan_assumed
    Kfan = MN/(omegaN/omega0)^2
    H = J*omega0^2/(2S)
    (; S, Uphase, n0, sN, sK, MN, MP, MM, Rs, Rr,
       Xls=2pi*f*Lsigma, Xlr=2pi*f*Lsigma, Xm=2pi*f*Lm,
       C, J, Kfan, H)
end
ssp_motor_and_load()

# 指定脚本和现有模型所在的目录。
ssp_project_dir = @__DIR__
model_name = "softstarter"
ssp_model_path = joinpath(ssp_project_dir,"$model_name.engee")
if !any(m -> string(m.name) == model_name, engee.get_all_models())
    engee.load(ssp_model_path)
end
# 计算模型会创建一个新的simout；不保存模型文件。
engee.run(model_name)
# 我们检查图表和旁路的所有信号都被记录下来.
ssp_required = ["$model_name/speed","$model_name/current",
    "$model_name/Set-Reset Flip-Flop.Q","$model_name/firing_angle",
    "$model_name/power","$model_name/pf"]
ssp_missing = setdiff(ssp_required,collect(keys(simout)))
isempty(ssp_missing) || error("Нет сигналов в simout: " * join(ssp_missing,", "))
Dict(k => size(collect(simout[k]),1) for k in ssp_required)

using Statistics, Plots
# 我们仅从上次simout运行的结果读取所有时间序列。
ssp_speed_data = collect(simout["$model_name/speed"])
ssp_current_data = collect(simout["$model_name/current"])
ssp_bypass_data = collect(simout["$model_name/Set-Reset Flip-Flop.Q"])
ssp_power_data = collect(simout["$model_name/power"])
ssp_pf_data = collect(simout["$model_name/pf"])
ssp_t = Float64.(ssp_speed_data.time)
ssp_n = Float64.(ssp_speed_data.value)
ssp_i = Float64.(ssp_current_data.value)
ssp_b = Bool.(ssp_bypass_data.value)
# 前1→0表示旁路接触器闭合。
ssp_edge = findfirst(ssp_b[1:end-1] .& .!ssp_b[2:end])
ssp_edge === nothing && error("Не найдено замыкание bypass")
ssp_tb = ssp_t[ssp_edge+1]
ssp_before = (ssp_t .>= ssp_tb-0.2) .& (ssp_t .< ssp_tb)
ssp_after = (ssp_t .>= ssp_tb) .& (ssp_t .< ssp_tb+0.2)
ssp_steady = ssp_t .>= 12.0
# 我们比较了旁路前后相同的0.2s窗口。
ssp_stats = (; bypass_s=ssp_tb, speed_pu=mean(ssp_n[ssp_steady]),
    current_pu=mean(ssp_i[ssp_steady]),
    current_before_pu=mean(ssp_i[ssp_before]),
    current_after_pu=mean(ssp_i[ssp_after]),
    current_peak_before_pu=maximum(ssp_i[ssp_before]),
    current_peak_after_pu=maximum(ssp_i[ssp_after]),
    current_peak_start_pu=maximum(ssp_i),
    peak_change_percent=100*(maximum(ssp_i[ssp_after])/maximum(ssp_i[ssp_before])-1),
    mean_change_percent=100*(mean(ssp_i[ssp_after])/mean(ssp_i[ssp_before])-1),
    active_power_kw=mean(ssp_power_data.value[ssp_steady])*ssp_motor_and_load().S/1000,
    pf=mean(ssp_pf_data.value[ssp_steady]))

# 两个图表面板：速度和电流。
gr(format=:png)
ssp_p1 = plot(ssp_t,ssp_n; xlabel="时间，从", ylabel="速度，O.E.", label="转子", lw=1.5)
vline!(ssp_p1,[ssp_tb]; label="Bypass", ls=:dash)
ssp_p2 = plot(ssp_t,ssp_i; xlabel="时间，从", ylabel="电流，O.E.", label="Terminal current", lw=1.4)
vline!(ssp_p2,[ssp_tb]; label="Bypass", ls=:dash)
ssp_zoom = (ssp_t .>= ssp_tb-0.2) .& (ssp_t .<= ssp_tb+0.2)
ssp_p3 = plot(ssp_t[ssp_zoom],ssp_i[ssp_zoom];
    xlabel="时间，从", ylabel="电流，O.E.", label="±0.2s", lw=1.5)
vline!(ssp_p3,[ssp_tb]; label="Bypass", ls=:dash)
ssp_fig = plot(ssp_p1,ssp_p2, ssp_p3; layout=(3,1), size=(1000,900));

# 下一行将图形保存为png格式
savefig(ssp_fig,joinpath(ssp_project_dir,"softstarter_script_results.png"))
# ssp_fig

ssp_zoom = (ssp_t .>= ssp_tb-0.2) .& (ssp_t .<= ssp_tb+0.2)
ssp_p3 = plot(ssp_t[ssp_zoom],ssp_i[ssp_zoom];
    xlabel="时间，从", ylabel="电流，O.E.", label="±0.2s", lw=1.5)
vline!(ssp_p3,[ssp_tb]; label="Bypass", ls=:dash)