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

86 阅读 0 评论 57 点赞

我是靠谱客的博主等待棒棒糖，最近开发中收集的这篇文章主要介绍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);

如果觉得靠谱客网站的内容还不错，欢迎将靠谱客网站推荐给程序员好友。

本图文内容来源于网友提供，作为学习参考使用，或来自网络收集整理，版权属于原作者所有。