多项式运算及符号方程求解.pptVIP

* 多项式运算及 符号方程求解 * Matlab 多项式运算 Matlab 中多项式的表示方法 在 Matlab 中，n 次多项式是用一个长度为 n+1的向量来表示，缺少的幂次项系数为 0。 在 Matlab中表示为向量： 注：系数中的零不能省！ 例： 多项式的符号形式：poly2str poly2str([2,-1,0,3],x) * 多项式四则运算 多项式加减运算 Matlab 没有提供专门进行多项式加减运算的函数，事实上，多项式的加减就是其所对应的系数向量的加减运算 对于次数相同的多项式，可以直接对其系数向量进行加减运算； 如果两个多项式次数不同，则应该把低次多项式中系数不足的高次项用 0 补足，然后进行加减运算。 例： * 多项式四则运算 多项式乘法运算： k = conv(p,q) 例：计算多项式 和 的乘积 p=[2,-1,0,3]; q=[2,1]; k=conv(p,q); 多项式除法运算： [k,r] = deconv(p,q) 其中 k 返回的是多项式 p 除以 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]); [k2,d]=polyder([2,-1,0,3],[2,1]); 例：已知 ， ， 求 * 多项式的值 计算多项式在给定点的值 代数多项式求值 y = polyval(p,x): 计算多项式 p 在 x 点的值 注：若 x 是向量或矩阵，则采用数组运算 (点运算)！ p=[2,-1,0,3]; x=2; y=polyval(p,x) x=[-1, 2;-2,1]; y=polyval(p,x) 例：已知 ，分别取 x=2 和一个 2?2 矩阵， 求 p(x) 在 x 处的值 * 多项式的零点 x=roots(p)：若 p 是 n 次多项式，则输出是 p=0 的 n 个根组成的 n 维向量。 若已知多项式的全部零点，则可用 poly 函数给出该多项式 p=poly(x) p=[1,-6,11,-6]; x=roots(p) 例：已知 ,求 p(x) 的零点。 * 线性方程组求解 线性方程组求解 X=linsolve(A,b)： 求解线性方程组AX=B，返回特解X 例：解方程组 A=[1 2 -1; 1 0 1; 1 3 0]; b=[2;3;8]; x=linsolve(A,b) * Matlab 符号方程求解器 s=solve(f,v)：求方程关于指定自变量的解； s=solve(f)：求方程关于默认自变量的解。 f 可以是用字符串表示的方程，或符号表达式； 用字符串表示时 ，若f 中不含等号，则表示解方程 f=0。 solve 符号对象的建立：sym 和 syms 例： syms x y z x=sym(‘x’); y=sym(‘y’); z=sym(‘z’); == Matlab 符号方程求解器 * 例：解方程 syms x; f=x^3-6*x^2+11*x-6; s=solve(f,x) s=solve(x^3-6*x^2+11*x-6,x) s=solve(x^3-6*x^2+11*x-6=0,x) * Matlab 符号方程求解器 solve 也可以用来解方程组 solve( f1 , f2 , ... , fN , v1 , v2 , ... , vN) 求解由 f1 , f2 , ... , fN 确定的方程组关于 v1 , v2 , ... , vN 的解 例：解方程组 [x,y,z]=solve(x+2*y-z=27,x+z=3, ... x^2+3*y^2=28,x,y,z) 输出变量的顺序要书写正确！ * k = conv(p,q) [k,r] = deconv(p,q) k = polyder(p) k = polyder(p,q) [k,d] = polyder(p,q) y = polyval(p,x) x = roots(p) 多项式运算小结 多项式运算中， 使用的是多项式 系数向量， 不涉及符号计算！ poly2str(p,’x’) * roots(p)