给出上述问题的完整matlab代码

完整MATLAB代码如下：

% 建立数学模型
h = 300; % 飞行高度，单位：m
v = 300/3.6; % 飞行速度，单位：m/s
r = 0.2; % 物资半径，单位：m
m = 50; % 物资重量，单位：kg
g = 9.8; % 重力加速度，单位：m/s^2
k = 1/2*(1.2)*pi*(r^2); % 空气阻力系数，单位：N*s^2/m^4

% 风速及夹角设置
v_w = 5; % 风速，单位：m/s
theta = [0, 180, 90]; % 夹角，单位：度

% 计算投放距离
for i = 1:length(theta)
t = sqrt((2*h+r)/(g/2));
v1 = sqrt((m*g)/(4/3*pi*r^3*k)) * (1 - exp(-2*k*t/m));
x = v*t;
d = x - integral(@(t) v1, 0, t);
if theta(i) == 0
d_same = d - v_w*t;
elseif theta(i) == 180
d_opposite = d + v_w*t;
elseif theta(i) == 90
d_perpendicular = d;
end
end

% 输出结果
fprintf("无人机飞行方向与风向相同（夹角为0度）的投放距离为%.2f米\n", d_same);
fprintf("无人机飞行方向与风向相反（夹角为180度）的投放距离为%.2f米\n", d_opposite);
fprintf("无人机飞行方向与风向垂直（夹角为90度）的投放距离为%.2f米\n", d_perpendicular);