09-12-2018
Pure awk FFT.
Hi Corona688...
I just knew you would create a translation and get it working on awk. ;o)
You could remove the SIN_ANGLE and COS_ANGLE and put sin(ANGLE) and cos(ANGLE) directly into the two lines below them, eliminating two lines...
Thanks for the translation, I will see if I can eliminate the Python version in AudioScope.sh and attach your name to it...
Another thing is that there is no need to issue the imaginary values for this basic FFT. As long as the real values are actually, cropped or padded to power(s) of 2 then the imaginary values of 0.0 can be internally generated to the same array length as the real values, eliminating the IMAGSTR initialisation on the command line call.
But nice to see an awk version that I can steal... <wink>
Last edited by wisecracker; 09-14-2018 at 07:45 AM..
Reason: Add title...
10 More Discussions You Might Find Interesting
1. UNIX for Dummies Questions & Answers
Hello All,
Below is a simple script i worte to find the 208th char in a file. If the char = "C" then I re-direct the line to a file called change.txt. If it is not "C" then I re-direct it to a file called delete.txt.
My problem is I have a file 0f 500K lines. this script is very slow. I am... (4 Replies)
Discussion started by: eja
4 Replies
2. Shell Programming and Scripting
Any means of running ksh93 in a ksh88-mode?
Might sound odd, but I want/need to restrict U/Win-developed scripts to correspond to the ksh88 version on my Solaris environment(s). Thanks. (2 Replies)
Discussion started by: curleb
2 Replies
3. Shell Programming and Scripting
Hi all!
Does anybody know how can I check if any UNIX installation has implemented ksh88 or ksh93?
Thanks in advance.
Néstor. (3 Replies)
Discussion started by: Nestor
3 Replies
4. Shell Programming and Scripting
I have a text file that contains 4 million lines, each line contains 2 fields(colon as field separator). as shown:
123:444,555,666,777,888,345
233:5444,555,666,777,888,345
623:454,585,664,773,888,345
......
Here I have to split the second field(can be up to 40,000 fields) by comma into an... (14 Replies)
Discussion started by: kevintse
14 Replies
5. Shell Programming and Scripting
Hi,
I have a script that, basically, has two input files of this type:
file1
key1=value1_1_1
key2=value1_2_1
key4=value1_4_1
...
file2
key2=value2_2_1
key2=value2_2_2
key3=value2_3_1
key4=value2_4_1
...
My files are 10k lines big each (approx).
The keys are strings that don't... (7 Replies)
Discussion started by: fzd
7 Replies
6. Shell Programming and Scripting
Hi,
I have a large number of input files with two columns of numbers.
For example:
83 1453
99 3255
99 8482
99 7372
83 175
I only wish to retain lines where the numbers fullfil two requirements. E.g:
=83
1000<=<=2000
To do this I use the following... (10 Replies)
Discussion started by: s052866
10 Replies
7. UNIX for Advanced & Expert Users
Hi Experts,
I have several shell scripts that have been developed on a Linux box for korn ksh93.
If we want to run this software on an AIX 6.1 box that runs ksh88 by default can we just change the she-bang line to reference /bin/ksh93 which ships with AIX as its "enhanced shell" to ensure... (6 Replies)
Discussion started by: Keith Turley
6 Replies
8. OS X (Apple)
I don't know if anyone is interested but I have been meddling with FFT for the AMIGA.
(Sadly we AMIGAns don't have these luxuries through any scripting language.
Below is a Python snippet that uses the builtin 'cmath' module to work with the lowly
Python 2.0.1 for the AMIGA. It is part of a... (0 Replies)
Discussion started by: wisecracker
0 Replies
9. AIX
Hi everyone !
Im trying to know from wich version of AIX KSH93 is available ?
Internet tell me 6.x and 7.x AIX are available, bue what about 5.x ?
Is KSH93 available on AIX 5.x ?
Is it the same way to manipulate variables as KSH93 on 7.x ?
Thanks for your support and have a nice day ! (2 Replies)
Discussion started by: majinfrede
2 Replies
10. OS X (Apple)
Hi guys...
This is code that was originally designed to work on an upgraded AMIGA A1200 using Python 2.0.x.
Unfortunately it broke inside much later versions, NOT because of the print statement/function but other minor subtleties. So this is the final result tested on various machines including... (0 Replies)
Discussion started by: wisecracker
0 Replies
LEARN ABOUT OPENSOLARIS
mlib_signalfft_4
mlib_SignalFFT_4(3MLIB) mediaLib Library Functions mlib_SignalFFT_4(3MLIB)
NAME
mlib_SignalFFT_4, mlib_SignalFFT_4_S16_S16, mlib_SignalFFT_4_S16C_S16C, mlib_SignalFFT_4_S16C_S16, mlib_SignalFFT_4_S16, mlib_Sig-
nalFFT_4_S16C - signal Fast Fourier Transform (FFT)
SYNOPSIS
cc [ flag... ] file... -lmlib [ library... ]
#include <mlib.h>
mlib_status mlib_SignalFFT_4_S16_S16(mlib_s16 *dstr, mlib_s16 *dsti,
const mlib_s16 *srcr, const mlib_s16 *srci, mlib_s32 order,
mlib_s32 *scale);
mlib_status mlib_SignalFFT_4_S16C_S16C(mlib_s16 *dstc,
const mlib_s16 *srcc,
mlib_s32 order, mlib_s32 *scale);
mlib_status mlib_SignalFFT_4_S16C_S16(mlib_s16 *dstc,
const mlib_s16 *srcr,
mlib_s32 order, mlib_s32 *scale);
mlib_status mlib_SignalFFT_4_S16(mlib_s16 *srcdstr, mlib_s16 *srcdsti,
mlib_s32 order, mlib_s32 *scale);
mlib_status mlib_SignalFFT_4_S16C(mlib_s16 *srcdstc, mlib_s32 order,
mlib_s32 *scale);
DESCRIPTION
Each of the functions in this group performs Fast Fourier Transform (FFT).
The following equation is used for forward FFT:
1 N-1
dst[k] = ---- SUM {src[n] * exp(-j2*PI*n*k/N)}
C1 n=0
and the following equation is used for inverse FFT (IFFT):
1 N-1
dst[n] = ---- SUM {src[k] * exp(j2*PI*n*k/N)}
C2 k=0
where
k = 0, 1, ..., (N - 1)
n = 0, 1, ..., (N - 1)
N = 2**order
The signal FFT/IFFT functions can be categorized into four groups according to the ScaleMode in the function names in the following form:
mlib_Signal[FFT|IFFT]_ScaleMode_OutType_InType_OpMode()
mlib_Signal[FFT|IFFT]_ScaleMode_DataType_OpMode()
The scaling factors C1 and C2 used in the equations are defined as follows:
o For ScaleMode = 1, C1 = 1 and C2 = 2**order.
o For ScaleMode = 2, C1 = 2**order and C2 = 1.
o For ScaleMode = 3, C1 = C2 = 2**(order/2) when order is even, or C1 = 2**((order+1)/2) and C2 = 2**((order-1)/2) when order is
odd.
o For ScaleMode = 4, C1 = 2**P and C2 = 2**Q, where P and Q are adaptive scaling factors and are generated by the functions.
For functions with only real parts for the source signal, the imaginary parts are assumed to be all zero. For functions with only real
parts for the destination signal, the imaginary parts are discarded. The functions with only one data type in their names perform the oper-
ation in place.
PARAMETERS
Each function takes some of the following arguments:
dstr Destination signal array that contains the real parts.
dsti Destination signal array that contains the imaginary parts.
srcr Source signal array that contains the real parts.
srci Source signal array that contains the imaginary parts.
dstc Complex destination signal array. dstc[2*i] contains the real parts, and dstc[2*i+1] contains the imaginary parts.
srcc Complex source signal array. srcc[2*i] contains the real parts, and srcc[2*i+1] contains the imaginary parts.
srcdstr Source and destination signal array that contains the real parts.
srcdsti Source and destination signal array that contains the imaginary parts.
srcdstc Complex source and destination signal array. srcdstc[2*i] contains the real parts, and srcdstc[2*i+1] contains the imaginary
parts.
order Order of the transformation. The base-2 logarithm of the number of data samples.
scale Adaptive scaling factor.
RETURN VALUES
The function returns MLIB_SUCCESS if successful. Otherwise it returns MLIB_FAILURE.
ATTRIBUTES
See attributes(5) for descriptions of the following attributes:
+-----------------------------+-----------------------------+
| ATTRIBUTE TYPE | ATTRIBUTE VALUE |
+-----------------------------+-----------------------------+
|Interface Stability |Committed |
+-----------------------------+-----------------------------+
|MT-Level |MT-Safe |
+-----------------------------+-----------------------------+
SEE ALSO
mlib_SignalFFT_1(3MLIB), mlib_SignalFFT_2(3MLIB), mlib_SignalFFT_3(3MLIB), mlib_SignalIFFT_1(3MLIB), mlib_SignalIFFT_2(3MLIB), mlib_Sig-
nalIFFT_3(3MLIB), mlib_SignalIFFT_4(3MLIB), attributes(5)
SunOS 5.11 2 Mar 2007 mlib_SignalFFT_4(3MLIB)