Esercizi di elettrostatica 
Calcolo del campo elettrico prodotto da cariche isolate

Caso della carica singola
 
 

Caso di due cariche di segno uguale

Caso di due cariche di segno opposto

programma per il calcolo delle componenti della carica isolata

% Paolo Villoresi DEI Universita' di Padova 1-3-2001
% esercitazioni di FISICA GENERALE II per Ingegneria delle Telecomunicazioni

%calcolo delle componenti del campo elettrico nel piano prodotto da una carica isolata

[x,y]=meshgrid(-3:.1:3);
eps0=8.85e-12;
e=1.61e-19; 
q=1e8 *e;       %considero 10^8 cariche elementari poste in P0=(1.2,2.3)
x0=0.25;
y0=0.25;
delta=1e-2; %parametro per evitare la divergenza della funzione
k=1/(4* pi*eps0);
ex=k*q./((x-x0).^2+(y-y0).^2+delta).^1.5 .*(x-x0);
ey=k*q./((x-x0).^2+(y-y0).^2+delta).^1.5 .*(y-y0);

figure(1)
subplot(121)
%mesh(x,y,ex)
meshc(x,y,ex)
%contour(x,y,ex)
axis equal
xlabel('asse X')
ylabel('asse Y')
title('componente Ex')
grid
subplot(122)
%mesh(x,y,ey)
meshc(x,y,ey)
xlabel('asse X')
ylabel('asse Y')
title('componente Ey')
grid

figure(2)
quiver(x,y,ex,ey,6)
axis image

programma per il calcolo delle componenti di due carice dello stesso segno

% Paolo Villoresi DEI Universita' di Padova 1-3-2001
% esercitazioni di FISICA GENERALE II per Ingegneria delle Telecomunicazioni

%calcolo delle componenti del campo elettrico nel piano prodotto da due cariche 
%di eguale segno 

[x,y]=meshgrid(-4:.5:4);
delta=1e-2; %parametro per evitare la divergenza della funzione
eps0=8.85e-12;
e=1.61e-19;
q=1e8 *e;       %considero 10^8 cariche elementari poste in P0=(1.2,2.3)
x1=1.5;
y1=0;
x2=-1.5;
y2=0;
k=1/(4* pi*eps0);
ex1=k*q./((x-x1).^2+(y-y1).^2 +delta).^1.5 .*(x-x1);
ey1=k*q./((x-x1).^2+(y-y1).^2+delta).^1.5 .*(y-y1);
ex2=k*q./((x-x2).^2+(y-y2).^2+delta).^1.5 .*(x-x2);
ey2=k*q./((x-x2).^2+(y-y2).^2+delta).^1.5 .*(y-y2);

figure(1)
subplot(121)
%mesh(x,y,ex)
contour(x,y,ex1+ex2,100)
axis equal
xlabel('asse X')
ylabel('asse Y')
title('componente Ex')
grid
subplot(122)
%mesh(x,y,ey)
%surf(x,y,ey1+ey2)
contour(x,y,ey1+ey2,100)
axis equal
xlabel('asse X')
ylabel('asse Y')
title('componente Ey')
grid

figure(2)
quiver(x,y,ex1+ex2,ey1+ey2,8)
axis image

figure(3)
subplot(121)
%mesh(x,y,ex)
surf(x,y,ex1+ex2)
%axis equal
xlabel('asse X')
ylabel('asse Y')
title('componente Ex')
grid
subplot(122)
%mesh(x,y,ey)
%surf(x,y,ey1+ey2)
surf(x,y,ey1+ey2)
%axis equal
xlabel('asse X')
ylabel('asse Y')
title('componente Ey')
grid

programma per il calcolo delle componenti di due cariche di segno opposto

% Paolo Villoresi DEI Universita' di Padova 1-3-2001
% esercitazioni di FISICA GENERALE II per Ingegneria delle Telecomunicazioni

%calcolo delle componenti del campo elettrico nel piano prodotto da due cariche 
%di segno opposto

[x,y]=meshgrid(-4:.5:4);
delta=1e-2; %parametro per evitare la divergenza della funzione
eps0=8.85e-12;
e=1.61e-19;
q=1e8 *e;       %considero 10^8 cariche elementari poste in P0=(1.2,2.3)
x1=1.5;
y1=0;
x2=-1.5;
y2=0;
k=1/(4* pi*eps0);
ex1=k*q./((x-x1).^2+(y-y1).^2 +delta).^1.5 .*(x-x1);
ey1=k*q./((x-x1).^2+(y-y1).^2+delta).^1.5 .*(y-y1);
ex2=-k*q./((x-x2).^2+(y-y2).^2+delta).^1.5 .*(x-x2);
ey2=-k*q./((x-x2).^2+(y-y2).^2+delta).^1.5 .*(y-y2);


figure(1)
subplot(121)
%mesh(x,y,ex)
contour(x,y,ex1+ex2,100)
axis equal
xlabel('asse X')
ylabel('asse Y')
title('componente Ex')
grid
subplot(122)
%mesh(x,y,ey)
%surf(x,y,ey1+ey2)
contour(x,y,ey1+ey2,100)
axis equal
xlabel('asse X')
ylabel('asse Y')
title('componente Ey')
grid

figure(2)
quiver(x,y,ex1+ex2,ey1+ey2,6)
axis image

figure(3)
subplot(121)
%mesh(x,y,ex)
surf(x,y,ex1+ex2)
%axis equal
xlabel('asse X')
ylabel('asse Y')
title('componente Ex')
grid
subplot(122)
%mesh(x,y,ey)
%surf(x,y,ey1+ey2)
surf(x,y,ey1+ey2)
%axis equal
xlabel('asse X')
ylabel('asse Y')
title('componente Ey')
grid

Contatto:
dr. Paolo Villoresi,
Dipartimento di Elettronica e Informatica,
stanza 236, via Gradenigo 6a
tel. 049 827 7644 - fax. 049 827 7699
pavil@dei.unipd.it