请给我一个MATLAB的轴向切深矩阵程序

以下是一个简单的 MATLAB 轴向切深矩阵程序：

% 定义铣削参数
vc = 60; % 铣削速度，单位：m/min
ae = 0.15; % 切削宽度，单位：mm
fz = 0.25; % 进给量，单位：mm/z

% 定义材料参数
kc = 430; % 切削力系数，单位：N/mm^2

% 定义初始温度和刀具半径
T0 = 20; % 初始温度，单位：℃
r = 5; % 刀具半径，单位：mm

% 定义轴向切削深度矩阵
z_matrix = [0.1, 0.2, 0.3; 0.4, 0.5, 0.6; 0.7, 0.8, 0.9]; % 轴向切削深度矩阵，单位：mm

% 计算切削速度和进给速度
vs = vc * pi * r / 1000; % 切削速度，单位：m/s
f = fz * vs; % 进给速度，单位：mm/s

% 计算切削力和切削温度
kc_1 = kc * z_matrix / ae; % 实际切削力系数
F = kc_1 * ae * vs; % 切削力，单位：N
T = T0 + F.^2 / (2 * kc_1 * ae) + kc_1 * ae * vs^2 / (2 * pi * r * kc); % 切削温度，单位：℃

% 显示轴向切削深度矩阵和切削温度矩阵
disp('轴向切削深度矩阵：');
disp(z_matrix);
disp('切削温度矩阵：');
disp(T);

该程序计算了在不同轴向切削深度下的切削力和切削温度，并展示了轴向切削深度矩阵和切削温度矩阵。请注意，这只是一个简单的示例程序，实际的铣削过程模拟需要更复杂的参数和计算。此外，该程序也没有考虑到刀具磨损和切削力的动态变化等因素。