⑴ 求数值计算方法答案(韩旭里)复旦大学出版社,谢谢。在这两天给出最好,急急急急急!
《数值分析 中南大学韩旭里 126讲》网络网盘资源免费下载
链接: https://pan..com/s/1ath5aUEumr5ueV5d_GRa5Q
数值分析 中南大学 韩旭里 126讲|线性方程组的迭代解法(一).mp4|线性方程组的迭代解法(五).mp4|线性方程组的迭代解法(四).mp4|线性方程组的迭代解法(三).mp4|线性方程组的迭代解法(七).mp4|线性方程组的迭代解法(六).mp4|线性方程组的迭代解法(二).mp4|数值积分与数值微分(一).mp4|数值积分与数值微分(五).mp4|数值积分与数值微分(四).mp4|数值积分与数值微分(十一).mp4|数值积分与数值微分(十五).mp4|数值积分与数值微分(十四).mp4|数值积分与数值微分(十三).mp4
⑵ 求数值计算方法 第三版 李有法 朱建新 课后答案
数值计算方法如下:
1、有限元法:有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式。
借助于变分原理或加权余量法,将微分方程离散求解。采用不同的权函数和插值函数 形式,便构成不同的有限元方法。
在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元 上的近似解构成。
根据所采用的权函数和插值函数的不同 ,有限元方法也分为多种计算格式。从权函数的选择来说,有配置法、矩量法、最小二乘法和伽辽金法,从计算单元网格的形状来划分,有三角形网格、四边形网格和多边形网格,从插值函数的精度来划分,又分为线性插值函数和高次插值函数等。不同的组合 同样构成不同的有限元计算格式。
2、多重网格方法:多重网格方法通过在疏密不同的网格层上进行迭代,以平滑不同频率的误差分量。具有收敛速度快,精度高等优点。
多重网格法基本原理微分方程的误差分量可以分为两大类,一类是频率变化较缓慢的低频分量;另一类是频率高,摆动快的高频分量。
一般的迭代方法可以迅速地将摆动误差衰减,但对那些低频分量,迭代法的效果不是很显着。高频分量和低频分量是相对的,与网格尺度有关,在细网格上被视为低频的分量,在粗网格上可能为高频分量。
多重网格方法作为一种快速计算方法,迭代求解由偏微分方程组离散以后组成的代数方程组,其基本原理在于一定的网格最容易消除波长与网格步长相对应的误差分量。
该方法采用不同尺度的网格,不同疏密的网格消除不同波长的误差分量,首先在细网格上采用迭代法,当收敛速度变缓慢时暗示误差已经光滑,则转移到较粗的网格上消除与该层网格上相对应的较易消除的那些误差分量,这样逐层进行下去直到消除各种误差分量,再逐层返回到细网格上。
3、有限差分方法:有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。
有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。
对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑,可分为中心格式和逆风格式。考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。
构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式:
一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。
4、有限体积法:有限体积法(Finite Volume Method)又称为控制体积法。其基本思路是:将计算区域划分为一系列不重复的控制体积,并使每个网格点周围有一个控制体积;将待解的微分方程对每一个控制体积积分,便得出一组离散方程。其中的未知数是网格点上的因变量的数值。
为了求出控制体积的积分,必须假定值在网格点之间的变化规律,即假设值的分段的分布的分布剖面。从积分区域的选取方法看来,有限体积法属于加权剩余法中的子区域法;从未知解的近似方法看来,有限体积法属于采用局部近似的离散方法。简言之,子区域法属于有限体积发的基本方法。
有限体积法的基本思路易于理解,并能得出直接的物理解释。离散方程的物理意义,就是因变量在有限大小的控制体积中的守恒原理,如同微分方程表示因变量在无限小的控 制体积中的守恒原理一样。
限体积法得出的离散方程,要求因变量的积分守恒对任意一组控制体积都得到满足,对整个计算区域,自然也得到满足。这是有限体积法吸引人的优点。有一些离散方法,例如有限差分法,仅当网格极其细密时,离散方程才满足积分守恒。
而有限体积法即使在粗网格情况下,也显示出准确的积分守恒。就离散方法而言,有限体积法可视作有限单元法和有限差分法的中间物。有限单元法必须假定值在网格点之间的变化规律(既插值函数),并将其作为近似解。
有限差分法只考虑网格点上的数值而不考虑值在网格点之间如何变化。有限体积法只寻求的结点值 ,这与有限差分法相类似;但有限体积法在寻求控制体积的积分时,必须假定值在网格点之间的分布,这又与有限单元法相类似。
在有限体积法中,插值函数只用于计算控制体积的积分,得出离散方程之后,便可忘掉插值函数;如果需要的话,可以对微分方程 中不同的项采取不同的插值函数。
5、近似求解的误差估计方法:近似求解的误差估计方法共有三大类:单元余量法,通量投射法及外推法。
单元余量法广泛地用于以FEM离散的误差估计之中,它主要是估计精确算子的余量,而不是整套控制方程的全局误差。
这样就必须假定周围的单元误差并不相互耦合,误差计算采用逐节点算法进行。单元余量法的各种不同做法主要来自对单元误差方程的边界条件的不同处理办法。基于此,该方法能够有效处理局部的残余量,并能成功地用于网格优化程序。
通量投射法的基本原理来自一个很简单的事实:精确求解偏微分方程不可能有不连续的微分,而近似求解却可以存在微分的不连续,这样产生的误差即来自微分本身,即误差为系统的光滑求解与不光滑求解之差。该方法与单元余量法一样,对节点误差采用能量范数,故也能成功地用于网格优化程序。
单元余量法及通量投射法都局限于局部的误差计算(采用能量范数),误差方程的全局特性没有考虑。另外计算的可行性(指误差估计方程的计算时间应小于近似求解计算时间)不能在这两种方法中体现,因为获得的误差方程数量,阶数与流场控制方程相同。
外推是指采用后向数值误差估计思想由精确解推出近似解的误差值。各类文献中较多地采用Richardson外推方法来估计截断误差。无论是低阶还是高阶格式,随着网格的加密数值计算结果都会趋近于准确解。但由于计算机内存与计算时间的限制,实际上不能采用这种网格无限加密的办法。
6、多尺度计算方法:近年来发展的多尺度计算方法包括均匀化方法、非均匀化多尺度方法、以及小波数值均匀化方法、多尺度有限体积法、多尺度有限元法等。
该方法通过对单胞问题的求解,把细观尺度的信息映射到宏观尺度上,从而推导出宏观尺度上的均匀化等式,即可在宏观尺度上求解原问题。均匀化方法在很多科学和工程应用中取得了巨大成功,但这种方法建立在系数细观结构周期性假设的基础上,因此应用范围受到了很大限制。
鄂维南等提出的非均匀化多尺度方法,是构造多尺度计算方法的一般框架。该方法有两个重要的组成部分:基于宏观变量的整体宏观格式和由微观模型来估计缺少的宏观数据,多尺度问题的解通过这两部分共同得到。
该方法基于多分辨分析,在细尺度上建立原方程的离散算子,然后对离散算子进行小波变换,得到了大尺度上的数值均匀化算子。此方法在大尺度上解方程,大大地减小了计算时间。
该法在宏观尺度上进行网格剖分,然后通过在每个单元里求解细观尺度的方程(构造线性或者振荡的边界条件)来获得基函数。从而把细观尺度的信息反应到有限元法的基函数里,使宏观尺度的解包含了细观尺度的信息。但多尺度有限元方法在构造基函数时需要较大的计算量。
借助于变分原理或加权余量法,将微分方程离散求解。采用不同的权函数和插值函数 形式,便构成不同的有限元方法。
在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元 上的近似解构成。
根据所采用的权函数和插值函数的不同 ,有限元方法也分为多种计算格式。从权函数的选择来说,有配置法、矩量法、最小二乘法和伽辽金法,从计算单元网格的形状来划分,有三角形网格、四边形网格和多边形网格,从插值函数的精度来划分,又分为线性插值函数和高次插值函数等。不同的组合 同样构成不同的有限元计算格式。
2、多重网格方法:多重网格方法通过在疏密不同的网格层上进行迭代,以平滑不同频率的误差分量。具有收敛速度快,精度高等优点。
多重网格法基本原理微分方程的误差分量可以分为两大类,一类是频率变化较缓慢的低频分量;另一类是频率高,摆动快的高频分量。
一般的迭代方法可以迅速地将摆动误差衰减,但对那些低频分量,迭代法的效果不是很显着。高频分量和低频分量是相对的,与网格尺度有关,在细网格上被视为低频的分量,在粗网格上可能为高频分量。
多重网格方法作为一种快速计算方法,迭代求解由偏微分方程组离散以后组成的代数方程组,其基本原理在于一定的网格最容易消除波长与网格步长相对应的误差分量。
该方法采用不同尺度的网格,不同疏密的网格消除不同波长的误差分量,首先在细网格上采用迭代法,当收敛速度变缓慢时暗示误差已经光滑,则转移到较粗的网格上消除与该层网格上相对应的较易消除的那些误差分量,这样逐层进行下去直到消除各种误差分量,再逐层返回到细网格上。
3、有限差分方法:有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。
有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。
对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑,可分为中心格式和逆风格式。考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。
构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式:
一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。
4、有限体积法:有限体积法(Finite Volume Method)又称为控制体积法。其基本思路是:将计算区域划分为一系列不重复的控制体积,并使每个网格点周围有一个控制体积;将待解的微分方程对每一个控制体积积分,便得出一组离散方程。其中的未知数是网格点上的因变量的数值。
为了求出控制体积的积分,必须假定值在网格点之间的变化规律,即假设值的分段的分布的分布剖面。从积分区域的选取方法看来,有限体积法属于加权剩余法中的子区域法;从未知解的近似方法看来,有限体积法属于采用局部近似的离散方法。简言之,子区域法属于有限体积发的基本方法。
有限体积法的基本思路易于理解,并能得出直接的物理解释。离散方程的物理意义,就是因变量在有限大小的控制体积中的守恒原理,如同微分方程表示因变量在无限小的控 制体积中的守恒原理一样。
限体积法得出的离散方程,要求因变量的积分守恒对任意一组控制体积都得到满足,对整个计算区域,自然也得到满足。这是有限体积法吸引人的优点。有一些离散方法,例如有限差分法,仅当网格极其细密时,离散方程才满足积分守恒。
而有限体积法即使在粗网格情况下,也显示出准确的积分守恒。就离散方法而言,有限体积法可视作有限单元法和有限差分法的中间物。有限单元法必须假定值在网格点之间的变化规律(既插值函数),并将其作为近似解。
有限差分法只考虑网格点上的数值而不考虑值在网格点之间如何变化。有限体积法只寻求的结点值 ,这与有限差分法相类似;但有限体积法在寻求控制体积的积分时,必须假定值在网格点之间的分布,这又与有限单元法相类似。
在有限体积法中,插值函数只用于计算控制体积的积分,得出离散方程之后,便可忘掉插值函数;如果需要的话,可以对微分方程 中不同的项采取不同的插值函数。
5、近似求解的误差估计方法:近似求解的误差估计方法共有三大类:单元余量法,通量投射法及外推法。
单元余量法广泛地用于以FEM离散的误差估计之中,它主要是估计精确算子的余量,而不是整套控制方程的全局误差。
这样就必须假定周围的单元误差并不相互耦合,误差计算采用逐节点算法进行。单元余量法的各种不同做法主要来自对单元误差方程的边界条件的不同处理办法。基于此,该方法能够有效处理局部的残余量,并能成功地用于网格优化程序。
通量投射法的基本原理来自一个很简单的事实:精确求解偏微分方程不可能有不连续的微分,而近似求解却可以存在微分的不连续,这样产生的误差即来自微分本身,即误差为系统的光滑求解与不光滑求解之差。该方法与单元余量法一样,对节点误差采用能量范数,故也能成功地用于网格优化程序。
单元余量法及通量投射法都局限于局部的误差计算(采用能量范数),误差方程的全局特性没有考虑。另外计算的可行性(指误差估计方程的计算时间应小于近似求解计算时间)不能在这两种方法中体现,因为获得的误差方程数量,阶数与流场控制方程相同。
外推是指采用后向数值误差估计思想由精确解推出近似解的误差值。各类文献中较多地采用Richardson外推方法来估计截断误差。无论是低阶还是高阶格式,随着网格的加密数值计算结果都会趋近于准确解。但由于计算机内存与计算时间的限制,实际上不能采用这种网格无限加密的办法。
6、多尺度计算方法:近年来发展的多尺度计算方法包括均匀化方法、非均匀化多尺度方法、以及小波数值均匀化方法、多尺度有限体积法、多尺度有限元法等。
该方法通过对单胞问题的求解,把细观尺度的信息映射到宏观尺度上,从而推导出宏观尺度上的均匀化等式,即可在宏观尺度上求解原问题。均匀化方法在很多科学和工程应用中取得了巨大成功,但这种方法建立在系数细观结构周期性假设的基础上,因此应用范围受到了很大限制。
鄂维南等提出的非均匀化多尺度方法,是构造多尺度计算方法的一般框架。该方法有两个重要的组成部分:基于宏观变量的整体宏观格式和由微观模型来估计缺少的宏观数据,多尺度问题的解通过这两部分共同得到。
该方法基于多分辨分析,在细尺度上建立原方程的离散算子,然后对离散算子进行小波变换,得到了大尺度上的数值均匀化算子。此方法在大尺度上解方程,大大地减小了计算时间。
该法在宏观尺度上进行网格剖分,然后通过在每个单元里求解细观尺度的方程(构造线性或者振荡的边界条件)来获得基函数。从而把细观尺度的信息反应到有限元法的基函数里,使宏观尺度的解包含了细观尺度的信息。但多尺度有限元方法在构造基函数时需要较大的计算量。
⑶ 数值计算方法上机题编程,,,用c语言编程序,用牛顿迭代法求18的倒数,精度为0.0005,求大神解
用牛顿迭代法求方程(2*(X-4)+3)X-6=0的根。
其迭代公式为X2=X1-F(X1)/F'(X1)
F'(X1)为对方程求导。本题中P'(X1)=(6*x1-8)*x1-3;
编译显示正确,但一运行就死机,我已经死了3次了。(一开始还以为电脑的问题)
#include<iostream.h>
#include<math.h>
void main(void)
{float x1,x2=100;
do
{x1=x2;
x2=(float)x1-(((2*x1-4)*x1+3)*x1-6)/((6*x1-8)*x1-3);
}while(fabs(x2-x1)>pow(10,-5));
cout<<x2;
}
⑷ 计算方法引论题
x的相对误差表示等于dlnx;而lnx的绝对误差表示为 dlnx,所以lnx得绝对误差为a
⑸ 跪求《数值计算方法》第二版课后题答案,丁丽娟,北京理工大学出版社
陶宗旺76地俊星铁扇子宋清
研究方法
2.1 文献资料法;
通过对现有文献资料的查阅、分析和筛选,首先确定了参数统计、非参数统计、多元统计分析方法及数值计算方法四大类内容,每一类中再细分为若干种方法(见下页表一),并对它们逐一进行甄别、测试和数据验证。
2.2 面向对象的程序设计方法(OOP);
体育常用数据分析处理方法通常数据传输( 数据输入、数据输出)量大、计算过程有的简单有的繁杂、计算结果数据常常成批产生,采用面向对象的程序设计方法(OOP),充分应用可视化技术,将体育领域中最常用的一些数据分析处理方法开发为在Windows下运行的全中文界面的“傻瓜”型实用软件。
2.3 系统分析法
按照软件工程学的思想对系统作结构化分析(SA),建立开发文档,列出数据流图,最后利用Visual Basic编程技术开发、调试,完成软件后期制作。
3. 结果与分析
以VB为软件开发工具,筛选了体育训练、科研中最常用的参数统计、非参数统计、多元统计、数值计算方法等四大类共60多个,在Windows上平台开发为全中文界面“傻瓜”型多功能实用软件:可为运动训练中的有关数据作量化分析,可为体育科研人员提供一个分析处理数据的实用工具,也可为高校开设相关课程的课堂教学、学生上机实习作教学辅助软件。软件的特点是:将四大类实用方法封装在四个功能模块中,使不熟悉各种体育用数据处理方法的体育专业人员可以在电脑上应用参数统计、非参数统计、多元统计、数值计算方法解决问题。
3.1 参数统计模块
参数统计方法用来估计总体的某一参数(例如总体平均数、标准差等),或是检验总体参数是否不同。因此,需要明确样本所来自的总体的分布或对此分布做出假设,而总体分布的特征是通过总体参数来决定的。本模块包括了体育统计中具有数字特征、量化分析的一些概念和常用方法,在实际应用中由于很难掌握总体的全部情况,
只能根据样本计算出相应的数字特征值来估计它,评分方法、差异的显着性检验、相关分析、回归分析等都是本模块中的重要内容。
3.2 非参数统计模块
非参数统计方法适用于未知分布的资料,所以应用范围广、方法简便。体育活动中未知分布的资料很多,对于那些只分`等级、只排名次或只用二值逻辑(例如只有成功/失败、正确/错误、阳性/阴性等两种结果)表示的资料的分析与处理,常常使用非参数统计方法。本模块包括了体育统计中具有非数字特征、定性分析的一些概念和常用方法,如各种检验方法、作图法、相关分析、权重回归等。
3.3 多元统计分析模块
多元统计分析是研究分析多个因素(变量或指标)之间关系的统计方法,体育领域中应用广泛,模块中包含了9种共计17个常用的多元统计方法,是体育科研和教练员分析问题、处理数据的主要方法和手段。
3.4 数值计算方法模块
数值计算方法近年来开始应用于运动生物力学分析、体育系统仿真技术研究等,按照“针对实际问题→抽象数学模型→确定数值计算方法→程序设计→上机处理出结果”的模式,模块中包括了函数插值、曲线拟合、数据平滑等三类数值处理方法。
对于上述四个模块中每一种数值方法,软件中配备了“例题演示”(如图一)和详尽的“使用说明”(如图二);如果用户对所选用的方法不是很熟悉,那么可先浏览一下软件为该方法配备的例题演示,通过例题,用户可以了解该方法输入/输出哪些初始数据(如图三)、中间结果和最后结果(如图四)。如果用户希望了解所用的方法的初始数据如何操作?有无参数需现场输入等,可阅读相应的“使用说明”,它会详尽地告诉用户这一切。
针对体育科研和训练的特点,为使软件的板块结构清晰、数据流畅、每个数值处理方法自成一体,软件中使用了多文档界面(MDI,Multiple Document Interface)即多窗体结构,选择数值方法的主菜单由父窗体控制,每个方法各自为一个子窗体,子窗体被包含在父窗体中,父窗体为每个子窗体提供工作空间。针对每个子窗体上的某一种方法,分别设置了“使用说明”、“初始数据录入”、“数值计算”、“ 打印输出”、“清窗口”和“返回主菜单”等六个功能块,这样,用户在处理数据时,需要做哪项工作,只需用鼠标点击相应的按钮就可以了。
4. 结束语
数据处理分析方法目前已广泛应用到体育科研和训练的许多领域,随着计算机的进一步普及和软件开发技术的“平民化”,开发一些体育常用的数据处理分析方法实用软件是必要、可行的,它为计算机数值处理技术在体育领域内开辟了一个应用窗口;
体育训练、科研中量化模型的研究,计算机数值方法是量化分析的最有效工具,软件的开发研制将数据处理分析方法实用化,为获取准确的量化数据提供了一种简捷、快速、有效的手段;
软件中的部分内容从一九九六年开始在国内推广应用,在体育领域取得了较好的社会效益和经济效益。