Matlab Hogyan - Hidrodinamikai Rendszerek Tanszék
Matlab Hogyan - Hidrodinamikai Rendszerek Tanszék
Matlab Hogyan - Hidrodinamikai Rendszerek Tanszék
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.