请解释下面每条matlab指令的含义：q = 16e-19; 带电粒子电荷量m = 91e-31; 带电粒子质量E = 001; 电场强度B = 010; 磁感应强度v0 = 100; 初始速度r0 = 000; 初始位置theta = 30pi180; 夹角tspan = 010; 时间范围f = tY Y46; qmE + crossY46B;; 带电粒子的运动方程opts

1.6e-19：表示科学计数法中的1.6乘以10的负19次方，即1.6乘以10的负19次方的数值，即元素电荷的大小。 9.1e-31：表示科学计数法中的9.1乘以10的负31次方，即9.1乘以10的负31次方的数值，即电子质量的大小。 E = [0,0,1]：表示一个大小为1的3维向量，代表电场强度的方向在z轴正方向上。 B = [0,1,0]：表示一个大小为1的3维向量，代表磁感应强度的方向在y轴正方向上。 v0 = [1,0,0]：表示一个大小为1的3维向量，代表带电粒子的初始速度方向在x轴正方向上。 r0 = [0,0,0]：表示一个大小为1的3维向量，代表带电粒子的初始位置在原点。 theta = 30pi/180：将角度制的30度转换为弧度制，结果为pi/6，即30度对应的弧度值。 tspan = [0,10]：表示一个时间范围，从0到10秒。 f = @(t,Y) [Y(4:6); q/m(E + cross(Y(4:6),B));]：定义了一个匿名函数f，其中t为时间变量，Y为带电粒子的状态向量，包括位置和速度。函数返回值是一个大小为6的列向量，前三个元素表示速度，后三个元素表示加速度，其中加速度的计算使用了电场和磁场的叉积。 opts = odeset('RelTol',1e-6,'AbsTol',1e-6)：定义了一个选项结构体，包括相对误差和绝对误差的容忍度，用于设置求解器的参数。 [T,Y] = ode45(f,tspan,[r0,v0],opts)：使用ode45求解微分方程，其中f为微分方程的右侧，tspan为时间范围，[r0,v0]为初始状态向量，opts为选项结构体。函数返回值包括时间向量T和状态向量Y。 plot3(Y(:,1),Y(:,2),Y(:,3))：绘制带电粒子的轨迹，其中Y(:,1)表示位置向量的x坐标，Y(:,2)表示位置向量的y坐标，Y(:,3)表示位置向量的z坐标。函数plot3绘制三维图像。 xlabel('x')：设置x轴标签为x。 ylabel('y')：设置y轴标签为y。 zlabel('z')：设置z轴标签为z