1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22
| function [TL]=TLL(f,D)% f etant la frequence,D la longueur du cable
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
mu_0=4*pi*10.^7; %V.s/A.m
eps_0=8.854*10.^-12; %F/m
eps_r=4;
tang_delta_c=0.01;
r=15.6*1e-3;%m
v=3.6*1e-3;%m
ro=2.8*1e-8;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
L=mu_0*v/(2*r);
C=2*eps_0*eps_r*r/v;
Rr=sqrt(ro*pi*f*mu_0/(r.^2));
G=2*pi*f*C*tang_delta_c;
Z_0=sqrt((Rr+(1j*2*pi*f*L))/(G+1j*2*pi*f*C));
sigma=sqrt((Rr+(1j*2*pi*f*L))*(G+1j*2*pi*f*C));
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
TL(1,1)=cosh(sigma*D);
TL(1,2)=Z_0*sinh(sigma*D);
TL(2,1)=(1/Z_0)*sinh(sigma*D);
TL(2,2)=TL(1,1);
end |
Partager