n = 92 ;
Load('EHField_3D_5to15fs.mat') ;
XX (1) = 0 ;
YY(1)=y(40) ;
ZZ (1) = 332 * 10 ^ -9 ;
VX(1) = 0 ;
Vy(1) = 0 ;
VZ(1) = 10 ^ 6 ;
pour i = 2:n
Vitesse = vt(xx(i-1),yy(i-1),zz(i-1),vx(i-1),vy(i-1),vz(i-1),i-1) ;
VX(i) = speed(1,1) ;
Vy(i) = speed(1,2) ;
VZ(i) = speed(1,3) ;
position = xz(xx(i-1),yy(i-1),zz(i-1),vx(i-1),vy(i-1),vz(i-1),i-1) ;
mesure = position(1,1) ;
YY(i) = position(1,2) ;
ZZ(i) = position(1,3) ;
Plot3(XX,YY,ZZ)
Attends
fin
Et maintenant, je dois simuler de nombreux différents temps initial. J’avais écrit ce code
effacer tout ;
n = 92 ;
Load('EHField_3D_5to15fs.mat') ;
XX (1) = 0 ;
YY(1)=y(40) ;
ZZ (1) = 332 * 10 ^ -9 ;
VX(1) = 0 ;
Vy(1) = 0 ;
VZ(1) = 10 ^ 6 ;
pour a = 90
j = a: n
pour i = 2:n-2 a +
Vitesse = vt(xx(i-1),yy(i-1),zz(i-1),vx(i-1),vy(i-1),vz(i-1),j-1) ;
VX(i) = speed(1,1) ;
Vy(i) = speed(1,2) ;
VZ(i) = speed(1,3) ;
position = xz(xx(i-1),yy(i-1),zz(i-1),vx(i-1),vy(i-1),vz(i-1),j-1) ;
mesure = position(1,1) ;
YY(i) = position(1,2) ;
ZZ(i) = position(1,3) ;
Plot3(XX,YY,ZZ)
Attends
fin
fin
fin