How to draw multiple surface plots on the same mesh?

Question

Muhamed Sewidan 2020-10-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/609021-how-to-draw-multiple-surface-plots-on-the-same-mesh

评论： Muhamed Sewidan 2020-10-14

I use this code to draw a Gaussian normal distribution. I want to plot several distributions at equal distance from the original one (like the figure).

How can I do so?

lambda = 1;               % wavelength
E0 = 1;                   % max value of E-field amplitude
% ==================================================================
% prepare the coordinates mesh
% ----------------------------
zLin = linspace(-201,201,400) * lambda;     % z-coordinate (horizontal)
rLin = linspace(-65,65,100) * lambda;       % radius-coordinate (vertical)
[z,r] = meshgrid(zLin,rLin);                % create mesh
% ===================================================================
wFactor = 4;                 
    
w0 = wFactor * lambda;                  % minimal waist
zR = pi* w0^2 / lambda;                 % Rayleigh length
w = w0 * sqrt( 1 + (z / zR) .^2 );      % Beam radius at any arbitrary distance (equ. 3)
% Irradiance
eta = 377;                              % for free space
I0 = abs(E0)^2/eta;
I = I0 .* (w0 ./ w).^2 .* exp( -2*(r./w).^2 );
surf(z,r,I); shading interp

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

VBBV 2020-10-9

编辑：VBBV 2020-10-9

在 MATLAB Online 中打开

Use contourf function as

%if true
  % code
% end
v2 = 0:0.00009:0.0025;
contourf(z,r,I,v2,'edgecolor','none');shading interp;

请先登录，再进行评论。

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

Answer 1

Durganshu 2020-10-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/609021-how-to-draw-multiple-surface-plots-on-the-same-mesh#answer_509561

编辑：Durganshu 2020-10-9

在 MATLAB Online 中打开

Did you try hold on?

figure;
surf(z1,r1,I1);
hold on
zlin2= linspace(-101,101,300) * lambda; %whatever shifted values you want
rLin2 = linspace(-65,65,100) * lambda; %whatever shifted values you want
[z2,r2] = meshgrid(zLin2,rLin2);
% similarly define I2 and then plot
surf(z2,r2,I2);

9 个评论
显示 7更早的评论隐藏 7更早的评论

Muhamed Sewidan 2020-10-12

在 MATLAB Online 中打开

Never mind, thank you for your interest.

Should I change rlin2 also the same way I change zlin2?

I did as you said, but it didn't plot I2, it is just extended the mesh of the first plot. Please have a look at my code after adding your suggestion.

lambda = 1;               % wavelength
E0 = 1;                   % max value of E-field amplitude
% ==================================================================
% prepare the coordinates mesh
% ----------------------------
zLin = linspace(-201,201,400) * lambda;     % z-coordinate (horizontal)
rLin = linspace(-65,65,100) * lambda;       % radius-coordinate (vertical)
[z,r] = meshgrid(zLin,rLin);                % create mesh
% ===================================================================
% Applying the code:
% ------------------
wFactor = 4;                 
    
w0 = wFactor * lambda;                  % minimal waist
zR = pi* w0^2 / lambda;                 % Rayleigh length
w = w0 * sqrt( 1 + (z / zR) .^2 );      % Beam radius at any arbitrary distance (equ. 3)
% Irradiance
eta = 377;                              % for free space
I0 = abs(E0)^2/eta;
I = I0 .* (w0 ./ w).^2 .* exp( -2*(r./w).^2 );
surf(z,r,I)
shading interp
%% Second source:
hold on
zlin2= linspace(201,602,400) * lambda; %whatever shifted values you want
rLin2 = linspace(-65,65,100) * lambda; %whatever shifted values you want
[z2,r2] = meshgrid(zlin2,rLin2);
w2 = w0 * sqrt( 1 + (z2 / zR) .^2 );
% similarly define I2 and then plot
I2 = I0 .* (w0 ./ w2).^2 .* exp( -2*(r2./w2).^2 );
surf(z2,r2,I2);

Durganshu 2020-10-12

Again, your code is absolutely okay. Let's see the differences between I and I2:

I includes w and I2 includes w2.
I includes r and I2 includes r2.

Now, w2 and w depend on z2 and z. Looking at the values of z2 and z (that you get in the MATLAB workspace), it is clearly visible that there must be a continuos boundary between the curves I and I2. That is also visible in the curve.

Now, r2 and r are the same.

COCLUSION: I2 is just a natural extension of I. I2 starts from where the range of I ends. Basically, by writing this code, you have extended the range of I to 602. It could be written that range(I') = range(I) + range(I2). That's why you're getting that kind of plot.

My answers were to guide you and not give you the final answer. You never specified anything about the inputs and what exactly is going in the problem. The code is alright. You just need to be clear about what exactly curves do you want. If you want a different curve, then you need to put something else in the equation for I2.

How can I tell you if you need to change r2 or not? It depends on your original problem statement.

So, please check on the exact requirements and change the equation accordingly.

If you're still stuck somewhere, then post what kind of output/shifts you want with detailed explanation.

Durganshu 2020-10-14

编辑：Durganshu 2020-10-14

Yes, that should work. Did you try it? All you're supposed to do is tweaking the I2 function in order to find the appropriate fit.

Muhamed Sewidan 2020-10-14

Thank you I will try it out.

请先登录，再进行评论。

How to draw multiple surface plots on the same mesh?

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

采纳的回答

9 个评论
显示 7更早的评论隐藏 7更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

How to draw multiple surface plots on the same mesh?

1 个评论 显示 -1更早的评论隐藏 -1更早的评论

采纳的回答

9 个评论 显示 7更早的评论隐藏 7更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

9 个评论
显示 7更早的评论隐藏 7更早的评论