MAtrice ABCD-Transmission Line

Bonjour,
j'ai besoin de quelqu'un pour me corriger un script de calcul de la matrice ABCD pour une ligne de transmission:

Code:

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

Merci