How to convert symbolic expressions to transfer functions

281 次查看（过去 30 天）
Qian Feng
Qian Feng 2016-10-31
评论： Walter Roberson 2023-8-7
I am encountering the problem of converting a symbolic expression to become a transfer function. Specifically, the linear system I am dealing with contains a non-constant distributed delay term which requires performing an integration to obtain the corresponding transfer function. However, it seems that the integration operator int cannot be applied with tf variables directly.
On the other hand, if there is a way to convert symbolic expressions to transfer functions, then this problem can be easily handled in symbolic setting first.
Thanks a lot
  3 个评论
显示 1更早的评论
Star Strider
Star Strider 2016-11-1
I suggested a similar approach yesterday. It’s apparently not a polynomial.
Qian Feng
Qian Feng 2016-11-2
编辑：Walter Roberson 2016-11-2
Here is the code,
r = 1;
s = tf('s');
syms x
A4 = [-1 x; -1-x^3 -1];
Ap = int(A4*exp(x*s),x, -r, 0);
The reason why we have an integration there is because I am dealing with a distributed delay term in the time domain.
The problem is that it seems we cannot mix a tf variable with a symbolic variable here.
However, the aforementioned integration can be easily handle if s is a symbolic variable, which is the reason why I asked about how to transfer a symbolic entity into a tf one.

请先登录，再进行评论。

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

采纳的回答

Walter Roberson
Walter Roberson 2016-11-2
  12 个评论
显示 10更早的评论
Paul
Paul 2021-2-27
@Walter Roberson
Cool code. Siight mod to also handle the case when symExp is a constant.
syms s
symExp(s) = 5;
ExpFun = matlabFunction(symExp);
ExpFun = str2func(regexprep(func2str(ExpFun), '\.([/^\\*])', '$1'));
TF = tf(ExpFun(tf('s')));
TF
TF =
5
Static gain.

请先登录，再进行评论。

更多回答（2 个）

HyunSang Park
HyunSang Park 2018-5-28
If you're just trying to find the peak value of the bode magnitude plot, might I suggest avoid using tf altogether? the peak value is when d(G(jw))/dw = 0. You can easily find the derivative with syms, and the plug in the w to the original tf.
Murugan venkatesan
In order to analyze the bifurcation, the input impedance expression how to plot the bode graph..
  1 个评论
Walter Roberson
Walter Roberson 2023-8-7
@Murugan venkatesan
I do not understand how people can use your answer to convert symbolic expressions to transfer functions? Could you show how your solution could be used for the example symExp = (s+2)/(s^2+5*s+9); ?

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Symbolic Math Toolbox 的更多信息

标签

产品

Community Treasure Hunt

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

Start Hunting!

Translated by