Keď som videl takýto nukleárny komplex s odpaľovacími rampami, došlo mi, čo sa deje: https://hrubos.tech/blogy/content/images/20260607004444-nuke_base.jpeg Ak 360/8 je 45 stupňový uhol a sin45==cos45, simulujme toto: Majme drôtik, ktorý ohýbame priemerovaním alebo bez neho: https://hrubos.tech/simulacie/splash_pendulum2.mp4

Bez priemerovania simulácia kopania lopty do rieky ukazuje toto: Tento async: https://hrubos.tech/blogy/content/images/20260607004738-Sni%CC%81mka%20obrazovky%202026-06-07%20o%200.37.17.png ukazuje pri maximálnej energii maximálnu vlnku čľupnutia.

Teraz sa pozrime na sync: https://hrubos.tech/blogy/content/images/20260607004949-Sni%CC%81mka%20obrazovky%202026-06-07%20o%200.38.48.png

V synchrónnej polohe je pri minimárnom area vlny maximálne čľupnutie. Teda pri minimálnej energii spravím na drôtiku maximálnu škodu priemerovaným ohybom, a toto som chcel ukázať ^^^ A toto chce dosiahnuť odpaľovacie 8 uholníkové zariadenie ;)

Poznámka: Pozrite si moje merania na vojne pri 12 minútovke, priemeroval som rýchlosť na kruhovej antuke pri behaní :)

Edit2: Bola tu otázka: Prečo sú sync proti sebe? Treba si uvedomiť, že nielen špeciálny prípad sin45=cos45 platí, ale zároveň pre každý prípad platí: sin(x)=cos(x-pi/2)

Edit3: Prečo je f1=2f2. No ohýbaš drôtik, ohneš ho o 2 celé, to je koľko polovičných? No 4 polovičné. Teraz kopeš lopty do rieky, no kopneš 2 alebo 4 za 1 sekundu, toť celô. A teraz chceš vidieť čo urobí kyvadlo pre amplitúdy, tj max hodnoty a plochu v tom momente pod grafom, to je energia. Toť vše.

A teraz čerešnička na torte: Čo spravili Američani: Hľadali prvočíslo v uhle: Prečo? Aby vyžaroval menej no nie?: Pentagon:

Ateraz uvažujme ďalej: Platí podobná simulácia pre vyžarovanie aj pre násobky frekvencií/uhlov? Ukážem len toto:

Popísal som toto, čo si prajem nakresliť do GPT a on vypľul tento kód: Díky parťák:

using GLMakie
using Statistics

# -----------------------
# PARAMETRE
# -----------------------
const FPS = 30
const DURATION = 30
const NFRAMES = FPS * DURATION

const WIDTH = 12.0
const NX = 500
x = range(0, WIDTH, length=NX)

const scale = 2.5   # 🔥 zosilnenie vizuálu aj kyvadiel

# -----------------------
# MODEL VLN
# -----------------------
function splash_profile(x, t, f)

    x0 = 3.0
    impact_time = 1.0

    if t < impact_time
        return zeros(length(x))
    end

    τ = t - impact_time

    envelope = exp.(-0.8 .* abs.(x .- x0))
    wave = sin.(2π * f * τ .- 3 .* abs.(x .- x0))

    splash =
        0.65 *
        exp(-1.6 * τ) *
        exp.(-((x .- x0).^2) / 0.03)

    return envelope .* wave .* 0.12 .+ splash
end

function penetration_depth(t)
    t < 1.0 && return 0.0
    τ = t - 1.0
    return min(0.55, 0.55 * (1 - exp(-4τ)))
end

# -----------------------
# FIGURE LAYOUT (2x2 + text)
# -----------------------
fig = Figure(resolution=(1600, 900))
gl = fig[1, 1] = GridLayout()

axL = Axis(gl[1, 1], title="Splash 2 Hz")
axR = Axis(gl[1, 2], title="Splash 4 Hz")

axPL = Axis(gl[2, 1], title="Pendulum 2 Hz")
axPR = Axis(gl[2, 2], title="Pendulum 4 Hz")

ylims!(axL, -1.5, 1.5)
ylims!(axR, -1.5, 1.5)

hidespines!(axPL)
hidespines!(axPR)
hidedecorations!(axPL)
hidedecorations!(axPR)

# -----------------------
# OBSERVABLES
# -----------------------
waterL = Observable(zeros(NX))
waterR = Observable(zeros(NX))

lines!(axL, x, waterL)
lines!(axR, x, waterR)

info = Observable("")
Label(gl[3, 1:2], info)

# -----------------------
# KYVADLÁ
# -----------------------
L = 0.9

θL = Observable(0.0)
θR = Observable(0.0)

rodL = Observable([Point2f(0,0), Point2f(0,0)])
rodR = Observable([Point2f(0,0), Point2f(0,0)])

massL = Observable(Point2f(0,0))
massR = Observable(Point2f(0,0))

lines!(axPL, rodL, linewidth=3)
lines!(axPR, rodR, linewidth=3)

scatter!(axPL, massL, markersize=22, color=:red)
scatter!(axPR, massR, markersize=22, color=:red)

# -----------------------
# MAX + AREA
# -----------------------
maxL = Ref(0.0)
maxR = Ref(0.0)

# -----------------------
# ANIMÁCIA
# -----------------------
record(fig, "splash_pendulum2.mp4", 1:NFRAMES; framerate=FPS) do frame

    t = (frame - 1) / FPS

    # -----------------------
    # SPLASH (zosilnený)
    # -----------------------
    zL = scale .* splash_profile(x, t, 2.0)
    zR = scale .* splash_profile(x, t, 4.0)

    waterL[] = zL
    waterR[] = zR

    sL = maximum(zL)
    sR = maximum(zR)

    maxL[] = max(maxL[], sL)
    maxR[] = max(maxR[], sR)

    # -----------------------
    # AREA (integrál nad nulou)
    # -----------------------
    dx = x[2] - x[1]
    areaL = sum(max.(zL, 0.0)) * dx
    areaR = sum(max.(zR, 0.0)) * dx

    # -----------------------
    # KYVADLÁ (silné zosilnenie)
    # -----------------------
    ampL = scale * 3.0 * sL
    ampR = scale * 3.0 * sR

    θL[] = ampL * sin(2π * 2.0 * t)
    θR[] = ampR * sin(2π * 4.0 * t)

    x1 = L * sin(θL[])
    y1 = -L * cos(θL[])

    x2 = L * sin(θR[])
    y2 = -L * cos(θR[])

    rodL[] = [Point2f(0, 0), Point2f(x1, y1)]
    rodR[] = [Point2f(0, 0), Point2f(x2, y2)]

    massL[] = Point2f(x1, y1)
    massR[] = Point2f(x2, y2)

    # -----------------------
    # REALTIME INFO
    # -----------------------
    info[] =
        "2Hz splash = $(round(maxL[]*100,digits=1)) cm | " *
        "4Hz splash = $(round(maxR[]*100,digits=1)) cm\n" *
        "Area 2Hz = $(round(areaL*100,digits=2)) | " *
        "Area 4Hz = $(round(areaR*100,digits=2))"
end

# -----------------------
# FINÁLNY PRINT
# -----------------------
println("===================================")
println("SIMULATION COMPLETE")
println("===================================")

println("2 Hz max splash = ", round(maxL[] * 100, digits=2), " cm")
println("4 Hz max splash = ", round(maxR[] * 100, digits=2), " cm")

---