to get the output of this system with differential equations.

i got this systema represented by his transfer function (and the equations):

i want to get the value of the output (surounded with a black circle)

and i got this code:

function xdot=diagrama_bloques_nuria(t,x)

Ka=10;
Kb=20;
Km=10;
Ra=10;
Ks=0.75;
Kf=1;
Jm=0.01;
Ke=0.25;
A=5;
Vp=5;

Vr=Vp*square(2*pi*t)+5;
Va=(Vr-Kf*x(1, :))*Ka;

xdot = [
        (x(2, :)*Ke)-(x(1, :)*Ks);
        x(3, :);
       (Km/(Ra*Jm))*(Va-(Kb*x(3, :)));
       ];

%[t,x]=ode45('diagrama_bloques_nuria',[0,10],[0,0,0])
%plot(t, x(:, 1));

this code represent the system without the 1/A gain. But i don't know how to get the output with the 1/A gain.

Answer 1

[jose luis guillan suarez]

I solved it in this way

function xdot=diagrama_bloques_nuria(t,x)

Ka=10;
Kb=20;
Km=10;
Ra=10;
Ks=0.75;
Kf=1;
Jm=0.01;
Ke=0.25;
A=5;
Vp=5;

Vr=Vp*square(2*pi*t)+5;
Va=(Vr-Kf*x(1, :))*Ka;

xdot = [
        ((x(2, :)*Ke)-(x(1, :)*Ks))*1/A;
        x(3, :);
       (Km/(Ra*Jm))*(Va-(Kb*x(3, :)));
       ];

i multiplied the xdot(1) by 1/A.