کد متلب رسم خطوط جریان و میدان سرعت


تمرین1:
رسم دسته خطوط Ψ و میدان سرعت زیر:
q =(x^2-y^2)i ⃑ -(2xy)j ⃗

[x,y]=meshgrid(0:0.1:1,0:0.1:1)
u=inline('x.^2-y.^2','x','y');
v=inline('-2.*x.*y','x','y');
U=u(x,y);
V=v(x,y);
quiver(x,y,U,V,'k')
title 'velocity field'
xlabel 'x',ylabel 'y'
axis image

[x,y]=meshgrid(0:0.1:1,0:0.1:1)
u=inline('x.^2-y.^2','x','y');
v=inline('-2.*x.*y','x','y');
U=u(x,y);
V=v(x,y);
quiver(x,y,U,V,'k')
xlabel 'x',ylabel 'y'
axis image
hold on
psi=(x.^2).*y-(y.^3)./3;
contour(x,y,psi,10,'k')
title 'velocity field & streamlines'

[x,y]=meshgrid(0:0.1:1,0:0.1:1)
z=(x.^2).*y-(y.^3)./3
h= contour(x,y,z)
clabel(h)
xlabel 'x',ylabel 'y'
title 'psi function for different c'