How to normalize plot?

Question

Athira T Das 2023-5-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1957989-how-to-normalize-plot

评论： Torsten 2023-5-4

在 MATLAB Online 中打开

How to normalize the plot to one?

clc;clear all;close all;

syms x y

x0 = -.02:0.001:.02;lambda=532*10^-9;

y0 = -.02:0.001:.02;w0 = 0.002;

omegax=1;omegay=2;zr=pi*w0^2/lambda;

m = 2;

s=1;k=0.0002;Rx=zr;Ry=Rx;

[X0,Y0] = meshgrid(x0,y0);

f = @(x0,y0)sinh((omegax*x0+omegay*y0)/w0).*exp(-(x0.^2+y0.^2)/w0^2).*(x0+(1i*s.*y0)).^m.*exp((1i*k/2)*((x0^2/Rx)+(y0^2/Ry)));

I_numerical = real(abs(vpa(f)));

figure

fcontour(real(I_numerical), [-1 1 -1 1]*1E-2, 'Fill','on', 'MeshDensity',150)

colormap('hot');

colorbar('vert');

axis('equal')

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

请先登录，再进行评论。

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

Answer 1

Torsten 2023-5-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1957989-how-to-normalize-plot#answer_1229089

在 MATLAB Online 中打开

x0 = -.02:0.001:.02;lambda=532*10^-9;

y0 = -.02:0.001:.02;w0 = 0.002;

omegax=1;omegay=2;zr=pi*w0^2/lambda;

m = 2;

s=1;k=0.0002;Rx=zr;Ry=Rx;

[X0,Y0] = meshgrid(x0,y0);

f = @(x0,y0)sinh((omegax*x0+omegay*y0)/w0).*exp(-(x0.^2+y0.^2)/w0^2).*(x0+(1i*s.*y0)).^m.*exp((1i*k/2)*((x0^2/Rx)+(y0^2/Ry)));

Z0 = arrayfun(@(x,y)f(x,y),X0,Y0);

contourf(X0,Y0,real(abs(Z0))/max(real(abs(Z0(:)))))

colormap('hot');

colorbar('vert');

axis('equal')

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

Athira T Das 2023-5-4

在 MATLAB Online 中打开

I tried this code for my function and showing some error

clc;clear all;close all;
z=100;
y=linspace(-.1,.1,500);x=y;
[X0,Y0] = meshgrid(x,y);
m =1;
lambda=532*10^(-9);
wo = 0.02;
Omega = 30;
k = (2.*pi)./(lambda); 
Pho=0.05;sgm=1;
Gamma = ((1i.*k)./(2.*z)) + (1./(Pho.^2)) + (1./(wo.^2));
Gamma_star = subs(Gamma, 1i, -1i);
alpha = (Gamma_star) - 1./((Gamma.*Pho.^4));
C = (1./(4.*(lambda.^2).*(z.^2))).*((pi.^2)./(alpha.*Gamma)).*(1./(2.*1i.*sqrt(Gamma)).^m);
fx = ((Omega.^2)./(4.*Gamma)) + ((1i.*k.*x.*Omega)./(2.*z.*Gamma)) - ((k.^2.*x.^2)./(4.*z.^2.*Gamma));
fy = ((Omega.^2)./(4.*Gamma)) + ((1i.*k.*y.*Omega)./(2.*z.*Gamma)) - ((k.^2.*y.^2)./(4.*z.^2.*Gamma));
gx = ((Omega.^2)./(4.*Gamma)) - ((1i.*k.*x.*Omega)./(2.*z.*Gamma)) - ((k.^2.*x.^2)./(4.*z.^2.*Gamma));
gy = ((Omega.^2)./(4.*Gamma)) - ((1i.*k.*y.*Omega)./(2.*z.*Gamma)) - ((k.^2.*y.^2)./(4.*z.^2.*Gamma));
Bx_p = ((-1i.*k.*x)./(2.*z)) + ((1i.*k.*x)./(2.*z.*Gamma.*Pho.^2)) + (Omega)./(2.*Gamma.*Pho.^2) + (Omega./2);
Bx_n = ((-1i.*k.*x)./(2.*z)) + ((1i.*k.*x)./(2.*z.*Gamma.*Pho.^2)) + (Omega)./(2.*Gamma.*Pho.^2) - (Omega./2);
By_p = ((-1i.*k.*y)./(2.*z)) + ((1i.*k.*y)./(2.*z.*Gamma.*Pho.^2)) + (Omega)./(2.*Gamma.*Pho.^2) + (Omega./2);
By_n = ((-1i.*k.*y)./(2.*z)) + ((1i.*k.*y)./(2.*z.*Gamma.*Pho.^2)) + (Omega)./(2.*Gamma.*Pho.^2) - (Omega./2);
Kx_p = ((-1i.*k.*x)./(2.*z)) + ((1i.*k.*x)./(2.*z.*Gamma.*Pho.^2)) - (Omega)./(2.*Gamma.*Pho.^2) + (Omega./2);
Kx_n = ((-1i.*k.*x)./(2.*z)) + ((1i.*k.*x)./(2.*z.*Gamma.*Pho.^2)) - (Omega)./(2.*Gamma.*Pho.^2) - (Omega./2);
Ky_p = ((-1i.*k.*y)./(2.*z)) + ((1i.*k.*y)./(2.*z.*Gamma.*Pho.^2)) - (Omega)./(2.*Gamma.*Pho.^2) + (Omega./2);
Ky_n = ((-1i.*k.*y)./(2.*z)) + ((1i.*k.*y)./(2.*z.*Gamma.*Pho.^2)) - (Omega)./(2.*Gamma.*Pho.^2) - (Omega./2);
X = 0;
for l = 0:m
  L = (((1i).*sgm).^(m-l)).*nchoosek(m,l);
  sum_j = 0;
  for j = 0:m
      J= (((-1i).*sgm).^(m-j)).*nchoosek(m,j);  
      
      sum1 = 0.0;
      for p = 0:1/2:l/2
            faktor1_P = (((-1).^p).*factorial(l))./(gamma(p+1).*factorial(l-(2.*p))).*(((2.*1i)./(sqrt(Gamma))).^(l-(2.*p))).*exp(fx).*exp(fy);
            for s1 = 0:l-2*p
                faktor1_s1 =  nchoosek((l-(2.*p)),s1).*((1./(Pho.^2)).^s1).*((((1i.*k.*x)./(2.*z))+(Omega./2)).^(l-(2.*p)-s1)).*(1./(((2.*1i.*sqrt(alpha))).^(j+s1)));
                for q = 0:1/2:(m-l)/2
                    faktor1_Q = (((-1).^q).*factorial(m-l))./(gamma(q+1).*factorial(m-l-(2.*q))).*(((2.*1i)./(sqrt(Gamma))).^(m-l-(2.*q)));
                   sum_inner_1=0;
                    for s2 = 0:m-l-2*q
                        E1 = (exp(((Bx_p).^2)./alpha).*hermiteH(j+s1,((1i.*Bx_p)./(sqrt(alpha)))).*exp(((By_p).^2)./alpha).*hermiteH(m-j+s2,((1i.*By_p)./(sqrt(alpha))))) - (exp(((Bx_n).^2)./alpha).*hermiteH(j+s1,((1i.*Bx_n)./(sqrt(alpha)))).*exp(((By_n).^2)./alpha).*hermiteH(m-j+s2,((1i.*By_n)./(sqrt(alpha)))));
                        sum_inner_1 =sum_inner_1+ nchoosek((m-l-(2.*q)),s2).*((1./(Pho.^2)).^s2).*((((1i.*k.*y)./(2.*z))+(Omega./2)).^(m-l-(2.*q)-s2)).*(1./(((2.*1i.*sqrt(alpha))).^(m-j+s2))).*E1;
                    end
                end
            end
            sum1 = sum1 + faktor1_P.*faktor1_s1.*faktor1_Q.*sum_inner_1;
       end
   sum_j= sum_j + J.*(sum1);
  end
 X = X + (L.*sum_j);
end  
X = subs(real(C.*X));
I_numerical = real(abs(vpa(X)));
f = @(x,y) I_numerical;
Z0 = arrayfun(@(x1,y1)f(x1,y1),X0,Y0);
Error using arrayfun
Non-scalar in Uniform output, at index 1, output 1.
Set 'UniformOutput' to false.
contourf(X0,Y0,real(abs(Z0))/max(real(abs(Z0(:)))))
colormap('hot');
colorbar('vert');
axis('equal')