Wideband autonomous SDR analysis engine forked from sdr-visual-suite
您最多选择25个主题 主题必须以字母或数字开头,可以包含连字符 (-),并且长度不得超过35个字符

129 行
4.2KB

  1. package main
  2. import (
  3. "context"
  4. "encoding/json"
  5. "log"
  6. "os"
  7. "runtime/debug"
  8. "strings"
  9. "sync"
  10. "time"
  11. "sdr-wideband-suite/internal/config"
  12. "sdr-wideband-suite/internal/detector"
  13. "sdr-wideband-suite/internal/dsp"
  14. "sdr-wideband-suite/internal/recorder"
  15. )
  16. func runDSP(ctx context.Context, srcMgr *sourceManager, cfg config.Config, det *detector.Detector, window []float64, h *hub, eventFile *os.File, eventMu *sync.RWMutex, updates <-chan dspUpdate, gpuState *gpuStatus, rec *recorder.Manager, sigSnap *signalSnapshot, extractMgr *extractionManager) {
  17. defer func() {
  18. if r := recover(); r != nil {
  19. log.Printf("FATAL: runDSP goroutine panic: %v\n%s", r, debug.Stack())
  20. }
  21. }()
  22. rt := newDSPRuntime(cfg, det, window, gpuState)
  23. ticker := time.NewTicker(cfg.FrameInterval())
  24. defer ticker.Stop()
  25. logTicker := time.NewTicker(5 * time.Second)
  26. defer logTicker.Stop()
  27. enc := json.NewEncoder(eventFile)
  28. dcBlocker := dsp.NewDCBlocker(0.995)
  29. for {
  30. select {
  31. case <-ctx.Done():
  32. return
  33. case <-logTicker.C:
  34. st := srcMgr.Stats()
  35. log.Printf("stats: buf=%d drop=%d reset=%d last=%dms", st.BufferSamples, st.Dropped, st.Resets, st.LastSampleAgoMs)
  36. case upd := <-updates:
  37. rt.applyUpdate(upd, srcMgr, rec, gpuState)
  38. dcBlocker.Reset()
  39. ticker.Reset(rt.cfg.FrameInterval())
  40. case <-ticker.C:
  41. art, err := rt.captureSpectrum(srcMgr, rec, dcBlocker, gpuState)
  42. if err != nil {
  43. log.Printf("read IQ: %v", err)
  44. if strings.Contains(err.Error(), "timeout") {
  45. if err := srcMgr.Restart(rt.cfg); err != nil {
  46. log.Printf("restart failed: %v", err)
  47. }
  48. }
  49. continue
  50. }
  51. if !rt.gotSamples {
  52. log.Printf("received IQ samples")
  53. rt.gotSamples = true
  54. }
  55. finished := art.finished
  56. thresholds := art.thresholds
  57. noiseFloor := art.noiseFloor
  58. var displaySignals []detector.Signal
  59. if len(art.iq) > 0 {
  60. displaySignals = rt.refineSignals(art, extractMgr, rec)
  61. if rec != nil && len(displaySignals) > 0 && len(art.allIQ) > 0 {
  62. aqCfg := extractionConfig{firTaps: rt.cfg.Recorder.ExtractionTaps, bwMult: rt.cfg.Recorder.ExtractionBwMult}
  63. streamSnips, streamRates := extractForStreaming(extractMgr, art.allIQ, rt.cfg.SampleRate, rt.cfg.CenterHz, displaySignals, rt.streamPhaseState, rt.streamOverlap, aqCfg)
  64. items := make([]recorder.StreamFeedItem, 0, len(displaySignals))
  65. for j, ds := range displaySignals {
  66. if ds.ID == 0 || ds.Class == nil {
  67. continue
  68. }
  69. if j >= len(streamSnips) || len(streamSnips[j]) == 0 {
  70. continue
  71. }
  72. snipRate := rt.cfg.SampleRate
  73. if j < len(streamRates) && streamRates[j] > 0 {
  74. snipRate = streamRates[j]
  75. }
  76. items = append(items, recorder.StreamFeedItem{Signal: ds, Snippet: streamSnips[j], SnipRate: snipRate})
  77. }
  78. if len(items) > 0 {
  79. rec.FeedSnippets(items)
  80. }
  81. }
  82. rt.maintenance(displaySignals, rec)
  83. } else {
  84. displaySignals = rt.det.StableSignals()
  85. }
  86. if sigSnap != nil {
  87. sigSnap.set(displaySignals)
  88. }
  89. eventMu.Lock()
  90. for _, ev := range finished {
  91. _ = enc.Encode(ev)
  92. }
  93. eventMu.Unlock()
  94. if rec != nil && len(finished) > 0 {
  95. evCopy := make([]detector.Event, len(finished))
  96. copy(evCopy, finished)
  97. rec.OnEvents(evCopy)
  98. }
  99. var debugInfo *SpectrumDebug
  100. if len(thresholds) > 0 || len(displaySignals) > 0 || noiseFloor != 0 {
  101. scoreDebug := make([]map[string]any, 0, len(displaySignals))
  102. for _, s := range displaySignals {
  103. if s.Class == nil || len(s.Class.Scores) == 0 {
  104. scoreDebug = append(scoreDebug, map[string]any{"center_hz": s.CenterHz, "class": nil})
  105. continue
  106. }
  107. scores := make(map[string]float64, len(s.Class.Scores))
  108. for k, v := range s.Class.Scores {
  109. scores[string(k)] = v
  110. }
  111. scoreDebug = append(scoreDebug, map[string]any{
  112. "center_hz": s.CenterHz,
  113. "mod_type": s.Class.ModType,
  114. "confidence": s.Class.Confidence,
  115. "second_best": s.Class.SecondBest,
  116. "scores": scores,
  117. })
  118. }
  119. debugInfo = &SpectrumDebug{Thresholds: thresholds, NoiseFloor: noiseFloor, Scores: scoreDebug}
  120. }
  121. h.broadcast(SpectrumFrame{Timestamp: art.now.UnixMilli(), CenterHz: rt.cfg.CenterHz, SampleHz: rt.cfg.SampleRate, FFTSize: rt.cfg.FFTSize, Spectrum: art.spectrum, Signals: displaySignals, Debug: debugInfo})
  122. }
  123. }
  124. }