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

Community Tip - Visit the PTCooler (the community lounge) to get to know your fellow community members and check out some of Dale's Friday Humor posts! X

20-Turquoise

## RMS bug.

RMS bug around 2pi. P6 file.

1 ACCEPTED SOLUTION

Accepted Solutions
24-Ruby V
(To:ttokoro)

Note that you must use the old legacy symbolics (muPad). It will not work when using the new friCAS.

Worksheet in P6 format attached.

12 REPLIES 12
24-Ruby V
(To:ttokoro)

Mathcad and Prime never had much luck with integrals, neither symbolically nor numerical as in your example.

Defining the function in a different way like the two shown in the pic does not help much:

23-Emerald III
(To:ttokoro)

Not necessarily restricted to Prime 6, Prime 4 shows the same.

The value of TOL has some effect.

As usual, symbolics is right:

How is that in Prime?

Luc

24-Ruby V
(To:LucMeekes)

How is that in Prime?

Here you are!

I gave it a try with Prime 6 using the new and future symbolics (friCAS) and the old one (muPad). The new symbolic is not really convincing. The simplification of muPad is (even without using simplify) even more compact than your Maple result 😉

21-Topaz II
(To:ttokoro)

Hi,

you shouldn't square the function containing the Heaviside function but only the function defined in the window created with the step function. See the picture:

don't square U(t). However, the rms is still zero.....

23-Emerald III
(To:-MFra-)

Squaring the sine 'within' gives you a different function:

with a different symbolic result:

but numerically it's still flawed...

Luc

24-Ruby V
(To:LucMeekes)

I think you misunderstood what Francesco meant.

The idea was not to replace V(t) by his U(t) but rather replace V(t)^2 by his U(t).

I am not sure what the benefit of doing so would be, though.

BTW, here is more clearly the deficiency of the new symbolics (friCAS) in Prime when it comes to integration with the Heaviside function with the step inside the integration interval:

We always moaned about Prime would not be developed (or not fast enough). Replacing the very old symbolics engine muPad by a new, more modern one sure is a good a big development step. When we moaned in the past about development, have we forgot to state the direction in which we would like to see development taking place? We did not meant the direction of making things worse!!

23-Emerald III
(To:Werner_E)

Yes, I misunderstood.

After investigation: The intended method should not make a difference, because PHI only takes on values 0 or 1, and likewise PHI^2 only takes on values 0^2=0 or 1^2=1. Mathcad 11 symbolics still gives the same (correct) answers. Numerically it's still a mess.

Luc

20-Turquoise
(To:ttokoro)

Using 3D plot, the bug region is visible.

For Fourier analysis, if the signal is f(t)=-f(t+T/2) integral from 0 to pi is OK. Otherwise, we should integral from 0 to 2*pi. Therefore, this bug is sever problem.

24-Ruby V
(To:ttokoro)

Note that you must use the old legacy symbolics (muPad). It will not work when using the new friCAS.

Worksheet in P6 format attached.

21-Topaz II
(To:ttokoro)

Hi,

Mathcad 15 gives almost the same results. However, if we change the integration algorithm the result is not null; see photo:

23-Emerald I
(To:ttokoro)

Attached is a slightly different approach (Prime 6.0)

20-Turquoise
(To:Fred_Kohlhepp)

Thanks for using time for my question.

In electricel engineering RMS of sin(t) must 1/sqrt(2)=0.7. The denominator is fixed to T, that is a period of fundamental signal, in this case 2pi.

Announcements
Top Tags