Poissons Equation with Point source

2 次查看（过去 30 天）

Hannah Burmester 2024-2-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2079316-poissons-equation-with-point-source

编辑： Dyuman Joshi 2024-2-17

Hello everyone :)

in matlab´s Example "Poisson's Equation with Point Source and Adaptive Mesh Refinement" ist a function "circlef" mentioned, that create the point source by returning 1/area for the triangle that contains the origin and zero elsewhere. How do I create such a funtion? Meaning how can I define f differently on specific traingles in the mesh?

Thanks in advance

Hannah

2 个评论
显示无隐藏无

Dyuman Joshi 2024-2-7

"How do I create such a funtion?"

Use the description given as its definition.

"Meaning how can I define f differently on specific traingles in the mesh?"

Use if, elseif, else conditions.

Hannah Burmester 2024-2-7

Thanks you for your reply. Yes I thought of something like

if (0,0) in triangle

f = 1/area(triangle)

else

f = 0

end

But how can I translate that into code?

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Dyuman Joshi 2024-2-7

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2079316-poissons-equation-with-point-source#answer_1404456

在 MATLAB Online 中打开

You can use inpolygon to see if the origin is inside the triangle or not and polyarea to calculate the area of the triangle -

Note - inpolygon() also considers the points that lie ON the edges of the polygon (see the case below).

x = [-2 -1 1];

y = [-2 1 -1];

plot(polyshape(x,y))

hold on

plot(0, 0, '.r', 'MarkerSize', 10)

if inpolygon(0, 0, x, y)

f = 1./polyarea(x,y);

else

f = 0;

end

f = 0.2500

This can be condensed to -

F = inpolygon(0, 0, x, y)./polyarea(x,y)
F = 0.2500

5 个评论
显示 3更早的评论隐藏 3更早的评论

Hannah Burmester 2024-2-12

Thank you all for your replays! Editing the matlab function should work.

Dyuman Joshi 2024-2-12

编辑：Dyuman Joshi 2024-2-17

在 MATLAB Online 中打开

From a cursory glance, I'd say my solution is simpler than whatever is being done here -

type circlef.m
function f=circlef(p,t,u,time)

%point source at (0,0)

% Copyright 1994-2001 The MathWorks, Inc.

x=0;
y=0;

np=size(p,2);
nt=size(t,2);

[ar,t1,t2,t3]=pdetrg(p,t);

[t1,tn,t2,t3]=tri2grid(p,t,zeros(np,1),x,y);

f=zeros(1,nt);

if ~isnan(tn)
  f(tn)=1/ar(tn);
end

请先登录，再进行评论。

请先登录，再回答此问题。