第3章 MATLAB的数值计.pptVIP

第3章 MATLAB的数值计算 3．1 多项式的创建与运算 3．1．1 多项式的描述与创建 1．多项式的描述 在MATLAB环境下多项式的表达式 可以转换成向量的形式加以描述，即 向量最右边的元素表示多项式的0阶系数，向左数依次表示多项式的1阶系数，2阶系数，3阶系数…。也就是说，MATLAB按降幂形式排列多项式的系数，例如多项式表达式 ，在MATLAB环境下用向量描述为：[6 3 4 2 7]。 2．多项式的创建 创建多项式的方法有许多。直接输入法；特征多项式输入法 ；由根向量创建多项式 【例3-1】在MATLAB环境下采用直接输入法，创建多项式 程序设计为： p=[1 3 -2 5]; %多项式的系数向量 poly2sym(p) %生成多项式的表达式 运行结果： ans = x^3+3*x^2-2*x+5 3．1．2 多项式的运算 1．求多项式的值 求多项式的值的方法有二种，一种按数组运算规则计算，由函数polyval实现；另一种按矩阵的运算规则计算，由函数polyvalm来实现。 【例3-4】求多项式 在5、6和8处的值。 程序设计为： p=[4 -2 10]; x=[5 6 8]; polyval(p,x) 运行结果： ans = 100 142 250 2．求多项式的根 求多项式的根即多项式为零时的值，在MATLAB环境下求根时有两种方法：一种是直接调用根函数roots来求解多项式的根的方法；另外还有一种方法就是先求多项式的伴随矩阵，然后再求特征值的方法来求多项式的根。这里仅介绍前面一种方法。 【例3-6】求多项式 的根。 程序设计为： a=[1 5 -6 10 3 20]; r=roots(a) 运行结果： r = -6.2232 1.1158 + 1.1853i 1.1158 - 1.1853i -0.5042 + 0.9791i -0.5042 - 0.9791i 3．多项式的加减运算 多项式的加减运算是多项式对应元素的加、减运算，多项式的阶数可以不同，但在进行多项式定义时，应当在低阶多项式的前面补充0，使其阶数相等，否则，不能进行加减运算。 【例3-7】对两个多项式 和 进行加、减运算 程序设计为： a=[1 3 0 -2 3]; b=[0 0 1 3 4]; c=a+b 运行结果： c= 1 3 1 1 7 poly2sym(c) ans = x^4+3*x^3+x^2+x+7 4．多项式的乘、除运算 多项式的乘法运算由函数conv实现；多项式的除法运算由函数deconv实现。 【例3-8】计算两个多项式 和 的乘积。 程序设计为： a=[1 5 -2 1]; b=[1 5]; c=conv(a,b) 运行结果： c = 1 10 23 -9 5 poly2sym(c) ans = x^4+10*x^3+23*x^2-9*x+5 5．多项式的微积分 多项式的微分由函数polyder实现，多项式的积分由函数polyint实现。 【例3-10】求多项式 的一阶导数。 程序设计为： p=[3 -2 5 -10 6]; polyder(p) 运行结果： ans = 12 -6 10 -10 poly2sym(ans)