科学与工程计算：引言预案.ppt

科学与工程计算;需要的基础知识;参考书;1.注意掌握各种方法的基本原理 2.注意各种算法（差分、有限元）的构造手法 3.重视各种方法的性态研究（收敛性、稳定性） 4.做一定量的习题和上机实践题 5.注意与实际问题相联系;考核; 随着计算机软件硬件的不断更新和计算方法的迅速发展，科学计算与实验以及理论研究成为现代科学研究的三大主要手段。科学计算还能解决实验及理论无法解决的问题，并由此发现一些新的物理现象，加深人们对物理机理的理解和认识，促进科学的发展。;; 偏微分方程数值解 ;偏微分方程数值解研究的对象;偏微分方程数值解的内容;第一章 引论、准备知识;§1 微分方程;一、偏微分方程简介 ;2、发展史 (1) 十八世纪 (2) 十九世纪中期，三类偏微分方程： 3）十九世纪末到二十世纪初，其它方程： 高阶方程： KDV方程： 薛定谔方程（Schr?dinger equation）:;偏微分方程分为线性和非线性，;3、建立方程的步骤 当确定了研究的是哪一类物理量后， (1)、划分出一小块，考虑其与邻近部分的关系； (2)、根据物理学规律，表示出此关系； 如：牛顿运动定律、能量守恒定律、麦克斯韦方程等。 (3)、化简、整理，即得偏微分方程。;这里主要讨论三类方程，弦振动方程、泊松方程、热传导方程。此三类方程并不包括所有物理问题，如：量子力学中的薛定谔方程，KDV方程等。;例. 一维波动方程的建立;T2;这个方程不仅仅可以描述上述物理过程，也可以用来描述管道气体小扰动的传播、电报方程、均匀杆的纵振动、均匀薄膜的微小横振动（二维）、电磁方程（三维）等。;薄膜振动问题;热传导问题;;?;热传导方程; 满足方程的函数u( x, t)称为方程的解，对于微分方程，仅仅满足方程还不够。还需要初始条件、边界条件。;定解问题=泛定方程+定解条件;典型方程;1、温度------- ?热量传输 2、污染物密度---------?环境工程 3、密度分布--------?统计力学 4、价格---------- ?金融工程 5、………..;三、 数值求解微分方程的意义;为什么要研究数值求解方法;2）在很多情况下，寻找解析解不可能。现实问题中归结的微分方程不满足解析解的存在条件的比比皆是，方程中出现的有些函数连续性都无法保证，它们并不存在前述意义的解析解。于是，求数值解便成了在这种情况下解决问题的重要手段了。;四、数值求解方法概述;3）初始和边界条件处理 离散后系统若是一个递推公式，那它需要初值才能启动。若是一个方程组，那它所含的方程个数一般少于未知量的个数，要想求解还需要补充若干个方程。这些需要补充的初值和方程往往可以通过微分方程的初始条件和边界条件来得到，这就是初始和边界条件处理过程。 ;4）离散系统的性态研究 我们主要研究：这个系统是否可解，即解的存在性、唯一性问题；它与精确解的差距有多大，这个差距当区域剖分的尺寸趋于零时，是否也会趋于零，趋于零的速度多快，即解的收敛性和收敛速度问题；当外界对数据有所干扰时，所得的解是否会严重背离离散系统的固有的解，即解的稳定性问题。;1、模型误差 2、观测误差 3、截断误差 4、舍入误差;1）模型误差 在将实际问题归结为数学模型时需要对问题作一定的简化和假设，由此产生的误差叫模型误差。;3）截断误差 将数学模型离散化时由于舍弃一些次要的项而导致的模型问题真解与离散问题真解的误差称为截断误差。;* 主要研究截断误差。一般来说，一个可用方法的截断误差不应该超过模型误差和观测误差;数值求解微分方程过程示意;求解过程中产生的误差