cancel
Showing results for 
Search instead for 
Did you mean: 
cancel
Showing results for 
Search instead for 
Did you mean: 

Community Tip - Did you get an answer that solved your problem? Please mark it as an Accepted Solution so others with the same problem can find the answer easily. X

invlaplace in Mathcad 15

JohnWoodgate
6-Contributor

invlaplace in Mathcad 15

See the attached file. When I try invlaplace in Mathcad 15 with a simple s-function, it works, but with a slightly more complex function, it doesn't give a useful result:

I am a very infrequent user of Mathcad, so please bear with me if I ask stupid questions.

I do not understand the repeated 'invlaplace' and the 's,t' part of the result. I expected a function of t that I could plot after including the numerical value of omega. Setting t=35 and [omega] = 2*pi*t doesn't help.

I tried putting ', omega' after the first 'invlaplace',  but it produces no more useful result.

 

1 ACCEPTED SOLUTION

Accepted Solutions
LucMeekes
23-Emerald III
(To:JohnWoodgate)

Í don't wish to spoil the fun, but the expression is its own magnitude:

LM_20190608_Complex.png

Success!
Luc

View solution in original post

20 REPLIES 20

The result means that Prime is not able to find a symbolic solution for the transform.

The "invlaplace(...)" you see at the RHS is just a function which normally is use internally by the symbolics to derive the inverse laplace.

Writing "omega" after "invlaplace"  is definitely wrong as thsi would be interpreted as that the variable in the transform domain, the complex freqency, is named omega in your function.

Default is that the complex frequence is "s" and the variable in the time domain is "t". If this is the case you don't need to provide any further modifyers to "invlaplace" but you sure can write "invlaplace,s,t".

 

BTW, Wolram Alpha gives up on your function, too:

B.PNG

 

JohnWoodgate
6-Contributor
(To:Werner_E)

Thank you for the information.//Including omega was more or less a
desperate 'try anything', because omega is actually a constant, 54*pi in
the case under study.//I am quite surprised that the transform doesn't
work, because the expression is more or less a 'standard form' - the
frequency response of a sixth-order maximally-flat low-pass filter.  I
deliberately chose that as 'learning example'  because i know the
answer. I wanted to go on to more complex expressions for which I do not
know the answer. Maybe I will try them to see if they can be resolved.

That Wolfram gives up, too, is even more surprising.

LucMeekes
23-Emerald III
(To:JohnWoodgate)

The 'more or less standard' would be 1/sqrt(a^2+s^2), the inverse Laplace transform of which is J0(a*t).

Mathcad (that is, version 11, using Maple as symbolic engine) transfoms this pair back and forth.

LM_20190606_Laplace.png

As shown, you could rewrite the Laplace transfomed as 1/sqrt(1+(s/w)^2), provided that w is positive, but that's still power 2, far from power 12.

 

Success!
Luc

I meant 'standard result in electronics', not math. I've given up with
Mathcad. Matlab handles it, but I don't have it. I'm trying to
understand Octave now, but it's written (even the manual) for code
writers, not Windows users.
LucMeekes
23-Emerald III
(To:JohnWoodgate)

Your reaction suggests to me that you have an incorrect understanding of what Mathcad, and Matlab, can deliver.

I'm convinced that neither Mathcad, nor Matlab will give you the inverse Laplace tranform of the expression in s that you gave, based on the fact that Mathematica cannot give it and I didn't find it in my textbooks. 

On the other hand I know that both Mathcad and Matlab have functions that allow you to calculate the response in time on a signal through a filter that is described by your expression in s.

In Mathcad look up the butter(), cheby1(), cheby2(), and bessel() functions in combination with iirlow(), iirhigh(), iirpass() and iirstop(). Matlab has similar functions.

You may know that a certain expression in s can describe a specific filter. But you should not expect Mathcad, when requested to deliver the inverse laplace transform of a mathematical expression in s, to react: "hey, this is the expression for such and so filter, so let's give that.". Mathcad, and the symbolic processor of Matlab too, deal with mathematics, and mathematics only.

 

Success!
Luc

You are quite right. Invlaplace is not the right approach at all. What I want to do is to take an expression with s as the single variable and substitute omega*i  and then rationalize the result (i.e. convert the expression to an equivalent in real quantities) so that I can plot it. I think I have worked out teh syntax for using substitute, but I can't find anything about rationalization in the Help. Perhaps it's called something else

LucMeekes
23-Emerald III
(To:JohnWoodgate)

Maybe you are looking for this:

LM_20190607_Bode.png

 

Note that mathcad 15 supports logspace() or so to define a log array.

Success!

Luc

Hi,

Graph on the Gauss plane of the function of s of which you want the inverse Laplace transform (ω=1.5):

function.jpg

JohnWoodgate
6-Contributor
(To:-MFra-)

I don't know what a Gauss plane is, but in any case I now know that I
don't want an inverse Laplace transform. I need to replace s in the
expressions by omega*i and then rationalize, but I can't find anything
about rationalization in the Help.
LucMeekes
23-Emerald III
(To:JohnWoodgate)

To substitute s with omega*i is simple:

LM_20190607_Subst.png

Please explain 'rationalization'...

 

Luc

Hi,

In this file of mine (in fact there would be another four or five) you can find many calculations with Laplace transforms:

https://community.ptc.com/t5/PTC-Mathcad/5-Analog-Equivalent-Digital-Low-Pass-II-Order-Filter-pdf/m-p/451105#M175805.

JohnWoodgate
6-Contributor
(To:-MFra-)

I have already posted that I now realise that I don't want invlaplace or
laplace. I just want to substitute w*i (AKA j-omega) for s in a fraction
formed of binomials in s and then calculate the magnitude.  I was
originally confused by the way another software uses 'Laplace' as a keyword.
LucMeekes
23-Emerald III
(To:JohnWoodgate)

Like this?

LM_20190607_Complex.png

Success!
Luc

It could be the answer, but 's' should be 's/omega_c', where omega_c is a constant. The imaginary component should not be zero for all omega. I need to check with other expressions in s that seem to give wrong results with other software.

LucMeekes
23-Emerald III
(To:JohnWoodgate)

I assume your omega_c is a (real) constant. Dividing s by it does not make a big change, and it certainly will not help in getting an imaginary part for the total expression. Note that you raise s (or s/omega_c) to the 12th power. Since that power is an even number, and with s replaced with j*omega, you're effectively raising j to the 12th power. That will make it real.

If, however, you keep s as a complex number (with a complex and a real part, both non-zero) you will get a real and a complex part for the total expression.

 

Success!

Luc

Yes, I realise that. I need to investigate further.  The correct result
for the magnitude is /sqrt(1/(1+X//^12 //))/, where /X/
=/omega//omega_c, omega_c being a constant. Of course, being a
magnitude, this is a real quantity.  In this sort of expression for the
frequency response of an electrical filter, the exponent is always even,
being twice the 'order' of the filter, a whole number.
LucMeekes
23-Emerald III
(To:JohnWoodgate)

Í don't wish to spoil the fun, but the expression is its own magnitude:

LM_20190608_Complex.png

Success!
Luc

I have now recalled how to handle the expressions in s or w*i.  We can
close this thread now. Since the title is very misleading, the moderator
should consider deleting the thread.

With thanks to all who offered help.

LucMeekes
23-Emerald III
(To:JohnWoodgate)

I'm afraid that's not how it works.

I suggest you mark your last reaction as the answer to this topic, in that case it will be moved to the top, and people who are looking for info can still read what was going on.

 

Success!
Luc

I will try to do that. Meanwhile, for the avoidance of doubt, this
thread adds nothing useful for explaining the invlaplace keyword.
Top Tags