资源数据采集技术方案 公司名称 2011年7月二O一一年七月
目录 第 1 部分概述 (3) 1.1 项目概况 (3) 1.2 系统建设目标 (3) 1.3 建设的原则 (4) 1.3.1 建设原则 (4) 1.4 参考资料和标准 (5) 第 2 部分系统总体框架与技术路线 (5) 2.1 系统应用架构 (6) 2.2 系统层次架构 (6) 2.3 关键技术与路线 (7) 第 3 部分系统设计规范 (9) 第 4 部分系统详细设计 (9)
第 1 部分概述 1.1 项目概况 Internet已经发展成为当今世界上最大的信息库和全球范围内传播知识的主要渠道，站 点遍布全球的巨大信息服务网，为用户提供了一个极具价值的信息源。无论是个人的发展还 是企业竞争力的提升都越来越多地依赖对网上信息资源的利用。 现在是信息时代，信息是一种重要的资源，它在人们的生活和工作中起着重要的作用。 计算机和现代信息技术的迅速发展，使Internet成为人们传递信息的一个重要的桥梁。网络 的不断发展，伴随着大量信息的产生，如何在海量的信息源中查找搜集所需的信息资源成为 了我们今后建设在线预订类旅游网重要的组成部分。 因此，在当今高度信息化的社会里，信息的获取和信息的及时性。而Web数据采集可以通过一系列方法，依据用户兴趣，自动搜取网上特定种类的信息，去除无关数据和垃圾数据，筛选虚假数据和迟滞数据，过滤重复数据。直接将信息按照用户的要求呈现给用户。可 以大大减轻用户的信息过载和信息迷失。 1.2 系统建设目标 在线预订类旅游网是在线提供机票、酒店、旅游线路等旅游商品为主，涉及食、住、行、游、购、娱等多方面的综合资讯信息、全方位的旅行信息和预订服务的网站。 如果用户要搜集这一类网站的相关数据，通常的做法是人工浏览网站，查看最近更新的信息。然后再将之复制粘贴到Excel文档或已有资源系统中。这种做法不仅费时费力，而且 在查找的过程中可能还会遗漏，数据转移的过程中会出错。针对这种情况，在线预订类旅游网信息自动采集的系统可以实现数据采集的高效化和自动化。
第六章习题解答 2、利用梯形公式和Simpson 公式求积分2 1 ln xdx ? 的近似值，并估计两种方法计算值的最大 误差限。 解：①由梯形公式： 21ln 2 ()[()()][ln1ln 2]0.3466222 b a T f f a f b --= +=+=≈ 最大误差限 3''2 ()111 ()()0.0833******** T b a R f f ηη-=-=≤=≈ 其中，(1,2)η∈ ②由梯形公式： 13()[()4()()][ln14ln()ln 2]0.38586262 b a b a S f f a f f b -+= ++=++≈ 最大误差限 5(4)4()66 ()()0.0021288028802880 S b a R f f ηη-=-=≤≈， 其中，(1,2)η∈。 4、推导中点求积公式 3''()()()()() ()224 b a a b b a f x dx b a f f a b ξξ+-=-+<
粘弹性结构动力学分析中的一种数值方法 彭 凡 傅 衣 铭 (湖南大学工程力学系, 长沙 410082, 中国) 摘要：针对材料具积分型本构关系，及松弛模量为Prony 级数形式的粘弹性结构动力学问题，本文结合Newmark 方法与Taylor 方法，建立了计算该类问题的一种数值算法。且以简支梁为例，应用该方法具体地分析了考虑线性与非线性粘弹性时梁的强迫振动响应。 关键词：粘弹性 动力学 数值方法 响应 1 引言 随着人们对结构材料物理与力学性质了解的不断深入，以及新型材料的广泛应用，粘弹性结构的动力学研究受得了愈来愈多的重视，数值计算已成为一种主要的分析手段。文[1]基于Newmark 方法建立了粘弹性结构动力学响应的有限元法，但只涉及到线性问题，而且在每一计算步，卷积积分的计算量较大。桂洪斌等[2]提出将粘弹性结构的动力学方程进行Laplace 变换，然后在相域中求解问题，显然这种处理方式同样只适应于线性情况。当考虑几何，物理包括损伤等非线性因素时，粘弹性结构动力学的数值分析就变得十分复杂与困难了。文[3,4]通过将微分－积分型非线性动力学方程化成高阶的微分方程，最终由Runge-Kutta 法来获得数值解，但只有当材料为标准线性固体或Prony 级数取较少项数时，这种方法才比较容易实现。本文针对材料服从积分型本构关系，且松弛模量为Prony 级数形式的粘弹性结构动力学问题，建立了从时域内直接求解的数值算法，它是基于Newmark 方法与Taylor 方法而得出的。其中Taylor 方法为卷积积分的递规算法，能使计算量显著降低[5]。文中通过对粘弹性梁的受迫振动分析来说明方法的应用。 2 简支粘弹性梁受迫振动的动力学方程 考虑一简支梁，其跨度为L ，高为h ，中点受横向周期激励t H θsin 。设材料具非线性粘弹性，可由Leaderman [6]本构关系描叙，则有 00 () ()(())(())()t E t t E g t g d t τσεεττ τ?-=+?-? (1) 式中)0(0E E =，)(t E 为松弛函数，)(εg 为应变ε的非线性函数： 23()g εεβεγε=++ (2) 其中β与γ为常数。在小挠度情况下，梁的受迫振动方程为： ()3 23452202422 334522422 0(,)(,)(,) 1280()(,)(,) sin ()1280t w x t h w x t h w x t A E t x x x E t h w x h w x d H x L t t x x x ργτττγτδθτ?? ????????++ ????????????? ???-?????++=- ??-????????? ? （3） 式中A ,ρ分别为梁的质量密度及横截面面积，δ为Dirac 函数， 满足两个简支端条件，即,,(0,)(0,)(,)(,)0xx xx w t w t w L t w L t ====的挠度),(t x w 取为 1(,)()sin k k k x w x t f t L π∞ ==∑ （4） 为说明问题起见，式(4)中只考虑1=k 的项，且令)()(1t f t f =。将式(4)代入式(3)后，作Galerkin 积分，并记
《数据采集及分析》实验指导书 实验一采样定理 一、实验目的 熟悉信号采样过程，并通过本实验观察欠采样时信号频谱的混迭现象，了解采样前后信号频谱的变化，加深对采样定理的理解，掌握采样频率的确定方法。 二、实验原理 模拟信号经过(A/D) 变换转换为数字信号的过程称之为采样，信号采样后其频谱产生了周期延拓，每隔一个采样频率fs，重复出现一次。为保证采样后信号的频谱形状不失真，采样频率必须大于信号中最高频率成份的两倍，这称之为采样定理。 a) 正常采样b)欠采样 图1.1 采样信号的频混现象 需要注意的是，在对信号进行采样时，满足了采样定理，只能保证不发生频率混叠，对信号的频谱作逆傅立叶变换时，可以完全变换为原时域采样信号，而不能保证此时的采样信号能真实地反映原信号。工程实际中采样频率通常大于信号中最高频率成分的3到5倍。 三、实验仪器和设备 1. 计算机 n台 2. 实验软件 1套 四、实验步骤及内容 1. 启动计算机。 2. 启动实验软件。
图1.2 采样定理实验 3. . 点击"采样定理"实验中的"正弦波"按钮，产生正弦波信号，然后选择不同的采样抽取率，分析和观察信号的时域波形与频谱的变化。 4. 点击"采样定理"实验中的"方波"按钮，产生方波信号，然后选择不同的采样抽取率，分析和观察信号的时域波形与频谱的变化。 5. 点击"采样定理"实验中的"三角波"按钮，产生三角波信号，然后选择不同的采样抽取率，分析和观察信号的时域波形与频谱的变化。 五、实验报告要求 1. 简述实验目的和原理。 2. 按实验步骤附上相应的信号波形和频谱曲线，说明采样频率的变化对信号时域和频域特性的影响，总结实验得出的主要结论。 六、思考题 1.为什么在实际测量中采样频率通常要大于信号中最高频率成分的3到5倍？
结构动力学历年试题（简答题） 1.根据荷载随时间的变化规律，动力荷载可以划分为哪几类？每一类荷载包括哪几种，请 简述每一种荷载的特点。P2 2.通过与静力问题的对比，试说明结构动力计算的特点。P3 3.动力自由度数目计算类 4.什么叫有势力？它有何种性质。P14 5.广义力是标量还是矢量？它与广义坐标的乘积是哪个物理量的量纲？P16 6.什么是振型的正交性？它的成立条件是什么？P105 7.在研究结构的动力反应时，重力的影响如何考虑？这样处理的前提条件是什么？P32 8.对于一种逐步积分计算方法，其优劣性应从哪些方面加以判断？P132 9.在对结构动力反应进行计算的思路上，数值积分方法与精确积分方法的差异主要表现在 哪里？第五章课件 10.利用Rayleigh法求解得到的振型体系的基本振型和频率及高阶振型和频率与各自的精确 解相比有何特点？造成这种现象的原因何在？P209 11.根据荷载是否预先确定，动荷载可以分为哪两类？它们各自具有怎样的特点？P1 12.坐标耦联的产生与什么有关，与什么无关？P96 13.动力反应的数值分析方法是一种近似的计算分析方法，这种近似性表现在哪些方面？ P132及其课件 14.请给出度哈姆积分的物理意义？P81 15.结构地震反应分析的反应谱方法的基本原理是什么？P84总结 16.某人用逐步积分计算方法计算的结构位移，得到如下的位移时程的计算结果：。。。 17.按照是否需要联立求解耦联方程组，逐步积分法可以分为哪两类？这两类的优劣性应该 如何进行判断？P132 18.根据荷载随时间的变化规律，动力荷载可以划分为哪几类？每一类荷载又包括哪些类型， 每种类型请给出一种实例。P2 19.请分别给出自振频率与振型的物理意义？P103 20.振型叠加法的基本思想是什么？该方法的理论基础是什么？P111参考25题 21.在振型叠加法的求解过程中，只需要取有限项的低阶振型进行分析，即高阶振型的影响 可以不考虑，这样处理的物理基础是什么？P115 22.我们需要用数值积分方法求解一座大型的高坝结构的地震反应时程，动力自由度的总数 为25000个，我们如何缩短计算所耗费的机时？P103 23.什么是结构的动力自由度？动力自由度与静力自由度的区别何在？P11及卷子上答案 24.一台转动机械从启动到工作转速正好要经过系统的固有频率（又称为转子的临界转速）， 为减小共振，便于转子顺利通过临界转速，通常采用什么措施比较直接有效？简要说明理由。详解见卷子上答案 25.简述用振型叠加法求解多自由度体系动力响应的基本原理及使用条件分别是什么？若 振型叠加法不适用，可采用何种普遍适用的方法计算体系响应？详解见卷子上答案 26.振型函数边界条件。。。 27.集中质量和一致质量有限元的差异和优缺点，采用这两种有限元模型给出的自振频率与 实际结构自振频率相比有何种关系？P242及卷子上答案 28.人站在桥上可以感觉到桥面的震动，简述当车辆行驶在桥上和驶离桥面的主要振型特征 有何不同？ 29.简述用Duhamel积分法求体系动力响应的基本原理，以及积分表达式中的t和τ有何差
数据采集方法有哪些 数据采集数据采集（DAQ），是指从传感器和其它待测设备等模拟和数字被测单元中自动采集非电量或者电量信号，送到上位机中进行分析，处理。数据采集系统是结合基于计算机或者其他专用测试平台的测量软硬件产品来实现灵活的、用户自定义的测量系统。 数据采集，又称数据获取，是利用一种装置，从系统外部采集数据并输入到系统内部的一个接口。数据采集技术广泛应用在各个领域。比如摄像头，麦克风，都是数据采集工具。被采集数据是已被转换为电讯号的各种物理量，如温度、水位、风速、压力等，可以是模拟量，也可以是数字量。采集一般是采样方式，即隔一定时间（称采样周期）对同一点数据重复采集。采集的数据大多是瞬时值，也可是某段时间内的一个特征值。准确的数据测量是数据采集的基础。数据量测方法有接触式和非接触式，检测元件多种多样。不论哪种方法和元件，均以不影响被测对象状态和测量环境为前提，以保证数据的正确性。数据采集含义很广，包括对面状连续物理量的采集。在计算机辅助制图、测图、设计中，对图形或图像数字化过程也可称为数据采集，此时被采集的是几何量（或包括物理量，如灰度）数据。 在互联网行业快速发展的今天，数据采集已经被广泛应用于互联网及分布式领域，数据采集领域已经发生了重要的变化。首先，分布式控制应用场合中的智能数据采集系统在国内外已经取得了长足的发展。其次，总线兼容型数据采集插件的数量不断增大，与个人计算机兼容的数据采集系统的数量也在增加。国内外各种数据采集机先后问世，将数据采集带入了一个全新的时代。 现在谈论大数据已经没有新意了，形形色色的产品、平台和公司都贴满大数据标签，但大数据却并没有掀起预期飓风，甚至还被冠以“伪命题”污名。 本末倒置，数据采集才是大数据产业的基石。都在说大数据应用、大数据价值挖掘，却不想，没有数据何来应用、价值一说。就好比不开采石油，一味想得到汽油。当然，石油开采并不容易，各行各业包括政府部门的信息化建设都是封闭式进行，海量数据被封在不同
结构力学简答题 1、结构动力分析的目的：是确定结构在动力荷载作用下的内力和变形，并通过动力分析确定结构的动力特性。 1、动力荷载的类型：（1）是否随时间变化:静荷载和动荷载（2）是否已预先确定：确定性荷载和非确定性荷载（3）随时间变化的规律：周期荷载：简谐荷载和非简谐周期荷载；非周期荷载：冲击荷载和一般任意荷载。 2、结构动力计算的特点：（1）动力反应要计算全部时间点上的一系列解，比静力计算复杂且要消耗很多的计算时间。（2）由于动力反应中结构的位置随时间迅速变化，从而产生惯性力，惯性力对结构的反应又产生重要影响。 3、结构离散化的方法：集中质量法、广义坐标法、有限元法。本质是无限自由度问题转化为有限自由度的过程。 4、有限元法：（1）与广义坐标法相似，有限元法采用了形函数的概念，但不同于广义坐标法在全部体系上插值，而是采用了分片的插值，因此形函数的表达式可以相对简单。（2）与集中质量法相比，有限元法中的广义坐标也采用了真实的物理量，具有直接、直观的优点，与集中质量法相同。 5、广义坐标：能决定质点系几何位置的彼此独立的量。选择原则：解题方便。 6、动力自由度：结构体系在任意瞬时的一切可能的变形中，决定全部质量位置所需的独立参数的数目。动力自由度不完全取决于质点的数目，也与结构是否静定有关。静力自由度：确定体系在空间中的位置所需的独立参数的数目。前者是由于系统的弹性变形而引起的各质点的位移分量，后者是指结构中的刚体由于约束不足而产生的刚体位移。 7、有势力：（1）每一个力的大小和方向只决定于体系所有各质点的位置。（2）体系从某一位置到另一位置所做的功只决定于质点的始末位置，而与路径无关。（3）沿任何封闭路线所作的功为零。 8、实位移：如果位移不仅满足约束方程，而且满足运动方程和初始条件，则称为体系的实位移。可能位移：满足所有约束方程的位移称为体系的可能位移。虚位移：在某一固定时刻，体系在约束许可的情况下产生的任意组微小位移。三者关系：实位移即为体系的真实位移，它必是可能位移中的一员。虚位移与可能位移的区别在于虚位移是约束冻结后许可产生的微小位移。当对于约束方程中不显含时间的稳定约束体系中虚位移与可能位移相同时，实位移必与某一虚位移重合。 9、广义力：广义力是标量而非矢量，广义力与广义坐标的乘积具有功的量纲。 10、阻尼力：引起结构能量的耗散，使结构振幅逐渐变小的这种作用。 产生原因：（1）固体材料变形时的内摩擦，或材料快速应变引起的热耗散。 （2）结构连接部位的摩擦，结构与非结构之间的摩擦。（3）结构周围外部介质引起的阻尼。 11、四种建立运动方程方法的特点 （1）达朗贝尔原理：矢量方法，直观，建立了动平衡概念。 （2）虚位移原理：半矢量方法，可处理复杂分布质量和弹性问题。 （3）哈密顿原理：标量方法，表达简洁。 （4）拉格朗日方法：标量方法，运用面广。 13、进行结构动力分析计算时，重力的影响如何考虑？这样处理的前提条件是什么？ 如果重力在动荷载作用前被弹簧预先平衡，则在研究结构的动力反应时可以完全不考虑重力的影响。建立体系的运动方程，直接解出体系的动力解。若未被预先平衡，则需考虑重力的影响，应用叠加原理将动静问题分开计算，将结果相加即得到结构的真实反应。这样做的前提条件是结构是线弹性且处于小变形范围之内。重力问题的分析和动力问题的分析可以分开讨论。
数据采集统计分析方法 目的：为检验员检验数据收集提供方法 适用范围：本公司内部对产品进行检验从而得到检验数据，为管理评审提供依据。 可用以下方法做为参考 QC旧七种工具 排列图，因果图，散布图，直方图，控制图，检查表与分层法 QC新七种工具(略) 关联图，KJ法，系统图法，矩阵图法，矩阵数据解析法，过程决策程序图法（PDPC)和箭头图法。 数据统计分析方法-排列图 数据统计分析方法-排列图 排列图是由两个纵坐标，一个横坐标，若干个按高低顺序依次排列的长方形和一条累计百分比折线所组成 的，为寻找主要问题或主要原因所使用的图。 例1： 排列图的优点 排列图有以下优点： 直观，明了－－全世界品质管理界通用 用数据说明问题－－说服力强 用途广泛：品质管理/ 人员管理/ 治安管理 排列图的作图步骤 收集数据(某时间)
作缺陷项目统计表 绘制排列图 画横坐标(标出项目的等分刻度) 画左纵坐标(表示频数) 画直方图形(按每项的频数画) 画右纵坐标(表示累计百分比) 定点表数，写字 数据统计分析方法-因果图 何谓因果图： 对于结果（特性）与原因（要因）间或所期望之效果（特性）与对策的关系，以箭头连接，详细分析原因 或对策的一种图形称为因果图。 因果图为日本品管权威学者石川馨博士于1952年所发明，故又称为石川图，又因其形状似鱼骨，故也可称 其为鱼骨图，或特性要因图 作因果图的原则 采取由原因到结果的格式 通常从‘人，机，料，法，环’这五方面找原因 ‘4M1E’, Man, Machine, Material, Method, Environment 通常分三个层次：主干线、支干线、分支线 尽可能把所有的原因全部找出来列上 对少数的主要原因标上特殊的标志 写上绘制的日期、作者、有关说明等
计算流体力学常用数值方法简介 李志印 熊小辉 吴家鸣 (华南理工大学交通学院) 关键词 计算流体力学 数值计算 一 前 言 任何流体运动的动力学特征都是由质量守恒、动量守恒和能量守恒定律所确定的,这些基本定律可以由流体流动的控制方程组来描述。利用数值方法通过计算机求解描述流体运动的控制方程,揭示流体运动的物理规律,研究流体运动的时一空物理特征,这样的学科称为计算流体力学。 计算流体力学是一门由多领域交叉而形成的一门应用基础学科,它涉及流体力学理论、计算机技术、偏微分方程的数学理论、数值方法等学科。一般认为计算流体力学是从20世纪60年代中后期逐步发展起来的,大致经历了四个发展阶段:无粘性线性、无粘性非线性、雷诺平均的N-S方程以及完全的N-S方程。随着计算机技术、网络技术、计算方法和后处理技术的迅速发展,利用计算流体力学解决流动问题的能力越来越高,现在许多复杂的流动问题可以通过数值计算手段进行分析并给出相应的结果。 经过40年来的发展,计算流体力学己经成为一种有力的数值实验与设计手段,在许多工业领域如航天航空、汽车、船舶等部门解决了大量的工程设计实际问题,其中在航天航空领域所取得的成绩尤为显著。现在人们已经可以利用计算流体力学方法来设计飞机的外形,确定其气动载荷,从而有效地提高了设计效率,减少了风洞试验次数,大大地降低了设计成本。此外,计算流体力学也己经大量应用于大气、生态环境、车辆工程、船舶工程、传热以及工业中的化学反应等各个领域,显示了计算流体力学强大的生命力。 随着计算机技术的发展和所需要解决的工程问题的复杂性的增加,计算流体力学也己经发展成为以数值手段求解流体力学物理模型、分析其流动机理为主线,包括计算机技术、计算方法、网格技术和可视化后处理技术等多种技术的综合体。目前计算流体力学主要向二个方向发展:一方面是研究流动非定常稳定性以及湍流流动机理,开展高精度、高分辩率的计算方法和并行算法等的流动机理与算法研究;另一方面是将计算流体力学直接应用于模拟各种实际流动,解决工业生产中的各种问题。 二 计算流体力学常用数值方法 流体力学数值方法有很多种,其数学原理各不相同,但有二点是所有方法都具备的,即离散化和代数化。总的来说其基本思想是:将原来连续的求解区域划分成网格或单元子区
数值分析 报告 班级： 专业： 流水号： 学号： 姓名：
常用的插值方法 序言 在离散数据的基础上补插连续函数，使得这条连续曲线通过全部给定的离散数据点。插值是离散函数逼近的重要方法，利用它可通过函数在有限个点处的取值状况，估算出函数在其他点处的近似值。 早在6世纪，中国的刘焯已将等距二次插值用于天文计算。17世纪之后，牛顿、拉格朗日分别讨论了等距和非等距的一般插值公式。在近代，插值法仍然是数据处理和编制函数表的常用工具，又是数值积分、数值微分、非线性方程求根和微分方程数值解法的重要基础，许多求解计算公式都是以插值为基础导出的。 插值问题的提法是：假定区间[a，b〕上的实值函数f（x）在该区间上 n＋1 个互不相同点x 0，x 1 (x) n 处的值是f（x ），……f（x n ），要求估算f（x）在[a，b〕 中某点的值。其做法是：在事先选定的一个由简单函数构成的有n＋1个参数C ， C 1，……C n 的函数类Φ(C ，C 1 ，……C n )中求出满足条件P（x i ）＝f（x i ）（i＝0，1，…… n）的函数P(x)，并以P(x)作为f(x)的估值。此处f（x）称为被插值函数，x 0，x 1 ，……xn 称为插值结(节)点，Φ(C 0，C 1 ，……C n )称为插值函数类，上面等式称为插值条件， Φ（C 0，……C n ）中满足上式的函数称为插值函数，R（x）＝ f（x）－P（x）称为 插值余项。
求解这类问题，它有很多种插值法，其中以拉格朗日(Lagrange)插值和牛顿(Newton)插值为代表的多项式插值最有特点，常用的插值还有Hermit 插值，分段插值和样条插值。 一．拉格朗日插值 1.问题提出： 已知函数()y f x =在n+1个点01,, ,n x x x 上的函数值01,, ,n y y y ，求任意一点 x '的函数值()f x '。 说明：函数()y f x =可能是未知的；也可能是已知的，但它比较复杂，很难计算其函数值()f x '。 2.解决方法： 构造一个n 次代数多项式函数()n P x 来替代未知（或复杂）函数()y f x =，则 用()n P x '作为函数值()f x '的近似值。 设()2012n n n P x a a x a x a x =+++ +，构造()n P x 即是确定n+1个多项式的系数 012,,,,n a a a a 。 3.构造()n P x 的依据： 当多项式函数()n P x 也同时过已知的n+1个点时，我们可以认为多项式函数 ()n P x 逼近于原来的函数()f x 。根据这个条件，可以写出非齐次线性方程组： 20102000 201121112012n n n n n n n n n n a a x a x a x y a a x a x a x y a a x a x a x y ?+++ +=?++++=??? ?+++ +=? 其系数矩阵的行列式D 为范德萌行列式： () 200021110 2 111n n i j n i j n n n n x x x x x x D x x x x x ≥>≥= = -∏
第六章 习 题 1. 计算下列矩阵的1A ，2A ，A ∞三种范数。 （1）1101A ???=????，（2）312020116A ????=??????? . 2. 用Jacobi 方法和Gauss-Seidel 迭代求解方程组 1231231 238322041133631236x x x x x x x x x ?+=??+?=??++=? 要求取(0)(0,0,0)T x =计算到(5)x ，并分别与精确解(3,2,1)T x =比较。 3. 用Gauss-Seidel 迭代求解 12312312 35163621122x x x x x x x x x ??=??++=???+=?? 以(0)(1,1,1)T x =?为初值，当(1)() 310k k x x +?∞?<时，迭代终止。 4. 已知方程组121122,2,x x b tx x b +=?? +=? （1）写出解方程组的Jacobi 迭代矩阵，并讨论迭代收敛条件。 （2）写出解方程组的Gauss-Seidel 迭代矩阵，并讨论迭代收敛条件. 5. 设有系数矩阵 122111221A ?????=?????? ， 211111112B ?????=??????? ， 证明：（1）对于系数矩阵A ，Jacobi 迭代收敛，而Gauss-Seidel 迭代不收敛. （2）对于矩阵B ，. 6. 讨论方程组 112233302021212x b x b x b ?????????????=??????????????????? 用Jacobi 迭代和Gauss-Seidel 迭代的收敛性；如果都收敛，比较哪种方法收敛更快.
数值分析重点考察内容 第一章： 基本概念 第二章： Gauss消去法，Lu分解法 第三章： 题型：具体题＋证明，误差分析 三个主要迭代法，条件误差估计，范数的小证明 第四章： 掌握三种插值方法：拉格朗日，牛顿，厄尔米特，误差简单证明，构造复合函数 第五章： 最小二乘法计算 第六章： 梯形公式，辛普森（抛物线）公式，高斯公式三个重要公式，误差分析。 高斯求积公式的构造 第七章： 几种常用的迭代格式构造，收敛性证明。 第九章： 基本概念（收敛阶，收敛条件，收敛区域等）, 简单欧拉法。
第一章 误差 1. 科学计算中的误差来源有4个，分别是________，________，________，________。 2. 用Taylor 展开近似计算函数000()()'()()f x f x f x x x ≈+-，这里产生是什么误差？ 3. 0.7499作34 的近似值，是______位有效数字，65.380是舍入得到的近似值，有____几位有效数字，相对误差限为_______. 0.0032581是四舍五入得到的近似值，有_______位有效数字. 4. 改变下列表达式，使计算结果比较精确： （1） 11,||1121x x x x --++ （2） ||1x （3） 1cos ,0,|| 1.x x x x -≠ （4） sin sin ,αβαβ-≈ 5. 采用下列各式计算61)时，哪个计算效果最好？并说明理由。 （1） （2） 99-（3） 6(3- （4 6. 已知近似数*x 有4位有效数字，求其相对误差限。 上机实验题： 1、利用Taylor 展开公式计算 0!k x k x e k ∞ ==∑，编一段小程序，上机用单精度计算x e 的函数 值. 分别取 x =1，5，10，20，-1，-5，-10，-15，-20，观察所得结果是否合理，如不合理请分析原因并给出解决方法. 2、已知定积分10,0,1,2,,206 n n x I dx n x ==+? ,有如下的递推关系 1111100(6)61666 n n n n n x x x x I dx dx I x x n ---+-===++-?? 可建立两种等价的计算公式 (1) 1016,0.154n n I I I n -=-=取； (2) 12011),0.6n n I nI I n -=-=（取 来计算123419,,,,,I I I I I ，编程比较哪种计算的数值结果好，并给出理论分析。
5章 动力反应的数值计算 如果激励[作用力)(t p 或地面加速度)(t u g ]是随时间任意变化的，或者体系是非线性的，那么对单自由度体系的运动方程进行解析求解通常是不可能的。这类问题可以通过数值时间步进法对微分方程进行积分来处理。在应用力学广阔的学科领域中，有关各种类型微分方程数值求解方法的文献（包括几部著作中的主要章节）浩如烟海，这些文献包括这些方法的数学进展以及它们的精度、收敛性、稳定性和计算机实现等问题。 然而，本章仅对在单自由度体系动力反应分析中特别有用的很少几种方法进行简要介绍，这些介绍仅提供这些方法的基本概念和计算算法。尽管这些对许多实际问题和应用研究已经足够了，但是读者应该明白，有关这个主题存在大量的知识。 5.1 时间步进法 对于一个非弹性体系，欲采用数值求解的运动方程为 )(),(t p u u f u c u m s =++ 或者 )(t u m g - (5.1.1) 初始条件 )0(0u u = )0(0u u = 假定体系具有线性粘滞阻尼，不过，也可以考虑其他形式的阻尼（包括非线性阻尼），后面会明显看到这一点。然而由于缺乏阻尼信息．因此很少这样做，特别是在大振幅运动时。作用力)(t p 由一系列离散值给出: )(i i t p p = ,0=i 到N 。时间间隔 i i i t t t -=?+1 （5.1.2）
图5.1.1 时间步进法的记号 通常取为常数，尽管这不是必需的。在离散时刻i t （表示为i 时刻）确定反 应，单自由度体系的位移、速度和加速度分别为i u 、i u 和i u 。假定这些值是已知的，它们在i 时刻满足方程 i i s i i p f u c u m =++)( （5.1.3） 式中，i s f )(是i 时刻的抗力，对于线弹性体系，i i s ku f =)(，但是如果体系是非弹性的，那么它会依赖于i 时刻以前的位移时程和速度。将要介绍的数值方 法将使我们能够确定i +1时刻满足方程(5.1.1)的反应1+i u 、1+i u 和1+i u ，即在i +1时刻 1111)(++++=++i i s i i p f u c u m （5.1.4） 对于i =0，1，2，3，…，连续使用时间步进法，即可给出i =0，l ，2，3，… 所有瞬时所需的反应。已知的初始条件)0(0u u =)0(0u u =和提供了起动该方法的必要信息。 从i 时刻到i +1时刻的步进一般不是精确的方法，许多在数值上可以实现的近似方法是可能的。对于数值方法，有三个重要的要求：(1)收敛性一随着时间步长的减少，数值解应逼近精确解；(2)稳定性一在存在数值舍入误差的情况下，数值解应是稳定的；(3)精度一数值方法应提供与精确解足够接近的结果。这些重要的问题在本书中均作简要的讨论，全面的论述可在着重微分方程数值解法的书中找到。 本章介绍三种类型的时间步进法：(1)基于激励函数插值的方法；(2)基于速度和加速度有限差分表达的方法；(3)基于假设加速度变化的方法。前两类中各
一 1、结构动力学计算的特点？ （对比静力问题）○ 1动力反应要计算全部时间点上的一系列的解，比静力问题复杂要消耗更多的计算时间。○ 2与静力问题相比，由于动力反应中结构的位置随时间迅速变化，从而产生惯性力，惯性力对结构的反应又产生重要的影响。 2、结构动力学是研究什么的？包含什么内容？ 结构动力学：是研究结构体系的动力特性及其在动力荷载作用下的动力反应分析原理和 方法的一门理论和技术学科。 目的：在于为改善工程结构体系在动力环境中的安全性和可靠性提供坚实的理论基础。 二、 1、动力系数（有阻尼、无阻尼。简谐、半功率点法、位移计……） 2、动力系数和哪些因素有关 动力放大系数受阻尼比控制，Rd 曲线形状可以反映出阻尼比的影响。主要有两点：其一是峰值大小；其二是曲线的胖瘦。 3、动力系数在工程（隔震、调频减震）的应用 4、如何用动力系数测阻尼比 三、 1、阻尼 阻尼也称阻尼力，是引起结构能量的耗散，使结构振幅逐渐变小的作用。 阻尼的来源：1固体材料变形时的内摩擦，或材料快速反应引起的热耗散；2结构连接部位的摩擦；3结构周围外部介质引起的阻尼。 2.阻尼比常用的测量方法及其优缺点： （1）对数衰减率法：相邻振动峰值比的自然对数值称为对数衰减率。采用自由振动试验，测一阶振型的阻尼比较容易。测量高阶振型阻尼比的关键是能激发出按相应振型的自由振动。 （2） 共振放大法：采用强迫振动试验，通过共振得到（Rd ）max 由于静荷载下的位移较难确定，应用上存在一定的技术困难，但通过一定数学上的处理还是可以用的。（Ust 是零频 时的静位移，不容易测得。） （3） 半功率点（带宽）法：采用强迫振动试验，测出Rd-w/wn 图上振 幅值等于倍最大振幅的点，对应的长度的1/2即为阻尼比。不但能用于单自由度体系，也可以用于多自由度体系，对多自由度体系要求共振频率稀疏，即多个自振频率应相隔较远，保证在确定相应于某一自振频率的半功率点时不受相邻自振频率的影响。 3、等效粘滞阻尼比 ○1、粘性阻尼是一种理想化的阻尼，具有简单和便于分析计算的优点。○ 2工程中结构的阻尼源于多方面，其特点和数学描述更为复杂，这时可以将复杂的阻尼在一定的意义上等效 成粘性阻尼。○3一般采用基于能量等效的原则。○4阻尼耗散能量的大小可以用阻尼力的滞回曲线反映。 m st d u u R 0max 2)(21=≈ζn k k ln 21+≈y y n πξn a b f f f 2-=ζ
审计数据采集与分析技术 计算机审计的含义 ?计算机审计有3层含义： –面向数据的审计 –面向现行信息系统的审计 –对信息系统生命周期的审计 面向数据的计算机审计流程 ?审前调查：电子数据的组织、处理和存储 ?数据采集：审计接口、数据库访问技术、数据采集技术 ?数据清理、转换、验证、建立中间表 ?数据分析：数据分析技术、SQL、审计软件 ?审计取证 一、审前调查及电子数据的组织、处理和存储 1.审前调查的内容和方法 ?对组织结构调查 ?对计算机信息系统的调查 ?提出数据需求 2.电子数据的组织、处理和存储 电子数据处理的特点 ?存储介质改变 ?基于一定的数据处理平台，有一定的数据模型 ?数据表示编码化（各种编码） ?带来系统控制和数据安全性的新问题 ?审计线索改变 如何表示数据 ?数据类型与数据取值 –数据类型决定了取值范围与运算范围 ?数据模型 –数据模型是对现实世界数据特征的抽象 –它提供模型化数据和信息的工具
数据模型的2个层次 ? ?概念模型 –E-R模型的要素 ?实体：客观存在并可以相互区分的事物，用方框表示 ?属性：实体的特征或性质，用椭圆表示 ?联系：实体之间的联系，用菱形表示 ?数据模型 –关系模型 –层次模型 –网状模型 数据模型的3个要素 ?数据结构 –描述模型的静态特征 –是刻画数据模型最重要的方面 ?数据操作 –描述模型的动态特性 ?数据检索 ?数据更新（增加、删除、修改） ?约束条件 –一组完整性规则的集合 ?实体完整性 ?引用（参照）完整性 ?用户定义的完整性 关系模型 ?关系模型是目前最常用的一种数据模型 ?关系数据库采用关系模型作为数据的组织方式 ?关系模型建立在严格的关系代数基础之上 ?关系模型概念单一，用关系表示实体以及实体之间的联系?关系数据库的标准语言SQL是一种非过程化语言，使用方便关系模型的数据结构 ?关系 –一张二维表，每一列都不可再分 –表中的行、列次序并不重要 ?元组 –二维表中的每一行，相当于一条记录 ?属性 –二维表中的每一列，属性有名称与类型。 –属性不可再分，不允许重复 ?主键 –由表中的属性或属性组组成，用于唯一确定一条记录?域
数值分析习题集 （适合课程《数值方法A》和《数值方法B》） 长沙理工大学 第一章绪论 1.设x>0,x的相对误差为δ,求的误差. 2.设x的相对误差为2％,求的相对误差. 3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指 出它们是几位有效数字: 4.利用公式求下列各近似值的误差限: 其中均为第3题所给的数. 5.计算球体积要使相对误差限为1％,问度量半径R时允许的相对误差限是多少? 6.设按递推公式 ( n=1,2,…) 计算到.若取≈(五位有效数字),试问计算将有多大误差? 7.求方程的两个根,使它至少具有四位有效数字(≈. 8.当N充分大时,怎样求? 9.正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝? 10.设假定g是准确的,而对t的测量有±秒的误差,证明当t增加时S的绝对误差增加,而 相对误差却减小. 11.序列满足递推关系(n=1,2,…),若(三位有效数字),计算到时误差有多大?这个计算过程 稳定吗? 12.计算,取,利用下列等式计算,哪一个得到的结果最好? 13.,求f(30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一等价公式 计算,求对数时误差有多大? 14.试用消元法解方程组假定只用三位数计算,问结果是否可靠? 15.已知三角形面积其中c为弧度,,且测量a ,b ,c的误差分别为证明面积的误差满足 第二章插值法 1.根据定义的范德蒙行列式,令 证明是n次多项式,它的根是,且 .
2.当x= 1 , -1 , 2 时, f(x)= 0 , -3 , 4 ,求f(x)的二次插值多项式. 3. 4.给出cos x,0°≤x ≤90°的函数表,步长h =1′=(1/60)°,若函数表具有5位有效数 字,研究用线性插值求cos x 近似值时的总误差界. 5.设,k=0,1,2,3,求. 6.设为互异节点(j=0,1,…,n),求证: i) ii) 7.设且,求证 8.在上给出的等距节点函数表,若用二次插值求的近似值,要使截断误差不超过,问使用函 数表的步长应取多少? 9.若,求及. 10.如果是次多项式,记,证明的阶差分是次多项式,并且为正整数). 11.证明. 12.证明 13.证明 14.若有个不同实根,证明 15.证明阶均差有下列性质: i)若,则; ii)若,则. 16.,求及. 17.证明两点三次埃尔米特插值余项是 并由此求出分段三次埃尔米特插值的误差限. 18.求一个次数不高于4次的多项式,使它满足并由此求出分段三次埃尔米特插值的误差限. 19.试求出一个最高次数不高于4次的函数多项式,以便使它能够满足以下边界条件,,. 20.设,把分为等分,试构造一个台阶形的零次分段插值函数并证明当时,在上一致收敛到. 21.设,在上取,按等距节点求分段线性插值函数,计算各节点间中点处的与的值,并估计误 差. 22.求在上的分段线性插值函数,并估计误差. 23.求在上的分段埃尔米特插值,并估计误差. i) ii) 25.若,是三次样条函数,证明 i); ii)若,式中为插值节点,且,则. 26.编出计算三次样条函数系数及其在插值节点中点的值的程序框图(可用式的表达式). 第三章函数逼近与计算 1.(a)利用区间变换推出区间为的伯恩斯坦多项式. (b)对在上求1次和三次伯恩斯坦多项式并画出图形,并与相应的马克劳林级数部分和误
数据采集方法 市场调查中的数据采集方法：入户访问，拦截访问，邮寄问卷调查，留置问卷调查，电话调查，网络调查，深度访问法，座谈法等。 入户访问: 入户访问是指采用随机抽样方式抽取一定数量的家庭或单位，访问员到抽取出来的家庭或单位，访问员到抽取出来的家庭或单位中进行访问，直接与被访者接触，然后依照问卷活调查提纲进行面对面的直接提问，并记录下对方的答案的调查方式。 拦截访问: 拦截访问是目前十分流行的一种询问调研法，该方法的特点是调查者在某一特定的人群相对集中到呃公共场所现场拦截被调查者进行的访谈。 邮寄问卷调查: 邮寄问卷调查是调查者将设计好的问卷通过邮寄的方式送达被调查者手中，请他们按要求和规定时间填写问卷并寄回调查者，以此来获取信息的一种方法。 留置问卷调查: 留置问卷调查是指调研人员将调查问卷送到被调查者的手中，征得同意后对填写事项做出说明并与被调查者约定交返问卷的时间，调研人员姐约定时间再次登门收取问卷，并向被调查者致谢的整个收集信息的过程。 电话调查: 电话调查抽样的基本原理跟其他调查方法一样，但在操作上有其独立性。具体包括以下是哪个步骤：1.抽取样本户2.选择受访者3.选择替代样本 网络调查: 网络调查，也叫网上调查，是指在互联网上针对调查问题进行调查设计，收集资料及分析咨询等活动。 深度访问法: 深度访问法类似于记者采访，是一种无结构访问，指事先不拟定问卷，访问提纲或访问的标准程序，由访者与受访者就某些问题自由交谈，从交谈中获取信息，用以揭示对某一问题的潜在动机，态度和情感的资料采集方法，它最适合用于探索性调查。 座谈法: 座谈法也叫重点小组或焦点访谈法，就是采用小型座谈会的形式，挑选一组具有代表性的消费者或客户，在一个装有单面镜或录音录像设备的房间，在主持人的组织下，某个专题进行讨论，从而获得对有关问题的深入了解。 随着科技的发展技术的进步，市场调查中的数据采集方法也与以往有所不同，多样化的采集方法、遍布全国的采集网络和抽样系统、高标准的质量控制体系，已成为市场调查公司为客户提供有价值，准确的，及时的商业消费者信息和渠道信息。
第五章 结构动力学中的常用数值方法 5．1．结构动力响应的数值算法 ... . 0()(0)(0)M x c x kx F t x a x v ? ++=??=??=?? 当c 为比例阻尼、线性问题→模态叠加最常用。但当C 无法解耦，有非线性存在，有 冲击作用（激起高阶模态，此时模态叠加法中的高阶模态不可以忽略）。此时就要借助数值积分方法，在结构动力学问题中，有一类方法称为直接积分方法最为常用。所识直接是为模态叠加法相对照来说，模态叠加法在求解之前，需要对原方程进行解耦处理，而本节的方法不用作解耦的处理，直接求解。（由以力学，工程中的力学问题为主要研究对象的学者发展出来的） 中心差分法的解题步骤 1. 初始值计算 （1） 形成刚度矩阵K ，质量矩阵M 和阻尼矩阵C 。 （2） 定初始值0x ，. 0x ，.. 0x 。 （3） 选择时间步长t ?，使它满足cr t t ?