dsp.SineWave

Generate discrete sine wave

Description

The dsp.SineWave System object™ generates a real or complex, multichannel sinusoidal signal with independent amplitude, frequency, and phase in each output channel.

For both real and complex sinusoids, the Amplitude, Frequency, and PhaseOffset properties can be scalars or length-N vectors, where N is the number of channels in the output. When you specify at least one of these properties as a length-N vector, scalar values specified for the other properties are applied to each of the N channels.

To generate a discrete-time sinusoidal signal:

1. Create the dsp.SineWave object and set its properties.

2. Call the object with arguments, as if it were a function.

To learn more about how System objects work, see What Are System Objects? (MATLAB).

Creation

Syntax

sine = dsp.SineWave

sine = dsp.SineWave(Name,Value)

sine = dsp.SineWave(amp,freq,phase,Name,Value)

Description

sine = dsp.SineWave creates a sine wave object that generates a real-valued sinusoid with an amplitude of 1, a frequency of 100 Hz, and a phase offset of 0. By default, the sine wave object generates only one sample.

sine = dsp.SineWave(Name,Value) creates a sine wave object with each specified property set to the specified value. Enclose each property name in single quotes.

Example: sine = dsp.SineWave('Amplitude',2);

sine = dsp.SineWave(amp,freq,phase,Name,Value) creates a sine wave object with the Amplitude property set to amp, Frequency property set to freq, PhaseOffset property set to phase, and anyother specified properties set to the specified values.

Properties

Unless otherwise indicated, properties are nontunable, which means you cannot change their values after calling the object. Objects lock when you call them, and the release function unlocks them.

If a property is tunable, you can change its value at any time.

For more information on changing property values, see System Design in MATLAB Using System Objects (MATLAB).

Amplitude — Amplitude of sine wave

1 (default) | scalar | vector

Amplitude of the sine wave, specified as one of the following:

• scalar –– A scalar applies to all channels.

• vector –– A length-N vector contains the amplitudes of the sine waves in each of the N output channels. The vector length must be the same as that specified for the Frequency and PhaseOffset properties.

Tunable: Yes

Dependencies

This property is tunable only when you set Method to either 'Trigonometric function' or 'Differential'.

Data Types: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64

Frequency — Frequency of sine wave

100 (default) | scalar | vector

Frequency of the sine wave in Hz, specified as one of the following:

• scalar –– A scalar applies to all channels.

• vector –– A length-N vector contains the frequencies of the sine waves in each of the N output channels. The vector length must be the same as that specified for the Amplitude and PhaseOffset properties.

Data Types: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64

PhaseOffset — Phase offset of sine wave

0 (default) | scalar | vector

Phase offset of the sine wave in radians, specified as one of the following:

• scalar –– A scalar applies to all channels.

• vector –– A length-N vector contains the phase offsets of the sine waves in each of the N output channels. The vector length must be the same as that specified for the Amplitude and Frequency properties.

Data Types: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64

ComplexOutput — Flag that indicates whether waveform is real or complex

false (default) | true

Flag that indicates whether the waveform is real or complex, specified as either:

• false –– The waveform output is real.

• true –– The waveform output is complex.

Method — Method used to generate sinusoids

'Trigonometric function' (default) | 'Table lookup' | 'Differential'

Method used to generate sinusoids, specified as one of the following:

• 'Trigonometric function' –– The object computes the sinusoid by sampling the continuous-time function.

• 'Table lookup' –– The object precomputes the unique samples of every output sinusoid at the start of the simulation, and recalls the samples from memory as needed.

• 'Differential' –– The object uses an incremental algorithm. This algorithm computes the output samples based on the output values computed at the previous sample time and precomputed update terms.

TableOptimization — Optimize table of sine values for speed or memory

'Speed' (default) | 'Memory'

Optimize table of sine values for speed or memory, specified as either:

• 'Speed' –– The table contains k elements, where k is the number of input samples in one full period of the sine wave. The period of each sinusoid must be an integer multiple of 1/Fs, where Fs is the value of the SampleRate property value. That is, each element of the Frequency property must be of the form Fs/m, where m is an integer greater than 1.

• 'Memory' –– The table contains k/4 elements.

Dependencies

This property applies only when you set the Method property to 'Table lookup'.

SampleRate — Sample rate of output signal

1000 (default) | positive scalar

Sample rate of output signal in Hz, specified as a positive scalar.

Example: 44100

Example: 22050

SamplesPerFrame — Number of samples per frame

1 (default) | positive integer

Number of consecutive samples from each sinusoid to buffer into the output frame, specified as a positive integer.

Example: 1000

Example: 5000

Data Types: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64

OutputDataType — Data type of the sine wave output

'double' (default) | 'single' | 'Custom'

Data type of the sine wave output, specified as 'double', 'single', or 'Custom'.

Fixed-Point Properties

CustomOutputDataType — Output word and fraction lengths

numerictype([],16) (default) | numerictype([],32,30)

Output word and fraction lengths, specified as an autosigned numeric type with a word length of 16.

Example: numerictype([],32,30)

Example: numerictype([],16,15)

Dependencies

This property applies only when you set the Method property to 'Table lookup' and the OutputDataType property to 'Custom'.

Usage

Syntax

sineOut = sine()

Description

sineOut = sine() creates the sine wave output, sineOut.

Output Arguments

sineOut — Sine wave output

vector | matrix

Sine wave output, returned as a vector or matrix. The SamplesPerFrame property determines the number of rows in the output matrix. If the Frequency or the PhaseOffset property is a vector, the length of the vector determines the number of columns (channels) in the output matrix. If the Frequency or the PhaseOffset properties is a scalar, then the number of channels in the output matrix is 1.

The OutputDataType property sets the data type of the output.

Data Types: single | double | fi

Object Functions

To use an object function, specify the System object as the first input argument. For example, to release system resources of a System object named obj, use this syntax:

release(obj)

Common to All System Objects

step	Run System object algorithm
release	Release resources and allow changes to System object property values and input characteristics
reset	Reset internal states of System object

Examples

Generate a Sine Wave Signal

Note: If you are using R2016a or an earlier release, replace each call to the object with the equivalent step syntax. For example, obj(x) becomes step(obj,x).

Generate a sine wave with an amplitude of 2, frequency of 10 Hz, and an initial phase of 0.

sine1 = dsp.SineWave(2,10);
sine1.SamplesPerFrame = 1000;
y = sine1();
plot(y)

Generate two sine waves offset by a phase of pi/2 radians.

sine2 = dsp.SineWave;
sine2.Frequency = 10;
sine2.PhaseOffset = [0 pi/2];
sine2.SamplesPerFrame = 1000;
y = sine2();
plot(y)

More About

Sinusoid

A real-valued, discrete-time sinusoid is defined as:

$$y(n)=A\sin (2\pi fn+\varphi )$$

where A is the amplitude, f is the frequency in Hz, and φ is the initial phase, or phase offset, in radians.

A complex sinusoid is defined as:

$$y(n)=A{e}^{j(2\pi fn+\varphi )}$$

Algorithms

Trigonometric Function

The trigonometric function method computes the sinusoid in the ith channel, y_i, by sampling the continuous function

$$\begin{array}{ll}{y}_i=A_i\sin \left(2\pi f_{i}t+{\varphi }_i\right)\hfill & \text{(real)}\hfill \\ \begin{array}{l}\\ \text{or}\\ \end{array}\hfill & \hfill \\ y_i=A_i{e}^{j\left(2\pi f_{i}t+{\varphi }_i\right)}\hfill & \text{(complex)}\hfill \end{array}$$

with a period of T_s, where you specify T_s in the sample time.

At each sample time, the algorithm evaluates the sine function at the appropriate time value within the first cycle of the sinusoid. By constraining trigonometric evaluations to the first cycle of each sinusoid, the algorithm avoids the imprecision of computing the sine of very large numbers. This constraint also eliminates the possibility of discontinuity during extended operations, when an absolute time variable might overflow. This method therefore avoids the memory demands of the table lookup method at the expense of many more floating-point operations.

Table Lookup

The table lookup method precomputes the unique samples of every output sinusoid at the start of the simulation, and recalls the samples from memory as needed. Because a table of finite length can only be constructed when all output sequences repeat, the method requires that the period of every sinusoid in the output be evenly divisible by the sample period. That is, 1/(f_iT_s) = k_i must be an integer value for every channel i = 1, 2, ..., N.

When the algorithm optimizes the table of sine values for Speed, the table constructed for each channel contains k_i elements. When the optimization is for Memory, the table constructed for each channel contains k_i/4 elements.

For long output sequences, the table lookup method requires far fewer floating-point operations than any of the other methods. However, the method can demand considerably more memory, especially for high sample rates (long tables). This method is recommended for models that are intended to emulate or generate code for DSP hardware, which need to be optimized for execution speed.

Note

The lookup table for this object is constructed from double-precision floating-point values. When you use the Table Lookup computation mode, the maximum amount of precision you can achieve in your output is 53 bits. Setting the word length of the output data type to values greater than 53 bits does not improve the precision of your output.

Differential

The differential method uses an incremental algorithm. This algorithm computes the output samples based on the output values computed at the previous sample time (and precomputed update terms) by making use of the following identities.

$$\begin{array}{l}\sin \left(t+T_s\right)=\sin \left(t\right)\cos \left(T_s\right)+\cos \left(t\right)\sin \left(T_s\right)\\ \cos \left(t+T_s\right)=\cos \left(t\right)\cos \left(T_s\right)-\sin \left(t\right)\sin \left(T_s\right)\end{array}$$

The update equations for the sinusoid in the ith channel, y_i, can therefore be written in matrix form as

$$\left[\begin{array}{c}\sin \left\{2\pi f_i\left(t+T_s\right)+{\varphi }_i\right\}\\ \cos \left\{2\pi f_i\left(t+T_s\right)+{\varphi }_i\right\}\end{array}\right]=\left[\begin{array}{cc}\cos \left(2\pi f_i{T}_s\right)& \sin \left(2\pi f_i{T}_s\right)\\ -\sin \left(2\pi f_i{T}_s\right)& \cos \left(2\pi f_i{T}_s\right)\end{array}\right]\left[\begin{array}{c}\sin \left(2\pi f_{i}t+{\varphi }_i\right)\\ \cos \left(2\pi f_{i}t+{\varphi }_i\right)\end{array}\right]$$

where you specify T_s in the sample time. Since T_s is constant, the right-hand matrix is a constant and can be computed once at the start of the simulation. The value of A_isin[2πf_i(t+T_s)+ϕ_i] is then computed from the values of sin(2πf_it+ϕ_i) and cos(2πf_it+ϕ_i) by a simple matrix multiplication at each time step.

This mode offers reduced computational load, but is subject to drift over time due to cumulative quantization error. Because the method is not contingent on an absolute time value, there is no danger of discontinuity during extended operations, when an absolute time variable might overflow.

Extended Capabilities

C/C++ Code Generation

Generate C and C++ code using MATLAB® Coder™.

Usage notes and limitations:

• This object has no tunable properties for code generation.

• See System Objects in MATLAB Code Generation (MATLAB Coder).