Pt = 1;
ht = 50;
hr = 2;
f = 9e8;
c = 3e8;
lambda = c/f;
R = -1;
Gl = 1;
%Gr = 1;
Gr = [1,0.316,0.1,0.01];
for i=1:4
dc = 4*ht*hr/lambda;
d = 1:100000;
dBd = log10(d);
h1 = ht + hr;
h2 = ht - hr;
delta_dis = sqrt(h1^2 + d.^2) - sqrt(h2^2 + d.^2);
delta_phs = 2*pi*delta_dis/lambda;
Pr = Pt*(lambda/(4*pi))^2*(abs(sqrt(Gl)./sqrt(h2^2 + d.^2) + R*sqrt(Gr(i))*exp(-j*delta_phs)./sqrt(h1^2 + d.^2))).^2;
dBPr = 10*log10(Pr)-10*log10(Pr(1));
switch i
case 1
plot(dBd,dBPr,'b','LineWidth',2);
case 2
plot(dBd,dBPr,'r','LineWidth',2);
case 3
plot(dBd,dBPr,'k','LineWidth',2);
case 4
plot(dBd,dBPr,'m','LineWidth',2);
end
hold on
end
grid on
plot([log10(dc) log10(dc)],[-120 40],':r','LineWidth',2);
plot([log10(ht) log10(ht)],[-120 40],'--r','LineWidth',2);
legend('Gr = 1','Gr = 0.316','Gr = 0.1','Gr = 0.01','dc','ht')
评论