scirs2-signal 0.1.0-rc.2

Signal processing module for SciRS2 (scirs2-signal)
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76

use scirs2_signal::window::*;
use std::error::Error;

#[allow(dead_code)]
fn main() -> Result<(), Box<dyn Error>> {
    // Create windows of length 64
    let length = 64;

    // Generate different window types
    let hamming_win = hamming(length, true)?;
    let hann_win = hann(length, true)?;
    let blackman_win = blackman(length, true)?;
    let flattop_win = flattop(length, true)?;
    let kaiser_win = kaiser(length, 8.6, true)?;
    let boxcar_win = boxcar(length, true)?;
    let tukey_win = tukey(length, 0.5, true)?;

    // Print statistics about each window
    println!("Window comparison (length = {}):\n", length);

    println!("Hamming window:");
    print_stats(&hamming_win);

    println!("\nHann window:");
    print_stats(&hann_win);

    println!("\nBlackman window:");
    print_stats(&blackman_win);

    println!("\nFlat top window:");
    print_stats(&flattop_win);

    println!("\nKaiser window (beta=8.6):");
    print_stats(&kaiser_win);

    println!("\nBoxcar window:");
    print_stats(&boxcar_win);

    println!("\nTukey window (alpha=0.5):");
    print_stats(&tukey_win);

    Ok(())
}

#[allow(dead_code)]
fn print_stats(window: &[f64]) {
    // Calculate some basic statistics
    let min = window.iter().fold(f64::INFINITY, |a, &b| a.min(b));
    let max = window.iter().fold(f64::NEG_INFINITY, |a, &b| a.max(b));
    let sum: f64 = window.iter().sum();
    let mean = sum / window.len() as f64;

    // Calculate energy (sum of squares)
    let energy: f64 = window.iter().map(|&x| x * x).sum();

    println!("  Min value: {:.6}", min);
    println!("  Max value: {:.6}", max);
    println!("  Mean value: {:.6}", mean);
    println!("  Sum: {:.6}", sum);
    println!("  Energy: {:.6}", energy);

    // Print first and last few values
    println!(
        "  First 5 values: {:.6}, {:.6}, {:.6}, {:.6}, {:.6}",
        window[0], window[1], window[2], window[3], window[4]
    );
    let len = window.len();
    println!(
        "  Last 5 values: {:.6}, {:.6}, {:.6}, {:.6}, {:.6}",
        window[len - 5],
        window[len - 4],
        window[len - 3],
        window[len - 2],
        window[len - 1]
    );
}