1function g = loss_COST231_large_cities_urban(d, hb, hm, f) 2 3 %% 4 %% function g = loss_COST231_large_cities_urban(d, c) 5 %% 6 %% returns the loss at d meters for f frequency and mobile height m and 7 %% base station height of hb 8 9 assert(isscalar(f)); 10 assert(f > 0); 11 12 Fhm = 3.2*(log10(11.75*hm))^2; 13 C = 3; 14 g = zeros(size(d)); 15 g(find(d > 0)) = 46.3 + (33.9*log10(f)) - (13.82*log10(hb)) + (44.9-(6.55*log10(hb))).*log10(d) - Fhm + C; 16 17 g(find(d <= 0)) = 1; 18 19 20