how can i plot radiation pattern in cartesian coordinate in matlab ?

Question

8167 2013-10-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/90085-how-can-i-plot-radiation-pattern-in-cartesian-coordinate-in-matlab

评论： Abdul Lawal 2016-1-28

采纳的回答： sixwwwwww

nn

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

sixwwwwww 2013-10-13

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/90085-how-can-i-plot-radiation-pattern-in-cartesian-coordinate-in-matlab#answer_99512

在 MATLAB Online 中打开

Dear Naveedp, here is the code you can use for plotting:

theta = -90:90;
theta_rad = degtorad(theta);
r = 1;
elevation = 0;
[x,y,z] = sph2cart(theta_rad,elevation,r);
E_theta = 1j * 30 * exp(-1j * 2 * pi) * sin(5 * pi * sin(theta)) ./ (5 * pi * sin(theta));
E_phi = 1j * 30 * exp(-1j * 2 * pi) * cos(theta) .* sin(3 * pi * sin(theta)) ./ (3 * pi * sin(theta));
figure, subplot(2,1,1), plot3(x, y, abs(E_theta)), xlabel('x'), ylabel('y'), zlabel('z'), title('E(theta) amplitude')
        subplot(2,1,2), plot3(x, y, angle(E_theta)), xlabel('x'), ylabel('y'), zlabel('z'), title('E(theta) phase')
figure, subplot(2,1,1), plot3(x, y, abs(E_phi)), xlabel('x'), ylabel('y'), zlabel('z'), title('E(phi) amplitude')
        subplot(2,1,2), plot3(x, y, angle(E_phi)), xlabel('x'), ylabel('y'), zlabel('z'), title('E(phi) phase')

Since, I did't see use of angle phi so I assumed it to be 0 and assumed value of radius to be 1. I hope it will help you

12 个评论
显示 10更早的评论隐藏 10更早的评论

sixwwwwww 2013-10-13

在 MATLAB Online 中打开

Dear Navdeep, here is your code:

r = 1;
theta_degree = -90:90;
theta = degtorad(theta_degree);
theta = theta(theta ~= 0);
phi_degree = [0 90];
phi = degtorad(phi_degree);
for i = 1:length(phi)    
    if phi_degree(i) == 0
        Y = 0;
        sinYY = 1;
    else
        Y = 5 * pi * sin(theta) * sin(phi(i));
        sinYY = sin(Y) ./ Y;
    end
    if phi_degree(i) == 90
        X = 0;
        sinXX = 1;
    else
        X = 3 * pi * sin(theta) * cos(phi(i));
        sinXX = sin(X) ./ X;
    end
    E_theta = (1j / r) *  30 * exp (-1j * 2 * pi) * sin(phi(i)) * sinXX * sinYY;
    E_phi = (1j / r) * 30 * exp(-1j * 2 * pi) * cos(theta) * cos(phi(i)) .* sinXX  .* sinYY;
    [x,y,z] = sph2cart(theta,phi(i),r);
    figure, plot3(x, y, abs(E_theta)), xlabel('x'), ylabel('y'), zlabel('E(theta) amplitude'), title(strcat('E(theta) amplitude for phi =',...
        num2str(phi_degree(i))))
    figure, plot3(x, y, angle(E_theta)), xlabel('x'), ylabel('y'), zlabel('E(theta) phase'), title(strcat('E(theta) amplitude for phi =',...
        num2str(phi_degree(i))))
    figure, plot3(x, y, abs(E_phi)), xlabel('x'), ylabel('y'), zlabel('E(phi) amplitude'), title(strcat('E(phi) amplitude for phi =',...
        num2str(phi_degree(i))))
    figure, plot3(x, y, angle(E_phi)), xlabel('x'), ylabel('y'), zlabel('E(phi) phase'), title(strcat('E(phi) amplitude for phi =',...
        num2str(phi_degree(i))))
end

If you can't understand something, feel free to ask. Also if you like my answer then accept the answer to make others to find the solution as well if they are looking for similar question. Good luck!

sixwwwwww 2013-10-13

在 MATLAB Online 中打开

Dear Navdeep, now see this code, it will give you exact output as you need:

r = 1;
theta_degree = -90:0.01:90;
theta = degtorad(theta_degree);
theta = theta(theta ~= 0);
phi_degree = [0 90];
phi = degtorad(phi_degree);
for i = 1:length(phi)    
    if phi_degree(i) == 0
        Y = 0;
        sinYY = 1;
    else
        Y = 5 * pi * sin(theta) * sin(phi(i));
        sinYY = sin(Y) ./ Y;
    end
    if phi_degree(i) == 90
        X = 0;
        sinXX = 1;
    else
        X = 3 * pi * sin(theta) * cos(phi(i));
        sinXX = sin(X) ./ X;
    end
    E_theta = (1j / r) *  30 * exp (-1j * 2 * pi) * sin(phi(i)) * sinXX * sinYY;
    E_phi = (1j / r) * 30 * exp(-1j * 2 * pi) * cos(theta) * cos(phi(i)) .* sinXX  .* sinYY;
    [x,y,z] = sph2cart(theta,phi(i),r);
    figure, plot3(x, y, abs(E_theta)), xlabel('x'), ylabel('y'), zlabel('E(theta) amplitude'), title(strcat('E(theta) amplitude for phi =',...
        num2str(phi_degree(i)))), view([90 0])
    figure, plot3(x, y, abs(E_phi)), xlabel('x'), ylabel('y'), zlabel('E(phi) amplitude'), title(strcat('E(phi) amplitude for phi =',...
        num2str(phi_degree(i)))), view([90 0])
end

Feel free to ask if you can't understand anything

sixwwwwww 2013-10-13

在 MATLAB Online 中打开

I did but it doesn't look good. But maybe you can also try by yourself. Here is the code:

r = 1;
theta_degree = -90:0.01:90;
theta = degtorad(theta_degree);
theta = theta(theta ~= 0);
phi_degree = [0 90];
phi = degtorad(phi_degree);
for i = 1:length(phi)    
    if phi_degree(i) == 0
        Y = 0;
        sinYY = 1;
    else
        Y = 5 * pi * sin(theta) * sin(phi(i));
        sinYY = sin(Y) ./ Y;
    end
    if phi_degree(i) == 90
        X = 0;
        sinXX = 1;
    else
        X = 3 * pi * sin(theta) * cos(phi(i));
        sinXX = sin(X) ./ X;
    end
    E_theta = (1j / r) *  30 * exp (-1j * 2 * pi) * sin(phi(i)) * sinXX * sinYY;
    E_phi = (1j / r) * 30 * exp(-1j * 2 * pi) * cos(theta) * cos(phi(i)) .* sinXX  .* sinYY;
    E = E_theta + E_phi;
    E_dB = 10 * log(E);
    [x,y,z] = sph2cart(theta,phi(i),r);
    figure, plot3(x, y, abs(E)), xlabel('x'), ylabel('y'), zlabel('Radiation Intensity'), title(strcat('Radiation Intensity at phi=',...
        num2str(phi_degree(i)))), view([90 0])
    figure, plot3(x, y, abs(E_dB)), xlabel('x'), ylabel('y'), zlabel('Radiation Intensity [dB]'), title(strcat('Radiation Intensity at phi=',...
        num2str(phi_degree(i)))), view([90 0])
end

You will get 4 figures. second and fourth figures are in dB scales

sixwwwwww 2013-10-13

在 MATLAB Online 中打开

Here is your code:

r = 1;
theta_degree = -90:0.01:90;
theta_degree = theta_degree(theta_degree ~= 0);
theta = degtorad(theta_degree);
phi_degree = [0 90];
phi = degtorad(phi_degree);
for i = 1:length(phi)    
    if phi_degree(i) == 0
        Y = 0;
        sinYY = 1;
    else
        Y = 5 * pi * sin(theta) * sin(phi(i));
        sinYY = sin(Y) ./ Y;
    end
    if phi_degree(i) == 90
        X = 0;
        sinXX = 1;
    else
        X = 3 * pi * sin(theta) * cos(phi(i));
        sinXX = sin(X) ./ X;
    end
    E_theta = (1j / r) *  30 * exp (-1j * 2 * pi) * sin(phi(i)) * sinXX * sinYY;
    E_phi = (1j / r) * 30 * exp(-1j * 2 * pi) * cos(theta) * cos(phi(i)) .* sinXX  .* sinYY;
    E = E_theta + E_phi;
    E = E / max(E);
    E_dB = 10 * log(E);
    figure, plot(theta_degree, abs(E)), xlabel('Angle [theta]'), ylabel('Normalized Radiation Intensity'), title(strcat('Radiation Intensity at phi=',...
        num2str(phi_degree(i))))
    figure, plot(theta_degree, abs(E_dB)), xlabel('Angle [theta]'), ylabel('Normalized Radiation Intensity [dB]'),...
        title(strcat('Radiation Intensity at phi=', num2str(phi_degree(i))))
end

I forgot to ask you. Is it your homework? Tell me truly hahhahha

sixwwwwww 2013-10-13

在 MATLAB Online 中打开

hahahhaha. It's not fair you should have told me before and now I should get marks instead of you. lolx. You can add grid the following way:

figure, plot(theta_degree, abs(E)), xlabel('Angle [theta]'), ylabel('Normalized Radiation Intensity'), title(strcat('Radiation Intensity at phi=',...
        num2str(phi_degree(i)))), grid on

Abdul Lawal 2016-1-28

Hey sixwwwwww, I need to plot the two-way radiation pattern of a planar array for my SAR antenna for both elevation and azimuth. The figures are necessary as input to my Noise Equivalent Sigma Zero (NESZ) computation. Can you please help?

请先登录，再进行评论。