多项式方程求解.ppt

Matlab多项式运算 与代数方程求解器  Matlab多项式运算 口 Matlab中多项式的表示方法 ◆在 Matlab中,n次多项式是用一个长度为n+1的向量来 表示,缺少的幂次项系数为0。 p(x)=anx"+am-x+.+ax+a 在 Matlab中表示为向量:Ian,n1,3…,a1,aol ←>[2,-1,0,3 注:系数中的零不能省! 多项式的符号形式:poly2sym >>po1y2sym([2;-1,0,3])  多项式四则运算 口多项式加减运算 Matlab没有提供专门进行多项式加减运算的函数,事实 上,多项式的加减就是其所对应的系数向量的加减运算 ●对于次数相同的多项式,可以直接对其系数向量进行 加减运算; ●如果两个多项式次数不同,则应该把低次多项式中系 数不足的高次项用0补足,然后进行加减运算。 例:P1=2x3 3 2,-1,0,3] P2=2x+1 L0,O,2,1 P1+P2=2x3-x2+2x+4 「2,-1,2,4  多项式四则运算 口多项式乘法运算:k=conv(P,q 例:计算多项式2x3-x2+3和2x+1的乘积 >>p=[2,-1,0,3]; >>q=[2,1] >> k=conv(p, q) 口多项式除法运算:[k,x]= decon(P,q】 其中k返回的是多项式p除以q的商,x是余式。 [k, r]=decon(p, q)<==> p=conv(q,k)+r  多项式的求导 日po1yder k= polder(p):多项式p的导数; k= polder(P,q):p*g的导数; [k,d]= polder(p,q):p/的导数,k是分子,d是分母 例:已知p(x)=2x3-x2+3,q(x)=2x+1, 求p',(p·q),(p/q) >>k1= polder([2,-1,0,3]) >>k2= polymer([2,-1,0,3],[2,1]) >>[k2,d]=po1yder([2,-1,0,3],[2,1])