RIS_Models/ArrayRCS_RandomPhaseShift.m at main · mitiebin/RIS_Models · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
clear all;close all;clc;
%%
deltaTheta = 0.1;

C = -2*1i;

lambda = 1;

%%
a = 0.1;
b = 0.1;

d = 0.5;

numUnits = 100;
unitsPos = (0:1:numUnits-1)*d;


%%
thetaInDeg = (-90:deltaTheta:90)';
thetaIn = thetaInDeg/180*pi;

thetaOutDeg = (-90:deltaTheta:90)';
thetaOut = thetaOutDeg/180*pi;

[thetaInDegMesh, thetaOutDegMesh] = meshgrid(thetaInDeg, thetaOutDeg);
[thetaInMesh, thetaOutMesh] = meshgrid(thetaIn, thetaOut);

thetaIn = reshape(thetaInMesh, [], 1);
thetaOut = reshape(thetaOutMesh, [], 1);

%%
Sa = sinc(b./lambda*(sin(thetaOut)+sin(thetaIn)));

RCS = numUnits*4*pi*(abs(C))^2*(a*b)^2.*(cos(thetaIn)).^2.*Sa;

RCS = reshape(RCS,size(thetaInMesh));

%%
figure
mesh(thetaInDegMesh, thetaOutDegMesh, RCS)

xlabel('$\theta^i$ (degree)','interpreter','latex')
ylabel('$\theta^s$ (degree)','interpreter','latex')
zlabel('$E [ \sigma(\theta^s; \theta^i) ]$','interpreter','latex')
axis tight
rotate3d on
% axis equal
% view(45, 30);
% view(0, 90);
box on

xh = get(gca,'XLabel'); % Handle of the x label
set(xh, 'Units', 'Normalized')
pos = get(xh, 'Position');
set(xh, 'Position',pos.*[1.1,-0.5,1],'Rotation',15)
yh = get(gca,'YLabel'); % Handle of the y label
set(yh, 'Units', 'Normalized')
pos = get(yh, 'Position');
set(yh, 'Position',pos.*[0.7,-0.6,1],'Rotation',-25)

exportgraphics(gcf, 'ArrayRCS_RandomPhaseShift.pdf');