
GCD of Polynomials and Polynomials raised to some Power including Fractional Power 1.0 File ID: 82530 


 GCD of Polynomials and Polynomials raised to some Power including Fractional Power 1.0 License: Shareware File Size: 92.2 KB Downloads: 8
Submit Rating: 



GCD of Polynomials and Polynomials raised to some Power including Fractional Power 1.0 Description 

Description: 15th July, 2005 : Poly_POWER.m is now corrected !
So, for most reasonable cases of Multiple Roots including Multiple Real Roots, this Programme should now work.
For eg, Poly_POWER works successfully for : { x^6 + (12+18j)*x^5 + (75180j)*x^4 + (920+180j)*x^3 + (1785+1800j)*x^2 + (7323582j )*x + (2035+828j) } ^ 0.5 The answer is "approx": { ( x^3 + (6 + 9j) * x^2 + (15  36j) * x + (46 + 9j) } and for Liouville's Constant based Polys, like : ( { x^6  75*x^3  190*x + 21 } ^ 3 } ^ 0.3333
********************
Functional Description of Poly_GCD :  If we need to verify the fact that a Polynomial has multiple roots iff (if and only if) it has a common factor with it's derivative, we need two things :
a) A function to compute the GCD of 2 Polynomials. Since I could not find a Standard Matlab function for this, I created this function : Poly_GCD.m : GCD = Poly_GCD ( sx_poly, rx_poly ) For eg, if : sx_poly = x^9  3x^8 + 0x^7 + 2x^6 + 6x^5 + 0x^4  4x^3  6x^2  3x ? 1 and rx_poly = 9x^8  24x^7 + 0x^6 + 12x^5 + 30x^4 + 0x^3  12x^2  12x ? 3, then their GCD = x^6  2x^5  x^4 + 3x^2 + 2x + 1
b) A function which can raise a Polynomial to some power in order to simulate an overall polynomial with multiple roots. Since I could not find a Standard Matlab function for this, I created Poly_POWER.m : P = Poly_POWER ( poly, n )
In addition to computing the GCD of the 2 input Polynomials : sx_Poly and rx_Poly, Poly_GCD.m and Poly_GCD_Main.m also find suitable polynomials : Nx & nx and cx & dx such that : check_GCD_Orig = Nx * sx_Poly + nx * rx_Poly check_GCD_Used = cx * sx_Poly_Used + dx * rx_Poly_Used (Internally used vars)
While finding HCF of 2 Polynomials, as the degrees of the polynomials increase, the accumulated FP errors increase, and can only be "salvaged" to some extent. The trick lies in being able to devise a way to detect the correct Zero limit on the Remainder, and thereby, stop after the correct number of HCF divisions. This requires lots of experiments to find what limits to apply for the definition of 0 (zero) ; this has been done with some complex empirical logics obtained by trial and error.
Convergence also depends upon whether we convert the input Polys to monic, or the intermediate computed Poly also to monic, or we do not convert at all. These 3 combinations create a whole lot of extra logics, almost 30 % of the development effort ! We choose the best combination. When compared to my very early submission of this zip file, the earlier Poly_GCD.m is now renamed as Poly_GCD_Main.m, and Poly_GCD.m is now a "top level" function.
I think that all these functions : Poly_GCD, Poly_POWER, CMPLX_GCD, Ch_Rem_Thr_Poly.m, Ch_Rem_Thr_Int.m, Gen_Primes_Eq_2_Sqs etc are essential as commonly required functions, and can be placed in the dir : ...matlabspecfun
Should generally work for R14, R13 and R12.
License: Shareware Related: apply, definition, Complex, Limits, experiments, divisions, requires, empirical, logics, depends O/S:BSD, Linux, Solaris, Mac OS X File Size: 92.2 KB Downloads: 8


More Similar Code 

chatFX combines the power of Flash and XML to produce an intuitive, multiplatform, highly customizable chat client anyone can install and use.
XSQL combines the power of XML and SQL to provide a language and database independent means to store and retrieve SQL queries and their results.
Sometimes you get a list of lists and want to swap rows and columns, i.e. transpose the list. Yet, what if the rows have different lengths? Here is some advice you might find useful in such situations.
BTL is a template language that combines power of JSTL and XSLT to produce documents in XML, HTML, XHTML, XSLFO, PDF or other formats, based on the JavaBean input.
The purpose of project is libarary that is foundation for importers from Rational Rose file formats to other formats and importers to some specific formats like nsuml xmi.
It is application of of Simulink Block to Cardiac PDE VI1 system of 2 non linear coupled equations of PDE and try to design solution with dicretized space and time and issue in 1D.I associate a Doc file to introduction theory.Just see this page...
These algorithms are an implementation of the iterative amplitude adapted Fourier Transform method of Schreiber and Schmitz to make surrogate time series. These surrogate time series or fields were made to have the Power spectrum and the amplitude...
Simios is a collection of libraries and tools to ease development on any CLI framwork, focusing in mono and MS.NET. All libraries are 100% managed and are beeing developed by the Simios.org community with some support by other Mono communities as...
SVGMap is an open source web tool to help you in the visualisation of experimental data associated to some graphical representation.
More at :
1. Project website http://bg.upf.edu/svgmap
2. Reference :
DSPack is a set of Components and class to write Multimedia Applications using MS Direct Show and DirectX technologies. DSPack 2.2 is designed to work with DirectX 9 on Win9X, ME, 2000, and Windows XP operating systems. DSPack 2.2 is designed to... 
User Review for GCD of Polynomials and Polynomials raised to some Power including Fractional Power 
