利用matlab进行EIT正问题有限元求解

对EIT正问题中节点电势采用有限元方法进行求解，为加深求解方法。

%总刚度矩阵
function ASTIF=ASSEMBLE（NPOINNELEMNJIEZHICOORDLNODS）

ASTIF（1:NPOIN1:NPOIN）=0;%生成特定大小的总刚度矩阵，并赋值为0.总刚度矩阵，为节点数*节点数，代表了每一个节点与另一个节点间的关系矩阵
% 组装总刚度矩阵
for i=1:NELEM
    % 计算单刚
    [ESTIFSMATX]=FORMESTIF（iNJIEZHICOORDLNODS）;
    nt=LNODS（i:）; % 临时向量 nt用来记录当前单元的结点号，三个节点
    for j=1:3
        for k=1:3   %两个节点互相交错行列所对应的节点间的关系。
            ASTIF（nt（j）nt（k））=...
            ASTIF（nt（j）nt（k））+...%叠加刚度矩阵
            ESTIF（jk）;
            % 利用结点位移编号对应关系将单元刚度分块叠加到总刚矩阵
        end
    end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        660  2018-05-09 15:14  fem-eit\ASSEMBLE.m

     文件       8016  2018-05-09 15:09  fem-eit\data.txt

     文件       1149  2018-05-09 15:13  fem-eit\FORMESTIF.m

     文件        924  2018-05-09 15:09  fem-eit\gaus.m

     文件       1494  2018-05-09 15:12  fem-eit\main1.m

     目录          0  2018-05-09 15:17  fem-eit

----------- ---------  ---------- -----  ----

                12243                    6