无源定位 matlab,关于无源定位TDOA算法GDOP分析的问题

有一段程序，不是我编写的，用MATLAB仿真总是出现 下面这句话，我知道是矩阵下标不对了，可是怎么改呢，求指教

”Subscripted assignment dimension mismatch.

Error in Untitled3 (line 67)

miu(k,l) = cov(delta_r(k,i,j),delta_r(l,i,j))/(dd_r(k)*dd_r(l));“

代码如下：

L = 10; % 站间距

delta_tdoa = 1e-8; % 时差测量误差

delta_s = 0; % 站址测量误差

c = 3e5; % 光速(公里)

alpha = 30; % 站间夹角

d_r = delta_tdoa*c; % 距离差误差

% 计算站址

% s0

x0 = 0;

y0 = 0;

z0 = 0;

% s1

x1 = L*cos(30*pi/180);

y1 = L*sin(30*pi/180);

z1 = 0.1;

% s2

x2 = L*cos(150*pi/180);

y2 = L*sin(150*pi/180);

z2 = 0.2;

% s3

x3 = 0;

y3 = -L;

z3 = 0.3;

% he

xd = [x0,x1,x2,x3];

yd = [y0,y1,y2,y3];