数值分析实验 第二章
- 格式:ppt
- 大小:461.00 KB
- 文档页数:19


第二章有限元法的基本原理1. 引言有限元法(Finite Element Method, FEM)是一种数值分析方法,用于解决工程和科学领域中的复杂物理问题。
它通过将连续的物理领域离散化成许多小元素,通过求解代表元素之间关系的离散方程来近似解决原问题。
本章将介绍有限元法的基本原理。
2. 有限元法的基本思想有限元法的基本思想是将复杂的问题分割成更小的、易于处理的部分,通过求解这些部分的解,并通过它们之间的关系来得到整体解。
在有限元法中,将连续问题离散化为有限元模型,分为以下几个步骤:2.1 建立几何模型首先,根据实际问题建立几何模型。
几何模型可以是二维或三维的,通常使用节点和单元表示。
节点表示模型中的离散点,单元表示连接节点的几何形状。
2.2 确定节点自由度每个节点都有与之关联的自由度,它们是用来表示节点状态的参数。
常见的自由度有位移、温度等。
2.3 建立单元和节点之间的关系根据单元类型和节点连接关系,建立单元与节点之间的关系。
通常,一个单元由若干个节点组成。
2.4 建立元素刚度矩阵根据单元类型和材料参数,建立元素刚度矩阵。
2.5 建立整体刚度矩阵利用单元刚度矩阵和节点关系,建立整体刚度矩阵。
整体刚度矩阵由元素刚度矩阵按照节点自由度的排列组成。
2.6 施加边界条件和载荷根据实际问题,施加边界条件和载荷。
边界条件可以是位移、力或温度等。
2.7 求解方程通过将边界条件和载荷应用于整体刚度矩阵,可以得到未知节点的解答。
3. 有限元法的优缺点3.1 优点•适用于复杂几何形状和复杂边界条件的问题。
有限元法可以通过将问题离散化为小元素来逼近实际几何形状和边界条件。
•高精度的数值解。
有限元法通过增加节点数量和使用高阶元素可以得到更精确的数值解。
•灵活性。
有限元法可以灵活地处理不同类型的物理问题,例如结构力学、热传导、电磁场等。
3.2 缺点•需要大量的计算资源。
有限元法需要求解大型稀疏矩阵,这导致了计算资源的要求较高。
第一章 绪论1. *x = n 21k a a a .010⨯±,如果|*x -x|≤0.5n k 10-⨯(这里n 是使此式成立的最大正整数),则称*x 为x 的具有n 位有效数字的近似值。
2.定理:设x 的近似值*x 有(1-1)的表示式: (1)如果*x 有n 位有效数字,则n 1110a 21|x ||x x |-**⨯≤-(2)如果n 1110)1a (21|x ||x x |-**⨯+≤-,则*x 至少有n 位有效数字。
第二章 非线性方程根求解 1. (零点存在定理)如果f(x)在[a,b]上连续,使f(a)⋅f(b)<0,则必存在α∈(a,b),使f(α)=0。
2.二分法的误差: |1k 1k k k 2ab |x x ||x x +-*-=-≤- 3. 局部收敛性:设α是f(x)=0的根,若存在α的一个邻域∆,当迭代初值属于∆时,迭代法得到的序列{k x }收敛到α,则称该迭代法关于根α具有局部收敛性。
4. 收敛速度:设i x 为第i 次迭代值,α是f(x)=0的根,令α-=εi i x ,且假设迭代收敛,即α=∞→i i x lim 。
若存在实数P ≥1,使 c ||||limpi 1i i =εε+∞→≠0 ,则称此方法关于根α具有P阶收敛速度。
C 称为渐近误差常数,渐近误差常数C 与f(x)有关。
C ≠0保证了P 的唯一性。
对于特殊的函数,C 可能为零,此时,由这个函数针对此方法迭代产生的序列收敛得更快。
一般情况下,P 越大,收敛就越快。
当P=1时,我们称为线性收敛。
P>1,称为超线性收敛。
P=2,称为平方收敛。
5.牛顿迭代法:)x (f )x (f x x k k k 1k '-=+定理3:如果方程f(x)=0的根α是单根,且在α的某领域内f(x)具有二阶的连续导数,则Newton 迭代法必是局部收敛的 且 )(f 2)(f lim2i1i i α'α''-=εε+∞→(即具有二阶收敛速度)定理4:如果α是方程f(x)=0的r 重根(r>1),且f(x)在α的某邻域内具有r 阶连续导数,则Newton 法具有局部收敛性,且具有线性收敛速度。
数值分析实验(2014,9,16~10,28)信计1201班,人数34人数学系机房数值分析计算实习报告册专业__________________学号_______________姓名_______________2014~2015年第一学期实验一数值计算的工具Matlab1. 解释下MATLABS序的输出结果程序:t=0.1n=1:10e=n/10-n*te 的结果:0 0 -5.5511e-017 0 0-1.1102e-016 -1.1102e-016 0 0 02. 下面MATLABS序的的功能是什么?程序:x=1;while 1+x>1,x=x/2,pause(0.02),e nd用迭代法求出x=x/2,的最小值x=1;while x+x>x,x=2*x,pause(0.02),e nd用迭代法求出x=2*x,的值,使得2x>Xx=1;while x+x>x,x=x/2,pause(0.02),e nd用迭代法求出x=x/2,的最小值,使得2x>X3. 考虑下面二次代数方程的求解问题2ax bx c = 0公式x=电上4ac是熟知的,与之等价地有_____________________________ ,对于2a-b ■ b -4aca =1,b =100000000,c =1,应当如何选择算法。
b ~4ac计算,因为b与b2— 4ac相近,两个相加减不宜应该用2a u做分母3 5 74. 函数sin(x)有幂级数展开sin x = x - x - - ■■3! 5! 7!利用幕级数计算sinx的MATLAB程序为fun cti on s=powers in(x)s=0;t=x;n=1;while s+t~=s;s=s+t ;t=-x A2/ ((n+1)*(n+2) ) *t ;n=n+2 ;endt仁cputime;pause(10);t2=cputime;t0=t2-t1(a) 解释上述程序的终止准则。
数值分析课程第五版课后习题答案(李庆扬等)数值分析课程第五版课后习题答案(李庆扬等)第一章:数值分析导论1. 解答:数值分析是一门研究如何使用计算机来解决数学问题的学科。
它包括了从数学理论到计算实现的一系列技术。
数值分析的目标是通过近似的方式求解数学问题,其结果可能不是完全精确的,但是能够满足工程或科学应用的要求。
2. 解答:数值分析在实际应用中起着重要的作用。
它可以用于求解复杂的数学方程、计算机模拟及建模、数据的统计分析等等。
数值分析是科学计算和工程计算的基础,对许多领域都有着广泛的应用,如物理学、经济学、生物学等。
3. 解答:数值方法指的是使用数值计算的方式来求解数学问题。
与解析方法相比,数值方法一般更加灵活和高效,可以处理一些复杂的数学问题。
数值方法主要包括了数值逼近、插值、数值积分、数值微分、线性方程组的求解、非线性方程的求根等。
4. 解答:计算误差是指数值计算结果与精确解之间的差异。
在数值计算中,由于计算机的有限精度以及数值计算方法本身的近似性等因素,都会导致计算误差的产生。
计算误差可以分为截断误差和舍入误差两种。
第二章:数值误差分析1. 解答:绝对误差是指实际值与精确值之间的差异。
例如,对于一个计算出的数值近似解x和精确解x_0,其绝对误差为| x - x_0 |。
绝对误差可以衡量数值近似解的精确程度,通常被用作评估数值计算方法的好坏。
2. 解答:相对误差是指绝对误差与精确解之间的比值。
对于一个计算出的数值近似解x和精确解x_0,其相对误差为| (x - x_0) / x_0 |。
相对误差可以衡量数值近似解相对于精确解的精确度,常用于评估数值计算方法的收敛速度。
3. 解答:舍入误差是由于计算机的有限精度而引起的误差。
计算机中使用的浮点数系统只能表示有限的小数位数,因此在进行数值计算过程中,舍入误差不可避免地会产生。
舍入误差会导致计算结果与精确结果之间存在差异。
4. 解答:误差限度是指对于给定的数值计算问题,所能容忍的误差范围。