3D pattern plotting function seems to be wrong !

3 次查看（过去 30 天）

Krishan Kumar Tiwari 2023-10-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2036756-3d-pattern-plotting-function-seems-to-be-wrong

编辑： Krishan Kumar Tiwari 2023-11-25

I wrote the following function for plotting 3D pattern of an array with position vectors in a matrix named ``RIS'' and the beamforming weights in W, there seems to be an error which I am unable to locate, any comments ?

function [G] = Full_3D_Pattern (W, RIS)
theta = -pi/2 : pi/2e2 : pi/2; phi = theta; %angle steps of 1 degree for plotting
G = zeros ( length(theta), length(phi) );
A = zeros ( length(W),     length(theta), length(phi) );
for k = 1:length(phi)
    for tta = 1:length(theta)
        ux = cos(phi(k))*sin(theta(tta)); 
        uy = cos(phi(k))*cos(theta(tta)); 
        uz = sin (phi(k));
        u = [ux; uy; uz];
        A (:,tta,k) = exp (-1i*pi*u'*RIS); % RIS steering vector for tta, k
        G (tta,k)   = abs(W' * A (:,tta,k)) * cos(theta(tta))*cos(phi(k)) * 2;
    end
end
figure('Name','3D Pattern')
set(0,'DefaultAxesFontSize',20,'DefaultTextFontSize',20);
mesh(theta,phi,abs(G).^2);xlim([-pi/2 pi/2]);ylim([-pi/2 pi/2]); %%%% radian space
xlabel('Phi','FontSize',22); ylabel('Theta','FontSize',22); colorbar
end

6 个评论
显示 4更早的评论隐藏 4更早的评论

Star Strider 2023-10-22

在 MATLAB Online 中打开

The code runs without error. Is there a problem with the result?

LD1 = load('RIS_Full.mat');

RIS = LD1.RIS_Full;

LD2 = load('W.mat');

W = LD2.W;

G = Full_3D_Pattern (W, RIS);

function [G] = Full_3D_Pattern (W, RIS)

theta = -pi/2 : pi/2e2 : pi/2; phi = theta; %angle steps of 1 degree for plotting

G = zeros ( length(theta), length(phi) );

A = zeros ( length(W), length(theta), length(phi) );

for k = 1:length(phi)

for tta = 1:length(theta)

ux = cos(phi(k))*sin(theta(tta));

uy = cos(phi(k))*cos(theta(tta));

uz = sin (phi(k));

u = [ux; uy; uz];

A (:,tta,k) = exp (-1i*pi*u'*RIS); % RIS steering vector for tta, k

G (tta,k) = abs(W' * A (:,tta,k)) * cos(theta(tta))*cos(phi(k)) * 2;

end

figure('Name','3D Pattern')

set(0,'DefaultAxesFontSize',20,'DefaultTextFontSize',20);

mesh(theta,phi,abs(G).^2);xlim([-pi/2 pi/2]);ylim([-pi/2 pi/2]); %%%% radian space

xlabel('Phi','FontSize',22); ylabel('Theta','FontSize',22); colorbar

end

Krishan Kumar Tiwari 2023-10-22

编辑：Krishan Kumar Tiwari 2023-10-22

Yes, there is a problem, if you give el = az = 0, then you get a fan / blade beam, when you steer it, it gets curved -- just steering causes the curvature -- it should not !!

Krishan Kumar Tiwari 2023-10-22