17.04.2014 Views

Matlab Hogyan - Hidrodinamikai Rendszerek Tanszék

Matlab Hogyan - Hidrodinamikai Rendszerek Tanszék

Matlab Hogyan - Hidrodinamikai Rendszerek Tanszék

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

9. FEJEZET. KÖZÖNSÉGES DIFFERENCIÁLEGYENLETEK (ODE) MEGOLDÁSA 37<br />

hold on;<br />

end;<br />

hold off;<br />

title([’Adatok: dp/dx=’,num2str(P),’ \nu=’,num2str(nu)],’FontSize’,14);<br />

xlabel(’aramlasi sebesseg, v’,’FontSize’,14);<br />

ylabel(’keresztemtszet, x’,’FontSize’,14);<br />

%---------------------------------------------<br />

function dydx=mat4ode(x,y)<br />

global P nu<br />

dydx=[ y(2); P/nu];<br />

%----------------------------------------------<br />

function res=mat4bc(ya,yb)<br />

global v0<br />

res=[ ya(1); yb(1)-v0];<br />

%---------------------------------------------<br />

function yinit=mat4init(x)<br />

global v0<br />

yinit = [ v0*x^2; 2*v0*x ];<br />

couette.eps<br />

9.4. ábra. Sebességeloszlás Couette áramlás esetén nyomásgradienssel, τ csúsztatófeszültség szerint színezve.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!