m = 0.11 # 冰雹质量 g = 9.8 # 重力加速度 k = 0.1 # 阻力与速度成正比时的比例常数 c = 0.01 # 阻力与速度及其二次方的线性组合时的比例常数 h = 0.01 # 时间步长 v = 0 # 初始速度为0 # 欧拉法迭代 for i in range(int(20/h)): v = v + h*(m*g - k*v - c*v**2)/m # 绘制图像 import matplotlib.pyplot as plt import numpy as np t = np.linspace(0, 20, int(20/h)+1) plt.plot(t, v_list) plt.show()将这段代码妆化为matlab程序

m = 0.11; % 冰雹质量
g = 9.8; % 重力加速度
k = 0.1; % 阻力与速度成正比时的比例常数
c = 0.01; % 阻力与速度及其二次方的线性组合时的比例常数
h = 0.01; % 时间步长
v = 0; % 初始速度为0

% 欧拉法迭代
for i = 1:(20/h)
v = v + h*(m*g - k*v - c*v^2)/m;
end

% 绘制图像
t = linspace(0, 20, (20/h)+1);
plot(t, v_list);
xlabel('时间');
ylabel('速度');
title('速度-时间图像');