Я пытаюсь решить переменную коэффициента уравнения Пуассона в 2-й прямоугольной области с граничными условиями Дирихле.Все работает нормально, указывая коэффициент c в качестве дескриптора функции, например:
C=@(location,state)abs(cos(location.y/2)+sin(location.x/2));
F=C;
specifyCoefficients(model,'m',0,...
'd',0,...
'c',C,...
'a',0,...
'f',F);
В конце я интерполирую решение на сетке сетки:
results = solvepde(model);
Nx=51;
Ny=21;
xq = linspace(a,b,Nx);
yq = linspace(c,d,Ny);
[X,Y] = meshgrid(xq,yq);
Z = interpolateSolution(results,X,Y);
Z = reshape(Z,size(X));
Теперь у меня есть толькоЗначения C в узлах этой сетки сетки.Есть ли способ поместить эту матрицу в функцию для использования в качестве ввода для specifyCoefficients
?