I have the following grid in the beginning of my code:
a = 25, % This will be the radius of the circle
r = linspace(0,sqrt(2)*2*a,501); phi = linspace(0,2*pi,501);[R,Phi] = meshgrid(r,phi);X = R.*cos(Phi);Y = R.*sin(Phi);
And I plot all my figures by using the code segment:
figure()h = pcolor(X,Y,abs(u_sum).^2); % u_sum is an electric field distribution
set(h,'EdgeColor','none');box ongrid onaxis equalxlim([-1,1]*r(end)/sqrt(2));ylim([-1,1]*r(end)/sqrt(2));colormap(hot);colorbar;xlabel('x');ylabel('y');title(['Title'])hold onrectangle('Position',[-1,-1,2,2]*a,'Curvature',1,... 'LineWidth',1,'EdgeColor','w','FaceColor','None',... 'LineStyle','--')hold offset(gca,'Layer','top')
The output of this code is:
Now, I want to draw a circle full of ones with radius "a" and place it on the same grid. Thus, I want the circle at the center of the above picture to be filled with ones, and zeros at rest. How can I do that using pcolor and the same grid?
Alternatively, how can I implement a Gaussian shape for the same grid using pcolor?
Best Answer