Scomponendo il numeratore e polinomi denominatore nelle loro parti pari e dispari

Question

Qui v'è una continua funzione di trasferimento tempo (G(s)) in forma di:

G(s) = N(s)/D(s); 
G(s) = (s^3+4s^2-s+1)/(s^5+2s^4+32s^3+14s^2-4s+50)  (1) 

e (s = j*w) dove w = frequency symbol.

Ora, come è possibile per decomporre il numeratore e il denominatore polinomi di Eq. (1) nelle loro parti pari e dispari e di ottenere il G(jw) come (Utilizzo di Matlab):

Answer 1

Si potrebbe probabilmente prendere le parti reale e immaginaria, dopo la sostituzione con s=j*w. Tuttavia, si può effettivamente selezionare le parti pari e dispari dei vostri polinomi:

% G(s) = N(s)/D(s); 

syms s; 
N = s^3+4*s^2-s+1; 

p = sym2poly(N); 

%// do this in fewer lines: 
%{ 
/* 
if mod(length(p),2)==0 %// then first index is odd 
    imin_o = 1; %// for odd part 
    imin_e = 2; %// for even part 
else 
    imin_o = 2; %// for odd part 
    imin_e = 1; %// for even part 
end 
*/ 
%} 
imin_o = mod(length(p),2) + 1; 
imin_e = 2 - mod(length(p),2); 

% odd part of numerator 
p_o = zeros(size(p)); 
p_o(imin_o:2:end) = p(imin_o:2:end); 
% even part of numerator 
p_e = zeros(size(p)); 
p_e(imin_e:2:end) = p(imin_e:2:end); 

% restore 
N_o = poly2sym(p_o,s); 
N_e = poly2sym(p_e,s); 

e lo stesso per il denominatore.