Diferencia entre revisiones de «Level sets of a scalar field»
De MateWiki
| Línea 47: | Línea 47: | ||
[[Categoría:TC20/21]] | [[Categoría:TC20/21]] | ||
[[Categoría:TC21/22]] | [[Categoría:TC21/22]] | ||
| − | [[Categoría: | + | [[Categoría:TC23/24]] |
[[Categoría:Informática]] | [[Categoría:Informática]] | ||
Revisión del 18:39 10 nov 2023
We show how to visualize level sets of scalar fields on plane regions with Octave UPM. We focus on the example in Dibujar un sólido 2-D, i.e. the rectangle [math] [-1/2,1/2]\times [0,2][/math] and the scalar field [math] f(x,y)=-\log (0.1+\sqrt{x^2+y^2})[/math]. We follow the steps:
- We introduce a sampling of the two segments with a suitable step
- With meshgrid command we define two matrixes with the x and y coordenates of the mesh points
- Compute the scalar field in the grid points.
- We use the contour command to draw the field and adjunst the axis. We see the picture from the top.
1 MATLAB code
x=-0.5:0.1:0.5; % sampling of the interval [-1/2,1/2]
y=0:0.1:2; % sampling of the interval [0,2]
[xx,yy]=meshgrid(x,y); % matrixes of x and y coordinates
figure(1)
f=-log(0.1+sqrt(xx.^2+yy.^2)); % The scalar field
contour(xx,yy,f) % Draw the level sets
hold on % We draw the boundary below
plot(x,x-x,'k','linewidth',1);
plot(x,2+x-x,'k','linewidth',1);
plot(-0.5+y-y,y,'k','linewidth',1);
plot(0.5+y-y,y,'k','linewidth',1);
axis([-2,2,-1,3]) % select region for drawing
view(2) % See the pisture from the top
colorbar % Include colorbar
2 To go further
Mesh of a parametrized 2-D solid
Visualization of a scalar field in a solid