% 大城市
function pathloss = hata_big_city(distance, frequency, hte, hre)
		for k = 1 : length(distance),
   pathloss(k) = 0;
   L = 69.55 + 26.16*log10(frequency) - 13.82*log10(hte) + (44.9 - 6.55*log10(hte))*log10(distance(k));
			% 大城市修正因子
   if frequency >= 300
     a = 3.2*(log10(11.75*hre)*log10(11.75*hre)) - 4.97;
   elseif frequency <= 300
     a = 8.29*(log10(1.54*hre)*log10(1.54*hre)) - 1.1;
   end
   L = L - a;
   pathloss(k) = L;
	end