Defining an improper function in Simulink

19 次查看（过去 30 天）
Michael
Michael 2024-1-18
编辑： Sam Chak 2024-1-18
I want to implement this transfer function model related to viscous friction in simulink : Ms^2 + ws which means it does not have a denominator thus cannot be defined in simulink since it has a greater value of numerator than the denominator meaning it is an improper function, are there any workarounds that i can implements so that i can define this transfer function? Thank you very much.
  1 个评论
Sam Chak
Sam Chak 2024-1-18
编辑：Sam Chak 2024-1-18
Ran in:
Hi @Michael,
Before incorporating @Andy Bartlett's practical suggestion to approximate the improper transfer function, are you attempting to cancelling out one pole at the origin and another at ? Additionally, "Permodelan Penyesuaian" is loosely translated to mean Adaptive Modeling. Could this possibly refer to an Adaptive Control System?
Update: Check if this is want you want
If you're attempting to model the motion of a mass M on a horizontal surface with a viscous coefficient W and an applied force F, the governing equation is likely given by
.
Taking the Laplace transform, you obtain
.
Thus, rearranging the equation to get the transfer function
.
M = 2;
W = 3;
G = tf(1, [M W 0])
G = 1 ----------- 2 s^2 + 3 s Continuous-time transfer function.
%% If X is measurable and is fed back to form a closed-loop, then
Gcl = feedback(G, 1)
Gcl = 1 --------------- 2 s^2 + 3 s + 1 Continuous-time transfer function.
%% Step response
step(Gcl), grid on

请先登录，再进行评论。

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

回答（2 个）

Andy Bartlett
Andy Bartlett 2024-1-18
The Laplace operator 1/s is an integrator and s is a derivative.
You can break apart your equation and implement it with a gain block and a derivative block.
But modeling things with derivatives is generally not great, super sensitive to noise.
A better practice is to model things using integrators.
I suggest looking at what you are trying to model from a higher level perspecitve. For that higher level view, see if you can re-arrange the full set of equations to use integrators instead of derivatives.
  2 个评论
Michael
Michael 2024-1-18
Thank you so much for the answer, any idea on how to implement it? i already tried using the residue method mentioned in other post to break apart the equation but numerator ended up having zero value.
Andy Bartlett
Andy Bartlett 2024-1-18
>> Implement it?
You mean the equation M*s^2 + W*s ?
Put two derivative blocks in series.
In parallel feed the output of the first derivative to a Gain of W
feed the output of the second derivative to a Gain of M
Another practical hack is to add "fast poles" to the transfer function.
p = value just big enough, maybe 10 rad/s, 100 rad/s, ...
So transfer function becomes
(M*s^2 + W*s) / ( (1/p)*s + 1 )^2
p big enough to allow accurate modeling of the low frequency dynamics that matter
but not so big that too much high frequency noise is allowed into the model.

请先登录，再进行评论。

Aquatris
Aquatris 2024-1-18
编辑：Aquatris 2024-1-18
Ran in:
Firstly, I do not think you understand the viscous friction model. Ms2+Ws is probably not what you think it is.
To answer your question, one workaround is to put a low pass filter with a large enough bandwidth to have negligible effects. This would still give you same dynamics at 'low' frequency which is what ou are interested. This is also one of the reason why in most PID implementations, the derivative term is implemented with a filter. Example;
m = 10;w = 5;
s = tf('s');
x = m*s^2+w*s; % actual tf we are interested in
x2 = x/(s/1000+1)^2; % tf to be used in simulation that would give same dynamic response
% 1000 was chosen arbitrarily. ideally it should be infinity. Too large of
% a number would result slow simulation and too small would change the
% dynamics too much so you should choose it with care. generaly 10 times
% the interested bandwidth should work fine
bode(x,x2)

类别

在 Help Center 和 File Exchange 中查找有关 General Applications 的更多信息

标签

产品


版本

R2021b

Community Treasure Hunt

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

Start Hunting!

Translated by