答案:```matlab
% 定义时间跨度和初始条件
tspan = [0 10];
y0 = [0; 0; 10]; % 初始位置、速度和抛射角度
% 定义微分方程
function dydt = projectile(t, y)
g = 9.81; % 重力加速度
dydt = [y(2); -g - y(1)^2]; % 速度和加速度
end
% 使用ode45求解微分方程
[t, y] = ode45(@projectile, tspan, y0);
% 绘制位置随时间变化的图像
plot(t, y(:,1));
xlabel('Time (s)');
ylabel('Position (m)');
title('Projectile Motion with Air Resistance');
grid on;
```