Exact value of a result
14 views (last 30 days)
Show older comments
Hello, I am trying to get the exact values of the following lines :
b=[-1, 3];
a=[1, -3, -2];
[r,p,k]=residuez(b,a);
right now, I obtain something like : p=[3.5615;-0.5615] and r=[-0.1361;-0.8638] and k=1.
Is there a format that would allow me to get me to see it displayed like : p=[(3-sqrt(17))/2 ; (3+sqrt(17))/2] and same for r?
thanks
0 Comments
Answers (2)
David Hill
on 29 Nov 2022
residuez requires doubles as input, but p is just the roots of polynomial a which you could use symbolic to produce.
a=sym([1, -3, -2]);
p=roots(a)
0 Comments
Walter Roberson
on 29 Nov 2022
No, the calculations are inherently carried out in floating point.
If you were to follow along in the source code you could strip out the checks for floating point, and you could convert the numeric roots() call to use poly2sym() and then solve for symbolic roots. But after that the code uses filter() to run an impulse through, and emulating that symbolically might take a bit of work.
0 Comments
See Also
Categories
Find more on Numbers and Precision in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!