Probleme fundamentale privind propagarea VHF si UHF

Probleme fundamentale privind propagarea VHF si UHF

Propagarea in spatiul liber:

Se vor realiza simularile programelor si se traseaza graficul pierderilor functie de distanta, luand frecventa f ca parametru si respectiv graficul pierderilor de propagare functie de frecventa, luand distanta ca parametru .

%Lb - Pierderile de propagare in spatiul liber si deasupra suprafetelor

plane, functie de distanta

cu frecventa parametru (reprezentare liniara pe axa Ox)

pause on

h=figure;

%scrsz = get(0,'ScreenSize');

%set(h,'Position',[scrsz(1)*2 scrsz(2)*2 scrsz(3)* . 9 scrsz(4)* . 9]);

GT=1; GR=1; hT=10; hR=1;

c=3*(10^8);

d= . 3: . 2:160;

f=100; %MHz

color=['b' 'g' 'r' 'c' 'm'

Lb=-32 . 44-20*log10(f)-20*log10(d/1000);

%Lbr=10*log10(GT) + 10*log10(GR) + 20*log10(hT) + 20*log10(hR)-40*log10(d);

Lbr=real(10*log10(4*GR*GT) + 20*log10(c)-20*log10(4*pi*f*(10^6)*d) + 20*log10(sin(2*pi*hT*hR*f*(10^6) . /(c*d))));

h1=plot (d,Lb,'LineWidth'

title('Reprezentare liniara pe axa Ox'

xlabel('d[m]'

ylabel('LB[dB]'

set(gcf, 'name' 'Variatia pierderilor de propagare in cazul reflexiei pe suprafete plane, comparativ cu propagarea in spatiul liber functie de distanta la frecventa de 100MHz' 'numbertitle' 'off'

set(h1,'Color',color(1));

grid; hold on

h2=plot (d,Lbr,'LineWidth'

set(h2,'Color',color(3));

%grid; hold on

legend('Propagare in spatiu liber' 'Propagare deasupra suprafetelor reflectante plane'

uiwait(h);

%Lb - Pierderile de propagare in spatiul liber si deasupra suprafetelor

plane, functie de distanta

cu frecventa parametru (reprezentare logaritmica pe axa Ox)

pause on

h=figure;

%scrsz = get(0,'ScreenSize');

%set(h,'Position',[scrsz(1)*2 scrsz(2)*2 scrsz(3)* . 9 scrsz(4)* . 9]);

GT=1; GR=1; hT=10; hR=1;

c=3*(10^8); %m/s

d= . 1: . 01:1000;

f=100; %MHz

color=['b' 'g' 'r' 'c' 'm'

Lb=-32 . 44-20*log10(f)-20*log10(d/1000);

%Lbr=10*log10(GT) + 10*log10(GR) + 20*log10(hT) + 20*log10(hR)-40*log10(d);

Lbr=real(10*log10(4*GR*GT) + 20*log10(c)-20*log10(4*pi*f*(10^6)*d) + 20*log10(sin((2*pi*hT*hR*f*(10^6)) . /(c*d))));

h1=semilogx(d,Lb,'LineWidth'

title('Reprezentare logaritmica pe axa Ox'

xlabel('d[m]'

ylabel('LB[dB]'

set(gcf, 'name' 'Variatia pierderilor de propagare in cazul reflexiei pe suprafete plane, comparativ cu propagarea in spatiul liber functie de distanta la frecventa de 100MHz' 'numbertitle' 'off'

set(h1,'Color',color(1));

grid; hold on

h2=semilogx(d,Lbr,'LineWidth'

set(h2,'Color',color(3));

%grid; hold on

legend('Propagare in spatiu liber' 'Propagare deasupra suprafetelor reflectante plane'

uiwait(h);

%Lb - Pierderile de propagare in proximitatea emitatorului functie

de distanta in cazul reflexiei pe suprafete plane

cu frecventa parametru (reprezentare liniara pe axa Ox)

pause on

h=figure;

%scrsz = get(0,'ScreenSize');

%set(h,'Position',[scrsz(1)*2 scrsz(2)*2 scrsz(3)* . 9 scrsz(4)* . 9]);

GT=1; GR=1; hT=10; hR=1;

c=3*(10^8); %m/s

d=0 . 1: . 1:25;

f=[1 4 12 100 1000];

f = f*10^6;

color=['b' 'g' 'r' 'c' 'm'

for y=1:5

Lbr=real(10*log10(4*GR*GT) + 20*log10(c)-20*log10(4*pi*f(y)*d) + 20*log10(sin((2*pi*hT*hR*f(y)) . /(c*d))));

h1=plot(d,Lbr,'LineWidth'

title('Reprezentare liniara pe axa Ox'

xlabel('d[m]'

ylabel('LB[dB]'

set(gcf, 'name' 'Variatia pierderilor de propagare in proximitatea emitatorului functie de distanta in cazul reflexiei pe suprafete plane cu frecventa parametru' 'numbertitle' 'off'

set(h1,'Color',color(y));

grid; hold on

end

legend('1MHz' '4MHz' '12MHz' '100MHz' '1GHz'

uiwait(h);

%Lb - Pierderile de propagare in proximitatea emitatorului functie

de distanta in cazul reflexiei pe suprafete plane

cu frecventa parametru (reprezentare logaritmica pe axa Ox)

pause on

h=figure;

%scrsz = get(0,'ScreenSize');

%set(h,'Position',[scrsz(1)*2 scrsz(2)*2 scrsz(3)* . 9 scrsz(4)* . 9]);

GT=1; GR=1; hT=10; hR=1;

c=3*(10^8); %m/s

d=0 . 5: . 01:25;

f=[1 4 12 100 1000];

f = f*10^6;

color=['b' 'g' 'r' 'c' 'm'

for y=1:5

Lbr=real(10*log10(4*GR*GT) + 20*log10(c)-20*log10(4*pi*f(y)*d) + 20*log10(sin((2*pi*hT*hR*f(y)) . /(c*d))));

h1=semilogx(d,Lbr,'LineWidth'

title('Reprezentare logaritmica pe axa Ox'

xlabel('d[m]'

ylabel('LB[dB]'

set(gcf, 'name' 'Variatia pierderilor de propagare in proximitatea emitatorului functie de distanta in cazul reflexiei pe suprafete plane cu frecventa parametru' 'numbertitle' 'off'

set(h1,'Color',color(y));

grid; hold on

end

legend('1MHz' '4MHz' '12MHz' '100MHz' '1GHz'

uiwait(h);

%Lb - Pierderile de propagare in proximitatea emitatorului functie

de distanta in cazul reflexiei pe suprafete plane

cu frecventa parametru (reprezentare liniara pe axa Ox)

pause on

h=figure;

%scrsz = get(0,'ScreenSize');

%set(h,'Position',[scrsz(1)*2 scrsz(2)*2 scrsz(3)* . 9 scrsz(4)* . 9]);

GT=1; GR=1; hT=10; hR=1;

c=3*(10^8); %m/s

d=0 . 1: . 1:5000;

f=[1 4 12 100 1000];

f = f*10^6;

color=['b' 'g' 'r' 'c' 'm'

for y=1:5

Lbr=real(10*log10(4*GR*GT) + 20*log10(c)-20*log10(4*pi*f(y)*d) + 20*log10(sin((2*pi*hT*hR*f(y)) . /(c*d))));

h1=plot(d,Lbr,'LineWidth'

title('Reprezentare liniara pe axa Ox'

xlabel('d[m]'

ylabel('LB[dB]'

set(gcf, 'name' 'Variatia pierderilor de propagare in proximitatea emitatorului functie de distanta in cazul reflexiei pe suprafete plane cu frecventa parametru' 'numbertitle' 'off'

set(h1,'Color',color(y));

grid; hold on

end

legend('1MHz' '4MHz' '12MHz' '100MHz' '1GHz'

uiwait(h);

%Lb - Pierderile de propagare in proximitatea emitatorului functie

de distanta in cazul reflexiei pe suprafete plane

cu frecventa parametru (reprezentare logaritmica pe axa Ox)

pause on

h=figure;

%scrsz = get(0,'ScreenSize');

%set(h,'Position',[scrsz(1)*2 scrsz(2)*2 scrsz(3)* . 9 scrsz(4)* . 9]);

GT=1; GR=1; hT=10; hR=1;

c=3*(10^8); %m/s

d=1: . 01:2000;

f=[1 4 12 100 1000];

f = f*10^6;

color=['b' 'g' 'r' 'c' 'm'

for y=1:5

Lbr=real(10*log10(4*GR*GT) + 20*log10(c)-20*log10(4*pi*f(y)*d) + 20*log10(sin((2*pi*hT*hR*f(y)) . /(c*d))));

h1=semilogx(d,Lbr,'LineWidth'

title('Reprezentare logaritmica pe axa Ox'

xlabel('d[m]'

ylabel('LB[dB]'

set(gcf, 'name' 'Variatia pierderilor de propagare in proximitatea emitatorului functie de distanta in cazul reflexiei pe suprafete plane cu frecventa parametru' 'numbertitle' 'off'

set(h1,'Color',color(y));

grid; hold on

end

legend('1MHz' '4MHz' '12MHz' '100MHz' '1GHz'

uiwait(h);

Variatia pierderilor de propagare in cazul reflexiei pe suprafete plane, comparativ cu propagarea in spatiul liber functie de distanta la frecventa de 100 MHz .

   

Pierderile de propagare in proximitatea emitatorului functie de distanta in cazul reflexiei pe suprafete plane cu frecventa parametru .

   

Pierderile de propagare in proximitatea emitatorului functie de distanta in cazul reflexiei pe suprafete plane cu frecventa parametru .

   

Coeficientul de reflexie al pamantului:

%Roh - Variatia modulului coeficientului de reflexie al pamantului

in cazul polarizarii orizontale

pause on

h=figure;

epsr=15;

eps0=8 . 85*10^(-12);

sigm=12*10^(-3);

psi=0:0 . 001:(pi/2);

f=[1 4 12 100 1000];

f = f*10^6;

color=['b' 'g' 'r' 'c' 'm'

for y=1:5,

hi=sigm/(2*pi*f(y)*eps0);

Roh1=abs( (sin(psi)-sqrt(((epsr-j*hi)-(cos(psi)) . ^2))));

Roh2=abs( (sin(psi) + sqrt((epsr-j*hi-((cos(psi)) . ^2))) ));

Roh=(Roh1 . /Roh2);

title('Polarizare orizontala'

xlabel('Psi[rad]'

ylabel('|rhoH|'

%subplot(2,1,1);

h1=plot (psi,Roh,'LineWidth'

set(h1,'Color',color(y));

set(gcf, 'name' 'Variatia modulului coeficientului de reflexie al pamantului' 'numbertitle' 'off'

grid; hold on

end

legend('1MHz' '4MHz' '12MHz' '100MHz' '1GHz'

uiwait(h);

%Rov - Variatia modulului coeficientului de reflexie al pamantului

in cazul polarizarii verticale

pause on

h=figure;

epsr=15;

eps0=8 . 85*10^(-12);

sigm=12*10^(-3);

psi=0:0 . 001:(pi/2);

f=[1 4 12 100 1000];

f = f*10^6;

color=['b' 'g' 'r' 'c' 'm'

for y=1:5,

hi=sigm/(2*pi*f(y)*eps0);

Rov1=abs( (epsr-j*hi)*sin(psi)-sqrt((epsr-j*hi)-cos(psi) . ^2));

Rov2=abs( (epsr-j*hi)*sin(psi) + sqrt((epsr-j*hi)-cos(psi) . ^2));

Rov=(Rov1 . /Rov2);   

title('Polarizare verticala'

xlabel('Psi[rad]'

ylabel('|rhoV|'

%subplot(2,1,1);

h1=plot (psi,Rov,'LineWidth'

set(h1,'Color',color(y));

set(gcf, 'name' 'Variatia modulului coeficientului de reflexie al pamantului' 'numbertitle' 'off'

grid; hold on

end

legend('1MHz' '4MHz' '12MHz' '100MHz' '1GHz'

uiwait(h);

%Roh - Variatia fazei coeficientului de reflexie al pamantului

in cazul polarizarii orizontale

pause on

h=figure;

epsr=15;

eps0=8 . 85*10^(-12);

sigm=12*10^(-3);

psi=0:0 . 001:(pi/2);

f=[1 4 12 100 1000];

f = f*10^6;

color=['b' 'g' 'r' 'c' 'm'

for y=1:5,

hi=sigm/(2*pi*f(y)*eps0);

Roh1=( (sin(psi)-sqrt(((epsr-j*hi)-(cos(psi)) . ^2))));

Roh2=( (sin(psi) + sqrt((epsr-j*hi-((cos(psi)) . ^2))) ));

Roh=angle(Roh1 . /Roh2);

title('Polarizare orizontala'

xlabel('Psi[rad]'

ylabel('arg(rhoH)'

%subplot(2,1,1);

h1=plot (psi,Roh,'LineWidth'

set(h1,'Color',color(y));

set(gcf, 'name' 'Variatia fazei coeficientului de reflexie al pamantului' 'numbertitle' 'off'

grid; hold on

end

legend('1MHz' '4MHz' '12MHz' '100MHz' '1GHz'

uiwait(h);

%Rov - Variatia fazei coeficientului de reflexie al pamantului

in cazul polarizarii verticale

pause on

h=figure;

epsr=15;

eps0=8 . 85*10^(-12);

sigm=12*10^(-3);

psi=0 . 001:0 . 001:(pi/2); % nu se porneste din origine!

f=[1 4 12 100 1000];

f = f*10^6;

color=['b' 'g' 'r' 'c' 'm'

for y=1:5,

hi=sigm/(2*pi*f(y)*eps0);

Rov1= (epsr-j*hi)*sin(psi)-sqrt((epsr-j*hi)-cos(psi) . ^2);

Rov2= (epsr-j*hi)*sin(psi) + sqrt((epsr-j*hi)-cos(psi) . ^2);

Rov= angle(Rov1 . /Rov2);

title('Polarizare verticala'

xlabel('Psi[rad]'

ylabel('arg(rhoV)'

%subplot(2,1,1);

h1=plot (psi,Rov,'LineWidth'

set(h1,'Color',color(y));

set(gcf, 'name' 'Variatia fazei coeficientului de reflexie al pamantului' 'numbertitle' 'off'

grid; hold on

end

legend('1MHz' '4MHz' '12MHz' '100MHz' '1GHz'

uiwait(h);

Variatia modulului coeficientului de reflexie al pamantului:

   

Variatia fazei coeficientului de reflexie al pamantului:

   

Pierderile de propagare in spatiul liber si deasupra suprafetelor plane, functie de distanta cu frecventa parametru