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

Community Tip - New to the community? Learn how to post a question and get help from PTC and industry experts! X

RMS bug.

ttokoro
20-Turquoise

RMS bug.

 

RMS bug around 2pi. P6 file.

image.png

1 ACCEPTED SOLUTION

Accepted Solutions

How about this approach?

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

clipboard_image_1.png

Worksheet in P6 format attached.

View solution in original post

12 REPLIES 12

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:

clipboard_image_0.png

LucMeekes
23-Emerald III
(To:ttokoro)

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

Mathcad 11 too:

clipboard_image_0.png

The value of TOL has some effect.

As usual, symbolics is right:

clipboard_image_1.png

How is that in Prime?

 

Luc


clipboard_image_1.png

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 😉

clipboard_image_0.png

-MFra-
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:

ttokoro bug.jpg 

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

LucMeekes
23-Emerald III
(To:-MFra-)

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

clipboard_image_0.png

with a different symbolic result:

clipboard_image_1.png

but numerically it's still flawed...

clipboard_image_2.png

Luc

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:

clipboard_image_1.png

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!!

LucMeekes
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

ttokoro
20-Turquoise
(To:ttokoro)

image.pngUsing 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.

 

How about this approach?

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

clipboard_image_1.png

Worksheet in P6 format attached.

-MFra-
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:

ttokoro bug +.jpg

Attached is a slightly different approach (Prime 6.0)

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.

Top Tags