pr1

clc
clear all
close all

% Пункт1
%Pt=52;
%Gt=25;
%Gr=15;

Pt=158.5;
Gt=17.8;
Gr=5.62;

f=[0.003*10^9 0.030*10^9 0.100*10^9 0.500*10^9 1*10^9 5*10^9 15*10^9 25*10^9];
r=1000:1:500000;
c=3*10^8;
figure(100)
for i=1:8
Pr(i,:)=(Pt.*Gt.*Gr.*(c./f(i)).^2)./((4.*pi.*r).^2);

Pr(i,:)=10*log10(Pr(i,:));

plot(r,Pr(i,:))
hold on
end
grid
title('График зависимости принимаемой мощности от расстояния')
legend('f1', 'f2','f3','f4','f5','f6','f7','f8')
hold off

% Пункт2
d=1:50:200000;
sigma=1;
d0=1;
n=3;
nsv=2;
hi=real(random('Normal',0,sigma,1,4000));


Pl=(10*2*log10(1)+hi)-10.*n.*log10(d./d0)+hi;
Plsv=(10*2*log10(1))-10.*n.*log10(d./d0);

figure(101)
semilogx(d,Pl)
hold on
semilogx(d,Plsv)
grid
title('График зависимости потерь от расстояния')

% Пункт3
sigma2=0.2;
N=100000000;
step=0.1;
n=0:0.1:3;

x=randn(1,N);
y=randn(1,N);

qqq=sqrt(sigma2.*(x.^2+y.^2));

h=hist(qqq,n);

www=h/(step*sum(h));

ReleyT=n./sigma2.*exp(-n.^2./(2.*sigma2));

figure(102)
plot(n,www)
hold on
plot(n,ReleyT)
grid
title('График плотности вероятности распределения Рэлея')
legend('Смоделир.','Теор.')
hold off

% Пункт4
sigma3=1;
m1=[0 1 2 4 8];
m2=0;
v=0:0.1:10;
figure(103)
for i=m1
s=sqrt(i+m2);
%Z(i,:)=v*s/sigma3;
%J= besseli(0,v*s/sigma3);

RaisT=(v / sigma3) .* exp(-(v.^2 + s^2) / (2*sigma3)) .* besseli(0, v * s / sigma3);

plot(v,RaisT)
grid
hold on
end
title('График плотности вероятности распределения Райса')