[经济学]第5章 数值计算基础.pptVIP

第5章 数值计算基础 本章目标 掌握多项式的构造和运算方法 掌握解线性方程的方法 能够使用常用的几种数值分析方法进行一般的数值问题求解 主要内容 5.1 多项式 5.2 线性代数 5.3 数值分析 5.4 函数极值和零点 5.5 插值和拟合 5.1 多项式 5.1.1 创建多项式(P47) 对多项式 P(x) = a0xn+a1xn-1+…+an-1x+an 在MATLAB中，多项式用行向量表示，向量中的元素按多项式降幂排列： P=[ a0 a1 … an-1 an] 例: X^4 - 34 X^3 -80X^2 +1可以表示为向量 [1 -34 -80 0 1] 多项式的创建方法： 1、系数矢量的直接输入法 在命令窗口中直接输入多项式的系数矢量，利用poly2sym将多项式由系数矢量形式转变为符号形式。 poly2sym([3 -6 7 0 -9]) ans = 3*x^4-6*x^3+7*x^2-9 例：p1=[3 2 0 4 5] p2=[6 3 -4 6] p3=[2 4 5 0 7 0 8] 2、特征多项式输入法 由矩阵的特征多项式系数建立多项式，由函数poly实现。 n阶方阵的特征多项式系数的矢量是n+1阶的，而且系数的第一个元素必须是1。 b=[3 5 6;-2 1 3;0 3 -2]; p=poly(b) p = 1.0000 -2.0000 -4.0000 89.0000 poly2sym(p) ans = x^3-2*x^2-4*x+89 5.1.2 多项式运算 MATLAB中在对多项式进行加减运算时,参加运算的多项式应该具有相同的阶次,如果阶次不同,低阶的多项式必须用零填补至高阶多项式的阶次。 1、多项式的加减法 a=[1 3 5 7]; b=[2 4 6 8]; c=a+b c = 3 7 11 15 d=[1 3 4 6 -5]; e=[0 c]+d e = 1 6 11 17 10 e=c+d ??? Error using == plus Matrix dimensions must agree. 2、求多项式的值 求多项式的值有两种算法 按数组运算规则计算，对应函数polyval 按矩阵运算规则计算，对应函数polyvalm 函数polyval的调用格式为： y=polyval(p,x)：求多项式p在x点的值，x也可以是一数组，表示求多项式p在各点的值。 p=[3 -6 7 0 -9]; polyval(p,[3 5 7]) ans = 135 1291 5479 函数polyvalm的调用格式 y=polyvalm(p,x)：求多项式p对于矩阵x的值，要求矩阵x必须是方阵，x如果是一标量，求得的值与函数polyval相同。 例1：求多项式4x2+2x-1对于向量[2 3 1 4] 的值。 例2：求多项式2x2+3x+1对于向量[2 1 0] 的值。 p=[4 2 -1]; polyval(p ,[2 3 1 4]) ans = 19 41 5 71 p=[2 3 1]; polyval(p ,[2 1 0]) ans = 15 6 1 3、求多项式的根 两种方法： 直接调用求根函数roots 先把多项式转化为伴随矩阵，然后求其特征值 p=[3 -6 7 0 -9]; r=roots(p) r = 0.6975 + 1.4641i 0.6975 - 1.4641i 1.4126 -0.8075 s=compan(p) s = 2.0000 -2.3333 0 3.0000 1.0000 0 0 0 0 1.0000 0 0 0 0 1.0000 0 r=eig(s) r = 0.6975 + 1.4641i 0.6975 - 1.4641i 1.4126 -0.8075 4、多项式的乘除运算 多项式的乘法由函数conv实现，除法由函数deconv实现。 a=[3 -4 6 0 2]; b=[1 2 -5 3]; c=conv(a,b) c =