MATLAB 指令解析：带电粒子运动模拟

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 的三维向量，代表电场强度的方向在 z 轴正方向上。

'B = [0, 1, 0]'：表示一个大小为 1 的三维向量，代表磁感应强度的方向在 y 轴正方向上。

'v0 = [1, 0, 0]'：表示一个大小为 1 的三维向量，代表带电粒子的初始速度方向在 x 轴正方向上。

'r0 = [0, 0, 0]'：表示一个大小为 1 的三维向量，代表带电粒子的初始位置在原点。

'theta = 30 * pi / 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。