At the first glance, you can replace
numerator=0;
for k = 1:n_trc
numerator = numerator + (H(k,i) - H_avg)*(P(k,j) - P_avg);
end
by
numerator = (H(1:n_trc,i) - H_avg)'*(P(1:n_trc,j) - P_avg);
as well as
denom_H = 0;
denom_P = 0;
for k = 1:n_trc
denom_H = denom_H + (H(k,i) - H_avg)^2;
denom_P = denom_P + (P(k,j) - P_avg)^2;
end
by
denom_H = (H(1:n_trc,i) - H_avg)'*(H(1:n_trc,i) - H_avg);
denom_P = (P(1:n_trc,j) - P_avg)'*(P(1:n_trc,j) - P_avg);
You can use ":" instead of "1:n_trc" if the corresponding array sizes exactly equal to n_trc.
You should also preallocate "raw" before the loop in order to avoid its sequential resizing: raw=zeros(...,...).
H_avg and P_avg are also recommended to compute before the loop as vectors, otherwise you compute the same values several times repeatedly.
Probably, fufther optimization is also possible.
Not all variables have clear sense (H_C_sum, P_C_sum, etc).
