所谓有限元就是利用有限维解决无限维问题,既然提到有限维空间,我们总能找到一组基函数的线性组合来表示这个空间内的任一元素。

对于分片一次多项式空间内的基函数该是怎样的呢?下面就来讲述基函数构造问题。

对于区间[x1,xn+1]上进行剖分(均匀或者非均匀剖分),然后构造如下图的分段函数:

其中Ti(x)为分段基函数,分别在区间[xi-1,xi]和[xi,xi+1]上都为一次函数.

其中,基函数分段函数表达形式如下:

有兴趣的同学可以证明这组基函数是线性无关的,并且可以线性张成这个有限维分片一次多项式空间,不妨记这个空间为Vh,无限维的分片一次多项式空间V,可以看出Vh是V的子空间。对于任意uh(x)在Vh空间中,uh(x)都可以被{Ti(x)}线性表示,不妨记为如下:

其中ui是常数.我们利用上述基函数来解决如下泊松问题:

解上述线性方程组即可.为了求解方便,我们采用均匀剖分的情况,当然非均匀的采用遍历法即可。假设xi+1-xi=h(固定常数)。

现在的问题就是如何算这些内积,因为基函数是分段的,所以必须依赖于它们各自的区间进行。

注意,一个基函数最多与左右相邻的基函数有关系,其它的基函数在此区间上都为0,所以不予考虑;这样就能够把上述线性方程组求解出来,对于不同的n逼近效果不同,如果n(程序里我用的N表示)取得足够大,那么逼近效果将会非常好。下面我给出如下泊松方程的逼近情况:

逼近效果如下图:

可以看出当N=20的情况下逼近效果已经非常好了,这也说明有限元的做法确实可行,那么逼近过程中误差的收敛的情况如何呢?收敛阶是多少呢?如果学过有限元逼近理论应该知道,对于分片一次多项式逼近,其误差收敛阶在2阶左右(不同的误差度量,误差结果不一样,收敛阶也不完全相同,但都在2阶附近),感兴趣的同学可以自己去编程测试.

至此,分片一次多项式空间的逼近理论就暂时讲到这。不懂得请留言评论……

有限元方法数学理论之分片一次多项式空间逼近理论相关推荐

  1. 有限元方法数学理论之分片二次多项式空间逼近理论

    在一个月之前我已经对有限元方法数学理论之分片一次多项式空间逼近理论进行了讲解,由于中间有点事情所以对二次多项式情况没有进行解说,下面开始正文. 其实思想跟一次多项式的情况一样,只不过这里在一个小区间上 ...

  2. 【信息科学技术与创新】数据压缩的理论方法与现实意义 信息论 压缩编码 通信的数学理论 Huffman编码 LZ算法 虚幻引擎与数据压缩

    数据压缩的理论方法与现实意义 摘要 首先通过信息论引出数据压缩编码的理论方法 接着结合目前技术发展分析压缩编码的现实意义 最后总结思考未来通信与存储的压缩方法 Navigator 数据压缩的理论方法与 ...

  3. 神经网络的数学方法及应用,神经网络的数学理论

    神经网络算法原理 4.2.1 概述人工神经网络的研究与计算机的研究几乎是同步发展的. 1943年心理学家McCulloch和数学家Pitts合作提出了形式神经元的数学模型,20世纪50年代末,Rose ...

  4. 二维三角元有限元方法matlab,有限元C++编程实践.doc

    有限元C++编程实践.doc 基于有限元算法的编程实践学号:2011043010031 姓名:廖校毅电子科技大学 物理电子学院[摘要]本文就有限元算法在电磁场分析中的应用展开研究与实践,从电磁场的 M ...

  5. pde中微元分析法的主要思想_有限元方法的核心思想

    有限元法(Finite Element Method)是基于近代计算机的快速发展而发展起来的一种近似数值方法,用来解决力学,数学中的带有特定边界条件的偏微分方程问题(PDE).而这些偏微分方程是工程实 ...

  6. 菲尔兹奖数学家丘成桐:人工智能中的数学理论尚无很大突破

    来源:机器人 10月17日,在2019中关村论坛上,菲尔兹奖首位华人获得者.美国国家科学院院士.哈佛大学教授丘成桐发表主旨演讲.他呼吁国家应重视基础科学和数学的发展.一个国家的强大和长治久安,离不开强 ...

  7. 深入理解张正友相机标定法:数学理论详细推导

    最近在项目中需要在激光雷达(Lidar)和相机(Camera)之间进行标定,即需要标定出相机内参和外参,使用的标定方法是张正友标定法,这里给出其数学理论推导过程. 论文原文:<A Flexibl ...

  8. 计算机科学数学理论浅谈 (转载)

    [来源:TSTC文档中心] 计算机自从其诞生之日起,它的主要任务就是进行各种各样的科学计算.文档处理,数据处理,图像处理,硬件设计, 软件设计等等,都可以抽象为两大类:数值计算与非数值计算.作为研究计 ...

  9. 从文本分类问题中的特征词选择算法追踪如何将数学知识,数学理论迁移到实际工程中去...

    博文转载请注明作者和出处(作者:finallyliuyu :出处博客园) 附:<卡方特征词选择算法> <DF特征词选择算法> 一.数学背景 将数学知识.数学理论以及数学思想迁移 ...

  10. 一文读懂APS系统的核心算法和数学理论

    APS系统拥有以数学模型为基础的最强算法,可以满足更高难度.更广泛的需求.但它不是可以直观理解的方法,真正的难点不在于把算法变成软件,而是在于寻找和证明这个算法.因此,各APS公司对其核心算法和数学理 ...

最新文章

  1. Dialog 带白色的边的处理方法
  2. k8s与CICD--将drone部署到kubernetes中,实现agent动态收缩
  3. 光电编码器的原理及应用场合_旋转式光电编码器工作原理及在视觉检测中的使用...
  4. mysql mysqldb_mysql模块mysqldb
  5. 苦逼的草根站长屡战屡败,屡败屡战是因为爱网站
  6. k210 hc-05 智慧课堂 人脸识别签到
  7. 微信群机器人管理系统源码
  8. 转载:医疗保险,公积金、养老、生育、工伤、失业保险
  9. python Excel xlsx file; not supported
  10. 每周分享第二期:人工智能对未来工作的影响
  11. 【JavaScript】实现移动小精灵
  12. Matplotlib绘图库入门(一):pyplot绘图基础
  13. [NameError]: name ‘F’ is not defined
  14. 如何在CTEX上对论文进行排版
  15. Spring 增强处理Advice
  16. ZLiTE Lightweight Clients for Shielded Zcash Transactions using Trusted Execution
  17. VMware 连接ISO镜像时,切换可移动设备时遇到错误: 对磁盘“ide1:0”执行的连接控制操作失败?
  18. 什么是胃肠神经官能症,会有哪些症状出现呢?
  19. 阿里云AMD服务器计算型实例规格族c6a实例CPU网络收发包云盘IOPS详解
  20. malloc函数java_malloc函数具体解释

热门文章

  1. 2019DTCC大会感受
  2. mysql如何讲ipv6改成ipv4_IPv4 和 IPv6 地址如何转换?
  3. 新代系统9服务器警报,新代系统数控铣加工中心。新代系统龙门加工中心故障与报警...
  4. 黑苹果 服务器系统安装,黑苹果原版安装从零开始之 :MBR转GPT和新建EFI分区详细教程 _ 黑苹果乐园...
  5. 【模块学习】关于磁力计HMC5883L型号探索
  6. 【原创】MapGIS K9 三维二次开发入门
  7. 金字塔固定资产管理软件_目前常用固定资产管理软件有哪些?
  8. 大疆无人机实现目标定位
  9. 蘑菇云matlab程序,蘑菇云刷机精灵怎么用?刷机视频详细图文教程
  10. python 赚钱 知乎_2020年,小红书、知乎与B站谁能赚钱?