Transfer function with non-zero initial value

Dammi87

New member
Joined
Sep 27, 2012
Messages
3
Hi there

This is my first time trying this webpage :)

I have this equation: Eq1.gif , with x(0) = E
Here x(t) is the output, D is the input

By taking Laplace transform, first with D = 0 and x(0) = E, then D = D and x(0) = 0 I get these two equations:

Eq5.gif and EQE2.gif So the final answer for X(S) would be the sum of these two or EQe3.gif

Now I know that x(t) = G(s)u(t) .. how can I find G(s) ? (I'm assuming I haven't done anything wrong up till now, unless you can find a error somewhere?)

This is a part of an assignment of Parametric Signal Identification which I've done, but I can't verify my solution unless I can derive the transfer function (Or state space model) to put into SimuLink.
I at least know that the input is a step and the output looks like the output for a first-order system.

You would be of great help if you can help me with this, it's really embarrassing to go through really hard PSI methods, and then get stuck on the small matter of a simple Transfer Function!
 

Attachments

  • Eq1.gif
    Eq1.gif
    953 bytes · Views: 9
  • Eq2.gif
    Eq2.gif
    983 bytes · Views: 9
  • Eq3.gif
    Eq3.gif
    851 bytes · Views: 9
  • Eq4.gif
    Eq4.gif
    1,001 bytes · Views: 9
  • Eq2.gif
    Eq2.gif
    1.1 KB · Views: 9
Just to give you guys some more insight into what I'm doing.

The differential equation actually looks like this: EQE4.gif , where p,c,C,W,T and Q are constants. What we are doing is heating tissue with the power Q and describing the heat in that tissue at time t --> x(t)
Now, the initial temperature of the tissue is of course 37°C.

What we are supposed to do is to estimate the parameter W according to a given dataset. I've done that and the DISCRETE results I get are the dots, on this figure. The actual data is the solid line

untitled.jpg

And this is an okay estimation for W I think. I use Backward Euler on x'(t) and convert the equation to discrete time and some more magic to be able to iterate my calculations until I get a good enough estimation of W.
Now to test this I wanted to make a simulink model with all these parameters and see if I'd get similar results.

Hopefully you have some more magic up your sleeve :)
 
Top