求多项式除法.doc

求多项式除法 设多项式被除式系数为，除式系数为构造矩阵 由普通带余除法可知在计算机上求多项式除法可以如下算法： 第一行减去第二行（第一行的第一个元素/第二行的第一个元素），得： 将第一行向前移一位最后补0得： 返回（1）执行，一直进行到第一行与对于的元素及以后的元素都为0，即： 则余项系数为，商为n-m次多项式，其系数为第二行所乘的那些数，… 设方程f(x)=0在区间(a，b)内有t(tl)个根，如何将(a，b)分隔为t个小区间，使每个小区间内只有一个根? 取(a，b)的中点c=(a+b)/2，求v(a)，v(c)，v(b)，(v(a)为多项式f(x)在a处的变号数，其余同理)，若v(a)-v(c)=0，则区间(a，c)内无根，取(c，b)作为新的区间，其内根的个数仍为t，将它们保存起来.否则若v(c)-v(b)=0则区间(c，b)内无根，取(a,c)作为新的区间，其内根的个数仍为t，将它们保存起来.否则必有v(a)v(c)v(b)则区间(a，c)，(c，b)内均有实根，记为(a,c)内实根个数，为(c,b)内实根个数，将(a,c)，，(c,b)，保存起来。 若保存下来的所有区间内都只有一个根，则分隔结束，否则从中取出一个其内根的个数大 于1的区间重复上述步骤. 用STURM定理将包含所有实根的区间分隔为若干只含一个根的小区间，并求各小区间 内方程近似根的程序的算法如下: (l)输入方程f(x)=0的次数n，系数，求f‘ (x)的系数， 求f(x)/f‘(x)的余项，将前次相除的除式作下次相除的被除式，将前次相除的余式乘-1保留 下来并作下次相除的除式，一直进行到余式为零即结束，这样求得f(x)的sturm组:. (2)若为零次多项式执行下一步，否则求的商F(x)，再求F(x)的sturm 组. (3)求包含所有实根的区间(a，b)，用上面求出的sturm组求出(a，b)内根的个数v(a)- v(b). (4)若(a，b)内根的个数为零则程序结束，否则将(a，b)分隔为几个小区间，使每一个小区间内 只有一个根. (5)对每一小区间()，取中点。，若，以作为的一个近似根，否则若则根在内以作为新小区间，否则根在内，以作为新的小区间，对新的小区间重复上述步骤，直到。 分别用六个函数来实现上述各部分功能. 用qsturmz(Pl,P2,n)求方程的sturm组. ②用fsturm:(Pl,P4,n)求F(x)=f(x)/fs(x). ③用v(P1,P2,x,n)求方程在x处的变号数. ④用qgdj(P1,P2)求方程根的界. ⑤用dxsv(P1,P2,x)求方程在x处的函数值. ⑥用gdgl(P1,P2,b,n)将(a,b)分隔为几个小区间使每个小区间内只有一个根.