[计算机软件及应用]第四讲 数值计算基础.pptVIP

第三章 数值计算基础  计算机与信息学院 有效使用数字计算机求数学问题近似解的方法与过程，以及由相关理论构成的学科。数值计算主要研究如何利用计算机更好的解决各种数学问题，包括连续系统离散化和离散形方程的求解，并考虑误差、收敛性和稳定性等问题。 从数学类型分，数值运算的研究领域包括数值逼近、数值微分和数值积分、数值代数、最优化方法、常微分方程数值解法、积分方程数值解法、偏微分方程数值解法、计算几何、计算概率统计等。 随着计算机的广泛应用和发展，许多计算领域的问题，如计算物理、计算力学、计算化学、计算经济学等都可归结为数值计算问题。 重要特征 数值计算的结果是离散的，并且一定有误差，这是数值计算方法区别与解析法的主要特征。 注重计算的稳定性。控制误差的增长势头，保证计算过程稳定是数值计算方法的核心任务之一。 注重快捷的计算速度和高计算精度是数值计算的重要特征。 数值计算主要是运用MATLAB这个数学软件来解决实际的问题 数值计算主要是运用有限逼近的的思想来进行误差运算 3.1 多项式 3.1.1 创建多项式 对于多项式, 可以用它的系数矢量表示： 这样，在MATLAB中，将多项式问题转化为矢量问题。 创建多项式的方法有以下几种。 1.???? 系数矢量的直接输入法 在MATLAB的命令窗口直接输入多项式的系数矢量，然后利用转换函数ploy2sym将多项式由系数矢量形式转换为符号形式。 【例3-1】 输入系数矢量，创建多项式。 在命令窗口输入系数矢量，并转换为多项式。 poly2sym([1 -4 3 2]) 2.???? 特征多项式输入法 创建多项式的另一种方法是由矩阵的特征多项式取得，由函数ploy实现。 需要说明的是，n阶方阵的特征多项式系数矢量一定是n+1阶的，同时特征多项式系数矢量的第一个元素必须为1。 【例3-2】 求矩阵 的特征多项式系数，并转换为多项式。 在命令窗口输入矩阵，求特征多项式系数并转换为多项式形式。 a=[1 2 3;4 5 6;7 8 0]; p=poly(a) poly2sym(p) 3.???? 由根矢量创建多项式 由给定的根矢量也可以创建多项式，同样由函数ploy实现。 由给定的根矢量创建多项式应注意两点： （1） 如果希望创建实系数多项式，根矢量的复数根必须共轭成对。 （2） 含复数根的根矢量所创建的多项式系数矢量的系数中，有可能带有很小（在截断误差数量级）的虚部，此时可采用取实部的命令（real）把虚部滤掉。 【例3－3】 根据根矢量[-0.5 –0.3+0.4i –0.3+0.4i]创建多项式。 在命令窗口输入根矢量，并创建多项式： r=[-0.5 -0.3+0.4i -0.3-0.4i]; p=poly(r) pr=real(p) ppr=poly2sym(pr) 3.1.2 多项式运算 1.???? 求多项式的值 求多项式的值有两种算法。一种按数组运算规则计算，对应的函数为ployval；另一种按矩阵的运算规则计算，对应的函数为ployvalm。 函数ployval的调用格式为： ?? y=ployval(p,x) 求多项式p在x点的值，x也可以是数组，表示求多项式p在各点的值。 y=ployvcal(p,x,[],mu) 在这种用法中，用 代 替x，求多项式p在x点的值，其中 mu=[ ], 。 [y,delta]=ployval(p,x,S) 和[y,delta]=ployval(p,x,S,mu) 使用可选的输出结构S产生误差估计 ，其他方法同上。 【例3-4】 求多项式 在5、7和9处的值。 在命令窗口输入多项式系数，并求值； p = [3 2 1]; polyval(p,[5 7 9]) 函数ployvalm的调用格式为： y=ployvalm(p,x) 求多项式p对于矩阵x的值，要求矩阵x必须是方阵，x如果是一标量，求得的值与函数ployval相同。 【例3-5】求多项式 对于矩阵 [2 5 ;7 9]及标量5的值。 在命令窗口输入多项式系数，并求值： p = [3 2 1]; polyvalm(p,[2 5;7 9]) polyvalm(p,5) 2.???? 求多项式的根 MATLAB进行多项式的求根运算时，可以有两种方法，一种是直接调用求根函数roots;另一种是先把多项式转化为伴随矩阵，然后再求其特征值。 在matlab中约定，多项式系数用行矢量表示，一组根用列矢量表示。 【例3-6】 分别用两种方法求多项式