gauspuls: Gaussian-modulated sinusoidal RF pulse

Description Usage Arguments Details Value Note Author(s) Examples

Description

Generate a Gaussian modulated sinusoidal pulse sampled at times t.

Usage

1
gauspuls(t, fc = 1000, bw = 0.5)

Arguments

t

Vector of time values at which the unit-amplitude Gaussian RF pulse is calculated.

fc

Center frequency of the Gaussian-modulated sinusoidal pulses, specified as a real positive scalar expressed in Hz. Default: 1000

bw

Fractional bandwidth of the Gaussian-modulated sinusoidal pulses, specified as a real positive scalar.

Details

yi = gauspuls(t,fc,bw) returns a unit-amplitude Gaussian-modulated sinusoidal RF pulse at the times indicated in array t, with a center frequency fc in hertz and a fractional bandwidth bw

Value

Inphase Gaussian-modulated sinusoidal pulse, returned as a vector of unit amplitude at the times indicated by the time vector t.

Note

The Octave gauspuls function differs from the function in the Matlab signal package. The Octave function does not use the bwr (fractional bandwidth reference level), the tpe (trailing pulse envelope level), and the 'cutoff' input arguments. In addition, it only returns the inphase Gaussian pulse, not the quadrature pulse nor the signal envelope and cutoff time.

Author(s)

Original Matlab/Octave code Copyright (C) 2007 Sylvain Pelissier, Copyright (C) 2018-2019 Mike Miller. Port to R by Geert van Boxtel [email protected].

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
fs <- 11025    # arbitrary sample rate
t <- seq(-10, 10, 1/fs)
yi1 <- gauspuls(t, 0.1, 1)
yi2 <- gauspuls(t, 0.1, 2)
plot(t, yi1, type="l", xlab = "Time", ylab = "Amplitude")
lines(t, yi2, col = "red")

fs <- 11025  # arbitrary sample rate
f0 <- 100    # pulse train sample rate
x <- pulstran (seq(0, 4/f0, 1/fs), seq(0, 4/f0, 1/f0), "gauspuls")
plot (0:(length(x)-1) * 1000/fs, x, type="l",
      xlab = "Time (ms)", ylab = "Amplitude",
      main = "Gaussian pulse train at 10 ms intervals")

gjmvanboxtel/gsignal documentation built on Jan. 19, 2020, 3:52 a.m.