code: def harmonicradar_cw(df, cols, fs, fc):
      for col in cols:
          ttxt = f"CW: {fc} Hz"
          # %% input
          t = df[col]
          tx = np.sin(2 * np.pi * fc * t)
          _, Pxx = signal.welch(tx, fs)
          # %% diode
          d = signal.square(2 * np.pi * fc * t)
          d[d < 0] = 0.0
          # %% output of diode
          rx = tx * d
          df[col + "_HARRAD"] = rx.values
      return df
  
imports: import pandas
  import numpy
  import statsmodels.api
  from scipy import signal, integrate
  from pykalman import UnscentedKalmanFilter
  from tsaug import *
  from fbprophet import Prophet
  import pylab
  from seasonal.periodogram import periodogram
label_desc: Use image node_id 19 for calling a global function with example usage: harmonicradar_cw(df, cols, fs, fc) and returns: df