Bonjour,
J'ai besoin d'un programme en marlab qui permet l'analyse temporel des ligne de transmissions par FDTD, j'ai essayé mais mon programme ne me donne pas les bon résultat.
Pouvez vous m'aider SVP?
Bonjour,
J'ai besoin d'un programme en marlab qui permet l'analyse temporel des ligne de transmissions par FDTD, j'ai essayé mais mon programme ne me donne pas les bon résultat.
Pouvez vous m'aider SVP?
Bonjour,
montres-nous ton programme, on pourra peut-être trouver ce qui ne va pas.
Pour une bonne utilisation des balises code c'est ici!
Petit guide du voyageur MATLABien : Le forum La faq Les tutoriels Les sources
La nature est un livre écrit en langage mathématique. Galilée.
bonjour,
Voici mon code:
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 clear all close all clc L=0.25*10^-6; C=10^-10; Lg=400; vit=3*10^8; dz=10; dt=5*10^-8; Rs=50;%[50 0;0 50]; RL=50;%[50 0;0 50]; Nx=40; itmax=100; a1=1/(dz/dt*Rs*C+1); a2=dz/dt*Rs*C-1; a3=dt/(dz*L); a4=1/(dz/dt*RL*C+1); a5=dz/dt*RL*C-1; a6=dt/(dz*C); V(1:Nx+1)=0; %INITIALISATION I(1:Nx)=0; for it=2:itmax; if it<5 Vs(it)=7.5*(it-1); else if (it>=5 & it<=25) Vs(it)=30; else if (it>25 & it<30) Vs(it)=-7.5*(it-1)+217.5; else Vs(it)=0.0; end end end end for it=2:itmax V(1)=a1*(a2*V(1)-2*Rs*I(1)+Vs(it)+Vs(it-1)); V(Nx+1)=a4*(a5*V(Nx+1)+2*RL*I(Nx)); for k=2:Nx V(k)=V(k)-a6*(I(k)-I(k-1)); end for k=1:Nx I(k)=I(k)-a3*(V(k+1)-V(k)); end end subplot(2,1,1),plot(Vs),title('Vs') subplot(2,1,2),plot(1:Nx+1,V),title('V')
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 clear all; %%%%%%%%les données%% l=1; deltx=0.005; ndx=l/deltx; itmax=250; Rs=50; Rl=50; L=0.25*10^-6; C=10^-10; v=1/sqrt(L*C); deltt=deltx/(v); a=((deltx/deltt)*Rs*C+1)^-1; b=((deltx/deltt)*Rs*C-1); d=((deltx/deltt)*Rl*C+1)^-1; m=((deltx/deltt)*Rl*C-1); %condition initial V(1:ndx+1)=0; I(1:ndx)=0; V1(1)=0; Vcharge(1)=0; %sourse for t=1:itmax if (t*deltt<0.1*10^-9) Vs(t)=10^10*(t-1)*deltt; else Vs(t)=1; end end for t=2:itmax V(1)=a*[b*V(1)-2*Rs*I(1)+(Vs(t)+Vs(t-1))]; V(ndx+1)=d*[m*V(ndx+1)+2*Rl*I(ndx)]; V1(t)=V(1) Vcharge(t)=V(ndx+1) for k=2:ndx V(k)=[V(k)-((deltt/(C*deltx))*(I(k)-I(k-1)))]; end for k=1:ndx I(k)=[I(k)-((deltt/(L*deltx))*(V(k+1)-V(k)))]; end end %affichages des résultats 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');
Vous avez un bloqueur de publicités installé.
Le Club Developpez.com n'affiche que des publicités IT, discrètes et non intrusives.
Afin que nous puissions continuer à vous fournir gratuitement du contenu de qualité, merci de nous soutenir en désactivant votre bloqueur de publicités sur Developpez.com.
Partager