pr4

clc
close all
clear all

EF=1.5;
f=1.5*10^9;
c=3e8;
lambda=c/f;

d=lambda/2;

thetaindegrees=-90:0.5:90;
theta=thetaindegrees*(pi/180);
theta0indegrees=[0 20 40 70];
theta0=theta0indegrees.*(pi/180);
M=16;
m=1:M;
EP=cos(theta).^(EF/2);
x=(m-0.5*(M+1))*d;
for ii=1:length(theta0)
for i=1:length(x)
AFMassiv(i,:)=1*exp(j*(2*pi./lambda.*x(i).*sin(theta)-2*pi./lambda.*x(i).*sin(theta0(ii))));
end
AF=sum(AFMassiv);
%s1=1*exp(j*(2*pi./lambda.*x(1).*sin(theta)-2*pi./lambda.*x(1).*sin(theta0)));
%F =cos(theta).^(EF/2).*s1;
Fall(ii,:) =EP.*AF;
%Fdb=20*log10(F);
Fallnorm(ii,:)=Fall(ii,:)/max(Fall(ii,:));
Falldb(ii,:)=20*log10(Fallnorm(ii,:));
AFnorm=AF/max(AF);
AFdb(ii,:)=20*log10(AFnorm);
EPnorm=EP/max(EP);
EPdb=20*log10(EPnorm);

figure(200)
plot(thetaindegrees,Falldb(ii,:))
xlabel('theta'), ylabel('F, дБ')
hold on
end
legend('угол ск.=0','угол ск.=20','угол ск.=40','угол ск.=70')
grid
ylim([-40 5])
figure(100)
plot(thetaindegrees,Falldb(1,:),thetaindegrees,EPdb,thetaindegrees,AFdb(1,:))
grid
title('График норм. ДН в дБ')
xlabel('theta'), ylabel('F, дБ')
legend('ДН итоговой ФАР', 'ДН одиночного элемента', 'ДН множителя решетки')
ylim([-40 5])


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

hamm=hamming(length(Falldb(1,:)));
Fallwindowed=Fall(1,:).*hamm';
Fallwindowednorm=Fallwindowed/max(Fallwindowed);
Falldbwindowed=20*log10(Fallwindowednorm);
figure(300)
plot(thetaindegrees,Falldb(1,:),thetaindegrees,Falldbwindowed(1,:))
grid
title('График норм. ДН в дБ')
xlabel('theta'), ylabel('F, дБ')
legend('ДН итоговой ФАР', 'ДН итоговой ФАР с окном')
ylim([-40 5])