## Simple symmetrical wave RMS DEMO using awk.

# 1
Simple symmetrical wave RMS DEMO using awk.

Hi everyone.
This is DEMO code for getting the RMS value of a symmetrical sampled waveform. It works a treat except see below.
I have used awk as the prime mover as CygWin does not havd 'bc' or 'dc' yet and this is going into AudioScope.sh soon.

However, I hit an anomaly.
In the code below parts B, C, D and G are more than acceptable and E was as expected and also within limits acceptable.
However I never expected this, assume the waveform starts at centreline 0V...
View part F in the code and using just 4 samples with the input waveform being an isosceles triangle then a stepped squarewave becomes the sample.
(The results on editing each waveform are in the code.)
As far as I am concerned there is no way of the sample knowing that the original signal was
either a sinewave or a triangle wave but the RMS result is that of a sinewave.
I can't find anything on a sample size of 4 on a triangle waveform and and its inversion so this anomaly stays.

Bazza.
 wisecracker View Public Profile for wisecracker Find all posts by wisecracker
# 2
It looks to me like you have a granularity problem. Four dots as a sample can fit a lot of functions, with intermediate datapoints from increased sampling rates you could then use fft to get good sine waves, for example.

You can code fft.c then call it using the system() function in awk. I think your awk supports that call.

A google for 'fft.c' got me 10 first page hits that looked respectably good. There are sites with an array of algorithmic variants to meet differing needs.

fft is gold standard for fitting data sampling output.
 jim mcnamara View Public Profile for jim mcnamara Find all posts by jim mcnamara
# 3
Hi JM...
Apologies for any typos, and hope it is lucid enough.
Quote:
Originally Posted by jim mcnamara
It looks to me like you have a granularity problem. Four dots as a sample can fit a lot of functions, with intermediate datapoints from increased sampling rates you could then use fft to get good sine waves, for example.

You can code fft.c then call it using the system() function in awk. I think your awk supports that call.

A google for 'fft.c' got me 10 first page hits that looked respectably good. There are sites with an array of algorithmic variants to meet differing needs.

fft is gold standard for fitting data sampling output.
Apologies for the delay, I have not been around for a while. My triangle wave is more or less hypothetical as sampling a real time random triangle wave that......
1) starts at the 0 point,
2) is initially positive going,
3) has 4 samples, or any number of the same 4 samples,
......is close to impossible.
But a file CAN be so easily be created as shown in the single cycle version in the code.
Although I was aware of the errors that creep in when creating the RMS values I was NOT aware of this one off bizarre situation, (and could find nothing on the WWW about it), so after creating the simple code I made a point of there being a reference to this anomaly.

As the RMS value(s) are near the exact value(s) for that of a sinewave hence my point.

In the days of Debian 6 and PCLinuxOS 2009, (my two main OSes in those days), when '/dev/dsp' was available, recreating the waveforms through '/dev/dsp' would result in a triangle wave as the sound system interpolation was Linear, this was also true for some Windows machines of bygone days. I am not sure what interpolations are used in current OSes and sound systems but it sure looks trigonometric to me, so I would guess that a sinewave would be recreated as a result.

I am happy with the idea of the RMS code at this point but looking more towards True RMS but would probably need to rethink my approach.

But to re-iterate for others, be aware of this anomaly.

LBNL, for any FFT work needed, I would do it through Python and SciPy as the heavy lifting is done already, see previous uploads of AudioScope.sh for the AF Spectrum display...

Bazza...

Last edited by wisecracker; 09-06-2017 at 12:18 PM.. Reason: Modify for error...
 wisecracker View Public Profile for wisecracker Find all posts by wisecracker
# 4
Your sample rate is way too low. It's just like trying to mathematically approximate the volume of a teapot by calculating the volume of a cube. Far too approximate to be useful!

When I plot your sine wave data, it looks awfully triangular. The crucial area going through the origin, accounting for the vast majority of the waveform, is a straight line from 38 to 217 -- and being symmetrical, can't be anything but a straight line. There's not enough points to represent anything else.

Mathematically speaking, triangle waves and square waves include higher frequencies than their base frequencies. Without a high enough sample rate to catch them, this information is thrown out like a brutal low pass filter. (Or worse, included as aliasing artifacts.)

Your waveforms should not be 8 points per wave. 16 at bare minimum, probably, but why not 1024? Then you'd see reliable differences, all the time.

Last edited by Corona688; 09-06-2017 at 07:42 PM..
 Corona688 View Public Profile for Corona688 Visit Corona688's homepage! Find all posts by Corona688

## Porting Rogue Wave to Linux

I am challenged with porting an old application from Solaris to Red Hat. The application uses Rogue Wave and I am searching for a Red Hat implementation. Your help is appreciated!