Discussion :

[Lotfi slimani]

salut,
s'il vous plait,j'ai besoin d'un programme MATLAB qui traite une ligne de transmission multifilaire(au moins 3 fils) avec la méthode FDTD.

j'ai fais un programme mais,il ne fonctionne pas,je sais pas ou l'erreur

s'il vous plait aidez moi

[FLB]

[Discussion liée : http://www.developpez.net/forums/d10...missions-fdtd/]
Bonjour,
Afin que l'on puisse t'aider montre nous ce que tu as fait - même si ça ne marche pas, et le cas échéant montre nous le message d'erreur complet que t'envoie Matlab.

[Lotfi slimani]

voila mon programme ,c'est le cas sans perte
je veux le cas avec et sans perte si possible
merci d'avance

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
%"""""""""""""""" FDTD POUR UNE LIGNE DE TRANSMISSION """"""""""""""""""
%-------------------------------------------------------
%"""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""
clear all;
  %DONNEES
l=0.1;
deltx=0.005;
ndx=l/deltx;
itmax=250;
Rs=50;
Rl=50;
L(1:ndx,1:ndx)=0.25*10^-6;
C(1:ndx,1:ndx)=10^-10;
v=inv(sqrt(L*C));
U(1:ndx,1:ndx)=1;
deltt=deltx*inv(v);
a=inv(deltx*inv(deltt)*Rs*C+U);
b=(deltx*inv(deltt)*Rs*C-U);
d=inv(deltx*inv(deltt)*Rl*C+U);
m=(deltx*inv(deltt)*Rl*C-U);
%COND INIT
V(1:ndx+1,1:ndx+1)=0;
I(1:ndx+1,1:ndx+1)=0;
V1(1)=0;
Vcharge(1)=0;
%SRC
for t=1:itmax
    T(1:ndx,1:ndx)=0.1*10^-9;
    Vs=0;
    if (t*deltt<T)
       Vs=(t-1)*deltt*10^10;
 
        else
            Vs=U;
    end
end
for t=2:itmax
 
   V(1,1)=a(1,1)*[b(1,1)*V(1,1)-2*Rs*I(1,1)+(Vs(t)+Vs(t-1))];
    V(ndx+1,ndx+1)=d(ndx,ndx)*[m(ndx,ndx)*V(ndx+1,ndx+1)+2*Rl*I(ndx,ndx)];
    V1(t)=V(1,1);
    Vcharge(t)=V(ndx+1,ndx+1);
    %COMMANDE DE TENSION ET COURANT
    for k=2:ndx
        for h=2:ndx
        V(k,h)=[V(k,h)-((deltt*inv(deltx*C))*(I(k,h)-I(k-1,h-1)))];   
        end
    end
    for k=1:ndx
        for h=1:ndx
        I(k,h)=[I(k,h)-((deltt/deltx)*inv(L(k,h)))*(V(k+1,h+1)-V(k,h))];
        end
    end
end
 
%AFFICHAGE DES RESULTATS
subplot(311)
plot(10^-6:itmax,Vs,'b')
title('source');
xlabel('t en  ns');
ylabel('s(t)');
subplot(312)
plot(1:itmax,V1,'g')
title('La tension');
xlabel('ttt');
ylabel('V1');
subplot(313)
plot(1:itmax,Vcharge,'m')
title('la tension');
xlabel('t');
ylabel('Vcharge');