如何将for i=1:LL向量化计算

将`for i=1:LL`向量化计算需要使用MATLAB的向量化计算功能。具体步骤如下：

1. 将`aa`、`ab`、`ac`、`ad`等数组改为二维数组，例如`aa(i,j,k)`。

2. 使用`meshgrid`函数生成`x`和`t`的网格矩阵，例如`[X,T] = meshgrid(x,t)`。

3. 计算`sqrt(rrrr + (z - x).(z - x))`和`sqrt(rrrr + (z + x).(z + x))`，并保存为新的变量，例如`temp1=sqrt(rrrr + (z - X).*(z - X))`和`temp2=sqrt(rrrr + (z + X).*(z + X))`。

4. 将`aa`、`ab`、`ac`、`ad`等数组中的循环计算转化为向量计算。例如：

aa = integral(@(x)0.25*exp(v*Xt*0.5 / Rap).*exp(-v * temp1 / Rap).*erfc((temp1 - v * T)*0.5 / sqrt(Rap*T))./sqrt(temp1)/(2 * 3.1415926*Rlamd),a,b);

在这个版本的代码中，`temp1`和`T`都是二维数组，`erfc`函数也是向量化的。这样就可以避免使用`for`循环，而是使用向量化的方式进行计算。

5. 最后，将计算结果保存到`Tj`数组中，例如`Tj(:,:,k,ii,a1,b1)=reshape(sum(sum(T,1),2),[LL 1])`。在这个版本的代码中，使用了`sum`函数进行求和，并使用`reshape`函数将结果转换为二维数组。

通过向量化计算，可以避免使用`for`循环，使代码更加简洁高效。