Message Boards Message Boards

GROUPS:

Create custom component with custom formula or transfer function?

Posted 1 month ago
404 Views
|
10 Replies
|
1 Total Likes
|

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

10 Replies

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: https://reference.wolfram.com/language/ref/CreateSystemModel.html

Thanks, Ankit

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

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

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

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. enter image description here

Thanks Lakshmi

Posted 1 month ago

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

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

Lakshmi,

To be clear, The model looks like this:

enter image description here

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

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

Thank you very much Neil Lakshmi

Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard

Group Abstract Group Abstract