基于MATLAB的有限元法求解EIT技术的正问题

有限元法（Finite Element Method）是基于!近代计算机的快速发展而发展起来的一种近似数值方法, 用来解决力学，数学中的带有特定边界条件的偏微分方程问题（PDE）。而这些偏微分方程是工程实践中常见的固体力学和流体力学问题的基础。有限元和计算机发展共同构成了现代计算力学（Computational Mechanics）的基础。有限元法的核心思想是“数值近似”和“离散化”，所以它在历史上的发展也是围绕着这两个点进行的。
作者：菲兹
链接：https://www.zhihu.com/question/27696855/answer/43989189
常常需要求解各类微分方程，而许多微分方程的解析解一般很难得到，使用有限元法将微分方程离散化后，可以编制程序，使用计算机辅助求解。有限元法在早期是以变分原理为基础发展起来的，所以它广泛地应用于以拉普拉斯方程和泊松方程所描述的各类物理场中（这类场与泛函的极值问题有着紧密的联系）。
求解EIT技术的正问题背景：圆形被测物，16个电极排列。

%利用网格划分工具箱将圆形分为多个三角形，并标记节点坐标与三角形单元编号。保存为txt文件为tri_in1.txt。
网格划分工具箱请看链接
这里写链接内容
然后读取文本的信息，%总结点数,单元数,受约束点数%COORD结点坐标数组 LNODS 单元定义数组通过ASSEMBLE函数ASTIF=ASSEMBLE(NPOIN,NELEM,YOUNG,POISS,THICK,COORD,LNODS)调用单元刚度矩阵函数[ESTIF,SMATX]=FORMESTIF(i,COORD,LNODS)%i为当前所计算的单元号获得刚度矩阵的集成。

%变量说明2008.11.10
%NPOIN NELEM NVFIX
%总结点数,单元数,受约束点数
%COORD结点坐标数组 LNODS 单元定义数组 %FORCE?结点力数组(n,3) n:受力结点数,(n,1):作用结点号,(n,2):x方向,(n,3):y方向
%FIXED?约束信息数组(n,3) n:受约束结点数, (n,1):约束结点号，%(n,2)与(n,3):约束点x方向和y方向的约束情况,受约束为1,否则为0
%ASTIF总刚度矩阵,ASLOD总荷载向量
% format short e %设定输出类型
clear
FP1=fopen('tri_in1.txt','rt');
%打开数据文件 FP1数据文件指针
%读入初始数据
NPOIN=fscanf(FP1,'%d',1);
NELEM=fscanf(FP1,'%d',1);NFORCE=fscanf(FP1,'%d',1);NVFIX=fscanf(FP1,'%d',1);
YOUNG=fscanf(FP1,'%e',1);POISS=fscanf(FP1,'%f',1);THICK=fscanf(FP1,'%f',1);
LNODS=fscanf(FP1,'%f',[3,NELEM])';%单元定义数组
COORD=fscanf(FP1,'%f',[2,NPOIN])';%结点坐标数组
FORCE=fscanf(FP1,'%d',[3,NFORCE])';%结点力数组
FIXED=fscanf(FP1,'%d',[3,inf])';%?约束数组
ASTIF=ASSEMBLE(NPOIN,NELEM,YOUNG,POISS,THICK,COORD,LNODS);%?生成总刚
m=fix(NPOIN/4);%m输入电流节点
n=fix(NPOIN/3);c=fix(NPOIN/2);%n输出电流节点
ASTIF(1:c-1,c)=0;
ASTIF(c+1,c)=0;
ASTIF(c,c)=1;%c参考电势为0
b=zeros(NPOIN,1);
b(m,1)=1;
b(n,1)=-1;
[Q1,Q2,Q3,Q4]=gaus(ASTIF,b);
scatter(COORD(:,1),COORD(:,2),5,Q4);%散点图
[X,Y,Z]=griddata(COORD(:,1),COORD(:,2),Q4,linspace(-1,1,100)',linspace(-1,1,100),'v4');
figure(1)
contourf(X,Y,Z,20); %等值线图
% figure(2)
% pcolor(X,Y,Z);shading interp%伪彩色图
figure(3)
surf(X,Y,Z);%三维曲面
figure(4)
scatter(COORD(:,1),COORD(:,2),5,Q4);%散点图
fclose(FP1);%?关闭文件

tri_in1.txt文档可以自己通过网格划分获取。
图像显示：

具体代码请留言！互相交流，或者点击链接
这里写链接内容

基于MATLAB的有限元法求解EIT技术的正问题相关推荐

【故障检测问题】基于matlab免疫算法求解故障检测问题【含Matlab源码 196期】
一.获取代码方式获取代码方式1: 完整代码已上传我的资源:[故障检测问题]基于matlab免疫算法求解故障检测问题[含Matlab源码 196期] 获取代码方式2: 通过订阅紫极神光博客付费专栏,凭 ...
10种基于MATLAB的方程组求解方法
线性方程组的求解包括直接法和迭代法,其中迭代法包括传统的高斯消元法,最速下降法,牛顿法,雅克比迭代法,共轭梯度法,以及智能启发式算法求解法和神经网络学习算法,传统算法可以相互组合改进,智能仿生启发式算 ...
matlab龙格库塔法求通解,基于matlab及龙格库塔法求解布拉修斯方程.doc
基于matlab及龙格库塔法求解布拉修斯方程 Runge-Kutta法求解布拉修斯解摘要薄剪切层方程主要有三种解法,即相似解,非相似条件下对偏微分方程组的数值解和近似解.布拉修斯解是布拉修斯于19 ...
【MVO TSP】基于matlab灰狼算法求解旅行商问题【含Matlab源码 1327期】
一.获取代码方式获取代码方式1: 完整代码已上传我的资源:[TSP]基于matlab灰狼算法求解旅行商问题[含Matlab源码 1327期] 点击上面蓝色字体,直接付费下载,即可. 获取代码方式2: ...
基于MATLAB的运动模糊图像恢复技术
基于MATLAB的运动模糊图像恢复技术摘要:MATLAB是当今流行的科学计算软件,它具有很强的数据处理能力.在其图像处理工具箱中有四个图像复原函数,本文就这些函数的算法原理.运用和恢复处理效果结合实 ...
【背包问题】基于matlab禁忌搜索算法求解背包问题【含Matlab源码 373期】
一.获取代码方式获取代码方式1: 完整代码已上传我的资源:[背包问题]基于matlab禁忌搜索算法求解背包问题[含Matlab源码 373期] 获取代码方式2: 通过订阅紫极神光博客付费专栏,凭支付 ...
【优化求解】基于matlab禁忌搜索算法求解函数极值问题【含Matlab源码 1204期】
⛄一.获取代码方式获取代码方式1: 完整代码已上传我的资源: [优化求解]基于matlab禁忌搜索算法求解函数极值问题[含Matlab源码 1204期] 点击上面蓝色字体,直接付费下载,即可. 获取 ...
【TS TSP】基于matlab禁忌搜索算法求解31城市旅行商问题【含Matlab源码 1143期】
一.获取代码方式获取代码方式1: 完整代码已上传我的资源:[TSP]基于matlab禁忌搜索算法求解31城市旅行商问题[含Matlab源码 1143期] 点击上面蓝色字体,直接付费下载,即可. 获取 ...
【单目标优化求解】基于matlab黑猩猩算法求解单目标问题【含Matlab源码 1413期】
一.获取代码方式获取代码方式1: 完整代码已上传我的资源:[单目标优化求解]基于matlab黑猩猩算法求解单目标问题[含Matlab源码 1413期] 点击上面蓝色字体,直接付费下载,即可. 获取代 ...

基于MATLAB的有限元法求解EIT技术的正问题

基于MATLAB的有限元法求解EIT技术的正问题相关推荐

最新文章

热门文章