结构动力响应数值计算方法简介

数值计算方法学习指导书内容简介

数值计算方法学习指导书内容简介数值计算方法学习指导书内容简介《数字信号处理学习指导》是浙江省高等教育重点建设教材、应用型本科规划教材《数字信号处理》(唐向宏主编，浙江大学出版社出版，以下简称教材)的配套学习指导书，内容包括学习要求、例题分析、教材习题解答、自测练习以及计算机仿真实验等。学习指导书紧扣教材内容，通过例题讲解，分析各章节的学习重点、难点以及需要理解、掌握和灵活运用的基本概念、基本原理和基本方法。全书共有66例例题分析、121题题解、2套自测练习和6个mat1ab计算机仿真实验。数值计算方法学习指导书目录绪论第1章离散时间信号与系统 1.1 学习要点 1.2 例题 1.3 教材习题解答第2章离散系统的变换域分析与系统结构 2.1 学习要点 2.2 例题 2.3 教材习题解答第3章离散时间傅里叶变换

3.1 学习要点 3.2 例题 3.3 教材习题解答第4章快速傅里叶变换 4.1 学习要点 4.2 例题 4.3 教材习题解答第5章无限长单位冲激响应(iir)数字滤波器的设计5.1 学习要点 5.2 例题 5.3 教材习题解答第6章有限长单位冲激响应(fir)数字滤波器的设计6.1 学习要点 6.2 例题 6.3 教材习题解答第7章数字信号处理中的有限字长效应 7.1 学习要点 7.2 例题 7.3 教材习题解答第8章自测题 8.1 自测题(1)及参考答案 8.2 自测题(2)及参考答案第9章基于matlab的上机实验指导 9.1 常见离散信号的matlab产生和图形显示

9.2 信号的卷积、离散时间系统的响应 9.3 离散傅立叶变换 9.4 离散系统的频率响应分析和零、极点分布 9.5 iir滤波器的设计 9.6 fir滤波器的设计数值计算方法学习指导书内容文摘第1章离散时间信号与系统 1.1 学习要点本章主要介绍离散时间信号与离散时间系统的基本概念，着重阐述离散时间信号的表示、运算，离散时间系统的性质和表示方法以及连续时间信号的抽样等。本章内容基本上是“信号与系统”中已经建立的离散时间信号与系统概念的复习。因此，作为重点学习内容，在概念上需要明白本章在整个数字信号处理中的地位，巩固和深化有关概念，注意承前启后，加强葙关概念的联系，进一步提高运用概念解题的能力。学习本章需要解决以下一些问题： (1)信号如何分类。 (2)如何判断一个离散系统的线性、因果性和稳定性。 (3)线性时不变系统(lti)与线性卷积的关系如何。 (4)如何选择一个数字化系统的抽样频率。 (5)如何从抽样后的信号恢复原始信号。因此，在学习本章内容时，应以离散时间信号的表示、离散时间系统及离散时间信号的产生为主线进行展开。信号的离散时间的表示主要涉及序列运算(重点是卷积和)、常用序列、如何判

数值分析综述-《数值分析与算法》徐士良

第2章矩阵与线性代数方程组一般的线性代数方程组，A非奇异可根据Cramer法则求解方程唯一解但是它的计算量很大。高斯消元法的算法时间复杂度是O(n3)，可以解一系列的线性方程；所占数据空间符合原地工作的原则。但是算法对数值计算不稳定（当分母为0或很小时）。可以用在计算机中来解决数千条等式及未知数。不过，如果有过百万条等式时，这个算法会十分费时。解决高斯法中的不稳定性，在每次归一化前增加选主元（列选主元、全选主元）过程。但是列选主元法仍不稳定，不适求解大规模线性代数方程组。全选主元的高斯消去法，则在复杂度降低的同时能够避免舍入误差，保证数值稳定性。高斯-约当消去法算法产生出来的矩阵是一个简化行梯阵式，而不是高斯消元法中的行梯阵式。相比起高斯消元法，此算法的效率比较低，却可把方程组的解用矩阵一次过表示出来。线性代数方程组的迭代解法简单迭代法：迭代格式发散但迭代值序列不一定发散，但收敛格式收敛，迭代值序列收敛于方程组的准确解与选取迭代初值无关。雅可比迭代法: 计算公式简单，且计算过程中原始矩阵A始终不变，比较容易并行计算。但是收敛速度较慢，而且占据的存储空间较大，所以工程中一般不直接用雅克比迭代法，而用其改进方法。高斯-赛德尔迭代法：较上面的迭代复杂，但是矩阵的条件相对宽松。松弛法：需要根据经验去调整，收敛速度依赖松弛参数的选择，收敛条件的要求更宽松。共轭梯度法：是介于最速下降法与牛顿法之间的一个方法，它仅需利用一阶导数信息，但克服了最速下降法收敛慢的缺点，又避免了牛顿法需要存储和计算Hesse矩阵并求逆的缺点，在各种优化算法中，共轭梯度法是非常重要的一种。其优点是所需存储量小，具有步收敛性，稳定性高，而且不需要任何外来参数。共轭梯度法不仅是解决大型线性方程组最有用的方法之一，也是解大型非线性最优化最有效的算法之一。第3章矩阵特征值乘幂法计算绝对值最大的特征值：其收敛速度受限于最大与次大特征值比值绝对值的大小，实际应用中采用加速技术。求对称特征值的雅克比方法96：每进行一次选装变换钱都需要在飞对角线的元素中选取绝对值最大的元素，很费时间，雅克比过关法对此做了改进。 QR方法求一般实矩阵的全部特征值98下100下：重复多次进行QR分解费时，计算工作量很大。一般先进行相似变换然后进行QR分解。但是这样仍然收敛速度慢，一般是线性收敛。实际应用中使用双重步QR变换将带原点的QR算法中相邻两步合并一步，加速收敛避免复数运算。第4章非线性方程与方程组二分法：每次运算后，区间长度减少一半，是线形收敛。优点是简单，但是不能计算复根和重根。简单迭代法：直接的方法从原方程中隐含的求出x，从而确定迭代函数 (x)，这种迭代法收敛速度较慢，迭代次数多。埃特金迭代法113中：对简单迭代进行改进，使在其不满足收敛条件下迭代过程也收敛，在其收敛时加快收敛速度，减少迭代次数降低时间复杂度。牛顿迭代法：其最大优点是在方程f(x) = 0的单根附近具有平方收敛，收敛速度快。而且该法还可以用来求方程的重根、复根。缺点：初值的选择会影响收敛结果。牛顿下山法：保证函数值稳定下降，且有牛顿法的收敛速度。

移动载荷作用下连续梁的动力响应分析

第八届全国振动理论及应用学术会议论文集，上海，2003年11月移动载荷作用下连续梁的动力响应分析钟卫洲1, 2，罗景润1，高芳清3，徐友钜1 （1.中国工程物理研究院结构力学研究所，绵阳 621900；2.中国工程物理研究院研究生部，绵阳 621900; 3.西南交通大学振动与强度实验室，成都 610031）摘要：本文以磁悬浮交通轮轨接触车桥动力行为研究为背景，把车辆对桥梁的动力作用简化为一个稳态力和一个低频扰动力，把连续钢桥梁简化为伯努力—欧拉梁，建立了车辆过桥的力学模型和振动微分方程，运用模态分析法得到了该微分方程的解析解，分析了连续桥梁频率方程、模态表达式以及低阶模态。援引德国TR06和连续钢梁的参数对不同速度的移动荷载下连续钢梁的动力响应进行计算分析，给出了相应条件下连续梁的动挠度曲线（w-t图和w-x图），并分析了桥梁的动力响应特征。本文的研究为评定桥梁在高速车辆作用下的稳定性和安全性提供了参考。关键词：连续梁；模态分析；动力响应；动挠度 Dynamic Response Analysis of Continuous Beam Under Moving Load ZHONG Wei-zhou 1, 2, LUO Jing-run 1, GAO Fang-qing3, XU You-ju 1 (1.Institute of Structural Mechanics of CAEP, Mianyang 621900; 2.Graduated School of CAEP, Mianyang 621900; https://www.360docs.net/doc/e012763674.html,boratory of Vibration and Intensity of SWJTU, Chengdu 610031) Abstract: This paper is based on the background of the study of the dynamic behavior between maglev vehicle and guideway. The moving force exerting on the bridge is simplified as a steady force and a pulsating force with low frequency. The continuous steel beam is taken as Bernoulli-Euler beam, then the corresponding force model and vibrating equation of the bridge is established. The modal analysis method is applied to solve the equation of vibration. Frequency equation, analytical solution of mode of the beam and the lower modes are analysed. By quoting the data of TR06 of German, the dynamic response of continuous beam is obtained under moving vehicle at several typical speeds. The results of this paper can be taken as reference to assess security and stability of a bridge under moving load.

数值计算方法课程报告

课程报告课程名称______《数值计算》 __ 学生学院_____机电工程学院___ 专业班级_____微电子（1）班____ 学号________ 学生姓名_______________ 指导教师_____ ________ XXXX年XX月XX日

姓名：线学号：订装专业：学院：广东工业大学考试试卷( A ) 课程名称: 数值计算试卷满分100 分考试时间: 2015 年 12 月 26 日(第 17 周星期六) 题号一二三四五六七八九十总分评卷得分评卷签名复核得分复核签名 “数值计算”考试要求 “数值计算”考试以开卷形式进行。在“数值计算”课程考试日（2015 年12 月 19 日，第 12 周星期五）考试时间，在考试教室领取试题，在 2015 年12 月 26 日（第 17 周星期六）进行答辩。不参加答辩者将取消考试成绩。 “数值计算”考试结果要求独立在计算机上完成，可使用Matlab或 C 程序编程实现。考试结果将以报告书形式提交，内容包括对问题描述、计算程序以及算例、计算结果、分析组成。计算程序要求具有通用性，能够处理异常情况，可以输入问题、算法参数、算例及初始值，在计算过程中显示当前计算状态、计算完成后显示计算结果。上述内容将作为试卷成绩的主要评定依据。特别提醒，不得使用教师在讲课和实验时的范例作为考试结果。报告书具体格式参考毕业设计手册。以考生学号命名的文件夹存放程序及报告书电子版，以班级为单位刻录在一张光盘中，与打印版报告书一起由班长和学习委员一起上交任课教师。数值计算课程总成绩将由试卷成绩（70%）、平时成绩（30%）组成。

基于 MATLAB 实现对结构动力响应的几种算法的验证

基于 MATLAB 实现对结构动力响应的几种算法的验证 1. 算例首先，本文给出一算例，结构在外力谐振荷载 P (t ) = P 0 sin θt 作用下，分别利用理论解法，杜哈梅积分， Wilson-θ 法求出该结构的位移时程反应。其中： m = 3.5×103 kg ， P 0 = 1.0×104 N ， k =1.3584515×107 ，ξ=0.05 ，θ=52.3s ?1 ，ω=62.3s ?1 ， ? D ω= ω 1-ξ2 =62.222 ，初始位移、速度v (0) = 0 ，v (0) = 0 ； 2. 算法验证 2.1 理论解法运动方程为： mv+cv+kv=0P sin t ?由线性代数解出其理论解为: ]cos 2)sin -[(] 4)-[(t] sin ) -(2cos [2] 4)-[(t] sin ) 0()0(cos )0([2 2222222022222 2 2 2 22 t t m P t m P e v v t v e v D D D t D D D t θξωθθθωθωξθωωωθωωξωξωθωξθωθωωξωωξωξω-++-+ ?++++=-- 由于初始位移v(0) =0 ，v(0) =0 ；则： ]cos 2)sin -[(] 4)-[(t] sin ) -(2cos [2] 4)-[(222 222220 22222222220 t t m P t m P e v D D D t θξωθθθωθωξθωωωθωωξωξωθωξθωθξω-++ -+ ?+=- v(t ) =-3.115t e ? 1.05269898×410-?[6.230cos62.222t ?18.106sin 62.222t] +2.012808757 6 10-??[1146sin 52.3t ?325.829cos52.3t] 可以用 MA TLAB 进行编程分析，画图位移时程图，详细程序见附录。 2.2Wilson-θ法 Wilson-θ法是Wilson 于1966年基于线性加速度法的基础上提出一种无条件收敛的计算方法。该方法假定在θt ?时程步长内，体系的加速度反应按线性变化。对于地震持续时间内

数值计算方法教学大纲

《数值计算方法》教学大纲课程编号：MI3321048 课程名称：数值计算方法英文名称：Numerical and Computational Methods 学时： 30 学分：2 课程类型：任选课程性质：任选课适用专业：微电子学先修课程：高等数学,线性代数集成电路设计与集成系统开课学期：Y3开课院系：微电子学院一、课程的教学目标与任务目标：学习数值计算的基本理论和方法，掌握求解工程或物理中数学问题的数值计算基本方法。任务：掌握数值计算的基本概念和基本原理，基本算法，培养数值计算能力。二、本课程与其它课程的联系和分工本课程以高等数学，线性代数，高级语言编程作为先修课程，为求解复杂数学方程的数值解打下良好基础。三、课程内容及基本要求 (一) 引论（2学时）具体内容：数值计算方法的内容和意义，误差产生的原因和误差的传播，误差的基本概念，算法的稳定性与收敛性。 1.基本要求（1）了解算法基本概念。（2）了解误差基本概念，了解误差分析基本意义。 2.重点、难点重点：误差产生的原因和误差的传播。难点：算法的稳定性与收敛性。 3.说明：使学生建立工程中和计算中的数值误差概念。 (二) 函数插值与最小二乘拟合（8学时）具体内容：插值概念，拉格朗日插值，牛顿插值，分段插值，曲线拟合的最小二乘法。 1.基本要求（1）了解插值概念。（2）熟练掌握拉格朗日插值公式，会用余项估计误差。（3）掌握牛顿插值公式。（4）掌握分段低次插值的意义及方法。

（5）掌握曲线拟合的最小二乘法。 2.重点、难点重点：拉格朗日插值, 余项，最小二乘法。难点：拉格朗日插值, 余项。 3.说明：插值与拟合是数值计算中的常用方法，也是后续学习内容的基础。 (三) 第三章数值积分与微分（5学时）具体内容：数值求积的基本思想，代数精度的概念，划分节点求积公式（梯形辛普生及其复化求积公式），高斯求积公式，数值微分。 1.基本要求（1）了解数值求积的基本思想，代数精度的概念。（2）熟练掌握梯形，辛普生及其复化求积公式。（3）掌握高斯求积公式的用法。（4）掌握几个数值微分计算公式。 2.重点、难点重点：数值求积基本思想，等距节点求积公式，梯形法，辛普生法，数值微分。难点：数值求积和数值微分。 3.说明：积分和微分的数值计算，是进一步的各种数值计算的基础。 (四) 常微分方程数值解法（5学时）具体内容：尤拉法与改进尤拉法，梯形方法，龙格—库塔法，收敛性与稳定性。 1.基本要求（1）掌握数值求解一阶方程的尤拉法，改进尤拉法，梯形法及龙格—库塔法。（2）了解局部截断误差，方法阶等基本概念。（3）了解收敛性与稳定性问题及其影响因素。 2.重点、难点重点：尤拉法,龙格－库塔法，收敛性与稳定性。难点：收敛性与稳定性问题。 3.说明：该内容是常用的几种常微分方程数值计算方法，是工程计算的重要基础。 (五) 方程求根的迭代法（4学时）具体内容：二分法，解一元方程的迭代法，牛顿法，弦截法。 1.基本要求（1）了解方程求根的对分法和迭代法的求解过程。（2）熟练掌握牛顿法。（3）掌握弦截法。 2.重点、难点重点：迭代法，牛顿法。

数值分析作业思考题汇总

￥数值分析思考题1 1、讨论绝对误差（限）、相对误差（限）与有效数字之间的关系。 2、相对误差在什么情况下可以用下式代替 3、查阅何谓问题的“病态性”，并区分与“数值稳定性”的不同点。 4、取，计算，下列方法中哪种最好为什么（1）(3 3-，（2）(2 7-，（3） ()3 1 3+ ，（4） ()6 1 1 ，（5）99- , 数值实验数值实验综述：线性代数方程组的解法是一切科学计算的基础与核心问题。求解方法大致可分为直接法和迭代法两大类。直接法——指在没有舍入误差的情况下经过有限次运算可求得方程组的精确解的方法，因此也称为精确法。当系数矩阵是方的、稠密的、无任何特殊结构的中小规模线性方程组时，Gauss消去法是目前最基本和常用的方法。如若系数矩阵具有某种特殊形式，则为了尽可能地减少计算量与存储量，需采用其他专门的方法来求解。 Gauss消去等同于矩阵的三角分解，但它存在潜在的不稳定性，故需要选主元素。对正定对称矩阵，采用平方根方法无需选主元。方程组的性态与方程组的条件数有关，对于病态的方程组必须采用特殊的方法进行求解。数值计算方法上机题目1 1、实验1. 病态问题实验目的：算法有“优”与“劣”之分，问题也有“好”和“坏”之别。所谓坏问题就是问题本身的解对数据变化的比较敏感，反之属于好问题。希望读者通过本实验对此有一个初步的体会。数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。 $ r e x x e x x ** * ** - == 141 . ≈)61

数值分析实验报告1

实验一误差分析实验（病态问题）实验目的：算法有“优”与“劣”之分，问题也有“好”与“坏”之别。对数值方法的研究而言，所谓坏问题就是问题本身对扰动敏感者，反之属于好问题。通过本实验可获得一个初步体会。数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。问题提出：考虑一个高次的代数多项式 )1.1() ()20()2)(1()(20 1∏=-=---=k k x x x x x p 显然该多项式的全部根为1,2,…,20共计20个，且每个根都是单重的。现考虑该多项式的一个扰动 )2.1(0 )(19=+x x p ε 其中ε是一个非常小的数。这相当于是对（）中19x 的系数作一个小的扰动。我们希望比较（）和（）根的差别，从而分析方程（）的解对扰动的敏感性。实验内容：为了实现方便，我们先介绍两个Matlab 函数：“roots ”和“poly ”。 roots(a)u = 其中若变量a 存储n+1维的向量，则该函数的输出u 为一个n 维的向量。设a 的元素依次为121,,,+n a a a ，则输出u 的各分量是多项式方程 01121=+++++-n n n n a x a x a x a 的全部根；而函数 poly(v)b =

的输出b 是一个n+1维变量，它是以n 维变量v 的各分量为根的多项式的系数。可见“roots ”和“poly ”是两个互逆的运算函数。 ;000000001.0=ess );21,1(zeros ve = ;)2(ess ve = ))20:1((ve poly roots + 上述简单的Matlab 程序便得到（）的全部根，程序中的“ess ”即是（）中的ε。实验要求：（1）选择充分小的ess ，反复进行上述实验，记录结果的变化并分析它们。如果扰动项的系数ε很小，我们自然感觉（）和（）的解应当相差很小。计算中你有什么出乎意料的发现表明有些解关于如此的扰动敏感性如何（2）将方程（）中的扰动项改成18x ε或其它形式，实验中又有怎样的现象出现（3）（选作部分）请从理论上分析产生这一问题的根源。注意我们可以将方程（）写成展开的形式， ) 3.1(0 ),(1920=+-= x x x p αα 同时将方程的解x 看成是系数α的函数，考察方程的某个解关于α的扰动是否敏感，与研究它关于α的导数的大小有何关系为什么你发现了什么现象，哪些根关于α的变化更敏感思考题一：（上述实验的改进）在上述实验中我们会发现用roots 函数求解多项式方程的精度不高，为此你可以考虑用符号函数solve 来提高解的精确度，这需要用到将多项式转换为符号多项式的函数poly2sym,函数的具体使用方法可参考Matlab 的帮助。

数值分析综述报告

淮阴工学院《数值分析》考试 ──基于Matlab的方法综合应用报告班级：金融1121 姓名：姚婷婷学号：1124104129 成绩：数理学院 2014年6月7日

《数值分析》课程综述报告前言：数值分析也称计算方法，它与计算工具的发展密切相关。数值分析是一门为科学计算提供必需的理论基础和有效、实用方法的数学课程，它的任务是研究求解各类数学问题的数值方法和有关的理论。正文：第一章近似计算与误差分析 1、产生误差的原因：①模型误差；②观测误差；③截断误差；④舍入误差。 2、四则运算的误差： ①加减法运算 ()()()****x y x y δδδ±=+ ②乘法运算 ()()() ****** *** ******xy x y xy xy xy x y x y y y x x x y x y y x δδδ-=-+-≤-+-?=+ ③ 除法运算： ()()() () () ***** ******* * * ** * * ** * *2 ** x x xy x y y y yy xy x y x y x y yy x x y y y x yy x y y x x y y δδ δ--=-+-=-+-= +?? ?≈ ??? 3、科学表示法、有效数字、近似值的精度任何一个实数都可以表示成如下的形式：其中：是正整数，是整数，如果是数的近似值并且则称该近似值具有位有效数字（significant digit ）。

此时，该近似值的相对误差为另一方面，若已知 ()() *111 1021n r x a δ-≤ + 那么， ()()***1112110.10 211 102 r m n n m n x x x x a a a a δ----≤?=+≤ 即：*x 至少有n 位有效数字。例： 3.141592653589793...π= 取其近似值如下： x*=3.14 x * =3.14159 x*=3.1415 x*=3.141 **213 100.314 110.0016...0.005101022 x x π--=?-=<=?=? **516 100.314159 110.0000026...0.00000510102 2 x x π--=?-=<=?=? **314 100.31415 110.000092...0.0001101022 x x π--=?-=<

用有限元方法进行摩托车动力响应分析报告

用有限元方法进行摩托车动力响应分析文>>月辉史春涛骞郝志勇摘要本文采用有限元方法对某125型骑式摩托车进行了动力响应分析。文章首先建立了摩托车整车的有限元模型，并利用该模型进行摩托车整车的动态特性计算，取得了和实验模态分析一致的结果。而后分析了摩托车在发动机激励和路面不平度激励下的整车动力学响应特性，得出了具有工程参考价值的结论。关键词摩托车应力有限元法本文采用有限元方法研究了摩托车整车结构的动态特性，并进行了在各种激励作用下的动力响应分析，得到了发动机车架的应力场，可用于进一步的摩托车强度分析。 1、摩托车有限元模型的建立摩托车有限元模型如图1所示。摩托车的车架结构大多是由各种截面形状的梁组合而成的空间框架结构，而且其截面尺寸，包括直径、壁厚，与构件长度相比很小，因此选用空间的直梁或者曲梁单元来离散车架结构，而车架的一些板件和加强盘可以采用空间板元模拟，各种梁单元的截面力学特性可用有限元程序的前处理模块或CAD软件计算。摩托车的发动机具有较大质量，同时也具有很大刚度。考虑到发动机在车体结构中所起的作用及变形小的特点，将发动机简化为若干个板单元，这些板的总质量应与发动机的质量相同。然后，根据发动机与车架的实际连接方式，将由这些板单元模拟的发动机与车架组装到一起。摩托车的减振器主要作用是支撑车体并缓和振动与冲击。考虑到减振器的结构与作用，简化后减振器的模型在受到载荷时应具有较大的轴向位移，同时又要有较大抗弯刚度。本文把减振器简化为一种梁单元和弹簧阻尼单元的综合体——轴向刚度由弹簧阻尼单元提供，而抗弯刚度由梁单元提供。摩托车车轮主要由轮胎和轮辋组成，其中轮胎直接与路面接触，与摩托车悬挂共同缓和摩托车行驶时所受到的冲击，并协助减振，轮辋是固定轮胎的骨架，它与轮胎共同承受作用在车轮上的负荷。轮辋可以采用若干个梁单元模拟，轮胎

《数值计算方法》上机实验报告

《数值计算方法》上机实验报告华北电力大学实验名称数值il?算方法》上机实验课程名称数值计算方法专业班级：电力实08学生姓名：李超然学号:200801001008 成绩: 指导教师:郝育黔老师实验日期：2010年04月华北电力大学实验报告数值计算方法上机实验报吿一. 各算法的算法原理及计算机程序框图1、牛顿法求解非线性方程 *对于非线性方程，若已知根的一个近似值，将在处展开成一阶 xxfx ()0, fx ()xkk 泰勒公式 "f 0 / 2 八八,fxfxfxxxxx 0 0 0 0 0 kkkk2! 忽略高次项，有，fxfxfxxx 0 ()()()，，, kkk 右端是直线方程，用这个直线方程来近似非线性方程。将非线性方程的 **根代入，即fx ()0, X ，* fxfxxx 0 0 0 0, ,, kkk fx 0 fx 0 0,

解出 fX 0 *k XX,, k' fx 0 k 水将右端取为，则是比更接近于的近似值，即xxxxk, Ik, Ik fx ()k 八XX, Ikk* fx()k 这就是牛顿迭代公式。 ,2,计算机程序框图:，见, ,3,输入变量、输出变量说明: X输入变量:迭代初值，迭代精度，迭代最大次数，\0 输出变量:当前迭代次数，当前迭代值xkl ,4,具体算例及求解结果: 2/16 华北电力大学实验报吿开始读入 l>k /fx()0?,0 fx 0 Oxx，，01* fx ()0 XX,,，?10 kk, ,1，kN, ?xx, 10 输出迭代输出X输出奇异标志1失败标志

,3,输入变量、输出变量说明: 结束例：导出计算的牛顿迭代公式，并il ?算。（课本P39例2-16） 115cc （0）, 求解结果: 10. 750000 10.723837 10. 723805 10. 723805 2、列主元素消去法求解线性方程组，1,算法原理: 高斯消去法是利用现行方程组初等变换中的一种变换，即用一个不为零的数乘 -个方程后加只另一个方程，使方程组变成同解的上三角方程组，然后再自下而上对上三角 3/16 华北电力大学实验报告方程组求解。列选主元是当高斯消元到第步时，从列的以下（包括）的各元素中选出绝 aakkkkkk 对值最大的，然后通过行交换将其交换到的位置上。交换系数矩阵中的两行（包括常ekk 数项），只相当于两个方程的位置交换了，因此，列选主元不影响求解的结 ,2,计算机程序框图:，见下页, 输入变量:系数矩阵元素，常向量元素baiji 输出变量:解向量元素bbb,,12n

动载荷作用下的结构动力响应分析

在动载荷作用下的结构动力响应分析一、问题的突出 1.求钢结构的动力学响应 2.结构及载荷数据立柱间距长180，宽120；顶距地面180 E=200GPa,v=0.3,p=7800kg/m3 A=100kg, w=20π 二、有限单元模型建立 1.单元类型选择以及截面性质定义 FINISH /CLE /PREP7 ET,1,BEAM188 mp,ex,1,200e9 mp,nuxy,1,0.3 mp,dens,1,7800 !截面1 sectype,1,beam,I secoffset,cent secdata,0.12,0.12,0.1212,0.00605, 0.00605,0.0039 !截面2 sectype,2,beam,I secoffset,cent secdata,0.06535,0.06535,0.0805,0.00465,0.00465,0.00285 !截面3 sectype,3,beam, HREC secoffset,cent secdata,0.06,0.06,0.0025,0.0025,0.0025,0.0025

选择Y方向为立柱方向，通过由下向上建模即由关键点到线。k,1, k,2,1.8 k,3,1.8,,1.2 k,4,,,1.2 k,5,,1.5 k,6,1.8,1.5 k,7,1.8,1.5,1.2 k,8,,1.5,1.2 l,1,5 l,2,6 l,3,7 l,4,8 k,9,1.8/2,1.8,1.2/2 l,9,5 l,9,6 l,9,7 l,9,8 l,5,8 l,8,7 l,7,6 l,6,5 3.对于不同的杆赋予不同的截面性质，并进行网格划分 lsel,s,,,1,4,1 lplot lesize,all,0.04 secnum,1 lmesh,all lsel,s,,,9,12 lplot lesize,all,0.04 secnum,2 lmesh,all lsel,s,,,5,8,1 lplot lesize,all,0.04 secnum,3 lmesh,all

导数的数值计算方法[文献综述]

毕业论文文献综述信息与计算科学导数的数值计算方法一、前言部分导数概念的产生有着直觉的起源，与曲线的切线和运动质点的速度有密切的关系．导数用于描述函数变化率，刻画函数的因变量随自变量变化的快慢程度．比如说，物理上考虑功随时间的变化率（称为功率），化学上考虑反应物的量对时间的变化率（称为反应速度），经济学上考虑生产某种产品的成本随产量的变化率（称为边际成本）等等，这些变化率在数学上都可用导数表示．导数由于其应用的广泛性，为我们解决所学过的有关函数问题提供了一般性的方法，导数是研究函数的切线、单调性、极值与最值等问题的有力工具；运用它可以简捷地解决一些实际问题，导数的概念是用来研究函数在一点及其附近的局部性质的精确工具，而对于函数在某点附近的性质还可以应用另一种方法来研究，就是通过最为简单的线性函数来逼近，这就是微分的方法．微分学是数学分析的重要组成部分，微分中值定理作为微分学的核心，是沟通导数和函数值之间的桥梁， Rolle 中值定理， Lagrange 中值定理， Cauchy 中值定理， Taylor 公式是微分学的基本定理，统称为微分学的中值定理，这四个定理作为微分学的基本定理，是研究函数形态的有力工具 ] 1[．在微分学中，函数的导数是通过极限定义的，但当函数用表格给出时，就不可用定义来求其导数，只能用近似方法求数值导数] 2[．最简单的数值微分公式是用差商近似地代替微商，常见的有 [3] ． ()()() 'f x h f x f x h +-≈ ， ()()() 'f x f x h f x h --≈， ()()() '2f x h f x h f x h +--≈ ．需要注意的是微分是非常敏感的问题，数据的微小扰动会使结果产生很大的变化] 4[．

数值分析实验报告模板

数值分析实验报告模板篇一：数值分析实验报告(一)(完整) 数值分析实验报告 1 2 3 4 5 篇二：数值分析实验报告实验报告一题目：非线性方程求解摘要：非线性方程的解析解通常很难给出，因此线性方程的数值解法就尤为重要。本实验采用两种常见的求解方法二分法和Newton法及改进的Newton法。利用二分法求解给定非线性方程的根，在给定的范围内，假设f(x,y)在[a,b]上连续，f(a)xf(b) 直接影响迭代的次数甚至迭代的收敛与发散。即若x0 偏离所求根较远，Newton法可能发散的结论。并且本实验中还利用利用改进的Newton法求解同样的方程，且将结果与Newton法的结果比较分析。前言：（目的和意义）掌握二分法与Newton法的基本原理和应用。掌握二分法的原理，验证二分法，在选对有根区间的前提下，必是收

敛，但精度不够。熟悉Matlab语言编程，学习编程要点。体会Newton使用时的优点，和局部收敛性，而在初值选取不当时，会发散。数学原理：对于一个非线性方程的数值解法很多。在此介绍两种最常见的方法：二分法和Newton法。对于二分法，其数学实质就是说对于给定的待求解的方程f(x)，其在[a,b]上连续，f(a)f(b) Newton法通常预先要给出一个猜测初值x0，然后根据其迭代公式xk?1?xk?f(xk) f'(xk) 产生逼近解x*的迭代数列{xk}，这就是Newton法的思想。当x0接近x*时收敛很快，但是当x0选择不好时，可能会发散，因此初值的选取很重要。另外，若将该迭代公式改进为 xk?1?xk?rf(xk) 'f(xk) 其中r为要求的方程的根的重数，这就是改进的Newton 法，当求解已知重数的方程的根时，在同种条件下其收敛速度要比Newton法快的多。程序设计：本实验采用Matlab的M文件编写。其中待求解的方程写成function的方式，如下 function y=f(x);

动力响应理论

第2章动力响应理论 2.1引言机柜结构动力响应的计算机仿真分析是以设备动力响应理论为基础的，是进行设备结构动力响应研究的一种有效手段。论文中主要研究设备动力响应两个方面的内容：设备结构固有特性分析和结构在地震波作用下的响应分析。固有特性分析可以得到结构的固有频率和固有振型，是进行响应分析的基础；地震波响应分析将得到设备响应的时间历程变化。在使用有限元工具对结构进行建模、分析之前必须掌握结构动力响应的理论和相关的有限元基本原理。因此，本章重点叙述了与设备结构动力响应相关的机械振动学理论及其有限元仿真技术。 2.2结构动力响应分析相关理论 2.2.1结构固有特性分析理论机柜设备结构的固有特性包括固有频率和振型，是响应分析的基础。通过进行结构的固有特性分析可以使设计有效地避开结构的共振频率。机柜设备是一个复杂振动系统，在理论分析过程中，常常可以把机柜设备简化为多自由度集中参数系统。一般，多自由度系统的自由振动方程可以写成如下形式： {}... []()[]{()}[]{()}{0}M x t C x t K x t ++= (21)a - 式中：[]M ， []C 和[]K 分别为系统的质量矩阵、阻尼矩阵和刚度矩阵;()x t 、 .()x t 、.. ()x t 分别为系统的位移列向量、速度列向量和加速度列向量。而多自由系统的无阻尼自由振动方程可以写成如下形式: {}.. []()[]{()}{0}M x t K x t += (21)b - 通常系统的自由振动是简谐振动，所以可以假设式 (21)b -的解为： {()}{}sin x t X pt = (22)- 式中：{}X 为系统的振幅列向量;p 为系统的自由振动频率。将(22)-代入(21)b -，就可以得到系统的振型方程，其具体形式如下： 2[][]{}{0}K M X p -= (23)- 可以看到，式(23)-是一个齐次线性方程组，根据线性代数知识，它具有非零解的充分必要条件为系数矩阵的行列式为零，亦即有下式成立。

数值计算方法教学大纲(本)

数值计算方法教学大纲(本) 本着“崇术重用、服务地方”的办学理念和我校“高素质应用型人才”的培养目标，特制定了适合我校工科专业本科生的新教学大纲。一、课程计划课程名称：数值计算方法Numerical Calculation Method 课程定位：数学基础课开课单位：理学院课程类型：专业选修课开设学期：第七学期讲授学时：共15周，每周4学时，共60学时学时安排：课堂教学40学时+实验教学20学时适用专业：计算机、电科、机械等工科专业本科生教学方式：讲授（多媒体为主）+上机考核方式：考试60%+上机实验30%+平时成绩10% 学分：3学分与其它课程的联系预修课程：线性代数、微积分、常微分方程、计算机高级语言等。后继课程：偏微分方程数值解及其它专业课程。二、课程介绍数值计算方法也称为数值分析，是研究用计算机求解各种数学问题的数值方法及其理论的一门学科。随着计算科学与技术的进步和发展，科学计算已经与理论研究、科学实验并列成为进行科学活动的三大基本手段，作为一门综合性的新科学，科学计算已经成为了人们进行科学活动必不可少的科学方法和工具。数值计算方法是科学计算的核心内容，它既有纯数学高度抽象性与严密科学性的特点，又有应用的广泛性与实际实验的高度技术性的特点，是一门与计算机使用密切结合的实用性很强的数学课程.主要介绍插值法、函数逼近与曲线拟合、线性方程组迭代解法、数值积分与数值微分、非线性方程组解法、常微分方程数值解以及矩阵特征值与特征向量数值计算，并特别加强实验环节的训练以提高学生动手能力。通过本课程的学习，不仅能使学生初步掌握数值计算方法的基本理论知识，了解算法设计及数学建模思想，而且能使学生具备一定的科学计算能力和分析与解决问题的能力，不仅为学习后继课程打下良好的理论基础，也为将来从事科学计算、计算机应用和科学研究等工作奠定必要的数学基础。科学计算是21世纪高层次人才知识结构中不可缺少的一部分，它潜移默化地影响着人们的思维方式和思想方法，并提升一个人的综合素质。

数值计算方法设计论文

课程设计（论文）题目: 三次样条插值问题学院: ___ 理学院 _ 专业: __ _ 数学与应用数学班级：数学08-2班学生姓名: 魏建波学生学号: 080524010219 指导教师：李文宇 2010年12月17日

课程设计任务书

目录摘要……………………………………………………………………… 一、前言………………………………………………………………… （一）Lagrange插值的起源和发展过程……………………………………… （二）本文所要达到的目的……………………………………………………… 二、插值函数…………………………………………………………… （一）函数插值的基本思想…………………………………………………… （二）Lagrange插值的构造方法……………………………………………… 三、MATLAB程序………………………………………………………… （一）Lagrange程序…………………………………………………………… （二）龙格程序………………………………………………………………… 四、理论证明…………………………………………………………… 五、综述……………………………………………………………………谢辞………………………………………………………………………参考文献…………………………………………………………………

摘要

前言要求：500字以上，宋体小四，行距20磅，主要内容写该算法的产生及发展、应用领域等。题目整体要求：报告页数，正文在8页以上字体：宋体小四（行距20磅）内容：1、理论依据 2、问题描述 3、问题分析 4、求解计算（程序） 5、结论注：（1）页码编号从正文页开始（2）标题可根据情况自己适当改动示例见下： 2判别…………………… 2.1 判……………… 2.1.1 判别……………… 所谓的判别分析，………………………………………………方法[3]。 2.1.2 判………………………… 常用的有四种判别方法：…………………………………………………步判别法[6]。 1. 马氏………………

数值分析论文

题目：论数值分析在数学建模中的应用学院: 机械自动化学院专业: 机械设计及理论学号: 学生姓名: 日期: 2011年12月5日

论数值分析在数学建模中的应用摘要为了满足科技发展对科学研究和工程技术人员用数学理论解决实际的能力的要求，讨论了数值分析在数学建模中的应用。数值分析不仅应用模型求解的过程中，它对模型的建立也具有较强的指导性。研究数值分析中插值拟合，解线性方程组，数值积分等方法在模型建立、求解以及误差分析中的应用，使数值分析作为一种工具更好的解决实际问题。关键词数值分析；数学建模；线性方程组；微分方程 the Application of Numerical Analysis in Methmetical Modeling Han Y u-tao 1 Bai Y ang 2 Tian Lu 2 Liu De-zheng 2 （1 College of Science ，Tianjin University of Commerce ，Tianjin ，300134 2 College of Science ，Tianjin University of Commerce ，Tianjin ，300134) Abstract In order to meet the technological scientific researchers who use mathematical theory to solve practical problems, the use of numerical analysis in mathematical modeling is discussed.Numerical analysis not only solve the model,but also relatively guide the model.Research on some numerical methods in numerical analysis which usually used in mathmetical modeling and error analysis will be a better way to solve practical problems. Key Words Numerical Analysis ；Mathematical Modeling; Linear Equations ；differential equation 1. 引言数值分析主要介绍现代科学计算中常用的数值计算方法及其基本原理，研究并解决数值问题的近似解，是数学理论与计算机和实际问题的有机结合[1]。随着科学技术的迅速发展，运用数学方法解决科学研究和工程技术领域中的实际问题，已经得到普遍重视。数学建模是数值分析联系实际的桥梁。在数学建模过程中，无论是模型的建立还是模型的求解都要用到数值分析课程中所涉及的算法,如插值方法、最小二乘法、拟合法等，那么如何在数学建模中正确的应用数值分析内容，就成了解决实际问题的关键。 2. 数值分析在模型建立中的应用在实际中，许多问题所研究的变量都是离散的形式，所建立的模型也是离散的。例如，对经济进行动态的分析时，一般总是根据一些计划的周期期末的指标值判断某经济计划执行的如何。有些实际问题即可建立连续模型，也可建立离散模型，但在研究中，并不能时时刻刻统计它，而是在某些特定时刻获得统计数据。例如，人口普查统计是一个时段的人口增长量，通过这个时段人口数量变化规律建立离散模型来预测未来人口。另一方面，对常见的微分方程、积分方程为了求解，往往需要将连续模型转化成离散模型。将连续模型转化成离散模型，最常用的方法就是建立差分方程。以非负整数k 表示时间，记k x 为变量x 在时刻k 的取值，则称k k k x x x -=?+1为k x 的一阶差分，称k k k k k x x x x x +-=??=?++1222)(为k x 的二阶差分。类似课求出k x 的n 阶差分k n x ?。由k ，k x ，及k x 的差分给出的方程称为差分方程[2]。例如在研究节食与运动模型时，发现人们往往采取节食与运动方式消耗体内存储的脂肪，引起体重下降，达到减肥目的。通常制定减肥计划以周为时间单位比较方便，所以采用差分方程模型进行讨论。记第k 周末体重为)(k w ，第k 周吸收热量为)(k c ，热量转换系数α，代谢消耗系数β，在不考虑运动情况下体重变化的模型