我正在尝试在Matlab中求解二维扩散方程。我的几何图形是一维和一维的二维网格(大小为1500X750),其中一维表示源(代码中的全局ProMat是1500X750的矩阵)。
现在,我正在使用Matlab偏微分方程(PDE)工具箱,并且模拟非常慢。问题似乎出在用于创建几何的decsg()
函数中。我想提高效率。
% Set the coefficients in a 3-D rectangular block geometry.
N1=size(ProMat,1);
N2=size(ProMat,2);
model = createpde();
R1 = [3 4 0 0 N1 N1 0 N2 N2 0]';
nm(1:2,1)='R1';
% The problem is here
g = decsg(R1,'R1',nm);
geometryFromEdges(model,g);
%fff = @(region,state)(region.x.^2);