clear all;
hb=200;
hm=3;
for d=[1 2 5 10 30 60 100]
f1=100:0.1:300;
f2=300:0.1:3000;
lb11=69.55+26.16*log10(f1)-13.82*log10(hb)-(8.29*(log10(1.54*hm).^2)-1.1)+((44.9-6.55*log10(hb))*log10(d));
lb12=69.55+26.16*log10(f2)-13.82*log10(hb)-(3.2*(log10(11.75*hm).^2)-4.97)+((44.9-6.55*log10(hb))*log10(d));
lb21=lb11-2*(log10(f1/28)).^2-5.4;
lb22=lb12-2*(log10(f2/28)).^2-5.4;
lb31=lb11-4.78*(log10(f1)).^2+18.33*log(f1)-40.98;
lb32=lb12-4.78*(log10(f2)).^2+18.33*log(f2)-40.98;
f=[f1 f2];
lb1=[lb11 lb12];
lb2=[lb21 lb22];
lb3=[lb31 lb32];
figure(1);
hold on;
plot(f,lb1,'r');
title('大城市');
xlabel('频率/MHz');
ylabel('损耗中值/dB');
grid;
figure(2);
hold on;
plot(f,lb2,'b');
title('郊区');
xlabel('频率/MHz');
ylabel('损耗中值/dB');
grid;
figure(3);
hold on;
plot(f,lb3,'g');
title('开阔区');
xlabel('频率/MHz');
ylabel('损耗中值/dB');
grid;
figure(4);
title('中、小城市' );
xlabel('频率/MHz');
ylabel('损耗中值/dB');
评论0