mirror of
git://projects.qi-hardware.com/ben-wpan.git
synced 2024-12-23 22:25:30 +02:00
bcd6e09189
- range: generate min/avg/max data from FFT files - vizr: plot the output of "range" - row: make "viz" keep the FFT output and use it to generate the summary
54 lines
870 B
Perl
Executable File
54 lines
870 B
Perl
Executable File
#!/usr/bin/perl
|
|
|
|
sub read
|
|
{
|
|
local ($n) = "$_[0].fft";
|
|
local (@f);
|
|
local $i = 0;
|
|
local $ym, $v;
|
|
|
|
open(FILE, $n) || die "$n: $!";
|
|
while (<FILE>) {
|
|
chop;
|
|
push(@f, $_);
|
|
$ym = $_ unless defined $ym && $ym > $_;
|
|
}
|
|
close FILE;
|
|
for (@f) {
|
|
$v = $_-$ym;
|
|
$min[$i] = $v unless defined $min[$i] && $min[$i] < $v;
|
|
$max[$i] = $v unless defined $max[$i] && $max[$i] > $v;
|
|
$sum[$i] += $v;
|
|
$i++;
|
|
}
|
|
}
|
|
|
|
|
|
sub usage
|
|
{
|
|
print STDERR "usage: $0\n";
|
|
exit(1);
|
|
}
|
|
|
|
|
|
&usage if @ARGV;
|
|
|
|
#
|
|
# @@@ for future study: use an x axis proportional to the center frequency,
|
|
# instead of an absolute frequency.
|
|
#
|
|
|
|
for ($f = 2405; $f <= 2480; $f += 5) {
|
|
&read($f);
|
|
&read($f+1);
|
|
$n += 2;
|
|
}
|
|
|
|
$mhz_per_bin = 100.0/16/@sum; # 100 MHz ADC, decimation 16
|
|
|
|
$, = " ";
|
|
for ($i = 0; $i != @sum; $i++) {
|
|
print $mhz_per_bin*($i-@sum/2)+0.5, "";
|
|
print $min[$i], $sum[$i]/$n, $max[$i], "\n";
|
|
}
|