matlab电子轨迹模拟程序

自己用matlab编写的电子在磁场中的轨迹模拟程序，利用了网格化的思想，给需要的同学以参考。

a=importdata（‘data.txt‘）;
Y=a.data（:2）;
Z=a.data（:3）;
BX=a.data（:4）;
plot3（YZBX）;
title（‘磁场分布‘）;
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

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      38701  2010-11-30 09:50  电子轨迹模拟\电子轨迹.jpg

     文件     224220  2010-11-30 09:51  电子轨迹模拟\磁场分布.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