反应扩散模型matlab程序

以下是基于Lengyel-Epstein模型的MATLAB程序，用于生成图案图像（2D）：

% 设置参数
a = 2.8e-4;
b = 5e-3;
tau = 5;
K = -0.005;

% 初始化矩阵
n = 200;
U = zeros(n,n);
V = zeros(n,n);

% 设置初始条件
U(:,:) = 1;
V(:,:) = 0.5;

% 迭代计算
for i = 1:1000
    % 计算Laplacian
    lapU = del2(U);
    lapV = del2(V);
    
    % 计算反应扩散方程
    dUdt = lapU + a*(1-U).*V.^2 - U + K;
    dVdt = lapV + b*(1-U).*V - V;
    
    % 更新矩阵
    U = U + tau*dUdt;
    V = V + tau*dVdt;
end

% 绘制图像
figure;
subplot(1,2,1);
imagesc(U);
title('U');
subplot(1,2,2);
imagesc(V);
title('V');

该程序使用Lengyel-Epstein模型生成图案图像，其中包含反应扩散方程。程序中的参数可以根据需要进行修改，例如可以调整初始条件和迭代次数来生成不同的图案。