Mercurial > dsp > rdsp
view src/lib/signal/sinusoid.rs @ 24:daa09aecdc40 default tip
lib: rename phase_shift to sample_shift
Instead of tracking the phase shift in signals, we should track the sample
shift (i.e., delay) of the signal. The sample shift is a constant
regardless of the frequency - unlike phase shift. This simplifies the
tracking considerably.
author | Josef 'Jeff' Sipek <jeffpc@josefsipek.net> |
---|---|
date | Thu, 15 Apr 2021 22:26:06 -0400 |
parents | 2d5ae243921c |
children |
line wrap: on
line source
/* * Copyright (c) 2020-2021 Josef 'Jeff' Sipek <jeffpc@josefsipek.net> * * Permission is hereby granted, free of charge, to any person obtaining a copy * of this software and associated documentation files (the "Software"), to deal * in the Software without restriction, including without limitation the rights * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell * copies of the Software, and to permit persons to whom the Software is * furnished to do so, subject to the following conditions: * * The above copyright notice and this permission notice shall be included in * all copies or substantial portions of the Software. * * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE * SOFTWARE. */ use std::f64::consts::PI; use std::ops::Mul; use crate::{IQSample, Scalar}; use super::Signal; use super::{Gain, Mixer}; pub struct Sinusoid { rate: u64, time: u64, exp: Vec<IQSample>, } impl Sinusoid { pub fn new(rate: u64, freq: i64) -> Sinusoid { let mut exp = Vec::<IQSample>::with_capacity(rate as usize); for i in 0..rate { let t = i as Scalar / rate as Scalar; let f = freq as Scalar; exp.push(IQSample::new(0., 2. * (PI as Scalar) * f * t).exp()); } Sinusoid { rate: rate, time: 0, exp: exp, } } } impl Signal for Sinusoid { fn sample_rate(&self) -> u64 { self.rate } fn sample_size(&self) -> usize { std::mem::size_of::<IQSample>() } fn sample_shift(&self) -> u64 { 0 } fn freq(&self) -> u64 { 0 /* the center frequency is always 0 */ } } impl Iterator for Sinusoid { type Item = IQSample; fn next(&mut self) -> Option<Self::Item> { let t = self.time as usize; self.time = (self.time + 1) % self.rate; Some(self.exp[t]) } } include!("math.rs"); signal_math!(Sinusoid);