blob: a4a3a5e25308a47c014170c2a827900a8fb477d5 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
|
function eps_debye = CalcDebyeMaterial(f, eps_r, kappa, eps_Delta, t_relax)
% eps_debye = CalcDebyeMaterial(f, eps_r, kappa, eps_Delta, t_relax)
%
% Calculate the Debye type dispersive material constant
%
% arguments:
% f: frequeny range of interest
% eps_r: eps_r infinity
% kappa: conductivity (losses)
% eps_Delta: (vector) delta of relative permitivity
% t_relax: (vector) relaxation time (losses)
%
% return:
% eps_debye: the complex relative permitivity
%
% See also: CalcLorentzMaterial
%
% CSXCAD matlab interface
% -----------------------
% author: Thorsten Liebig (2013)
EPS0 = 8.85418781762e-12;
eps_debye = ones(size(f))*eps_r - 1j*kappa./(2*pi*f)/EPS0;
for n=1:numel(eps_Delta)
eps_debye = eps_debye + eps_Delta(n)./(1+2j*pi*f*t_relax(n));
end
end
|