问在for循环中进行了一些迭代之后，Matlab的运行速度非常慢(几乎停滞)。

EN

ro1 = (180*3600)/pi;
ro = ro1^2;  
sigma_h = horizontal_prec_secs;
sigma_v = vertical_prec_secs;
error_ellip_dim = nan(5823,8153,3);
for i = 1:5823
    for j = 1:8153
        sigma_r = range_precision(i,j);
        h = horizontal_angle(i,j);
        v = vertical_angle(i,j);
        r = range(i,j)*1000;
        if (isnan(h))||(isnan(sigma_r))
           error_ellip_dim(i,j,:) = NaN;
        else
            EXX(1,1) = sigma_r^2*COSDH(i,j)^2*COSDV(i,j)^2 + (r^2*sigma_h^2*COSDV(i,j)^2*SINDH(i,j)^2)/ro + (r^2*sigma_v^2*COSDH(i,j)^2*SINDV(i,j)^2)/ro;
            EXX(1,2) = sigma_r^2*COSDH(i,j)*COSDV(i,j)^2*SINDH(i,j) - (r^2*sigma_h^2*COSDH(i,j)*COSDV(i,j)^2*SINDH(i,j))/ro + (r^2*sigma_v^2*COSDH(i,j)*SINDH(i,j)*SINDV(i,j)^2)/ro;
            EXX(1,3) = COSDH(i,j)*COSDV(i,j)*SINDV(i,j)*sigma_r^2 - (r^2*sigma_v^2*COSDH(i,j)*COSDV(i,j)*SINDV(i,j))/ro;
            EXX(2,1) = EXX(1,2); 
            EXX(2,2) = sigma_r^2*COSDV(i,j)^2*SINDH(i,j)^2 + (r^2*sigma_h^2*COSDH(i,j)^2*COSDV(i,j)^2)/ro + (r^2*sigma_v^2*SINDH(i,j)^2*SINDV(i,j)^2)/ro;
            EXX(2,3) = COSDV(i,j)*SINDH(i,j)*SINDV(i,j)*sigma_r^2 - (r^2*sigma_v^2*COSDV(i,j)*SINDH(i,j)*SINDV(i,j))/ro;
            EXX(3,1) = EXX(1,3);
            EXX(3,2) = EXX(2,3);
            EXX(3,3) = sigma_r^2*SINDV(i,j)^2 + (r^2*sigma_v^2*COSDV(i,j)^2)/ro;
            [eig_vec_mat,eig_val_mat] = eig(EXX);
            ellip_params = sqrt(diag(eig_val_mat));
            error_ellip_dim(i,j,:) = ellip_params';            
        end
    end
    i
    j
end