zcy
/
openEMS
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

documentation update GetField_TD2FD.m

pull/1/head
Sebastian Held 2010-07-15 11:46:36 +02:00
parent 7d8f12fc08
commit 1cc0bb7803
1 changed files with 9 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
matlab/GetField_TD2FD.m
View File

@ -1,18 +1,19 @@
function field_FD = GetField_TD2FD(field_TD, freq) function field_FD = GetField_TD2FD(field_TD, freq)
% function field_FD = GetField_TD2FD(field_TD, freq) % function field_FD = GetField_TD2FD(field_TD, freq)
% %
% transform time-domain field data into frequency domain % Transforms time-domain field data into the frequency domain
% Autocorrects the half-timestep offset of the H-field
% %
% example: % example:
% freq = linspace(0,1e9,100); %target frequency vector % freq = linspace(0,1e9,100); %target frequency vector (Hz)
% field = ReadHDF5FieldData('tmp/Ht.h5'); % field = ReadHDF5FieldData('tmp/Ht.h5');
% field_FD = GetField_TD2FD(field, freq); % field_FD = GetField_TD2FD(field, freq);
% %
% See also ReadHDF5FieldData
%
% openEMS matlab interface % openEMS matlab interface
% ----------------------- % -----------------------
% author: Thorsten Liebig % author: Thorsten Liebig
%
% See also ReadHDF5FieldData
t = field_TD.time; t = field_TD.time;
@ -28,6 +29,10 @@ for n=1:numTS
for nf = 1:numel(freq) for nf = 1:numel(freq)
f = freq(nf); f = freq(nf);
field_FD.values{nf} = field_FD.values{nf} + 2/numTS * field_TD.values{n}.*exp(-1i*2*pi*f*t(n)); field_FD.values{nf} = field_FD.values{nf} + 2/numTS * field_TD.values{n}.*exp(-1i*2*pi*f*t(n));
% t(n) is absolute time and therefore the half-timestep offset of
% the H-field is automatically compensated
% openEMS output: E-fields start at t=0
% openEMS output: H-fields start at t=delta_t/2
end end
end end