15. Metoda konÄnih elementov in reÅ¡evanje BDE in PDE v Matlabu
15. Metoda konÄnih elementov in reÅ¡evanje BDE in PDE v Matlabu
15. Metoda konÄnih elementov in reÅ¡evanje BDE in PDE v Matlabu
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Matlab implementacija za zgled 2<br />
function [c,f,s] = sisde(x,t,u,ux)<br />
c = [ 1; 1 ];<br />
f = ux/2;<br />
s = [ 1/(1+u(2)^2); 1/(1+u(1)^2) ];<br />
function [pa,qa,pb,qb] = sisbc(xl,ul,xr,ur,t)<br />
pa = [0; 0];<br />
qa = [1; 1];<br />
pb = [0; 0];<br />
qb = [1; 1];<br />
function y = sis<strong>in</strong>i(x)<br />
y = [ 1+0.5*cos(2*pi*x); 1-0.5*cos(2*pi*x)];<br />
m = 0;<br />
x = l<strong>in</strong>space(0,1,20);<br />
t = l<strong>in</strong>space(0,0.2,10);<br />
u = pdepe(m,@sispde,@sis<strong>in</strong>i,@sisbc,x,t);<br />
u1 = u(:,:,1); u2 = u(:,:,2);<br />
sublot(221); surf(x,t,u1);<br />
sublot(222); surf(x,t,u2);<br />
Bor Plestenjak - Numerična analiza 2005/06