a=importdata（‘data.txt‘）;
Y=a.data（:2）;
Z=a.data（:3）;
BX=a.data（:4）;
plot3（YZBX）;
title（‘磁場(chǎng)分布‘）;
xlabel（‘Y‘）;ylabel（‘Z‘）;zlabel（‘BX‘）;
axis（[-8801602000]）;
e=8;
m0=9.1*10^（-31）;
m=（0.511+8）*m0/0.511;
v0=2.99792458*10^8*sqrt（1-（0.511/（0.511+e））^2）;
q=1.6021917*10^（-19）;
B0=BX.*10^（-4）;
B=reshape（B0321321）;
r=[];
s=[];
cas=0;
i=1;
j=1;
y=0;
z=0;
cta=0;
y0=0;
z0=0;
while?z<=0.16&&y>=-0.08
????B1=（B（ij）+B（i+1j）+B（ij+1）+B（i+1j+1））/4;
????R=m*v0/（B1*q）;
????if?cas==0
????datz=sqrt（R^2-（R*cos（cta）-y0）^2）-R*sin（cta）;
????if?datz<=0.0005
????????z=z+datz;
????????y=y-y0;
????????i=i+1;
????????cas=1;
????????z0=datz;
????????cta=cta+2*asin（sqrt（y0^2+datz^2）/（2*R））;
????else
????????daty=sqrt（R^2-（R*sin（cta）+0.0005）^2）-R*cos（cta）+y0;
????????z=z

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件??????38701??2010-11-30?09:50??電子軌跡模擬\電子軌跡.jpg

?????文件?????224220??2010-11-30?09:51??電子軌跡模擬\磁場(chǎng)分布.jpg

?????文件???????1730??2010-11-30?09:46??電子軌跡模擬\guiji.m

?????文件????8449437??2010-11-12?20:06??電子軌跡模擬\data.txt

?????目錄??????????0??2010-11-30?09:55??電子軌跡模擬

-----------?---------??----------?-----??----

??????????????8714088????????????????????5