How to turn S11 to time domain by Matlab

Hi,everyone.I run HFSS to get S11 for my circuit,HFSS also provide time domain for S11.I want to obtain the same time domain response by Matlab.However the result is not even close.Can anyone help~Thx!

The input signal is an impulse and I perform the frequency domain simulation from 0~7.5GHz.

And both attached files are the result from HFSS.

Here's my code:

load S11_re.txt -ascii

load S11_im.txt -ascii

freq=1e+9*S11_re(:,1);

S11=S11_re(:,2)+1i*S11_im(:,2);

TDR=ifft(S11);

Fs=2*max(freq);

Ts=1/Fs;

N=numel(TDR);

tvec=(0:(N-1))*Ts;

plot(tvec,abs(TDR))

[Mathieu NOE]

the frequency domain data must contain the phase also (we need to compute the ifft of a complex valued transfer function, here you provide only the modulus)

[Guan Hao]

@Mathieu NOE Sorry,I'm not that familiar with the time domain signal.

Thanks for your help ! The phase data of frequency domain has been uploaded.

I've shared the complex value of S11 and also modified my code.

Answer 1

[Mathieu NOE]

hello again

I am mot sure to understand the shape of the impulse response from your S11_time.txt file :

why only two positive peaks ? from the Bode plot (first figure) I was expecting some kind of damped oscillations - alike what we get from the ifft (nb the symmetrical computation, including negative and positive frequency data)

also I was unsure what is the time unit in the S11_time.txt file ?

% freq domain (transfer function)
load S11_re.txt -ascii % freq + S11_re
load S11_im.txt -ascii % freq + S11_im
freq=1e+9*S11_re(:,1);
frf=S11_re(:,2)+1i*S11_im(:,2);

figure(1),
subplot(2,1,1),plot(freq,abs(frf))
xlabel('freq (Hz)')
ylabel('amplitude')
title('Impulse Response (IR)');
subplot(2,1,2),plot(freq,180/pi*angle(frf))
xlabel('freq (Hz)')
ylabel('phase(°)')


% IR (impulse response) obtained with ifft method
if mod(length(frf),2)==0 % iseven
    frf_sym = conj(frf(end:-1:2));
else
    frf_sym = conj(frf(end-1:-1:2));
end
    
TDR = real(ifft([frf; frf_sym])); % NB we need the negative and positive frequency complex FRF
TDR = TDR(1:101); % truncation is possible if TDR decays fast enough

Fs=2*max(freq);
Ts=1/Fs;
N=numel(TDR);
tvec=(0:(N-1))*Ts;

% compare to impulse response generated by your software
load S11_time.txt -ascii % time + IR
time = S11_time(:,1);
IR = S11_time(:,2);
time = time/1e9; % assuming time data was given in nanoseconds

figure(2),
plot(tvec,TDR./max(TDR),'b',time,IR./max(IR),'r')
xlabel('time (s)')
ylabel('amplitude')
title('Impulse Response (IR)');
legend('IR from re/im data','IR from file');

[Guan Hao]

@Mathieu NOE Thanks for your help!The S11_time file is the magnitude of the signal,so it returns two positive peak.I'll send another file.

By the way,the time unit is nanosecond.

[Mathieu NOE]

this time it's ok - how can you obtain either two positive peaks or pos / neg ?

results with provided file are matching , I simply removed the header lines so load will work fine (see attached txt file)

all the best

% freq domain (transfer function)
load S11_re.txt -ascii % freq + S11_re
load S11_im.txt -ascii % freq + S11_im
freq=1e+9*S11_re(:,1);
frf=S11_re(:,2)+1i*S11_im(:,2);

figure(1),
subplot(2,1,1),plot(freq,abs(frf))
xlabel('freq (Hz)')
ylabel('amplitude')
title('Impulse Response (IR)');
subplot(2,1,2),plot(freq,180/pi*angle(frf))
xlabel('freq (Hz)')
ylabel('phase(°)')


% IR (impulse response) obtained with ifft method
if mod(length(frf),2)==0 % iseven
    frf_sym = conj(frf(end:-1:2));
else
    frf_sym = conj(frf(end-1:-1:2));
end
    
TDR = real(ifft([frf; frf_sym])); % NB we need the negative and positive frequency complex FRF
TDR = TDR(1:101); % truncation is possible if TDR decays fast enough

Fs=2*max(freq);
Ts=1/Fs;
N=numel(TDR);
tvec=(0:(N-1))*Ts;

% compare to impulse response generated by your software
load S11_time.txt -ascii % time + IR
% remove the header lines otherwise load will throw an error message
% ============================================================================================
% ANSOFT                                                                              04/18/24
% Terminal S Parameter Plot 2                                                         00:46:05
% 
% --------------------------------------------------------------------------------------------
%         Time [ns]                   St(1,1) []
                                              
       
time = S11_time(:,1);
IR = S11_time(:,2);
time = time/1e9; % assuming time data was given in nanoseconds

figure(2),
plot(tvec,TDR./max(TDR),'b',time,IR./max(IR),'r')
xlabel('time (s)')
ylabel('amplitude')
title('Impulse Response (IR)');
legend('IR from re/im data','IR from file');

[Guan Hao]

@Mathieu NOE The software can display S11 in many ways,I sent the absolute value in former file,so it shows two positive peak.If I display the original data,it shows one postive peak and one negative peak.

Thank you for your help,I do learn a lot about how to do ifft.

By the way,I think the software inputs a rectangular pulse with the width of Ts(Time step) and the magnitude is 1/Ts.So maybe I should create a rectangular pulse and convolution with the impulse response?

But I got a question,how do I know where the rectangular pulse locates on the axis of time?

% create input rectangular pulse
t = 0 : Ts : 2e-9; 
Amp=1/Ts; % amplitude of the pulse
u=[zeros(1,numel(0:Ts:Ts)) Amp*ones(1,numel(Ts:Ts:2*Ts)) zeros(1,numel((2*Ts)+(Ts):Ts:2e-9))]; 
plot(u)
TDR=real(ifft([frf; frf_sym])); % NB we need the negative and positive frequency complex FRF
TDR=TDR(1:numel(u)); % truncation is possible if TDR decays fast enough
TDR1=conv(TDR,u)

[Mathieu NOE]

yes, you need to display the data with positive and negative peaks , otherwise we are comparing oranges and apples.

for the convolution with an impulse , an impulse will not modify your IR as a "true" impulse has falt spectrum (infinite)

as you create a rectangular implse with some samples delay, the effect is to delay your IR with a bit or smoothing (less ringing) as the spectrum of a rectangualr impulse has decaying energy as we go into higher frequencies.

% freq domain (transfer function)
load S11_re.txt -ascii % freq + S11_re
load S11_im.txt -ascii % freq + S11_im
freq=1e+9*S11_re(:,1);
frf=S11_re(:,2)+1i*S11_im(:,2);

figure(1),
subplot(2,1,1),plot(freq,abs(frf))
xlabel('freq (Hz)')
ylabel('amplitude')
title('Impulse Response (IR)');
subplot(2,1,2),plot(freq,180/pi*angle(frf))
xlabel('freq (Hz)')
ylabel('phase(°)')


% IR (impulse response) obtained with ifft method
if mod(length(frf),2)==0 % iseven
    frf_sym = conj(frf(end:-1:2));
else
    frf_sym = conj(frf(end-1:-1:2));
end
    
TDR = real(ifft([frf; frf_sym])); % NB we need the negative and positive frequency complex FRF
TDR = TDR(1:101); % truncation is possible if TDR decays fast enough

Fs=2*max(freq);
Ts=1/Fs;
N=numel(TDR);
tvec=(0:(N-1))*Ts;

% compare to impulse response generated by your software
load S11_time.txt -ascii % time + IR
% remove the header lines otherwise load will throw an error message
% ============================================================================================
% ANSOFT                                                                              04/18/24
% Terminal S Parameter Plot 2                                                         00:46:05
% 
% --------------------------------------------------------------------------------------------
%         Time [ns]                   St(1,1) []
                                              
       
time = S11_time(:,1);
IR = S11_time(:,2);
time = time/1e9; % assuming time data was given in nanoseconds

figure(2),
plot(tvec,TDR./max(TDR),'b',time,IR./max(IR),'r')
xlabel('time (s)')
ylabel('amplitude')
title('Impulse Response (IR)');
legend('IR from re/im data','IR from file');



% create input rectangular pulse
t = 0 : Ts : 2e-9; 
Amp=1/Ts; % amplitude of the pulse
u=[zeros(1,numel(0:Ts:Ts)) Amp*ones(1,numel(Ts:Ts:2*Ts)) zeros(1,numel((2*Ts)+(Ts):Ts:2e-9))]; 
figure(3),
% plot(u)
TDR=TDR(1:numel(u)); % truncation is possible if TDR decays fast enough
TDR1=conv(TDR,u)
tt = (0:numel(TDR1)-1)*Ts;
plot(tvec(1:numel(u)),TDR./max(TDR),'b',tt,TDR1./max(TDR1),'r')
xlabel('time (s)')
ylabel('amplitude')
title('Impulse Response (IR)');
legend('TDR','TDR1');

[Guan Hao]

@Mathieu NOE Was it nearly impossible to reconstruct the data produced by the software?

I saw some rectangular shape on the peak which was produced by the software.

[Mathieu NOE]

you notice that the IR lacks a bit of time resolution

we would need to acquire data at a faster rate (and your frequency data must also contain more high frequency data)

the rectangular shape is simply due to a "relatively" low sampling rate - so as you would obtain after downsampling a sharp peak , result is a coarse peak with rectangular shaped peak

[Guan Hao]

@Mathieu NOE I see.Thank you for your help sir!

[Mathieu NOE]

my pleasure !