Search
Code Directory
 ASP
 ASP.NET
 C/C++
 CFML
 CGI/PERL
 Delphi
 Development
 Flash
 HTML
 Java
 JavaScript
 Pascal
 PHP
 Python
 SQL
 Tools
 Visual Basic & VB.NET
 XML
New Code
HTMLPad 2018 15.5
WeBuilder 2018 15.5
SentiMask SDK Trial 2.0.0
Track Order For Magento 2 1.0.0
Calendar 365 For Dynamics CRM 4.0
Scimbo 1.2
Odoo Furnito Theme 1.0
.Net Assemblies/WSDL Import Utility for Delphi 6.0.4.0
Devart Excel Add-in for Zendesk 1.7
Excel Add-in for Marketo 1.7
GetOrgChart 2.5.3
AnyGantt JS Gantt Charts 8.6.0
AnyChart JS Charts and Dashboards 8.6.0
Multi Vendor Ecommerce Script & Apps 4.0
Magento 2 Inventory Management Extension 1.0
Top Code
HTMLPad 2018 15.5
MCU 8051 IDE 1.4.7
jQuery Slider Ken Burns 1.0
Extreme Injector 3.7
Azizi search engine script PHP 4.1.10
BedRock IDE 1.0
Email Grabber
Guest Book V. 2 V2
XML Products Showcase Slider 1.0
Free Pageant Contest Source Code 1.0
ActivePerl 5.8.4.810/5.6.1.638
Calendar 365 For Dynamics CRM 4.0
Gz Hotel Booking Script 2.0
Image Reconstruction Demo (GUI) 1.0
CardCheck COM DLL for Credit Card Processing 1.1
Top Rated
VisualNEO Web 2018.12.15
Paste phpSoftPro 1.4.1
Extreme Injector 3.7
Deals and Discounts Website Script 1.0.2
ADO.NET Provider for ExactTarget 1.0
Solid File System OS edition 5.1
Classified Ad Lister 1.0
Aglowsoft SQL Query Tools 8.2
Invoice Manager by PHPJabbers 3.0
ICPennyBid Penny Auction Script 4.0
PHP Review Script 1.0
ATN Resume Finder 2.0
ATN Site Builder 3.0
Availability Booking Calendar PHP 1.0
PHP GZ Blog Script 1.1
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
Download GCD of Polynomials and Polynomials raised to some Power including Fractional Power 1.0http://www.mathworks.comReport Error Link
License: Shareware
File Size: 92.2 KB
Downloads: 6
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 + (-75-180j)*x^4 + (920+180j)*x^3 + (-1785+1800j)*x^2 + (-732-3582j )*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: 6



More Similar Code

chatFX combines the power of Flash and XML to produce an intuitive, multi-platform, 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, XSL-FO, 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
- required fields
     

Please enter text on the image