1
$\begingroup$

I have a slowly varying sinusoidal system. The output is separated using methods that produce a real signal $x(t)$ and quadrature component $y(t)$ such that $x(t)$ represents the real physical value of the signal, $x^2 + y^2$ is instantaneous amplitude and $\mbox{atan}(y,x)$ is the instantaneous cosine argument with respect to a central frequency $\omega_c$. So you can think of it as $\omega_c + \phi(t)$ but those two pieces are bundled and the phase would have to be unwrapped to get at $\phi$.

I guess these qualities are unique enough to say that $z(t)=x+iy$ is the Analytic Signal representation of $x(t)$ for that frequency band and that $iy=\hat{x}$ is the Hilbert transform of $x$ in some bandpassed sense, but my knowledge is limited in this area.

Let's say this signal is the output of a linear and mostly time-invariant system with complex response $H(\omega)$. I am only interested in $H(\omega_c)$ at the widely spaced central frequencies and assume it is constant across the passed band around {\omega_c}. $H$ produces both amplitude and phase change but its characteristics are not known quantitatively. I can only probe it with the input and output I have, not probe it.

The idea I am exploring here is how to use the data I have to infer $H(\omega_c)$. If I take the complex output $Z$ as the Analytic Signal and input $W$ and do a complex regression $Z(\omega,t) \approx H(\omega) W(\omega,t) + \varepsilon(t)$ on the input using time to provide samples and yielding an estimate $\hat{H}$ is that an estimate of $H$? A good one? One concern I have is that $Z$ and $W$ likely have no negative frequency content, so I'm concerned that I'll be learning about the analytical analog of $H$ or something like that.

$\endgroup$

1 Answer 1

1
$\begingroup$

First to clarify- real signals include the positive and negative frequencies, so we would be looking at $H(-\omega)$ when we evaluate with real signals. If $H(\omega)$ itself was complex, then the resulting output even with a real input would be a complex output since the positive and negative frequencies would be processed with different responses. When $H(\omega)$ is real, then the resulting output for a real input would also be real.

What needs to be clarified in the question is if the output is the baseband equivalent signal (I assume so), in that the magnitude and phase represents the magnitude and phase changing versus time of some particular carrier for the passband signal. This way, given a real output as described in the first paragraph, we can by using the Analytic Signal representation, represent it as:

$$x(t)+j\hat{x}(t) = A(t)e^{j\phi(t)}e^{j\omega_c t}$$

Where $\hat{x}(t)$ is the Hilbert Transform of $x(t)$, and $\omega_c(t)$ is a carrier frequency of interest. The real waveform, $x(t)$ in this case would be $x(t)=A(t)\cos(\omega_c t + \phi(t)$) where $A(t)$ is the amplitude modulation and $\phi(t)$ is the phase modulation.

The "Analytic Signal" by definition has no negative frequency components so would be a single-sided (positive frequencies only) passband signal as represented above. What I suspect is the OP is looking at the complex baseband equivalent signal, which does have positive and negative frequency components, and is the frequency translated Analytic Signal as give by:

$$x(t) = A(t)e^{j\phi(t)}$$

In any event, what is desired here is a channel estimation where we know the input and output and from that we wish to estimate the channel response $H(\omega)$. If the channel is stationary and we only wish to know this response at discrete frequencies, then one option is to simply compare the input and output at each of these frequencies and from that derive the magnitude and phase for each given frequency. For more complicated cases where we wish to have a complete estimate of the channel over a given bandwidth, we can use the Wiener-Hopf equations for least squares channel estimate as detailed in this link.

Also related is this link showing the use of a frequency chirp (ramp) along with a Tukey windowed FFT for accurate frequency response measurement.

$\endgroup$
5
  • $\begingroup$ I'm happy to clarify, but I'm learning a lot fast. There is no carrier signal here. It is akin to taking the original real signal (x(t)) and decomposing it using a wavelet/windowed method such that the output isolates the band and has two outputs (x(t), y(t)) where x(t) is the bandlimited part of the signal and y(t) provides the additional info needed to get amplitude and phase. $\endgroup$
    – Eli S
    Commented Jan 7, 2023 at 18:17
  • $\begingroup$ @EliS If there is no "carrier", what does the phase represent to you? If the signal is real, then the phase in the analytic signal must be the phase of your signal relative to some particular reference sinusoid, is it not? (That is what I am referring to as a "carrier"). Sorry for the confusion. $\endgroup$ Commented Jan 7, 2023 at 19:28
  • $\begingroup$ I've clarified. I've learned a lot from your response, but I'm hoping to learn how to do the best I can do with the ingredients I have. Your comment "...simply compare the input and output at each of these frequencies" is the essence of what I am asking in my last paragraph. Do I just average out/in? Does it matter that both are transformed to analytical? Will the answer be an estimate of $H$ or its analytical transform? $\endgroup$
    – Eli S
    Commented Jan 8, 2023 at 18:41
  • $\begingroup$ //"If there is no "carrier", what does the phase represent to you?"// ---- well, $\omega_c$ is sorta arbitary and the time-varying phase $\phi(t)$ can absorb whatever is left over after setting $\omega_c$ to whatever you want. You could set it to zero. $\endgroup$ Commented Jan 8, 2023 at 18:59
  • $\begingroup$ @EliS I don't think I'll be able to provide further details here concisely to get you to where you need to be. It's going to turn into a lot of back and forth which is discouraged in this forum (your question isn't in the simple Q/A category). Feel free to email me boschen at loglin dot com and perhaps I can help you further directly. $\endgroup$ Commented Jan 9, 2023 at 3:58

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.