脚本中的函数定义必须出现在文件的结尾。 请将 "particle_acceleration" 函数定义后面的所有语句都移到第一个局部函数定义前面。

好的，以下是调整后的代码示例：

% 物理常数
q = 1.602e-19; % 电荷量
m = 9.109e-31; % 粒子质量
c = 299792458; % 光速

% 模拟参数
B = 1; % 磁场强度
E = 1; % 电场强度
theta = pi/4; % 粒子进入磁场的方向与磁场方向之间的夹角
v0 = 1e6; % 初始速度大小
tspan = [0, 1e-7]; % 时间范围
options = odeset('RelTol', 1e-6, 'AbsTol', 1e-9); % ODE求解器参数

function a = particle_acceleration(t, x, v, B, E, theta)
% 计算粒子加速度
% q = 1.602e-19; % 电荷量
% m = 9.109e-31; % 粒子质量
% c = 299792458; % 光速

% 计算洛伦兹力和电场力
F_B = q*cross(v, B*[0;0;1]);
F_E = q*E*[cos(theta); sin(theta); 0];

% 计算加速度
a = (F_B + F_E)/m;
end

% 情况1：电场强度和磁场强度都不为零
[T1, X1] = ode45(@(t, x) particle_acceleration(t, x, v, B*[0;0;1], E*[cos(theta); sin(theta); 0]), tspan, [0;0;0;v0*cos(theta);v0*sin(theta);0], options);

% 情况2：电场强度为零，磁场强度不为零
[T2, X2] = ode45(@(t, x) particle_acceleration(t, x, v, B*[0;0;1], [0;0;0]), tspan, [0;0;0;v0*cos(theta);v0*sin(theta);0], options);

% 情况3：电场强度不为零，磁场强度为零
[T3, X3] = ode45(@(t, x) particle_acceleration(t, x, v, [0;0;0], E*[cos(theta); sin(theta); 0]), tspan, [0;0;0;v0*cos(theta);v0*sin(theta);0], options);

% 绘制情况1的运动轨迹
figure;
plot3(X1(:,1), X1(:,2), X1(:,3));
xlabel('x');
ylabel('y');
zlabel('z');
title('情况1');

% 绘制情况2的运动轨迹
figure;
plot3(X2(:,1), X2(:,2), X2(:,3));
xlabel('x');
ylabel('y');
zlabel('z');
title('情况2');

% 绘制情况3的运动轨迹
figure;
plot3(X3(:,1), X3(:,2), X3(:,3));
xlabel('x');
ylabel('y');
zlabel('z');
title('情况3');

现在，"particle_acceleration" 函数定义已经出现在文件的开头，可以正确执行。