Create custom component with custom formula or transfer function?

Questions
, , , ,
1

How can I create a custom component with custom equation of transfer function? I have an equation for a low pass filter as shown below

L(t)= exp(-t/tau)/tau

Thanks

2

Hi Lakshmi,

You can check out the following videos in Wolfram U on creating custom components in SystemModeler. https://www.wolfram.com/wolfram-u/modeling-systems/

If you are interested in creating your component in Mathematica, see the documentation on CreateSystemModel here: CreateSystemModel: Create a System Model—Wolfram Documentation

Thanks,
Ankit

3

Lakshmi,

You can also implement your transfer function using a continuous “TransferFuntion” block. In your example, you would also need a FixedDelay to handle the transportation lag e(-t *s) in the numerator.

Regards,

Neil

4

Hi, Thank you.
I created this transfer function using CreateSystemModel

RTFV2= CreateSystemModel[TransferFunctionModel[{{{E^(-s/\[Tau])}}, \[Tau]}, s]]

but when I used it in system model it gives me an error that there is no scope for variable s.
Alternatively, I also tried to modify the in-built transfer function in modelica but was not sure how to use the exponent. What is the equation of y? Please let me know.
Lakshmi

5

Lakshmi,

You have a few issues.

  1. Your system is a time delayed system (it has a transportation lag). This system cannot be represented with a transfer function Model in modelica because it is not a polynomial transfer function.

  2. You can make the correct system as a transfer function of 1/s and a time delay (e^-(tau*s)). The time delay is tau in this example but see my point 3 below.

  3. you wrote e^-s/tau are you sure this is correct? In your case your time delay is 1/tau – did you intend that? or did you intend a time delay of tau seconds?

  4. You can construct your model in WSM or you can construct the model in WMA by “wiring” the blocks together to form the model described above using CreateSystemModel to create the model.

Regards,

Neil

6

Lakshmi,

To be clear, The model looks like this:

ScreenShot2019-05-22at9.44.59PM.png
ScreenShot2019-05-22at9.44.59PM.png760Ă—408 16.4 KB

and the code is

model example
  Modelica.Blocks.Continuous.TransferFunction transferFunction1(a = {1, 0}, b = {0, 1}) annotation(Placement(visible = true, transformation(origin = {-55, 0}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
  Modelica.Blocks.Nonlinear.FixedDelay fixedDelay1(delayTime = 1 / 1) annotation(Placement(visible = true, transformation(origin = {-5, 0}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
equation
  connect(transferFunction1.y, fixedDelay1.u) annotation(Line(visible = true, origin = {-30.5, 0}, points = {{-13.5, 0}, {13.5, 0}}, color = {1, 37, 163}));
  annotation(Diagram(coordinateSystem(extent = {{-150, -90}, {150, 90}}, preserveAspectRatio = true, initialScale = 0.1, grid = {5, 5})), Icon(coordinateSystem(extent = {{-100, -100}, {100, 100}}, preserveAspectRatio = true, initialScale = 0.1, grid = {10, 10}), graphics = {Rectangle(visible = true, lineColor = {0, 114, 195}, fillColor = {255, 255, 255}, extent = {{-100, -100}, {100, 100}}, radius = 25), Text(visible = true, textColor = {64, 64, 64}, extent = {{-150, 110}, {150, 150}}, textString = "%name")}));
end example;

The fixed delay is your exponent term.

Regards

7

Hi Neil,
Thanks for your response. I have 2 variables in the function. One is tau and other is “s”. What i am trying to build is a standard equation describing the impulse response of a low pass filter with time constant tau. I want to use this as my filter in system modeler.

Ltau(s)= [Exp(-s/?au)]/ ?au where
Ltau is a function of time (s) and tau is my time constant

Please find attached picture of the WMA script and output.

myfunction.png
myfunction.png233Ă—150 1.71 KB

Thanks
Lakshmi

8

Lakshmi,

I’m confused. You refer to this as your transfer function and use it that way in the “TransferFunction” command in the way you would specify a frequency domain expression. However, now you say it’s a function of time. Is your expression an impulse response in time and not a transfer function in frequency ?

Regards,

Neil

9

Hi,
Ltau(s) is a function of time, and tau is the time constant. Its the impulse response of an IIR filter. I recently started learning Mathematica and system modeler. I thought the appropriate component in blocks to do this would be a transfer function or a first order filter. Please let me know

Thanks
Lakshmi

10

Lakshmi,

You confused me because TransferFunctionModel takes a frequency domain expression and your expression looked like a frequency domain expression because time was called “s” (which is usually reserved for the complex variable in Laplace Transforms).

You should do this

exp = E^(-t/\[Tau])/\[Tau]
CreateSystemModel["myFilter", 
 TransferFunctionModel[Simplify[LaplaceTransform[exp, t, s]], s]]

Regards,

Neil

11

Thank you very much Neil
Lakshmi