MATLAB利用streamslice()函数绘制流线

1、第一，想要绘制的图形如下图：流场绕过半圆形物体时的流线图。

3、第三，在上述代码的基础上进行了改进，给出如下MATLAB代码，该段代码可以运行出流场绕过半癣嗡赧箬圆形物体时的流线。c造婷用痃lose all; clear all; clc[x,y] = meshgrid(-20:0.1:20,0:0.1:10);r = sqrt(x.*x+y.*y);costheta=x./r; sintheta=y./r; d=10./2./r; ur=-0.5*costheta.*(d.^3-d*3+2);utheta=-0.5*sintheta.*(d.^3-d*3+4);u=-ur.*costheta-utheta.*sintheta;v=-ur.*sintheta+utheta.*costheta;figure('Position',[100 100 800 300]);hold onstreamslice(x,y,u,v);fai = 0:0.1:1.1*pi; R=5; x=R*cos(fai);y=R*sin(fai);% 圆的参数方程fill(x,y,[1,1,1]); % 填充式画半圆plot(R*cos(fai),R*sin(fai),'k-','LineWidth',3); % 线条式画半圆axis([-20,20,0,10]);axis fill;box onset(gca,'TickDir','in','XTick',[-20:10:20],'YTick',[0:2:10])

5、第五，如果把代码中的fill(x,y,[1,1,1])注释掉或者删除掉，得到图像如下，就会发现流场穿过半圆形物体，这显然不合适，所以fill()这个用来填充颜色（这里填充的是白色）的函数很有用。