Extract num and den in symbolic expression.

The 'tf' function in Control Systems needs numerator and denometer in a vector of decreasing powers of s.

For example a numerator of a*s^2 + b*s + c needs to be num= [a b c];. I have been cutting and pasting the necessary coefficients

but to do this I have to pause my program. Is there a way to extract the appropriate vector(s) programatically?

Thanks,

Bruce Taylor

Answer 1

[Star Strider]

Use the sym2poly funciton:

syms a b c s 
num = 1*s^2 + 2*s + 3;
n = sym2poly(num)

n = 
     1     2     3

[Bruce Taylor]

This answers only part of the questuon.

given

sys=(a*s^2+b*s+c)/(d*s^3+e*s^2+f*s+g)

What I need is

num=[a b c];
den=[d e f g];

then I can do

sys=tf(num,den);

[Star Strider]

I thought you had already isolated the numerator and denominator.

Try this:

syms a b c d e f g s 
sys=(a*s^2+b*s+c)/(d*s^3+e*s^2+f*s+g)
sys = subs(sys,{a,b,c,d,e,f,g},{3,6,8,1,2,9,5})
[np,dp] = numden(sys)
n = sym2poly(np)
d = sym2poly(dp)

producing:

sys =
(a*s^2 + b*s + c)/(d*s^3 + e*s^2 + f*s + g)
 
sys =
(3*s^2 + 6*s + 8)/(s^3 + 2*s^2 + 9*s + 5)
 
np =
3*s^2 + 6*s + 8
 
dp =
s^3 + 2*s^2 + 9*s + 5
 
n =
     3     6     8

d =
     1     2     9     5

Note that this will only work with numeric coefficients. It will error with symbolic coefficients.

[Walter Roberson]

Star Strider is correct, MATLAB transfer functions, tf() and ss(), cannot support symbolic expressions.

Answer 2

[Walter Roberson]

https://gist.github.com/maruta/1035254