用matlab 2022将下列步骤编写为一段完整的代码： 1、存在一个名称为Shuju的EXCEL表，表中第1列是中心树的顺序编号，第2列是中心树的x坐标，第3列是中心树的y坐标，第4列是中心树的胸径数据，第5列是中心树的种类。 2、存在一个名称为zyshu的excel表,表中第1列是中心树的顺序编号，第2列是离中心树最近的第1棵树的编号，第3列是离中心树最近的第2棵树的编号，第4列是离中心树最近的第3棵树的编号，第5列是离中心树最近的第4棵树的编号。 3、根据中心树编号，在名称为Shuju的EXCEL表中，读取该编号对应的中心树的胸径数据，保存在名称为xiongj的excel表的第1列中。 4、根据名称为zyshu的excel表中，第2、3、4、5列中的编号，依次从该编号对应的名称为Shuju的EXCEL表中读取出该编号对应的胸径数据，保存在名称为xiongj的excel表的第2、3、4、5列中。 5、在名称为xiongj的excel表中，遍历每一行，并顺序编号，将编号依次保存在名称为daxiaobi的excel表的第1列中。依次比较名称为xiongj的excel表中每1行的第1列与第2、3、4、5列的数据值大小，当第1列的数值小于其他列的数据值时，变量赋值为1，否则赋值为0。比较完一行后，取变量赋值和的平均值，将值保存在名称为daxiaobi的excel表中每1行的第2列中。 6、保存名称为daxiaobi的excel表，运行后在命令行窗口输出名称为daxiaobi的excel表中数据。

以下是代码实现：

% 1. 读取名称为Shuju的EXCEL表
shuju = readtable('Shuju.xlsx');

% 2. 读取名称为zyshu的excel表
zyshu = readtable('zyshu.xlsx');

% 3. 根据中心树编号，读取胸径数据保存在名称为xiongj的excel表的第1列中
xiongj = table(shuju{:, 4}, 'VariableNames', {'xiongj'});

% 4. 依次读取第2、3、4、5列中的编号，读取对应的胸径数据保存在名称为xiongj的excel表的第2、3、4、5列中
for i = 2:5
col_name = strcat('tree', num2str(i-1));
xiongj.(col_name) = shuju{zyshu{:, i}, 4};
end

% 5. 遍历每一行、顺序编号，并将编号依次保存在名称为daxiaobi的excel表的第1列中
daxiaobi = table((1:size(xiongj, 1))', 'VariableNames', {'id'});

% 比较每一行的胸径数据大小，计算变量赋值和的平均值
for i = 1:size(xiongj, 1)
row_val = xiongj{i, 2:end};
n = length(row_val);
sum_val = sum(xiongj{i, 1} < row_val);
daxiaobi{i, 'daxiaobi'} = sum_val / n;
end

% 6. 保存daxiaobi表格
writetable(daxiaobi, 'daxiaobi.xlsx');

% 在命令行窗口输出daxiaobi表格
disp(daxiaobi);