matlab解坐标方程,matlab程序(解泊松方程)

求解泊松方程的

function Finite_element_tri(Imax)

% 用有限元法求解三角形形区域上的Possion方程

Jmax=2*Imax;

% 其中Imax Jmax分别表示x轴和y轴方向的网格数，其中Jmax等于Imax的两倍

% 定义一些全局变量

global ndm nel na

% ndm 总节点数

@

% nel 基元数

% na 表示活动节点数

V=0; J=0;X0=1/Imax;Y0=X0;%V=0为边界条件

domain_tri % 调用函数画求解区域

[X,Y,NN,NE]=setelm_tri(Imax,Jmax); % 给节点和三角形元素编号，并设定节点坐标

% 以下求解有限元方程的求系数矩阵

T=zeros(ndm,ndm);

for n=1:nel

…

n1=NE(1,n); n2=NE(2,n); n3=NE(3,n);%整体编号

s=abs((X(n2)-X(n1))*(Y(n3)-Y(n1))-(X(n3)-X(n1))*(Y(n2)-Y(n1)))/2;%三角形面积

for k=1:3

if n1<=na|n2<=na

T(n1,n2)=T(n1,n2)+((Y(n2)-Y(n3))*(Y(n3)-Y(n1))+(X(n3)-X(n2))*(X(n1)-X(n3)))/(4* s);