step response from an equation

Question

Juan Jose Ortiz Torres 2018-9-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/420131-step-response-from-an-equation

评论： Dimitris Kalogiros 2018-9-21

Hello guys.

I am makig a program to plot the step response of a transfer function. But my program starts defining the numerator and the denominator of the transfer function like the following:

if true
  syms s
  num_planta=1; %numerator of the transfer function
  den_planta=s*(s+2); %denominator of the transfer function
end

Then I make the fraction

if true
transfer_fnc=num_planta/den_planta
end

Then with that variable I make a code works, but the code does not plot anything. So I decided to do it possible. So I use some techniques but nothing work. I am gonna explain the techniques, so if you can guide me, I will be grateful.

1. Using step(transfer_fcn)

if true
  aux_transfer_fcn=transfer_fcn
  SS=tf('s')
  subs(transfer_fcn,s,SS)
  step(aux_transfer_fcn)
end

But matlab says that I cannot substitute a "symbolic" variable to a "tf" variable

2. Making with matrix

if true
   [num_transfer_fcn,den_transfer_fcn]=numden(transfer_fcn)
   coefs_num_transfer_fcn=coeffs(num_transfer_fcn,s,'all')
   coefs_den_transfer_fcn=coeffs(den_transfer_fcn,s,'all') %%so it will be [1, 2, 0]
   system=tf(coefs_num_transfer_fcn,coefs_den_transfer_fcn)
end

But matlab does no accept matrices like [1,2,3] only [1 2 3], so I was looking up how to convert but I failed lot of times.

Could you please guide me?

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

请先登录，再进行评论。

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

Answer 1

Dimitris Kalogiros 2018-9-21

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/420131-step-response-from-an-equation#answer_337764

在 MATLAB Online 中打开

Use this:

clear; clc; 
  close all;
    % sympolic math toolbox
    syms s, 
    num_planta=sym(1); %numerator of the transfer function
    den_planta=s*(s+2); %denominator of the transfer function
    transfer_fnc=num_planta/den_planta
    %%extract coeffs of num and denum
    nump=sym2poly(num_planta)
    denp=sym2poly(den_planta)
    %%create transfer function using control toolbox
    H=tf(nump, denp);
    step(H) % step function