Resolução de EDP com pdepe no Matlab

Prévia do material em texto

8 ) Resolva a EDP utilizando o comando pdepe do Matlab:
Condição inicial:
Condições de contorno:
	Solução: Inicialmente, precisamos escrever 3 “function”. Uma para a EDP, uma para a condição inicial e uma para as condições de contorno. Assim:
No script 
function ped
 
m = 0;
x = linspace(0,0.1,200);
t = linspace(0,10,50);
 
sol = pdepe(m,@pdex1pde,@pdex1ic,@pdex1bc,x,t);
Temperature = sol(:,:,1);
 
 %Perfil de Temperatura
 figure(1)
 plot(x,Temperature(end, :))
 title(strcat('Solução em t =', num2str(tend)))
 xlabel('Distância X')
 ylabel('Temperatura (ºC)')
 grid on
 
 %Variação da Temperatura com o tempo
 figure(2)
 plot(t,Temperature(:,1))
 title('Temperatura na Superfície')
 xlabel('Tempo(s)')
 ylabel('Temperatura(ºC)')
 grid on
 
 % --------------------------------------------------------------
 function [c,f,s] = pdex1pde(x,t,u,DuDx)
c = 10000*500;
f = 200*DuDx;
s = 0;
 % --------------------------------------------------------------
 function u0 = pdex1ic(x)
u0 = 0;
 % --------------------------------------------------------------
 function [pl,ql,pr,qr] = pdex1bc(xl,ul,xr,ur,t)
pl = 6;
ql = l;
pr = ur;
qr = 0;