Hi Shahid,
As per my understanding, you wish to seek clarification on how to get the expected result from the code provided while plotting the colormap.
After reviewing the code, consider changing the variable from “M(j,i)” to “M(i,j)” in the nested “for” loop as shown below:
for j=1:23
for i=1:23
M(i,j) = mod(-ko*(x_v(j)*(sind(theta_i)*cosd(phi_i)+sind(theta_r)*cosd(phi_r)) ...
+y_v(i)*(sind(theta_i)*sind(phi_i)+sind(theta_r)*sind(phi_r))),360);
end
end
I hope this resolves your questions as making the above change also ensures x cells and y cells are plotted on X and Y axis respectively.
I would also like to suggest some optimization for the code to perform the same tasks using “meshgrid” function. “meshgrid” function creates a 2-D grid of the required values on X and Y-axis and eliminates the need to use “for” loops. Consider the below example implementation:
% Define the x and y coordinates
x_v = 1e-3*(4.5 + 9*(0:22));
y_v = 1e-3*(4.5 + 9*(0:22));
% Create a grid of coordinates
[X, Y] = meshgrid(x_v, y_v);
% Compute the values of each cell
M = mod(-ko*(X.*(sind(theta_i)*cosd(phi_i)+sind(theta_r)*cosd(phi_r)) ...
+Y.*(sind(theta_i)*sind(phi_i)+sind(theta_r)*sind(phi_r))),360);
% Plot the result
imagesc(M)
colormap(jet)
colorbar
Please refer to below MathWorks documentation for more information on “meshgrid” function: https://www.mathworks.com/help/matlab/ref/meshgrid.html
I hope this helps!
