Transfer function to differential equation

barneygumble742
Messages
27
Reaction score
0

Homework Statement


Write the differential equation that is equivalent to the transfer function given below. Plot y(t). Assume that r(t) = 4t^{2}
Y(s) = 2s^{4}+3s^{3}+2s^{2}+s+1
R(s) = 2s^{5}+3s^{4}+2s^{3}+2s^{2}+4s+2


The transfer function is Y(s)/R(s).

Homework Equations



I'm a little lost on how to get started with this problem. Could anyone please help?

The Attempt at a Solution



Given r(t), I thought of converting it to LaPlace and then multiplying it with the numberator so I would be left with Y(s) = numerator / denominator. After that I'll have a mess that I don't think will factor without imaginary numbers. I'm thinking of using partial fraction expansion.
OR
I could have it in this form:
Y(s) [2s^{5}+3s^{4}+2s^{3}+2s^{2}+4s+2] = R(s) [2s^{4}+3s^{3}+2s^{2}+s+1]
and then convert each item to the time domain and then put it back in the transfer function form. However if I did this, then what about the final r(t) = 4t^{2} that's left over?

Thanks,
BG742
 
Physics news on Phys.org
I don't quite understand the point of making up a differential equation for which you have the solution (or really that part at all, is there more to the problem?), but as far as I can tell you have a transfer function in the Fourier domain, and you need to find the differential equation that it represents in the time domain. There is some response to a certain function given by the transfer function Y(s) = R(s)T(s), where T(s) is the transfer function. You know that y(t) = F(R(s)*T(s)), and you know that T(s) = Y(s)/R(s) with r(t) = 4t^2. So with all this you can figure out what y(t) is.

Whether you want to use Laplace Transforms or Fourier Transforms depends on the context of the problem. If this is a question about filters, which it seems like it is, then you should use Laplace. Actually dividing the two equations out will probably make the Laplace transform a lot easier.
 
The problem asks me to convert the transfer function to the time domain and plot it in terms of y(t). We have to use Laplace Transforms because of the nature of the course. Could you please explain what you meant by dividing the two equations out?
 
Actually go through the division of \frac{2s^4+3s^3++2s^2+s+1}{2s^5+3s^4+2s^3+2s^2+4s+2}. Matlab or Mathematica can do the job quick, but by hand won't take long either.

You might be over-analyzing the whole thing. The problems intent seems like it is to have you become more familiar with transfer functions and how the convolution theorem works.
 
Last edited:
I think it's worth doing because if you ever have to implement a transfer function in anything less transparent than MATLAB you'll need to know.

-factor the numerator and denominator, make sure there isn't anything you can cancel between the two (there's not in this case).

-Convert any complex conjugate pairs back into 2nd order polynomials.

x=r1=(a+bi)
x=r2=(a-bi)

(x-r1)(x-r2) = x2-(r1+r2)x+r1*r2 = x2+2ax+(a2+b2)

-Now we don't have to deal with imaginary numbers.

-Think of the system not as one transfer function but a series of small ones.

-you can factor any polynomial into a series multiplication of first and second order polynomials.

-A first or second order denominator in a transfer function is a low pass filter or an under damped oscillator respectively.

-A multiplication of transfer functions, in the time domain means the output from one, is the input for the next.

- it doesn't matter what order you multiply things in-> it doesn't matter what order you chain these sub-systems together in if all you're interested in is the original input and output variables, but it will make a difference on the hidden variables.

-it's not difficult to turn a low pass filter back into it's differential equation. 1/(Ts+1) => X'=(X-Xi)/T

X is the output,
Xi is the input,
T is the time constant

-an oscilator isn't dificult either: 1/(ms2+bs+k) => V'= F-bV-KX, X'=V

F is the input or 'external force'
X is the output 'mass location'
V is the velocity of the mass
b is the 'damping coeficient'
K is the 'spring constant'.

(of course the physical analogy depends on your system. it's also possible to have negative spring constants and damping factors here).

- multiplying by 's' is taking the derivative. multiplying by (3s+2) takes x and outputs (3x'+2x). To include the numerator all you have to do is take all the derivatives that will be necessary (up to the maximum power of 's'), weight them (with the coeficients in the polynomial), and add them back together. I prefer to put the numerator at the end, so that you never have to take the derivative (or in this case the 4th derivative) of a step input.
 
suki said:
-you can factor any polynomial into a series multiplication of first and second order polynomials.

-A first or second order denominator in a transfer function is a low pass filter or an under damped oscillator respectively.


actually now that I think a little more : you don't need to factor the denominator. You can get a differential equation directly from it using the same pattern as for the second order system. the max power of s in the denominator, put that many integrators in series, after each integrator put a negative feedback link, with a constant coefficient, to before the first integrator except for the coefficient on the highest power, the feedback coefficients go directly to the coefficients of the polynomial. and it'll be obvious:

Xin=2X(5)+3X(4)+2X(3)+2X''+4X'+2X

I just prefer to have it in terms of first and order systems, which I understand.
 
There are two things I don't understand about this problem. First, when finding the nth root of a number, there should in theory be n solutions. However, the formula produces n+1 roots. Here is how. The first root is simply ##\left(r\right)^{\left(\frac{1}{n}\right)}##. Then you multiply this first root by n additional expressions given by the formula, as you go through k=0,1,...n-1. So you end up with n+1 roots, which cannot be correct. Let me illustrate what I mean. For this...
Back
Top