UNIVERSIDAD TCNICA DE AMBATOFACULTAD DE SISTEMAS, ELECTRNICA E INDUSTRIALCarrera de Ingeniera Electrnica y ComunicacionesTEORIA II

INTEGRANTES: Cunalata WalterNIVEL: SextoFECHA: 20/05/2015TEMA: MATLAB CAMPO ELCTRICO Y MAGNTICO

AMBATO - ECUADOR

Cdigo campo elctrico

close all; clear; x = -1:0.1:1; y = -1:0.1:1; [xx,yy] = meshgrid(x,y); px = exp(-(xx.^2+yy.^2)); py = sin(xx.*yy); figure;quiver(x,y,px,py,2); xlabel('x');ylabel('y'); title('2-D vector field: {\bf P} = e^{-(x^2+y^2)}{\bf x} + sin(xy) {\bf y} '); pause; [pxx,pxy]=gradient(px,0.1,0.1); [pyx,pyy]=gradient(py,0.1,0.1); curl_p_z=(pyx-pxy); [m,n]=size(curl_p_z); curl_p_x=zeros(m,n); curl_p_y=zeros(m,n); figure;quiver3(xx,yy,zeros(m,n),curl_p_x,curl_p_y,curl_p_z,4,'r'); hold on; surf(xx,yy,zeros(m,n));shading flat;colormap pink; rotate3d; xlabel('x');ylabel('y'); zlabel('z'); title('Curl of 2-D vector field: \nabla \times {\bf P}'); pause; figure;surf(x,y,curl_p_z); rotate3d; xlabel('x');ylabel('y'); zlabel('|\nabla \times {\bf P}|'); title('Magnitude of Curl of 2-D vector field: |\nabla \times {\bf P}|');

Campo magntico

%Vector fields w and wo sources/sinksclose all; clear;x = -2:.2:2; y = -1:.2:1;[xx,yy] = meshgrid(x,y);v = xx.*exp(-xx.^2-yy.^2);[ex,ey] = gradient(v,.2,.2);figure(1);quiver(x,y,ex,ey,2);xlabel('x');ylabel('y');title('2-D vector field: {\bf E} = (1-2x^2)e^{-(x^2+y^2)} {\bf x} - 2xye^{-(x^2+y^2)} {\bf y}');pause;div_e=del2(v);figure(2);surf(x,y,div_e);shading interp;rotate3d;xlabel('x');ylabel('y');title('Divergence of 2-D vector field: \nabla \cdot {\bf E}');pause;x = -2:0.2:2; y = -2:0.2:2;[xx,yy] = meshgrid(x,y);bx = -yy./(xx.^2+yy.^2);by = xx./(xx.^2+yy.^2);figure;quiver(x,y,bx,by);xlabel('x');ylabel('y');title('2-D vector field: {\bf B} = -y /(x^2+y^2){\bf x} + x /(x^2+y^2) {\bf y} ');pause;[bxx,bxy]=gradient(bx,0.2,0.2);[byx,byy]=gradient(by,0.2,0.2);div_b=(bxx+byy);figure;surf(x,y,div_b);shading interp;rotate3d;xlabel('x');ylabel('y');title('Divergence of 2-D vector field: \nabla \cdot {\bf B}: {\bf Must be Zero!!!}');