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

// Simple test of elliptic filter design

use scirs2_signal::filter::ellip;

#[allow(dead_code)]
fn main() {
    println!("Testing Elliptic Filter Design");

    // Test different elliptic filter configurations
    let test_cases = vec![
        (3, 0.5, 30.0, 0.3, "lowpass"),
        (4, 0.5, 40.0, 0.3, "lowpass"),
        (5, 1.0, 35.0, 0.25, "lowpass"),
        (3, 0.5, 30.0, 0.4, "highpass"),
        (4, 1.0, 40.0, 0.35, "highpass"),
    ];

    for (i, (order, ripple, atten, cutoff, ftype)) in test_cases.iter().enumerate() {
        match ellip(*order, *ripple, *atten, *cutoff, *ftype) {
            Ok((b, a)) => {
                println!(
                    "\nTest {}: Order={}, Ripple={:.1}dB, Atten={:.0}dB, Type={}",
                    i + 1,
                    order,
                    ripple,
                    atten,
                    ftype
                );
                println!("  B coeffs: {:?}", format_coeffs(&b));
                println!("  A coeffs: {:?}", format_coeffs(&a));
                println!("  B length: {}, A length: {}", b.len(), a.len());
            }
            Err(e) => {
                println!("Test {}: Failed - {:?}", i + 1, e);
            }
        }
    }

    println!("\nElliptic filter design test completed!");
}

/// Format coefficients for display
#[allow(dead_code)]
fn format_coeffs(coeffs: &[f64]) -> String {
    coeffs
        .iter()
        .map(|&c| format!("{:.4}", c))
        .collect::<Vec<_>>()
        .join(", ")
}