Macaulay functions in MUPad

47 vues (au cours des 30 derniers jours)
Ed Davis
Ed Davis le 12 Nov 2014
In general the Macaulay function (<x-a>^n) is evaluated as 0 for x<a and x^n for x>a.
The integral is defined as int(<x-a>^n)=(<x-a>^(n+1))/(n+1)
I would like to implement functions of this type in MuPad to allow solving of beam deflection problems symbolically. However, I am not able to find a function that behaves as the Macaulay functions do. Any suggestions would be appreciated.
  3 commentaires
alejandro buitrago lopez
alejandro buitrago lopez le 20 Mai 2022
were any of you able to solve this, i' currently trying to program this so i can print the diagrams in matlab.
Nathan Hardenberg
Nathan Hardenberg le 3 Déc 2022
I did not solve the deflection problem, but I wrote a symbolic solver for shear force and bending moment. It is pretty well documented in the ReadMe file and you obviously can get some inspiration from the code.
I also found a Python implementation that should be able to symbolically calculate deflection, shear force etc. Having said this, I did not use it yet.

Connectez-vous pour commenter.

Réponses (1)

Nathan Hardenberg
Nathan Hardenberg le 29 Mai 2020
Not quite the solution you are looking for, but it is certainly possible to use the heaviside function. The heaviside function is also integratable and differentiable, so it is possible to use it instead.
<x - a>^0 is heaviside(x-a)
For functions with a higher exponent you have to write something like this
x^2 * heaviside(x-a)
this "switches" the function x^2 on (after a). Before it is 0
  2 commentaires
Walter Roberson
Walter Roberson le 3 Déc 2022
<x - a>^0 is heaviside(x-a)
No, the ^0 form is 1 everywhere, including at -inf and +inf and 0
a = 2;
x = linspace(-5,5);
y1 = (x-a).^0;
y2 = heaviside(x-a);
stairs(x, y1, '-*');
hold on
stairs(x, y2, ':+')
hold off
legend({'\^0', 'heaviside'}, 'Location', 'best')
ylim([-0.1 1.1])
Nathan Hardenberg
Nathan Hardenberg le 6 Déc 2022
You are right in the fact, that is not !
But here < and > are meant as macaulay brackets. With this notation my statement is correct. See the equation here.
It's probably better written as equation:
But note that the term is not defined at . In MATLAB is defined as at this position. It is sometimes also defined as 1.

Connectez-vous pour commenter.

Community Treasure Hunt

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

Start Hunting!

Translated by