Approximating Differential Equation on the Basis of Measured Values

Newsbeitragnews: snipped-for-privacy@i12g2000prf.googlegroups.com...

ftp://ftp.deltacompsys.com/public/NG/Mathcad%20-%20Type%201%20First%20Order%20System.pdf

Your data

K = 3 alpha = 2

K*alpha 1

----------- = ------------------------ s*(s+alpha) 0.16666s^2+0.33333s + 0)

I approximated in Page 2

formatting link
^^^^^^^^

0.1666572 y'' + 0.3333362 y' + 2.174717E-06 y = 1 0.1666572 y'' + 0.3333362 y' + 0 y = 1

For 2.174717E-06 ~ 0

The velocity is the green line (y') on the same page.

Any problems left?

Reply to
JCH
Loading thread data ...

Newsbeitragnews: snipped-for-privacy@i29g2000prf.googlegroups.com...

That looks good. The other answers don't.

Peter Nachtwey

Reply to
pnachtwey

Newsbeitragnews: snipped-for-privacy@i29g2000prf.googlegroups.com...

Velocity must be ok because y' is part of

0.1666572 y'' + 0.3333362 y' + 0 y = 1

formatting link

Reply to
JCH

Newsbeitragnews: snipped-for-privacy@i12g2000prf.googlegroups.com...

Newsbeitragnews: snipped-for-privacy@i29g2000prf.googlegroups.com...

Yes, the velocity is OK. However, your web site is still wrong. The steady state position for the equation on your web site is

1/2,174717E-06 or 459830 which is not right. You must integrate velocity to get position so you can control position. The PID gains must be calculated for position feedback. The PID gains for velocity control and the PID gains for position control are different even if the gain and time constant are the same. For starters a first order velocity control system doesn't need a derivative gain.

Peter Nachtwey

Reply to
pnachtwey

Newsbeitragnews: snipped-for-privacy@i12g2000prf.googlegroups.com...

Newsbeitragnews: snipped-for-privacy@i29g2000prf.googlegroups.com...

1 T(s) = ----------------------------- 0.1666572s^2 + 0.3333362s + 0

The level goes theoretically to infinity as any 'level' process transfer function does if e.g. more water [m^3/h] is pumped in then goes out of a tank.

In simulations you can ignore an error of 2.174717E-06

Reply to
JCH

That is your hint that you need to use two differential equations like I showed in my pdf file. In this case the control output determines the steady state velocity. Then the velocity must be integrated to get position.

No, you just must use two differential equations. Not one.

See this ftp://ftp.deltacompsys.com/public/NG/Mathcad%20-%20Ode%20Opt%20First%20Order%20no%20offset.pdf

Peter Nachtwey

Reply to
pnachtwey

PolyTech Forum website is not affiliated with any of the manufacturers or service providers discussed here. All logos and trade names are the property of their respective owners.