Finding Transfer Function from Step Response

A bit of context:

I'm in a stundent project at my universtiy. We're desingning and building an Mechanical ventilator for Corona-Patients. I'm having some trouble with the control loop of the system.

The ouput of the control loop is controlling a proportional valve. The purpose of this is to reach the desired pressure (Setpoint) in less than 3 seconds. I'm using a PID control, yet it's not tuned correctly. That's why I decided to go to basics, and get the transfer function from a step response, to then simulate the different constants in MATLAB.

The following step response was achieved by opening the proportional valve fully, and waited till the Setpoint was achieved. There's a 1 bar relative pressure in the system. Any pressure above 25mbar flows out, through a PEEP valve. The proplem with the system is that we can't overshot, because we cannot get rid of this pressure as the patient is inhaling.

If someone could help me, calculate the transfer function of the following step response, I would be really gratefull.

Thanks in advance.

Answer 1

[Ameer Hamza]

You can do it using the system identification app: https://www.mathworks.com/help/ident/ref/systemidentification-app.html. For this, you also need the toolbox: https://www.mathworks.com/help/ident/index.html

[Max Mai Tobon]

Thanks a lot. I'm trying to figure out how the idetification toolbox works

[Ameer Hamza]

I am glad to be of help! You can find many examples in the documentation.

[Max Mai Tobon]

The transfer function I got was. H (s) = 847.4 / ( s + 757.2)

I'm now trying different constants for my PID control. Tomorrow I'll test the results.

Thanks again

Answer 2

[Ameer Hamza]

Another approach can be something like this. Consider the transfer function of the system is

and the transfer function of input step function is

The transfer function of output will be

which in time-domain become

Here you have y(t). You can use fmincon() to estimate a and b.

[Max Mai Tobon]

Gonna try this after I try the ToolBox.

Thanks again

[Max Mai Tobon]

The code for this should be fairly straightforward, yet I'm having some trouble with the fmincon() function.

Could you help me out with it?

Thanks

[Ameer Hamza]

Something like this

t; % time values
y; % measured signal values

fun = @(p) 25*p(1)/p(2)*(1-p(1)*exp(-p(2)*t));
objFun = @(p) norm(fun(p)-y);
sol = fmincon(objFun, rand(2,1));
a = sol(1);
b = sol(2);

[Max Mai Tobon]

I'm getting the following error:

Unrecognized function or variable 'getIpOptions'.

Error in fmincon (line 846)

options = getIpOptions(options,sizes.nVar,mEq,flags.constr,defaultopt,10,0.01);

Error in PID_Pressure (line 11)

sol = fmincon(objFun, rand(2,1));

How can I get this to work?

Thanks again. Sorry for the inconvenience.

[Ameer Hamza]

Do you have the optimization toolbox installed? Can you show the output of

ver optim

[Max Mai Tobon]

Oh, that was the problem! I didn't know I had to download an extra Toolbox.

It's working now.

Thanks a lot for your help! Very appreciated.

[Ameer Hamza]

I am glad to be of help!