用matlab将下列文字编写为一段完整的代码： 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.xlsx中的数据
Shuju = readmatrix('Shuju.xlsx'); % 读取Shuju.xlsx中的数据
xiongj = Shuju(:, 4); % 获取中心树的胸径数据

% 2. 读取zyshu.xlsx中的数据
zyshu = readmatrix('zyshu.xlsx');

% 3. 根据中心树编号，在Shuju.xlsx中读取对应的胸径数据并保存到xiongj.xlsx中
xiongj_xls = zeros(size(zyshu, 1), 1); % 创建一个和zyshu一样大小的数组
for i = 1:size(zyshu, 1)
idx = zyshu(i, 1); % 获取中心树的编号
xiongj_xls(i) = xiongj(idx); % 获取中心树的胸径数据
end
xlswrite('xiongj.xlsx', xiongj_xls, 'Sheet1', 'A1'); % 保存到xiongj.xlsx中的第1列

% 4. 根据zyshu.xlsx中的编号，在Shuju.xlsx中读取对应的胸径数据并保存到xiongj.xlsx中
for i = 2:5 % 遍历zyshu.xlsx中的第2到第5列
xiongj_xls = zeros(size(zyshu, 1), 1); % 创建一个和zyshu一样大小的数组
for j = 1:size(zyshu, 1)
idx = zyshu(j, i); % 获取编号
xiongj_xls(j) = xiongj(idx); % 获取胸径数据
end
xlswrite('xiongj.xlsx', xiongj_xls, 'Sheet1', char('A'+i-1)); % 保存到xiongj.xlsx中对应的列
end

% 5. 遍历xiongj.xlsx中每一行，并依次比较大小
daxiaobi = zeros(size(xiongj, 1), 2); % 创建一个和xiongj一样大小的数组
for i = 1:size(xiongj, 1)
daxiaobi(i, 1) = i; % 保存顺序编号
sum = 0; % 计算变量赋值和
for j = 2:5 % 比较第2到第5列的数据
if xiongj(i) < xiongj(j)
sum = sum + 1;
end
end
daxiaobi(i, 2) = sum / 4; % 保存平均值
end
xlswrite('daxiaobi.xlsx', daxiaobi, 'Sheet1', 'A1:B1'); % 保存到daxiaobi.xlsx中

% 6. 输出daxiaobi.xlsx中的数据
daxiaobi = readmatrix('daxiaobi.xlsx'); % 读取daxiaobi.xlsx中的数据
disp(daxiaobi); % 输出到命令行窗口