请用matlab写一个画出绕圆柱无环量流动流线分布图的程序

以下是一个用MATLAB画出绕圆柱无环量流动流线分布图的程序示例：

% 设置流场参数
Uinf = 1; % 自由流速度
R = 1; % 圆柱半径
L = 10; % 流场长度
N = 50; % 离散化网格数
x = linspace(-L/2, L/2, N); % X方向网格坐标
y = linspace(-L/2, L/2, N); % Y方向网格坐标
[X, Y] = meshgrid(x, y); % 构建网格坐标矩阵

% 计算流场速度分量
r = sqrt(X.^2 + Y.^2); % 极径
theta = atan2(Y, X); % 极角
Vr = Uinf*(1 - R^2./r.^2); % 径向速度分量
Vt = -Uinf*R^2./r.^2.*sin(theta); % 切向速度分量
Vx = Vr.*cos(theta) - Vt.*sin(theta); % X方向速度分量
Vy = Vr.*sin(theta) + Vt.*cos(theta); % Y方向速度分量

% 画流线图
figure;
streamline(X, Y, Vx, Vy);
xlim([-L/2, L/2]);
ylim([-L/2, L/2]);
title('无环量圆柱绕流流线图');
xlabel('X方向');
ylabel('Y方向');

运行此程序，即可得到绕圆柱无环量流动流线分布图。需要注意的是，此程序只是一个简单的示例，实际的流场模拟需要更加复杂的处理和计算方法。