给出一份基于matlab的计算绕圆柱有环量流动的程序代码

以下是一个基于MATLAB的计算绕圆柱有环量流动的程序代码，可以通过修改相应参数来适应不同的情况：

clear all; clc;

% 给出参数
U = 1;      % 无量纲速度
D = 1;      % 圆柱直径
nu = 0.1;   % 运动粘性系数
L = 10;     % 计算区域边长
N = 100;    % 离散化网格数
dx = L/N;   % 网格大小

% 生成网格
x = linspace(-L/2, L/2, N+1);
x = x(1:N);
y = x;
[X, Y] = meshgrid(x, y);

% 初始化速度势和流函数
phi = zeros(N, N); 
psi = zeros(N, N);

% 定义圆柱的位置和半径
xc = 0; yc = 0; r = D/2;

% 计算速度势和流函数
for i = 1:N
    for j = 1:N
        if sqrt((X(i,j)-xc)^2 + (Y(i,j)-yc)^2) < r
            phi(i,j) = U*(r^2/2 - (X(i,j)-xc)^2/2 - (Y(i,j)-yc)^2/2);
            psi(i,j) = U*(Y(i,j)-yc);
        else
            phi(i,j) = U*(r^2/2)*log(sqrt((X(i,j)-xc)^2 + (Y(i,j)-yc)^2)/r);
            psi(i,j) = U*atan2(Y(i,j)-yc, X(i,j)-xc);
        end
    end
end

% 计算绕圆柱有环量
Gamma = -2*pi*r*U/nu;

% 输出结果
fprintf('绕圆柱的有环量为: %f\n', Gamma);