to get the output of this system with differential equations.

Question

jose luis guillan suarez 2018-6-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/404905-to-get-the-output-of-this-system-with-differential-equations

回答： jose luis guillan suarez 2018-6-9

采纳的回答： jose luis guillan suarez

在 MATLAB Online 中打开

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.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

jose luis guillan suarez 2018-6-9

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/404905-to-get-the-output-of-this-system-with-differential-equations#answer_323985

在 MATLAB Online 中打开

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, :)));
       ];