2MATLAB2014a多项式计算.pptVIP

应用举例 应用举例 应用举例 应用举例 * 多项式计算 Matlab多项式运算 在 Matlab 中，n 次多项式是用一个长度为 n+1的向量来表示，缺少的幂次项系数为0。例如： 在 Matlab中表示为相应的向量： 例： 注：系数中的零不能省！ 多项式四则运算 多项式加减运算：Matlab没有提供专门进行多项式加减运算的函数，事实上，多项式的加减就是其所对应的系数向量的加减运算。 例： 对于次数相同的多项式，可以直接对其系数向量进行加减运算； 如果两个多项式次数不同，则应该把低次多项式中系数不足的高次项用0补足，然后进行加减运算。 多项式乘法运算： k=conv(p,q) 例1 计算多项式 和 的乘积 p=[2,-1,0,3]; q=[2,1]; k=conv(p,q) 多项式除法运算：[k,r]=deconv(p,q) 其中 k 返回的是多项式 p 除以 q 的商，r 是余式。 [k,r]=deconv(p,q) p=conv(q,k)+r == p=[1 2 0 -5]; q=[1 -1 2]; [k,r]=deconv(p,q) conv(k,q)+r 多项式的导数：polyder k=polyder(p):多项式 p 的导数； k=polyder(p,q): p*q 的导数； [k,d]=polyder(p,q) p/q 的导数，k是分子，d是分母。 k1=polyder([2,-1,0,3]); k2=polyder([2,-1,0,3],[2,1]); [k,d]=polyder([2,-1,0,3],[2,1]) 例3 已知 求 多项式求值 代数多项式求值： y=polyval(p,x):计算多项式 p 在 x 点的值 注：若 x 是向量或矩阵，则采用数组运算（点运算）！ p=[2,-1,0,3]; x=2; polyval(p,x) x=[-1, 2;-2,1]; polyval(p,x) 例4 已知 ，分别取 x=2和一个2?2矩阵， 求 p(x) 在 x 处的值 多项式求根 p=[2,-1,0,3]; x=roots(p) x=roots(p)：若p是n次多项式，则输出x为包含p=0的n个根的n维向量。 例：已知 ，求p(x)的零点。 有理多项式的展开 有理多项式 num=[5,3,-2,7] den=[-4,0,8,3] [r,p,k]=residue(num,den) num = 5 3 -2 7 den = -4 0 8 3 r = -1.4167 -0.6653 1.3320 p = 1.5737 -1.1644 -0.4093 k = -1.2500 有理多项式展开的逆运算 [n,d]=residue(r,p,k) n = -1.2500 -0.7500 0.5000 -1.7500 d = 1.0000 -0.0000 -2.0000 -0.7500 非线性方程的根 Matlab 非线性方程的数值求解 fzero(f,x0)：求方程 f=0 在 x0 附近的根。 方程可能有多个根，但 fzero 只给出距离 x0 最近的一个 fzero 先找出一个包含 x0 的区间，使得 f 在这个区间两个端点上的函数值异号，然后再在这个区间内寻找方程 f=0 的根；如果找不到这样的区间，则返回 NaN。 x0 是一个标量，不能缺省 由于 fzero 是根据函数是否穿越横轴来决定零点，因此它无法确定函数曲线仅触及横轴但不穿越的零点，如 |sin(x)| 的所有零点。 fzero 的另外一种调用方式 fzero(f,[a,b]) 方程在 [a,b] 内可能有多个根，但 fzero 只给出一个 求方程 f=0 在 [a,b] 区间内的根。 f 不是方程！也不能使用符号表达式！ fzero(sin(x)-0.1*x,6) fzero(sin(x)-0.1*x,[2,6]) [x,f,h]=fsolve(‘F’,X0) F为字符串表示的函数； x返回F在x0附近的一个零点，f返回F在x的函数值；h返回值如果大于0，说明计算结果可靠。 非线性方程组求解 非线性方程组求解 例 解方程组 [x,y,h]=fsolve([4*x(1)-x(2)+exp(x(1))/10-1,-x(1)+4*x(2)+x(1).^2/8],[0,0]) 插值法又称