From 031a0c0fdf46d429a9e55ac1677e5932b9584b8b Mon Sep 17 00:00:00 2001 From: cw2007 Date: Sun, 5 Jul 2026 14:02:48 +0900 Subject: [PATCH] =?UTF-8?q?feat:=20=EB=8C=80=EC=8B=9C=EB=B3=B4=EB=93=9C=20?= =?UTF-8?q?=EC=98=88=EB=B3=B4=C2=B7=EB=8C=80=EC=9D=91=20=EC=A0=95=EB=B3=B4?= =?UTF-8?q?=20(=EC=98=88=EB=B0=A9=20=EC=8B=A0=ED=98=B8=20+=20=EC=A1=B0?= =?UTF-8?q?=EC=B9=98=20=EA=B0=80=EC=9D=B4=EB=93=9C)?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit 기존 측정 데이터만으로 "경보가 터지기 전"을 보여주고, 위험 시 무엇을 할지 안내한다. 새 센서·ML·DB 스키마 변경 없음(폐쇄망 정합). 신규 php/forecast.php — 순수 계산 함수(DB·출력 없음, 단위테스트 용이): - dew_point_c(Magnus 이슬점), linear_trend(최소자승 회귀), threshold_eta_min (임계 도달 예상 분), headroom_state(여유 3단계), forecast_metric/freshness, risk_items(대응 섹션용 정규화 위험 목록). tests/forecast_test.php — 31 검사(이슬점 알려진값·회귀·ETA·여유·신선도·위험 엣지). config.php — FORECAST_*(지평/최소점/근접여유/급변 기준) + RESPONSE_PLAYBOOKS (상태별 대응 체크리스트) + RESPONSE_CONTACTS. 전부 config.local.php 오버라이드. dashboard.php — 센서별 최근점으로 예보 계산 → - 카드 예보줄(상시): 임계 여유·추세 ETA·이슬점 결로·변화율(근접/초과만 강조). - 대응 필요 섹션(위험 시만): 조건·ETA + 상태별 권장 조치 체크리스트 + 담당자. 기존 offline 판정을 단일 진실원으로 사용, 누수 노드 중복 없음. app.css — .forecast-chip(near/over 상태 3중 인코딩)·.response-panel·.playbook· .contact (토큰 기반, 라이트/다크 자동). run_tests.sh 에 forecast_test 등록. 검증: forecast_test 31/0, php -l clean(7페이지+forecast), 헤드리스 렌더 확인. --- php/assets/app.css | 39 +++++ php/config.local.example.php | 11 ++ php/config.php | 31 ++++ php/dashboard.php | 131 +++++++++++++++++ php/forecast.php | 272 +++++++++++++++++++++++++++++++++++ run_tests.sh | 1 + tests/forecast_test.php | 100 +++++++++++++ 7 files changed, 585 insertions(+) create mode 100644 php/forecast.php create mode 100644 tests/forecast_test.php diff --git a/php/assets/app.css b/php/assets/app.css index ca8004f..464f188 100644 --- a/php/assets/app.css +++ b/php/assets/app.css @@ -321,6 +321,45 @@ tbody tr:hover td { background: var(--surface-2); } .refresh { float: right; font-size: .8125rem; color: var(--ink-2); } +/* ── 예보(카드 상시 예방 신호) ───────────────────────────────────────────── */ +.forecast { display: flex; flex-wrap: wrap; gap: var(--sp-2); margin-top: var(--sp-2); } +.forecast-chip { + display: inline-flex; align-items: center; gap: .35em; border-radius: var(--r-pill); + padding: 2px 9px; font-size: .75rem; font-weight: 700; line-height: 1.5; + background: var(--warn-bg); color: var(--warn-ink); +} +.forecast-chip::before { content: ""; width: .45em; height: .45em; border-radius: 50%; background: currentColor; flex: none; } +.forecast-chip.near { background: var(--warn-bg); color: var(--warn-ink); } +.forecast-chip.over { background: var(--leak-bg); color: var(--leak-ink); } +.forecast-note { font-size: .75rem; color: var(--ink-3); } + +/* ── 대응 필요 섹션(위험 시에만 노출: 조건 + 플레이북 + 담당자) ─────────────── */ +.response-panel { + background: var(--surface); border: 1px solid var(--line); border-left: 6px solid var(--leak); + border-radius: var(--r); padding: var(--sp-5); margin-bottom: var(--sp-5); +} +.response-panel h2 { margin: 0 0 var(--sp-3); display: flex; align-items: center; gap: var(--sp-2); } +.response-count { + font-size: .8125rem; font-weight: 700; color: var(--leak-ink); + background: var(--leak-bg); border-radius: var(--r-pill); padding: 2px 10px; +} +.response-item { border-top: 1px solid var(--line-2); padding: var(--sp-3) 0; } +.response-item:first-of-type { border-top: none; padding-top: 0; } +.response-head { display: flex; flex-wrap: wrap; align-items: baseline; gap: var(--sp-2); } +.response-head strong { font-size: .95rem; } +.response-detail { color: var(--ink-2); font-size: .875rem; } +.playbook { margin-top: var(--sp-2); background: var(--surface-2); border-radius: var(--r-sm); padding: var(--sp-3) var(--sp-4); } +.playbook-title { font-size: .8125rem; font-weight: 700; color: var(--ink-2); margin-bottom: 4px; } +.playbook-steps { margin: 0; padding-left: 1.4em; } +.playbook-steps li { font-size: .875rem; line-height: 1.7; } +.response-contacts { + display: flex; flex-wrap: wrap; align-items: center; gap: var(--sp-3); + margin-top: var(--sp-4); padding-top: var(--sp-3); border-top: 1px solid var(--line-2); +} +.contacts-label { font-size: .8125rem; font-weight: 700; color: var(--ink-2); } +.contact { font-size: .875rem; } +.contact a { font-weight: 700; } + /* ── 폼 / 로그인 / 카드형 페이지 ────────────────────────────────────────── */ .login-box { background: var(--surface); border: 1px solid var(--line); border-radius: var(--r-lg); diff --git a/php/config.local.example.php b/php/config.local.example.php index 7d0beb3..45e51d7 100644 --- a/php/config.local.example.php +++ b/php/config.local.example.php @@ -53,6 +53,17 @@ return [ 'METRIC_DEWPOINT_MARGIN_C' => 0, // 결로(이슬점) 경보 비활성(0 이면 판정 분기 스킵) 'METRIC_ALERT_COOLDOWN_SEC' => 1800, // 동일 종류 경보 재발송 억제(초, 기본 30분) + // 예보·대응(대시보드 예방 신호 + 조치 가이드). 전부 기본값으로 동작 — 필요 시만 override. + // FORECAST_ENABLED(true), FORECAST_WARN_MARGIN_TEMP_C(2.0), FORECAST_WARN_MARGIN_RH(8.0), + // FORECAST_HORIZON_MIN(180), FORECAST_RATE_TEMP_C_PER_10MIN(3.0) 등. config.php 참고. + // 대응 담당자(대시보드 '대응 필요' 섹션에 표시). 비우면 담당자 블록 생략. + 'RESPONSE_CONTACTS' => [ + // ['role' => '시설담당', 'name' => '홍길동', 'phone' => '010-0000-0000'], + // ['role' => '야간당직', 'name' => '김철수', 'phone' => '010-0000-0000'], + ], + // 상태별 대응 체크리스트는 config.php 의 기본 플레이북을 사용. 현장 절차로 바꾸려면 + // 'RESPONSE_PLAYBOOKS' => ['high_temp' => ['title'=>'고온 대응','steps'=>['...']], ...] 로 통째 override. + // 미복귀 경보 2차(야간) 에스컬레이션 수신자. 비우면 비활성(opt-in). 'SMS_RECIPIENTS_TIER2' => [], 'ESCALATE_AFTER_SEC' => 1800, // 경보 미복귀 30분 후 2차 통지 diff --git a/php/config.php b/php/config.php index b90631f..75d6a1a 100644 --- a/php/config.php +++ b/php/config.php @@ -152,6 +152,37 @@ define('LEAK_ALERT_COOLDOWN_SEC', (int)cfg('LEAK_ALERT_COOLDOWN_SEC', 1800)); define('METRIC_DEWPOINT_MARGIN_C', (float)cfg('METRIC_DEWPOINT_MARGIN_C', 3.0)); define('METRIC_DEWPOINT_HYSTERESIS_C', (float)cfg('METRIC_DEWPOINT_HYSTERESIS_C', 1.0)); +// ─── 예보·대응 (대시보드 예방 신호 + 조치 가이드) ───────────────────────────── +// 기존 측정 데이터만으로 "경보가 터지기 전"을 대시보드에 표시한다(php/forecast.php). +// 임계 여유·추세 ETA·이슬점 마진·변화율을 계산하고, 위험 시 상태별 대응 체크리스트를 +// 보여준다. DB 스키마 변경 없음(전부 기존 leak_sensor_metric/status 로 산출). +define('FORECAST_ENABLED', (bool)cfg('FORECAST_ENABLED', true)); +define('FORECAST_RECENT_POINTS', (int)cfg('FORECAST_RECENT_POINTS', 12)); // 회귀에 쓸 최근 측정 수(≈1시간) +define('FORECAST_MIN_POINTS', (int)cfg('FORECAST_MIN_POINTS', 6)); // ETA 산출 최소 점수 +define('FORECAST_MIN_R2', (float)cfg('FORECAST_MIN_R2', 0.5)); // 추세 신뢰(결정계수) 하한 +define('FORECAST_HORIZON_MIN', (int)cfg('FORECAST_HORIZON_MIN', 180)); // ETA 표시 지평(분) +define('FORECAST_WARN_MARGIN_TEMP_C', (float)cfg('FORECAST_WARN_MARGIN_TEMP_C', 2.0)); // 온도 여유 근접 기준(℃) +define('FORECAST_WARN_MARGIN_RH', (float)cfg('FORECAST_WARN_MARGIN_RH', 8.0)); // 습도 여유 근접 기준(%) +define('FORECAST_RATE_TEMP_C_PER_10MIN', (float)cfg('FORECAST_RATE_TEMP_C_PER_10MIN', 3.0)); // 급변 기준(℃/10분) +define('FORECAST_RATE_RH_PER_10MIN', (float)cfg('FORECAST_RATE_RH_PER_10MIN', 10.0)); // 급변 기준(%/10분) + +// 상태별 대응 체크리스트(플레이북). 조건 키 → ['title'=>, 'steps'=>[...]]. +// config.local.php 에서 'RESPONSE_PLAYBOOKS' 로 통째 오버라이드 가능(현장 절차 반영). +$default_playbooks = [ + 'leak' => ['title' => '누수 대응', 'steps' => ['현장 즉시 확인', '해당 구역 급수 밸브 차단', '누수 지점·범위 사진 확인', '전기 설비 침수 여부 안전 점검', '조치 후 복구 확인']], + 'high_temp' => ['title' => '고온 대응', 'steps' => ['냉방·환기 동작 확인', '실외기/흡배기·필터 막힘 점검', '발열 장비·부하 확인', '필요 시 임시 냉각']], + 'low_temp' => ['title' => '저온 대응', 'steps' => ['난방 동작 확인', '외기 유입·틈새 점검', '동결 위험 배관 확인']], + 'high_humidity' => ['title' => '고습 대응', 'steps' => ['제습기 동작 확인', '누수·침수 여부 확인', '환기 상태 점검']], + 'low_humidity' => ['title' => '저습 대응', 'steps' => ['가습 설비 확인', '정전기 민감 장비 점검']], + 'condensation' => ['title' => '결로 위험 대응', 'steps' => ['차가운 표면·배관 결로 확인', '제습·환기로 이슬점 마진 확보', '장비 표면 물기 점검']], + 'offline' => ['title' => '오프라인 대응', 'steps' => ['장비 전원·상태 표시등 확인', '유선/무선 네트워크 연결 확인', 'NTP 시각 동기 확인', '필요 시 장비 재부팅', '복구 후 수신 재개 확인']], +]; +define('RESPONSE_PLAYBOOKS', cfg('RESPONSE_PLAYBOOKS', $default_playbooks)); + +// 대응 담당자 연락처: [['role'=>'시설담당','name'=>'홍길동','phone'=>'010-...'], ...]. +// 비면(기본) 대응 섹션에서 연락처 블록을 생략한다. config.local.php 에서 채운다. +define('RESPONSE_CONTACTS', cfg('RESPONSE_CONTACTS', [])); + // ─── SMS 잔액(잔여 건수) 사전경보 ───────────────────────────────────────────── // Cafe24 응답의 잔여 건수가 이 값 이하로 내려가면 1회 '[SMS 잔액부족]' 경고를 보낸다 // (경보 채널이 잔액 소진으로 조용히 죽는 것을 사전 방지). 쿨다운으로 반복을 억제. diff --git a/php/dashboard.php b/php/dashboard.php index a355cb5..75abfbc 100644 --- a/php/dashboard.php +++ b/php/dashboard.php @@ -314,6 +314,70 @@ if ($alert_count > 0) { $system_title = '정상 감시 중'; $system_detail = '현재 온습도 경보와 중요 운영 경고가 없습니다.'; } + +// ── 예보·대응 신호 (forecast.php — 순수계산) ────────────────────────────── +// 기존 측정 데이터만으로 임계 여유·추세 ETA·이슬점·변화율·신선도를 산출한다. +// $statuses 는 모든 센서(누수 노드 포함)를 담으므로 한 번만 순회한다(중복 없음). +require_once __DIR__ . '/forecast.php'; +$forecasts = []; // sensor_id => forecast_metric() 결과 (카드 예보줄에 사용) +$risk_items_all = []; // 대응 섹션용 정규화 위험 목록 +if (FORECAST_ENABLED) { + $fc_cfg = [ + 'temp_high' => METRIC_TEMP_HIGH_C, + 'temp_low' => METRIC_TEMP_LOW_C, + 'rh_high' => METRIC_RH_HIGH, + 'rh_low' => (METRIC_RH_LOW >= 0 ? METRIC_RH_LOW : null), // -1 = 비활성 + 'dew_margin'=> METRIC_DEWPOINT_MARGIN_C, // 0 이하 = 결로 판정 스킵 + 'warn_temp' => FORECAST_WARN_MARGIN_TEMP_C, + 'warn_rh' => FORECAST_WARN_MARGIN_RH, + 'horizon_min'=> FORECAST_HORIZON_MIN, + 'min_points'=> FORECAST_MIN_POINTS, + 'min_r2' => FORECAST_MIN_R2, + 'rate_temp_10'=> FORECAST_RATE_TEMP_C_PER_10MIN, + 'rate_rh_10' => FORECAST_RATE_RH_PER_10MIN, + ]; + $now_ts = time(); + $recent_stmt = $db->prepare( + "SELECT temperature_c, humidity_percent, UNIX_TIMESTAMP(measured_at) AS t + FROM leak_sensor_metric WHERE sensor_id = :sid ORDER BY id DESC LIMIT :lim" + ); + foreach ($statuses as $s) { + $sid = (int)$s['sensor_id']; + $metric = $s['latest_metric'] ?? null; + $temp = ($metric && $metric['temperature_c'] !== null) ? (float)$metric['temperature_c'] : null; + $hum = ($metric && $metric['humidity_percent'] !== null) ? (float)$metric['humidity_percent'] : null; + + $recent = []; + if ($temp !== null || $hum !== null) { + $recent_stmt->bindValue(':sid', $sid, PDO::PARAM_INT); + $recent_stmt->bindValue(':lim', (int)FORECAST_RECENT_POINTS, PDO::PARAM_INT); + $recent_stmt->execute(); + foreach (array_reverse($recent_stmt->fetchAll()) as $rr) { // 오래된→최신 + $recent[] = [ + 't' => (int)$rr['t'], + 'temp' => $rr['temperature_c'] !== null ? (float)$rr['temperature_c'] : null, + 'hum' => $rr['humidity_percent'] !== null ? (float)$rr['humidity_percent'] : null, + ]; + } + } + + $fc = forecast_metric($temp, $hum, $recent, $fc_cfg); + $forecasts[$sid] = $fc; + + // 신선도 — offline 은 대시보드 기존 판정($s['is_offline'])을 단일 진실원으로 사용 + $fr = forecast_freshness(strtotime($s['last_seen'] ?? '') ?: null, $now_ts, HEARTBEAT_TIMEOUT_SEC); + $fr['offline'] = (bool)$s['is_offline']; + if ($fr['offline']) { $fr['imminent'] = false; } + + $risk_items_all = array_merge($risk_items_all, risk_items( + ['sensor_id' => $sid, 'name' => $s['sensor_name'], 'location' => $s['device_location'], + 'is_leak' => (int)($s['is_leak'] ?? 0)], + $fc, $fr + )); + } + // critical 먼저, 그다음 warn + usort($risk_items_all, fn($a, $b) => ($a['severity'] === 'critical' ? 0 : 1) <=> ($b['severity'] === 'critical' ? 0 : 1)); +} ?> @@ -417,6 +481,44 @@ if ($alert_count > 0) { + +
+

대응 필요 건

+ +
+
+ + + +
+ +
+
— 권장 조치
+
    +
    2. +
+
+ +
+ + +
+ 담당자 + + + + + +
+ +
+ +
@@ -564,6 +666,35 @@ if ($alert_count > 0) { (분 전)
+ $h['state'], 'text' => $h['label'] . ' 여유 ' . $h['margin'] . $h['unit']]; + } + foreach (['temp', 'hum'] as $vk) { + $e = $fc['trend'][$vk]['eta'] ?? null; + if ($e) { $chips[] = ['state' => 'near', 'text' => '↗ ~' . (int)$e['min'] . '분 뒤 ' . metric_status_label($e['key'])]; } + } + foreach ($fc['rate'] as $r) { + $chips[] = ['state' => 'near', 'text' => ($r['spike'] === 'rise' ? '급상승 ' : '급락 ') . ($r['per10'] > 0 ? '+' : '') . $r['per10'] . '/10분']; + } + if (!empty($fc['dew']['risk'])) { + $chips[] = ['state' => 'near', 'text' => '결로주의 (마진 ' . $fc['dew']['margin'] . '℃)']; + } + ?> +
+ + + + 예보: 임계 여유 충분 + +
+ diff --git a/php/forecast.php b/php/forecast.php new file mode 100644 index 0000000..12c844d --- /dev/null +++ b/php/forecast.php @@ -0,0 +1,272 @@ + 100.0) { + return null; + } + $a = 17.62; + $b = 243.12; + $gamma = log($rh / 100.0) + ($a * $t) / ($b + $t); + if ($a - $gamma == 0.0) { + return null; + } + return ($b * $gamma) / ($a - $gamma); +} + +/** + * 최소자승 선형회귀. points = [[x, y], ...] (x·y float). + * 반환 ['slope'=>Δy/Δx, 'intercept'=>, 'r2'=>0..1, 'n'=>] 또는 null(점<2, x 분산 0). + */ +function linear_trend(array $points): ?array +{ + $n = count($points); + if ($n < 2) { + return null; + } + $sx = $sy = $sxx = $sxy = $syy = 0.0; + foreach ($points as $p) { + $x = (float)$p[0]; + $y = (float)$p[1]; + $sx += $x; $sy += $y; $sxx += $x * $x; $sxy += $x * $y; $syy += $y * $y; + } + $denom = $n * $sxx - $sx * $sx; + if (abs($denom) < 1e-9) { + return null; // x 분산 0 (모두 같은 시각) → 기울기 미정 + } + $slope = ($n * $sxy - $sx * $sy) / $denom; + $intercept = ($sy - $slope * $sx) / $n; + + $ss_tot = $syy - ($sy * $sy) / $n; + $ss_res = 0.0; + foreach ($points as $p) { + $e = (float)$p[1] - ($intercept + $slope * (float)$p[0]); + $ss_res += $e * $e; + } + $r2 = $ss_tot > 1e-12 ? max(0.0, 1.0 - $ss_res / $ss_tot) : 1.0; + + return ['slope' => $slope, 'intercept' => $intercept, 'r2' => $r2, 'n' => $n]; +} + +/** + * 임계 도달 ETA(분). cur 현재값, thr 임계, slope_per_min 분당 변화, horizon_min 지평. + * 임계 방향으로 수렴 중이고 지평 내일 때만 양수 분, 아니면 null. + */ +function threshold_eta_min(float $cur, float $thr, float $slope_per_min, float $horizon_min): ?float +{ + if (abs($slope_per_min) < 1e-6) { + return null; // 정체 → 도달 안 함 + } + $eta = ($thr - $cur) / $slope_per_min; + if ($eta <= 0.0 || $eta > $horizon_min) { + return null; // 이미 넘음/멀어짐 or 지평 밖 + } + return $eta; +} + +/** + * 여유 상태 3단계. margin = 현재값과 임계 거리(임계 안 넘었으면 양수). + * over: 이미 임계 초과, near: warn_margin 이내 근접, safe: 여유. + */ +function headroom_state(float $margin, float $warn_margin): string +{ + if ($margin <= 0.0) { + return 'over'; + } + if ($margin <= $warn_margin) { + return 'near'; + } + return 'safe'; +} + +/** + * points 를 [분(minutes-from-first), value] 로 변환(회귀 x축). null 값은 제외. + * recent = [['t'=>epoch_sec, vkey=>?float], ...] 오래된→최신. + */ +function _series_minutes(array $recent, string $vkey): array +{ + $pts = []; + $t0 = null; + foreach ($recent as $r) { + $v = $r[$vkey] ?? null; + if ($v === null) { + continue; + } + $t = (int)($r['t'] ?? 0); + if ($t0 === null) { + $t0 = $t; + } + $pts[] = [($t - $t0) / 60.0, (float)$v]; + } + return $pts; +} + +/** + * 한 온습도 센서의 예보 계산. + * $temp/$hum: 현재값(없으면 null). $recent: [['t'=>epoch,'temp'=>,'hum'=>],...] 오래된→최신. + * $c: 튜닝/임계 배열 — temp_high, temp_low, rh_high, rh_low(각 null=비활성), + * dew_margin(>0 시 결로 판정), warn_temp, warn_rh, horizon_min, min_points, min_r2, + * rate_temp_10, rate_rh_10. + * 반환: ['headroom'=>[...], 'trend'=>[...], 'dew'=>..|null, 'rate'=>..]. + */ +function forecast_metric(?float $temp, ?float $hum, array $recent, array $c): array +{ + $out = ['headroom' => [], 'trend' => ['temp' => null, 'hum' => null], 'dew' => null, 'rate' => []]; + + // ── 임계 여유 ────────────────────────────────────────────── + $band = [ + ['key' => 'high_temp', 'label' => '고온', 'cur' => $temp, 'thr' => $c['temp_high'] ?? null, 'unit' => '℃', 'above' => true, 'warn' => $c['warn_temp'] ?? 2.0], + ['key' => 'low_temp', 'label' => '저온', 'cur' => $temp, 'thr' => $c['temp_low'] ?? null, 'unit' => '℃', 'above' => false, 'warn' => $c['warn_temp'] ?? 2.0], + ['key' => 'high_humidity', 'label' => '고습', 'cur' => $hum, 'thr' => $c['rh_high'] ?? null, 'unit' => '%', 'above' => true, 'warn' => $c['warn_rh'] ?? 8.0], + ['key' => 'low_humidity', 'label' => '저습', 'cur' => $hum, 'thr' => $c['rh_low'] ?? null, 'unit' => '%', 'above' => false, 'warn' => $c['warn_rh'] ?? 8.0], + ]; + foreach ($band as $b) { + if ($b['cur'] === null || $b['thr'] === null) { + continue; // 값 없거나 임계 비활성 + } + // above=고온/고습: margin = thr - cur (양수면 아직 여유). below=저온/저습: margin = cur - thr. + $margin = $b['above'] ? ($b['thr'] - $b['cur']) : ($b['cur'] - $b['thr']); + $out['headroom'][] = [ + 'key' => $b['key'], 'label' => $b['label'], 'thr' => $b['thr'], 'unit' => $b['unit'], + 'margin' => round($margin, 1), 'state' => headroom_state($margin, (float)$b['warn']), + ]; + } + + // ── 추세 ETA + 변화율 ────────────────────────────────────── + $minPts = (int)($c['min_points'] ?? 6); + $minR2 = (float)($c['min_r2'] ?? 0.5); + $horizon = (float)($c['horizon_min'] ?? 180.0); + foreach (['temp' => [$temp, $c['temp_high'] ?? null, $c['temp_low'] ?? null], + 'hum' => [$hum, $c['rh_high'] ?? null, $c['rh_low'] ?? null]] as $vk => $vv) { + [$cur, $thrHigh, $thrLow] = $vv; + if ($cur === null) { + continue; + } + $series = _series_minutes($recent, $vk); + if (count($series) < $minPts) { + continue; + } + $tr = linear_trend($series); + if ($tr === null) { + continue; + } + $slope = $tr['slope']; // 분당 + $dir = abs($slope) < 1e-4 ? 'flat' : ($slope > 0 ? 'up' : 'down'); + $eta = null; + if ($tr['r2'] >= $minR2) { + if ($slope > 0 && $thrHigh !== null) { + $m = threshold_eta_min($cur, (float)$thrHigh, $slope, $horizon); + if ($m !== null) { $eta = ['key' => $vk === 'temp' ? 'high_temp' : 'high_humidity', 'min' => (int)round($m)]; } + } elseif ($slope < 0 && $thrLow !== null) { + $m = threshold_eta_min($cur, (float)$thrLow, $slope, $horizon); + if ($m !== null) { $eta = ['key' => $vk === 'temp' ? 'low_temp' : 'low_humidity', 'min' => (int)round($m)]; } + } + } + $out['trend'][$vk] = ['slope_min' => round($slope, 3), 'dir' => $dir, 'r2' => round($tr['r2'], 2), 'eta' => $eta]; + + // ── 변화율 급변(최근 ~10분 환산) ── + $rateKey = $vk === 'temp' ? 'rate_temp_10' : 'rate_rh_10'; + $rateThr = (float)($c[$rateKey] ?? ($vk === 'temp' ? 3.0 : 10.0)); + $per10 = $slope * 10.0; + if (abs($per10) >= $rateThr) { + $out['rate'][$vk] = ['per10' => round($per10, 1), 'spike' => $slope > 0 ? 'rise' : 'drop']; + } + } + + // ── 이슬점 결로 마진 ─────────────────────────────────────── + $td = dew_point_c($temp, $hum); + if ($td !== null && $temp !== null) { + $margin = $temp - $td; + $dewCfg = (float)($c['dew_margin'] ?? 0.0); + $out['dew'] = [ + 'td' => round($td, 1), + 'margin' => round($margin, 1), + 'risk' => ($dewCfg > 0.0 && $margin <= $dewCfg), + 'enabled' => ($dewCfg > 0.0), + ]; + } + + return $out; +} + +/** + * 신선도/오프라인 임박. last_seen_ts=epoch(없으면 null), now=epoch, timeout=하트비트 초. + * warn_ratio: timeout 대비 이 비율 넘으면 임박(기본 0.5). + * 반환 ['elapsed_sec','elapsed_min','offline','imminent']. + */ +function forecast_freshness(?int $last_seen_ts, int $now, int $timeout, float $warn_ratio = 0.5): array +{ + if ($last_seen_ts === null || $timeout <= 0) { + return ['elapsed_sec' => null, 'elapsed_min' => null, 'offline' => false, 'imminent' => false]; + } + $elapsed = max(0, $now - $last_seen_ts); + return [ + 'elapsed_sec' => $elapsed, + 'elapsed_min' => (int)round($elapsed / 60), + 'offline' => $elapsed >= $timeout, + 'imminent' => ($elapsed < $timeout) && ($elapsed >= $timeout * $warn_ratio), + ]; +} + +/** + * 대응 섹션용 위험 항목 목록. sensor(sensor_id·name·location·is_leak) + forecast(metric) + freshness 를 받아 + * [['condition'=>, 'severity'=>'critical|warn', 'sensor_id'=>, 'name'=>, 'location'=>, 'detail'=>, 'eta_min'=>?], ...] 정규화. + * 이미 발생한 것(critical) + 임박한 것(warn) 모두 포함. 대시보드가 이 목록으로 플레이북/담당자를 붙인다. + */ +function risk_items(array $sensor, array $forecast, array $freshness): array +{ + $risks = []; + $meta = ['sensor_id' => $sensor['sensor_id'] ?? 0, 'name' => $sensor['name'] ?? '', 'location' => $sensor['location'] ?? '']; + + // 오프라인(critical) / 임박(warn) + if (!empty($freshness['offline'])) { + $risks[] = $meta + ['condition' => 'offline', 'severity' => 'critical', 'detail' => "마지막 수신 {$freshness['elapsed_min']}분 전", 'eta_min' => null]; + } elseif (!empty($freshness['imminent'])) { + $risks[] = $meta + ['condition' => 'offline', 'severity' => 'warn', 'detail' => "수신 지연 {$freshness['elapsed_min']}분 — 오프라인 임박", 'eta_min' => null]; + } + + // 누수(critical) + if (!empty($sensor['is_leak'])) { + $risks[] = $meta + ['condition' => 'leak', 'severity' => 'critical', 'detail' => '누수 감지', 'eta_min' => null]; + } + + // 임계 초과(critical) / 근접(warn) — headroom 기반 + foreach ($forecast['headroom'] ?? [] as $h) { + if ($h['state'] === 'over') { + $risks[] = $meta + ['condition' => $h['key'], 'severity' => 'critical', 'detail' => "{$h['label']} 임계 초과 (현재 여유 {$h['margin']}{$h['unit']})", 'eta_min' => null]; + } elseif ($h['state'] === 'near') { + $eta = null; + foreach (['temp', 'hum'] as $vk) { + $e = $forecast['trend'][$vk]['eta'] ?? null; + if ($e && $e['key'] === $h['key']) { $eta = $e['min']; break; } + } + $detail = "{$h['label']} 임계 근접 (여유 {$h['margin']}{$h['unit']})" . ($eta !== null ? " · 현 추세 ~{$eta}분 뒤 도달" : ''); + $risks[] = $meta + ['condition' => $h['key'], 'severity' => 'warn', 'detail' => $detail, 'eta_min' => $eta]; + } + } + + // 결로(warn) — config 활성 시 + if (!empty($forecast['dew']['risk'])) { + $risks[] = $meta + ['condition' => 'condensation', 'severity' => 'warn', 'detail' => "결로 위험 (이슬점 마진 {$forecast['dew']['margin']}℃)", 'eta_min' => null]; + } + + return $risks; +} diff --git a/run_tests.sh b/run_tests.sh index a6a7c89..2da75f4 100644 --- a/run_tests.sh +++ b/run_tests.sh @@ -28,6 +28,7 @@ run_php tests/raw_body_signature_test.php run_php tests/security_hardening_test.php run_php tests/threshold_sms_test.php run_php tests/chaos_drill.php # SMS 실패 주입 경로 드릴 +run_php tests/forecast_test.php # 예보(이슬점/회귀/ETA/여유/신선도/위험) 순수함수 echo "" echo "=== 2) 펌웨어 호스트 패리티 (C↔Python↔PHP) ===" diff --git a/tests/forecast_test.php b/tests/forecast_test.php new file mode 100644 index 0000000..46b77da --- /dev/null +++ b/tests/forecast_test.php @@ -0,0 +1,100 @@ + Td~13.9"); +ok(approx(dew_point_c(20.0, 100.0), 20.0, 0.2), "RH100 -> Td~T"); +ok(dew_point_c(20.0, 0.0) === null, "RH0 -> null"); +ok(dew_point_c(20.0, 150.0) === null, "RH>100 -> null"); +ok(dew_point_c(null, 50.0) === null, "temp null -> null"); + +echo "== linear_trend ==\n"; +$t = linear_trend([[0, 0], [1, 2], [2, 4]]); +ok($t !== null && approx($t['slope'], 2.0, 1e-6) && approx($t['r2'], 1.0, 1e-6), "완전선형 slope=2 r2=1"); +ok(linear_trend([[1, 5], [1, 7]]) === null, "동일 x -> null"); +ok(linear_trend([[0, 1]]) === null, "점<2 -> null"); +$t2 = linear_trend([[0, 5], [1, 5], [2, 5]]); +ok($t2 !== null && approx($t2['slope'], 0.0, 1e-6), "평탄 slope~0"); + +echo "== threshold_eta_min ==\n"; +ok(approx(threshold_eta_min(20.0, 27.0, 0.1, 180.0), 70.0, 0.5), "20->27 @0.1/min = 70분"); +ok(threshold_eta_min(20.0, 27.0, -0.1, 180.0) === null, "멀어지는 방향 -> null"); +ok(threshold_eta_min(20.0, 27.0, 0.01, 180.0) === null, "지평(700분>180) 밖 -> null"); +ok(threshold_eta_min(20.0, 27.0, 0.0, 180.0) === null, "정체 slope0 -> null"); + +echo "== headroom_state ==\n"; +ok(headroom_state(5.0, 2.0) === 'safe', "여유 5 warn2 -> safe"); +ok(headroom_state(1.0, 2.0) === 'near', "여유 1 warn2 -> near"); +ok(headroom_state(-1.0, 2.0) === 'over', "여유 -1 -> over"); + +echo "== forecast_metric ==\n"; +// 상승 온도 시리즈: 24.0 -> 26.5 (6점, 5분 간격) -> slope ~ 0.1/min (완전선형) +$recent = []; +for ($i = 0; $i < 6; $i++) { + $recent[] = ['t' => 1751000000 + $i * 300, 'temp' => 24.0 + $i * 0.5, 'hum' => 55.0]; +} +$cfg = ['temp_high' => 30.0, 'temp_low' => 10.0, 'rh_high' => 70.0, 'rh_low' => 20.0, + 'dew_margin' => 3.0, 'warn_temp' => 3.0, 'warn_rh' => 8.0, + 'horizon_min' => 180.0, 'min_points' => 6, 'min_r2' => 0.5, 'rate_temp_10' => 3.0, 'rate_rh_10' => 10.0]; +$fc = forecast_metric(26.5, 55.0, $recent, $cfg); +$hh = null; +foreach ($fc['headroom'] as $h) { if ($h['key'] === 'high_temp') $hh = $h; } +ok($hh !== null && approx($hh['margin'], 3.5, 0.01), "고온 여유 = 30-26.5 = 3.5"); +ok($hh !== null && $hh['state'] === 'safe', "여유 3.5 > warn 3.0 -> safe"); +ok($fc['trend']['temp'] !== null && $fc['trend']['temp']['dir'] === 'up', "온도 추세 상승"); +ok($fc['trend']['temp']['eta'] !== null && $fc['trend']['temp']['eta']['key'] === 'high_temp', "고온 ETA 산출(추세로 도달 예상)"); +ok($fc['dew'] !== null && $fc['dew']['enabled'] === true, "이슬점 마진 계산(활성)"); + +// 저온 임계 비활성(null) -> low_temp headroom 없음 +$fc2 = forecast_metric(15.0, 50.0, [], ['temp_high' => 30.0, 'temp_low' => null, 'rh_high' => 70.0, 'rh_low' => null]); +$hasLow = false; +foreach ($fc2['headroom'] as $h) { if ($h['key'] === 'low_temp') $hasLow = true; } +ok(!$hasLow, "저온 임계 null -> low_temp 여유 스킵"); + +// 누수전용(메트릭 없음) -> headroom 비고, dew null +$fc3 = forecast_metric(null, null, [], $cfg); +ok(count($fc3['headroom']) === 0 && $fc3['dew'] === null, "온습도 null -> 메트릭 신호 스킵"); + +echo "== forecast_freshness ==\n"; +$now = 1751000000; +ok(forecast_freshness($now - 100, $now, 1200)['offline'] === false, "100초 경과 -> 온라인"); +ok(forecast_freshness($now - 1300, $now, 1200)['offline'] === true, "1300>1200 -> 오프라인"); +$fr = forecast_freshness($now - 700, $now, 1200, 0.5); +ok($fr['imminent'] === true && $fr['offline'] === false, "700>=600(50%) 아직<1200 -> 임박"); +ok(forecast_freshness(null, $now, 1200)['offline'] === false, "last null -> offline false"); + +echo "== risk_items ==\n"; +$leakRisks = risk_items(['sensor_id' => 3, 'name' => '지하배관', 'is_leak' => 1], ['headroom' => []], ['offline' => false]); +ok(count($leakRisks) === 1 && $leakRisks[0]['condition'] === 'leak' && $leakRisks[0]['severity'] === 'critical', "누수 -> leak critical"); +$offRisks = risk_items(['sensor_id' => 1, 'name' => 'A'], ['headroom' => []], ['offline' => true, 'elapsed_min' => 25]); +ok(count($offRisks) === 1 && $offRisks[0]['condition'] === 'offline' && $offRisks[0]['severity'] === 'critical', "오프라인 -> offline critical"); +$overRisks = risk_items(['sensor_id' => 2, 'name' => 'B'], ['headroom' => [['key' => 'high_temp', 'label' => '고온', 'margin' => -1.2, 'unit' => '℃', 'state' => 'over']], 'trend' => []], ['offline' => false]); +ok(count($overRisks) === 1 && $overRisks[0]['severity'] === 'critical', "고온 초과 -> critical"); +$nearRisks = risk_items(['sensor_id' => 2, 'name' => 'B'], + ['headroom' => [['key' => 'high_temp', 'label' => '고온', 'margin' => 1.5, 'unit' => '℃', 'state' => 'near']], + 'trend' => ['temp' => ['eta' => ['key' => 'high_temp', 'min' => 40]], 'hum' => null]], ['offline' => false]); +ok(count($nearRisks) === 1 && $nearRisks[0]['severity'] === 'warn' && $nearRisks[0]['eta_min'] === 40, "고온 근접 -> warn + ETA 40"); + +echo "\n==== forecast_test: $tests 검사, $fails 실패 ====\n"; +exit($fails === 0 ? 0 : 1);