How can I obtain parameters K' and tau in transfer function G2?

18 views (last 30 days)
I have a transfer function and I need to approximate it for another which is simpler, a linear one. So I have:
U1=[K K*c];
D1=[1 p];
D11=[1 2*delta*w w^2];
D2=[tau 1];
G1 is the complex transfer function and G2 is the simple transfer function. The coefficients in G1 are known, but K' and tau in G2 are unknown and I have to find them according to the response I receive from y1's graph (y1 stands for the output due to an input, which is a step function, and to the transfer function G1). How can I determine K' and tau?
PS: For y1 I've used:
xlabel('Time (seconds)')

Accepted Answer

Christiaan on 11 Mar 2015
Dear Ricardo,
I wrote a matlab code for you (edited from this Mathworks example ) which could help you further:
clc;clear all;close all;
mass = 1; spring = 10; damping = 2;
num_orig = 1;
den_orig = [mass damping spring]; %order of tf is now 2
h_orig = tf(num_orig,den_orig) % this is the original transfer function
[sys,g] = balreal(h_orig); % balanced state-space realization
h_simpl = modred(sys,2,'MatchDC'); % simplified balanced state-space realization
[num_simpl,den_simpl] = ss2tf(h_simpl.a,h_simpl.b,h_simpl.c,h_simpl.d);
h_simpl = tf(num_simpl,den_simpl) % simplified transfer function
step(h_orig,'-r'); hold on
step(h_simpl,'-b'); hold off
legend('step response, original transfer function','step response, simplified transfer function','Location','northoutside')
If you modify if for your own code, a crititcal point is the second input argument of the modred function .
Good luck! Christiaan

More Answers (0)


Find more on Dynamic System Models in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!