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

compensate H-field half time-step offset

pull/1/head
Sebastian Held 2010-04-09 17:30:54 +02:00
parent 9523de1350
commit a08f602a55
1 changed files with 5 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
matlab/examples/Coax_partial_CylinderCoords.m
View File

@ -103,7 +103,11 @@ cd(savePath);
UI = ReadUI({'ut1_1','ut1_2','it1'},'tmp/'); UI = ReadUI({'ut1_1','ut1_2','it1'},'tmp/');
u_f = (UI.FD{1}.val + UI.FD{2}.val)/2; %averaging voltages to fit current u_f = (UI.FD{1}.val + UI.FD{2}.val)/2; %averaging voltages to fit current
i_f = UI.FD{3}.val / partial; i_f = UI.FD{3}.val / partial;
Z = u_f./i_f;
delta_t = 1.35929e-11; % time-step (s) FIXME will change, if mesh is changed!
i_f2 = i_f .* exp(-1i*2*pi*UI.FD{1}.f*delta_t/2); % compensate half time-step advance of H-field
Z = u_f./i_f2;
plot(UI.FD{1}.f,real(Z),'Linewidth',2); plot(UI.FD{1}.f,real(Z),'Linewidth',2);
hold on; hold on;
grid on; grid on;