Matlab：一维数组和其应用.ppt

一维数组相关函数的应用以线性拟合为例 四、一维数组在一元多项式运算中的应用 2、多项式的符号表示 练习 练习 11、部分分式展开（留数计算） 求函数在指定点处的留数 12、多项式的矩阵运算 p是多项式的系数向量，相当于用矩阵X代替多项式的变量来对矩阵进行计算，X必须为方阵。 0 -2.5 -4 -5.7 -3.5 -2 -1 2 3.5 4 7 7.5 9.9 10.9 0 0.2 0.6 1 1.3 1.6 1.7 1.8 1.9 2.2 2.3 2.5 2.6 2.9 y x 11.9 13.5 13 11.9 9 6.5 4 1.5 0 -2.5 -5 3.1 3.4 3.8 4.1 4.4 4.7 4.8 4.9 5 5.1 5.3 y x 多项式拟合时的注意事项 在进行曲线拟合时对多项式阶次的选择是任意的。 虽然高阶的多项式可以更准确地拟合数据（仅指该曲线与给定的数据之间的均方误差最小），但在进行曲线拟合时，并不需要采用太高阶的多项式，这主要基于以下原因： 1、越是高阶的多项式其数值特性越差，计算起来也越耗时； 2、随着多项式阶次的升高，拟合的曲线变的越来越不平滑，通常会出现用户不愿意看到的局部波形； 3、由于数据本身的近似性，因此在进行数据拟合时没有必要仅仅考虑使拟合的曲线无限接近数据点，而要在曲线的阶次合均方误差之间综合考虑，因为越是高阶的多项式在物理实现时越困难。 从数学原理上，n+1个数据点可以惟一定义一个n次曲线（或n阶多项式）。 3、向量的方向角与方向余弦，单位向量 a = [4,3,1] d = sqrt(sum(a.^2)) e0 = a./d 4、向量的线性运算 a = [4,3,1] b = [5,2,3] c = 1.6 直接利用matlab的数组运算规则即可进行运算 pab = a + b sab = a – b d = c*a 5、向量的内积（点积，点乘） dot(a,b) = sqrt(dot(a,a)) Matlab中提供了 dot() 函数实现向量的内积 6、两个向量之间的夹角 tmp = dot(a,b)/sqrt(dot(a,a))/sqrt(dot(b,b)) theta = acos(tmp) 7、向量的叉积（叉乘） cross(a,b) = Matlab中提供了 cross() 函数实现向量的内积 8、向量的混合积： dot(cross(a,b),c) 上机练习： 3.以A(1,2,3)，B(2,0,5)，C(4,2,-1)为顶点的三角形的面积。 4.求以A(0,0,2)，B(3,0,5)，C(1,1,0)，D(4,1,2)为顶点的四面体的体积。 　　借助matlab提供的函数，处理多项式是一件非常简单的事情，很容易对多项式进行积分、微分以及求根的操作。 一元多项式在代数中占有非常重要的地位。在实际应用中如对实验数据的插值、微商和曲线拟合等，都要大量用到多项式；在矩阵分析时，也要用到一元多项式的概念。多项式函数是形式最简单的函数，也是最容易计算的函数，从理论上讲，它可以表示绝大多数复杂函数。在许多计算机的计算和编程中，很多函数值如sin(x),cos(x)等的计算都是先将函数进行Tailor展开为多项式进行逼近计算的，并且都能达到很高的精度。 1、多项式的表示和创建 在matlab中，一个一元多项式用一个行向量来表示的，向量元素为多项式系数的降幂排列，其中最后一个元素代表多项式中的0幂项。 p = [1,-12,0,25,116] poly2sym() 把系数向量表示的多项式转为符号表示 poly2sym(p) 用默认字符’x’来表示多项式的变量 poly2sym(p,’v’) 用字符‘v’来表示多项式的便量 p = [1,-12,0,25,116] px = poly2sym(p) pt = poly2sym(p,’t’) 3、计算多项式的值 p = [1,-12,0,25,116]; x1 = 2.3; x2 = [-1,3,7]; y = polyval(p,x) p1 = polyval(p,x1); p2 = polyval(p,x2); 计算多项式p在x处的值，x可以是标量，向量或矩阵 4、求多项式的根 r = roots(p) 此函数返回多项式p的根组成的向量 　　多项式和多项式的根都是用向量表示的，为加以区别，多项式的根表示成列向量。 p = [1,-12,0,25,116] r = root(p) 　　对于系数为实数的多项式，如果其根出现复数，则复数必是成对出现的。 5、由根创建相应的多项式 r = [1;2;-1;3;7] p = poly(r) p = poly(r)　 r是一个向量，构造一个多项式，其根为r