Ejercicios Con Matlab-Teoría Electromagnética

EJERCICIOS
1. Ejercicios: graficar la funciones siguientes:

2
a)
f ( x , y ) = ( x 2 y 2 )
b)
f ( x , y )=
c)
f ( x , y ) =sen (| x|+| y|)
1
2
2
1+ x + y
1
2
2
x +y
d)
f ( x , y ) =e
e)
f ( x , y ) =sen x + sen y
f)
f ( x , y )=
sen x sen y
xy
Funcin seno(x)
1
0.8
0.6
0.4
0.2
2. x=-4:.01:4; y=sin(x); plot(x,y), grid, title('Funcin

seno(x)')
0
-0.2
-0.4
-0.6
-0.8
-1
-4
-3
-2
-1
1
0.9
0.8
0.7
0.6
1. x= -2.9:0.2: 2.9; y=exp(-x.*x); plot(x,y)
0.5
0.4
0.3
0.2
0.1
0
-3
-2
-1
15
10
5
2. [X,Y]=meshgrid(-2:.2:2,3:0.5:3);Z=X.^2.*Y;plot3(X,Y,Z)
0
-5
-10
-15
4
2
0
-2
-4
-2
-1
40
20
f ( x , y ) =x3 + y 3
3. [X,Y]=meshgrid(-2:.2:2,3:0.5:3);Z=X.^3.+Y.^3.;plot3(X,Y,Z)
-20
-40
4
2
0
-2
-4
-2
-1
x 10
5
z=e
x +y
3
2
4. X,Y]=meshgrid(-2:.2:2,3:0.5:3);Z=exp(X.^2.+Y.^2.);plot3(X,Y,Z)
1
0
4
2
0
-2
-4
-2
-1
z= x 2+ y 2
3
2
5. [X,Y]=meshgrid(-2:.2:2,3:0.5:3);Z=sqrt(X.^2.+Y.^2.);plot3(X,Y,Z)
1
0
4
2
0
-2
-4
-2
-1
4
3
6. [X,Y]=meshgrid(-2:.2:2,3:0.5:3);Z=sqrt(X.^2.+Y.^2.);plot3(X,Y,Z);mesh(X,
Y,Z);surf(X,Y,Z)
2
1
0
4
2
0
-2
-4
-2
-1
4
3
7. [X,Y]=meshgrid(-2:.2:2,-3:0.5:3);
Z=sqrt(X.^2.+Y.^2.);plot3(X,Y,Z);
mesh(X,Y,Z);surfc(X,Y,Z)
2
1
0
4
2
0
-2
-4
-2
-1
15
z=x 2+ y 2
10
8. [X,Y]=meshgrid(-2:.2:2,-3:0.5:3); Z=X.^2.+Y.^2.;
plot3(X,Y,Z);mesh(X,Y,Z);surf(X,Y,Z)
0
4
2
0
-2
-4
-1
-2
1.5
z=
3 y
2
x + y 2 +1
1
0.5
0
9. [X,Y]=meshgrid(-3:.2:3,-3:.2:3); Z=-3.*Y./
(X.^2.+Y.^2.+1.); plot3(X,Y,Z);
mesh(X,Y,Z);surf(X,Y,Z)
-0.5
-1
-1.5
4
-2
-4
-2
-1
-3
0.2
2
x y
z=xy e
0.1
10.[X,Y]=meshgrid(-3:.2:3,-3:.2:3); Z=X.*Y.*exp(-X.^2.-Y.^2.);
plot3(X,Y,Z);mesh(X,Y,Z);surfc(X,Y,Z)
0
-0.1
-0.2
4
2
0
-2
-4
f ( x , y ) =( x 2 +3 x2 )ex y
11.[X,Y]=meshgrid(-5:.2:5,-5:.2:5);
Z=(X.^2.+3.*Y.^2.)*exp(-X.^2.Y.^2.);plot3(X,Y,Z);
mesh(X,Y,Z);surfc(X,Y,Z)
-4
-2
800
600
400
200
0
5
5
0
0
-5
-5
1.5
1. [X,Y]=meshgrid(-1:0.5:1);
U=2*X;V=-X+Y;
quiver(X,Y,U,V)
0.5
-0.5
-1
-1.5
-1.5
-1
-0.5
0.5
1.5
1.5
2. [X,Y]=meshgrid(-1:0.5:1);
U=exp(X);V=-X+Y;
quiver(X,Y,U,V)
0.5
-0.5
-1
-1.5
-1
-0.5
0.5
1.5
1.5
3. [X,Y]=meshgrid(-1:0.5:1);
U=X;V=sqrt(X^2.+Y^2.);
quiver(X,Y,U,V)
0.5
-0.5
-1
-1.5
-2
-1.5
-1
-0.5
0.5
1.5
1
0.8
0.6
0.4
4. [X,Y]=meshgrid(-1:0.5:1);
Z=X+Y;[c,h]=contour(X,Y,Z);
clabel(c,h)
0.2
0
-0.2
-0.4
-0.6
-0.8
-1
-1
-0.8
-0.6
-0.4
-0.2
0.2
0.4
0.6
0.8
5. Determina y grafica los contornos de un campo

escalar y grafica una distribucin vectorial de un
campo de gradiente asociado.
% Gradient Example
% Draws z-Contours and Then
Plots Gradient Vector Field
>> v = -2:0.1:2;
[x,y] = meshgrid(v);
z = sqrt(x.^2 + y.^2);
[px,py] = gradient(z,.2,.2);
figure(1)
contour(v,v,z,'linewidth',2)
hold on
quiver(v,v,px,py,'linewidth',1)
hold off
title('Gradient Plot');
axis equal
6. Cdigo de onda plana:
[x,y]=meshgrid(5:0.05:15,-5:0.05:5);
lambda=1;
k=2*pi/lambda;
f=3e8/lambda;
w=2*pi*f;
t=linspace(0,60e-9,200);
%z=exp(i*k*abs(x+i*y));
z=exp(i*k*x);
for n=1:length(t)
surf(x,y,real(z*exp(-i*w*t(n))));
%view(ceil(90*n*1/length(t)),ceil(90*n*2/length(t)));
view(3)
zlim([-1.5 1.5])
xlim([min(min(x)) max(max(x))])
ylim([min(min(y)) max(max(y))])
shading interp
getframe();
end
7. f(x,y)=(x sen y, x cos y, y/3):
X=linspace (0,2*pi,40);
Y=linspace (-3 3, 20);
[x,y]=meshgrid (X,Y);
Surf (y.*sin (x), y.*cos (x), x./3);
View (-50, 50);
Axis ([-3 3 -3 3 0 2.*pi/3])
Gradient Plot
2
1.5
1
0.5
0
-0.5
-1
-1.5
-2
-2.5
-2
-1.5
-1
-0.5
0.5
1.5
2.5
1.5
1
0.5
0
-0.5
-1
-1.5
5
15
0
10
-5
x = cos(t), y = sin(t), z = t/(2 )
1.5
8. t=linspace(0,4*pi,1000);x=cos(t);y=sin(t);z=t/(2*pi);
plot3(x,y,z);syms t; x=cos(t);y=sin(t);z=t/(2*pi);
ezplot3(x,y,z,[0,4*pi])
1
0.5
0
1
0.5
0
-0.5
y
-1
-1
-0.5
x
0.5
9. t=linspace(0, 2*pi, 200);

a = 10; b = 1.0; c = 0.3;
x = b*cos(t);
y = b*sin(t);
z = c*cos(a*t);
plot3(x, y, z, 'k')
axis equal
0.2
0
-0.2
0.5
1
0.5
-0.5
-0.5
0.5
1. x=-8:0.5:8; y=-6:0.5:6; [X,Y]=meshgrid(x,y);

r=sqrt(X.^2+Y.^2); Z=sin(r)./r; surf(X,Y,Z)
-0.5
10
5
0
-5
-10
-10
-5
10
1
0.5
1. x=-8:0.5:8; y=-3:0.5:3; [X,Y]=meshgrid(x,y); Z=(X.^2Y.^2)./(X.^2+Y.^2); surf(X,Y,Z)
0
-0.5
-1
4
2
0
-2
-4
1. x=linspace(-20,20,100); y=linspace(-20,20,100);
z=linspace(-20,20,100); [X,Y,Z]=meshgrid(x,y,z);
%Se define Q(x,y,z)=x^2+y^2+z^2-9
Q= X.^2 +(Y.^2) +(Z.^2) -9; %Se representa la
superficie de nivel Q(x,y,z)=0
H=patch(isosurface(X,Y,Z,Q,0)); isonormals(X,Y,Z,Q,
H); %Para dibujar la superficie en color rojo
%e iluminada set(H, 'FaceColor', 'red', 'EdgeColor',
'none');
daspect([1 1 1])
view(3)
camlight
lighting phong
axis equal
xlabel('x'); ylabel('y'); zlabel('z')
-10
-5
10
1
0.9
0.8
0.7
0.6
2. zz = linspace(0, 2*pi, 26); [x, y, z] =

cylinder(1.1+sin(zz), 16); surf(x, y, z) axis of
0.5
0.4
0.3
0.2
0.1
0
Actividad Experimental # 1: Grafica varias funciones paramtricas:

a) r(u,t)={cos t (3+cos u), sen t (3+cos u), sen u}
[u v]=meshgrid(linspace(0,2*pi,40));
surf(cos(u).*(3+cos(v)),sin(u).*(3+cos(v)),sin(v));
view(-37.5,50);
axis([-4 4 -4 4 -3 3])
2
0
-2
4
2
4
2
-2
-2
-4
-4
b) r(u,t)={cos tcos u, sen t cos u, sen u}

U=linspace(0,2*pi-pi/4,40);
V=linspace(-pi/2,pi/2-pi/4,40);
[u v]=meshgrid(U,V);
surf(cos(u).*cos(v),sin(u).*cos(v),sin(v));
view(60,50);
axis([-1.2 1.2 -1.2 1.2 -1.2 1.2])
1
0.5
0
-0.5
-1
-1
-0.5
0
0.5
1
c) r(u,t)={sen t, cos t, u}
U=linspace(0,2*pi,40);V=linspace(-3,3,20);
surf(sin(u),cos(u),v);
-1
-0.5
0.5
view(-37.5,50);
axis([-2 2 -2 2 -3 3])
2
0
-2
2
1
2
1
-1
-1
-2
-2
d) r(u,t)={u sen t, ucos t, t/3}

U=linspace(0,2*pi,40);V=linspace(-3,3,20);
surf(v.*sin(u),v.*cos(u),u./3);
view(-50,50);
axis([-3.5 3.5 -3.5 3.5 0 2*pi/3])
2
1.5
1
0.5
0
2
0
-2
-3
-2
-1
Actividad Experimental # 2: Dibuja el campo vectorial gradiente y campo direccional de:

2
1.
z=f ( x , y )=x e x y , -1x1, -1y1

[x,y] = meshgrid(-2:.2:2, -2:.2:2);
z = x .* exp(-x.^2 - y.^2);
[px,py] = gradient(z,.2,.2);quiver(px,py)
25
20
15
10
2.
z=f ( x , y )= x2 + y 2 , -1x1, -1y1

[x,y] = meshgrid(-2:.2:2, -2:.2:2);
z = sqrt(x.^2 + y.^2);
[px,py] = gradient(z,.2,.2);quiver(px,py)
10
15
20
25
25
20
15
10
-5
-5
10
15
20
25
a) Representa el campo gradiente

b) Grafica su campo direccional
Actividad Experimental # 3: Utilizando Matlab, dibuja el
E ( r ) radial de:
Si: F(x, y, z) = x2 + y2 1=0, (cilindro circular de eje z y radio 1).

[X,Y] = meshgrid(cylinder(5,50));
U = X^2;
V = Y^2;
W = 0.;
subplot(1,2,1)
for z = [-1,0,1]
Z = z +0*X;
quiver3(X,Y,Z,U,V,W)
hold on
end
axis image
Actividad Experimental # 4: El campo de una onda plana uniforme es:
E ( r )=E o ( j ^x ^y + j z^ ) e j 2 (x+ z)
a)
b)
c)
d)
calcula su campo magntico correspondiente.

cul es su frecuencia?
grafique la forma de la onda en (a) con Matlab
determine la forma del vector de poynting
Operadores diferenciales
Para describir el comportamiento de una magnitud vectorial en cada punto del espacio, hay que usar el
concepto de campo vectorial, que es una funcin que asigna un vector a cada punto de R o R3. En la
prctica, expresamos un campo vectorial V en R2 de forma:
V(x,y)=u(x,y)i+v(x,y)j
Las funciones u y v, son funciones escalares llamadas componentes del campo V. para hacerse una
idea visual del campo vectorial, se dibujan vectores V(x,y) en forma de flechas, en puntos seleccionados
del dominio de V; un diagrama de este tipo es la grfica de campo.
La grfica de campo vectorial suministra informacin interesante sobre las propiedades de campo. Con
Matlab, los campos vectoriales se representan con la instruccin quiver. As, para representar campos
vectoriales en R2 se utiliza la instruccin:
>>quiver (x,y,u,v)
(o de forma ms simple: quiver (u,v)), que dibuja los vectores de componentes u y v en los puntos (x,y).
Este comando admite una versin ms completa para definir una escala, especifiacaciones de lnea y
relleno que puede ser consultada ejecutando la sentencia: help quiver.
Ejemplos:
1. Para obtener la grfica del campo vectorial F(x,y)=-yi+xj en el cuadrado [-1,1]x[-1,1], se sigue:
>>[x,y]=meshgrid(-1:.5:1);
>>u=-y;v=x;
>>quiver(u,v), axis square
>>legend(CAMPO CIRCULAR EN R^2)
6
CAMPO CIRCULAR EN R2
5
0
0
2. >> [x,y]=meshgrid(-1:.5:1);u=zeros(size(y));v=x;quiver(u,v)
7
6
5
4
3
2
1
0
-1
0.5
1.5
2.5
3.5
3. >> [x,y]=meshgrid(-1:.5:1);v=zeros(size(x));u=y;quiver(u,v)
4.5
5.5
5.5
5
4.5
4
3.5
3
2.5
2
1.5
1
0.5
-1
En la prctica, un campo vectorial V en R3 de componentes u,v,w se expresa en la forma:

V(x,y,z)=u(x,y,z)i+v(x,y,z)j+w(x,y,z)k
Para representar campos vectoriales en R3 se utiliza la instruccin:
>>quiver3(x,y,z,u,v,w)
Cuyo efecto es dibujar los vectores de componentes u, v, w en los puntos (x,y,z). Como en el caso
bidimensional, admite una versin ms completa para definir una escala, especificaciones de lnea y
relleno (help quiver 3).
Algunos ejemplos comunes de campos vectoriales en Fsica son los campos de velocidades, los
campos de fuerza elctrica y los campos gravitatorios. stos ltimos vienen definidos por la Ley de
Newton y son de la forma:
F ( x , y , z) =
Gm
3
2 2
(x i+ y j+ z k )
( x2 + y2 + z )
Donde m es la masa y G es la constante de gravitacin universal. Mediante el siguiente problema,

vemos que un campo gravitatorio siempre apunta hacia el origen.

Ejercicios Con Matlab-Teoría Electromagnética

Cargado por

Copyright:

Formatos disponibles

Ejercicios Con Matlab-Teoría Electromagnética

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Ejercicios Con Matlab-Teoría Electromagnética

Cargado por

Copyright:

Formatos disponibles

EJERCICIOS

1. Ejercicios: graficar la funciones siguientes:

f ( x , y ) =sen (| x|+| y|)

2. x=-4:.01:4; y=sin(x); plot(x,y), grid, title('Funcin

1. x= -2.9:0.2: 2.9; y=exp(-x.*x); plot(x,y)

5. Determina y grafica los contornos de un campo

x = cos(t), y = sin(t), z = t/(2 )

9. t=linspace(0, 2*pi, 200);

1. x=-8:0.5:8; y=-6:0.5:6; [X,Y]=meshgrid(x,y);

1. x=-8:0.5:8; y=-3:0.5:3; [X,Y]=meshgrid(x,y); Z=(X.^2Y.^2)./(X.^2+Y.^2); surf(X,Y,Z)

2. zz = linspace(0, 2*pi, 26); [x, y, z] =

Actividad Experimental # 1: Grafica varias funciones paramtricas:

b) r(u,t)={cos tcos u, sen t cos u, sen u}

d) r(u,t)={u sen t, ucos t, t/3}

Actividad Experimental # 2: Dibuja el campo vectorial gradiente y campo direccional de:

z=f ( x , y )=x e x y , -1x1, -1y1

z=f ( x , y )= x2 + y 2 , -1x1, -1y1

a) Representa el campo gradiente

Si: F(x, y, z) = x2 + y2 1=0, (cilindro circular de eje z y radio 1).

calcula su campo magntico correspondiente.

En la prctica, un campo vectorial V en R3 de componentes u,v,w se expresa en la forma:

Donde m es la masa y G es la constante de gravitacin universal. Mediante el siguiente problema,

También podría gustarte