Я пытаюсь решить двумерное уравнение диффузии в 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);