육군 비호복합 30mm 쌍열 대공포 탄막 형성(Barrage Firing) 분산 타원체 & 신궁 적외선 유도탄 듀얼 모드 표적 할당 1,000Hz 0-Alloc C# 사격통제 엔진

using System;
using System.Runtime.CompilerServices;
using System.Runtime.InteropServices;
using System.Runtime.Intrinsics;
using System.Runtime.Intrinsics.X86;

namespace ZeliDesk.Defense.AntiAirSystems
{
    public enum WeaponChoice : byte
    {
        None = 0,
        Twin30mmGun = 1,      // 30mm 쌍열 기관포 탄막 사격
        ChironMissile = 2,    // 신궁 적외선 유도 미사일 발사
        SimultaneousDual = 3  // 동시 복합 교전
    }

    [StructLayout(LayoutKind.Sequential, Pack = 16)]
    public struct FireControlSolution
    {
        public Vector256<float> PredictedInterceptPos; // [x, y, z, t_go, 0, 0, 0, 0]
        public float TargetDistanceM;                  // 표적 거리 (m)
        public float GunKillProbability;               // 30mm 탄막 격추 확률
        public float MissileLockProbability;           // 신궁 IR 락온 확률
        public WeaponChoice RecommendedWeapon;         // 최적 무장 선택
    }

    public sealed class BihoHybridFireControlEngine
    {
        private const float MuzzleVelocity = 1175.0f; // 30mm 포구초속 (m/s)
        private const float GunMaxEffectiveRange = 3000.0f; // 대공포 최대 유효사거리 (m)
        private const float MissileMaxRange = 6000.0f;      // 신궁 최대 사거리 (m)
        private const float MissileMinRange = 800.0f;       // 신궁 최소 사거리 (m)

        private const float AngularDispersionRad = 0.0025f; // 2.5 mrad 분산각
        private const float LethalRadiusM = 3.5f;           // 30mm 파편 치사 반경 (m)

        [MethodImpl(MethodImplOptions.AggressiveOptimization)]
        public FireControlSolution ComputeSolution(
            float targetX, float targetY, float targetZ,
            float targetVx, float targetVy, float targetVz,
            float targetIrIntensity, int gunBurstCount)
        {
            FireControlSolution solution = default;

            float dist = MathF.Sqrt(targetX * targetX + targetY * targetY + targetZ * targetZ);
            solution.TargetDistanceM = dist;

            // 1. 탄환 도달 시간 및 예측 조우점(PIP) 계산 (1차 반복 근사)
            float tGo = dist / MuzzleVelocity;
            float pipX = targetX + targetVx * tGo;
            float pipY = targetY + targetVy * tGo;
            float pipZ = targetZ + targetVz * tGo - 0.5f * 9.80665f * tGo * tGo;

            solution.PredictedInterceptPos = Vector256.Create(pipX, pipY, pipZ, tGo, 0, 0, 0, 0);

            // 2. 30mm 탄막 분산 타원체 반경 및 격추 확률 계산
            float sigmaCross = dist * AngularDispersionRad;
            float sigmaAlong = 12.0f; // 종방향 분산
            float lethalArea = MathF.PI * LethalRadiusM * LethalRadiusM;
            float ellipsoidVol = (4.0f / 3.0f) * MathF.PI * sigmaCross * sigmaCross * sigmaAlong;

            float singleShotHitProb = MathF.Min(0.45f, lethalArea / (MathF.PI * sigmaCross * sigmaCross + 1e-3f));
            // N발 점사 시 누적 격추 확률
            float cumulativeSurvival = MathF.Pow(1.0f - singleShotHitProb * 0.85f, (float)gunBurstCount);
            solution.GunKillProbability = (dist <= GunMaxEffectiveRange) ? (1.0f - cumulativeSurvival) : 0.0f;

            // 3. 신궁 미사일 적외선 탐색기 락온 확률 계산
            float snr = (targetIrIntensity * 10000.0f) / (dist * dist + 1.0f);
            float pLock = 1.0f / (1.0f + MathF.Exp(-1.2f * (snr - 2.5f)));
            solution.MissileLockProbability = (dist >= MissileMinRange && dist <= MissileMaxRange) ? pLock : 0.0f;

            // 4. TEWA 무장 할당 결정론적 판정
            if (dist > GunMaxEffectiveRange && dist <= MissileMaxRange)
            {
                solution.RecommendedWeapon = (solution.MissileLockProbability > 0.65f) 
                    ? WeaponChoice.ChironMissile : WeaponChoice.None;
            }
            else if (dist <= GunMaxEffectiveRange && dist >= 1200.0f)
            {
                if (solution.MissileLockProbability > 0.85f && solution.GunKillProbability < 0.75f)
                    solution.RecommendedWeapon = WeaponChoice.ChironMissile;
                else if (solution.GunKillProbability >= 0.80f)
                    solution.RecommendedWeapon = WeaponChoice.Twin30mmGun;
                else
                    solution.RecommendedWeapon = WeaponChoice.SimultaneousDual;
            }
            else if (dist < 1200.0f)
            {
                // 근거리 긴급 표적은 무조건 30mm 기관포 탄막 요격
                solution.RecommendedWeapon = WeaponChoice.Twin30mmGun;
            }

            return solution;
        }
    }
}