1
\$\begingroup\$

I asked this question and one of the answers suggested using gating technique for deembedding purposes. I've created a script to transform S-parameters to time domain.

S_ifft = np.fft.ifft(S11)
plt.plot(abs(S_ifft))
plt.show()

When I plot it, it looks like this

enter image description here

Well, that looks like gibberish to me. The board itself is merely a trace with two connectors. The measured data goes from 20 GHz to 28 GHz.

My question is what measurements do I need to take to get a good time domain plot so I can identify the signal of interest? Or how do I analyse this time domain signal to see the connectors?

\$\endgroup\$
  • \$\begingroup\$ Not so bad really. You may need to increase the n \$\endgroup\$ – johnnymopo Nov 17 '15 at 4:34
  • \$\begingroup\$ Are you sure that your data is formatted correctly for input to your ifft() routine? The usual numpy ifft() requires an input that includes a 0 frequency entry and also negative frequency components in the upper half of the input array. \$\endgroup\$ – The Photon Nov 17 '15 at 5:02
  • \$\begingroup\$ With the 20 to 28 GHz frequency information, this can be operated on as if it was -4 to 4GHz already. \$\endgroup\$ – johnnymopo Nov 17 '15 at 5:31
  • \$\begingroup\$ @johnnymopo, well, sure. But if the value at 27 GHz isn't the complex conjugate of the value at 21 GHz (for example), you're going to get unusual (unreal, even) results. \$\endgroup\$ – The Photon Nov 17 '15 at 5:39
  • \$\begingroup\$ That's true, but it is ok. We only care about the magnitude of time. You would get nothing extra by doing it another way \$\endgroup\$ – johnnymopo Nov 17 '15 at 6:25
2
\$\begingroup\$

Edited: Initial answer was too brief and visually not acceptable

Lets think about what you are showing in this image. You are taking some arbitrary set of frequencies and trying to show the time domain representation. As this is a stepped-set of frequencies, we can view the alias-free time as:

T alias-free = (N-1) / BW

For S11 or S22, this is round trip time, in the medium. That puts the effective "range" as:

R alias-free = c * (N-1) / (2 * sqrt(e_effective) * BW)

for 8 GHz bandwidth, 101 points, effective permittivity 2, the alias free range is ~ 1.33 meters. You will likely not alias in your analysis of your PCB.

As for interpreting the image, remember bounce-diagrams from a microwave class: All of the impedances will reflect back to the source, but they will also ring between each other, take a look at the following image from High Frequency Electronics and a paper I wrote on a stepped frequency radar. The radar is viewing a steel pipe covered with ABS plastic. Reflection 1 is from the ABS covering, reflection 2 is from the steel pipe, and if you'll note, there is a lower power right after 2 the exact same dt or dR as 1 to 2. This is a ringing reflection in the ABS plastic. Further more, at point 3 you see the radar glancing over the steel component and ringing from the backside of the ABS. The inset image shows time/range using the inverse DFT and the main part of the image uses a chirp-Z transform to "zoom" in on the detail of interest.

enter image description here

To sum it up, your time-domain S11 will show you the round trip electrical delay between the normalized phase center (your calibration) and the impedance mismatches but care must be taken to discriminate your intended "target" from ringing (or ghosting).

\$\endgroup\$
  • \$\begingroup\$ updated the answer with more detail \$\endgroup\$ – johnnymopo Nov 17 '15 at 13:51

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

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