Estudante

PTC 313
Eletromagnetismo

Trabalho computacional 1

Professor: Plinio

A. Listagem do programa utilizado.
Utilizamos para a resolução desse exercício programa o Matlab. Segue o código desenvolvido.

clear all close all a=0.11; b=0.07; g=0.03; c=0.05; d=b-0.03; h=(b-d)/2; esp=1;
V0=100;

eps=1.9*8.854187817620e-12; sigma=3.2e-3; passo=1e-3; N=a/passo;
M=b/passo;

ilim1 = cast((b-h-d)/passo, 'uint16'); ilim2 = cast((b-h)/passo, 'uint16'); jlim1 = cast( g/passo, 'uint16'); jlim2 = cast ((g+c)/passo, 'uint16'); V=zeros(M,N); for j= jlim1 : jlim2 for i= ilim1 : ilim2 V(i,j)=V0; end end % V(:,0)=0;
% V(:,N)=0;
% V(0,:)=0;
% V(M,:)=0; maxvar=V0; while(maxvar >= 1e-6*V0) maxvar=0; for i=2:M-1 for j=2:N-1 if ( i< ilim1 | i> ilim2 | j < jlim1 | j > jlim2 ) Vnew = ( V(i+1,j)+V(i,j+1)+V(i-1,j)+V(i,j-1) )/4; if( abs(Vnew- V(i,j)) > maxvar ) maxvar = abs(Vnew- V(i,j)) ; end V(i,j) = Vnew; end end end end %Calculo para o problema dual
Vdual = zeros(M,N); posicaoV0 = cast((b/2)/passo, 'uint16'); for j= 1 : jlim1 Vdual(posicaoV0,j)=V0; end for j = jlim2 : N Vdual(posicaoV0,j) = 0; end maxvar = V0; while(maxvar >= 1e-6*V0) maxvar=0; for i= 1:M for j=1:N if i ~= posicaoV0 if ( i= ilim2 || j = jlim2 ) if i == 1 && j == 1 Vnew = (2*Vdual(i+1,j) + 2*Vdual(i,j+1))/4; elseif i==1 && j==N Vnew = (2*Vdual(i,j-1) + 2*Vdual(i+1,j))/4; elseif i==M && j==1 Vnew = (2*Vdual(i-1,j) + 2*Vdual(i,j+1))/4; elseif i==M && j==N Vnew = (2*Vdual(i-1,j) + 2*Vdual(i,j-1))/4;