for i=1:nj T1b(i)=T2(1,ZK(1,1),ZK(1,2),i)+T2(2,ZK(1,1),ZK(1,2),i)+T2(3,ZK(1,1),ZK(1,2),i)+T2(4,ZK(1,1),ZK(1,2),i)+T0(i); T2b(i)=T2(1,ZK(2,1),ZK(2,2),i)+T2(2,ZK(2,1),ZK(2,2),i)+T2(3,ZK(2,1),ZK(2,2),i)+T2(4,ZK(2,1),ZK(2,2),i)+T0(i); T3b(i)=T2(1,ZK(3,1),ZK(3,2),i)+T2(2,ZK(3,1),ZK(3,2),i)+T2(3,ZK(3,1),ZK(3,2),i)+T2(4,ZK(3,1),ZK(3,2),i)+T0(i); T4b(i)=T2(1,ZK(4,1),ZK(4,2),i)+T2(2,ZK(4,1),ZK(4,2),i)+T2(3,ZK(4,1),ZK(4,2),i)+T2(4,ZK(4,1),ZK(4,2),i)+T0(i); end向量化代码

可以使用矩阵运算来向量化代码，具体实现如下：

T1b = sum(T2(1:4, ZK(1:4, 1), ZK(1:4, 2), :), 1) + T0;
T2b = sum(T2(1:4, ZK(5:8, 1), ZK(5:8, 2), :), 1) + T0;
T3b = sum(T2(1:4, ZK(9:12, 1), ZK(9:12, 2), :), 1) + T0;
T4b = sum(T2(1:4, ZK(13:16, 1), ZK(13:16, 2), :), 1) + T0;

其中，使用了 `sum` 函数对第1维进行求和，实现了对第1个for循环的向量化。