comm.CPMDemodulator

Demodulate signal using CPM method and Viterbi algorithm

expand all in page

Description

Thecomm.CPMDemodulatorSystem object™ demodulates an input signal that was modulated using the continuous phase modulation (CPM) method. The input is a baseband representation of the modulated signal. For more information about the demodulation and filtering applied, seeCPM Demodulation MethodandPulse Shape Filtering.

To demodulate a signal that was modulated using the CPM method:

Create thecomm.CPMDemodulatorobject and set its properties.
Call the object with arguments, as if it were a function.

To learn more about how System objects work, seeWhat Are System Objects?

Creation

Syntax

cpmdemod = comm.CPMDemodulator

cpmdemod = comm.CPMDemodulator(Name,Value)

cpmdemod = comm.CPMDemodulator(M,Name,Value)

Description

cpmdemod = comm.CPMDemodulatorcreates a demodulator System object to demodulate input CPM signals using the Viterbi algorithm.

example

cpmdemod = comm.CPMDemodulator(Name,Value)sets properties using one or more name-value arguments. For example,'SymbolMapping','Gray'specifies gray-ordered symbol mapping for the modulated symbols.

example

cpmdemod = comm.CPMDemodulator(M,Name,Value)sets theModulationOrderproperty toMand optional name-value arguments.

Properties

expand all

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

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

For more information on changing property values, seeSystem Design in MATLAB Using System Objects.

`ModulationOrder`—Modulation order
4(default) |power of two scalar

Modulation order, specified as a power-of-two scalar. The modulation order,M= 2^kspecifies the number of points in the signal constellation, wherekis a positive integer indicating the number of bits per symbol.

Data Types:double

`BitOutput`—Option to output data as bits
`0`or`false`(default) |`1`or`true`

Option to output data as bits, specified as a logical0(false) or1(true）.

Set this property tofalseto output data as integers.
Set this property totrueto output data as bits.

For more information, seeInteger-Valued and Binary-Valued Output Signals.

Data Types:logical

`SymbolMapping`—Symbol encoding
`'Binary'`(default) |`'Gray'`

Symbol encoding mapping of constellation bits, specified as'Binary'or'Gray'.

Set this property to'Binary'地图符号使用自然二进制编码排序.
Set this property to'Gray'to map symbols using Gray-coded ordering.

For more information, seeInteger-Valued and Binary-Valued Output Signals.

Dependencies

To enable this property, set theBitOutputproperty totrue.

`ModulationIndex`—Modulation index
`0.5`(default) |nonnegative scalar|column vector

Modulation index, specified as a nonnegative scalar or column vector. For more information, seeCPM Demodulation Method.

Data Types:double

`FrequencyPulse`—Type of frequency pulse shaping
`'Rectangular'`(default) |`'Raised Cosine'`|`'Spectral Raised Cosine'`|`'Gaussian'`|`'Tamed FM'`

Type of frequency pulse shaping used by the modulator to smooth the phase transitions of the modulated signal, specified as'Rectangular','Raised Cosine','Spectral Raised Cosine','Gaussian', or'Tamed FM'. For more information, seePulse Shape Filtering.

`MainLobeDuration`—Main lobe duration
`1`(default) |positive integer

Main lobe duration of the largest lobe in the spectral raised cosine pulse, specified as a positive integer representing the number of symbol intervals used by the demodulator to pulse-shape the modulated signal.

Dependencies

To enable this property, set theFrequencyPulseproperty to'Spectral Raised Cosine'.

Data Types:double

`RolloffFactor`—Roll-off factor
`0.2`(default) |scalar in the range [0, 1]

Roll-off factor of the spectral raised cosine pulse, specified as a scalar in the range [0, 1].

Dependencies

To enable this property, set theFrequencyPulseproperty to'Spectral Raised Cosine'.

Data Types:double

`BandwidthTimeProduct`—Product of bandwidth and symbol time of Gaussian pulse shape
`0.3`(default) |positive scalar

Product of the bandwidth and symbol time of the Gaussian pulse shape, specified as a positive scalar. UseBandwidthTimeProductto reduce the bandwidth, at the expense of increased intersymbol interference.

Dependencies

To enable this property, set theFrequencyPulseproperty to'Gaussian'.

Data Types:double

`PulseLength`—Length of frequency pulse shape
`1`(default) |positive integer

Length of the frequency pulse shape in symbol intervals, specified as a positive integer. For more information on the frequency pulse length, refer toLT在Pulse Shape Filtering.

Data Types:double

`SymbolPrehistory`—Symbol prehistory
`1`(default) |scalar|vector

Symbol prehistory, specified as scalar or vector with odd integer elements in the range [– (ModulationOrder– 1), (ModulationOrder– 1)]. This property defines the data symbols used by the modulator prior to the first call of the object in reverse chronological order. When you specify this property as a vector, the length must be (PulseLength– 1).

Data Types:double

`InitialPhaseOffset`—Initial phase offset
`0`(default) |scalar

Initial phase offset in radians of the modulated waveform, specified as a scalar.

Data Types:double

`SamplesPerSymbol`—Number of samples per input symbol
`8`(default) |positive integer

Number of samples per input symbol, specified as a positive integer. This property represents the number of samples input for each integer or binary word output. For all nonbinary schemes, as defined by the pulse shapes, this value must be greater than 1.

Data Types:double

`TracebackDepth`—Traceback depth for Viterbi algorithm
`16`(default) |positive integer

Traceback depth for the Viterbi algorithm, specified as a positive integer representing the number of trellis branches that the Viterbi algorithm uses to construct each traceback path. The value of this property is also the output delay and the number of zero symbols that precede the first meaningful demodulated symbol in the output. For more information, seeTraceback Depth and Output Delays.

Data Types:double

`OutputDataType`—Data type of output
`'double'`(default) |`'logical'`|`'int8'`|`'int16'`|`'int32'`

Data type of the output, specified as one of these values.

When you set theBitOutputproperty tofalse, you can set the output data type to'double','int8','int16', or'int32'.
When you set theBitOutputproperty totrue, you can set the output data type to'logical'or'double'.

Usage

Syntax

y = cpmdemod(x)

Description

y= cpmdemod(x)applies CPM demodulation method to the input signal and returns the demodulated signal.

Input Arguments

expand all

`x`—CPM-modulated signal
column vector

CPM-modulated signal, specified as a column vector with a length equal to an integer multiple of theSamplesPerSymbolproperty.

Data Types:double|single

Output Arguments

expand all

`y`— Output signal
column vector | matrix

Output signal, returned as a column vector or matrix. To specify whether the object outputs values as integers or bits, use theBitOutputproperty. To specify the output data type, use theOutputDataTypeproperty.

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 namedobj, use this syntax:

release(obj)

expand all

Common to All System Objects

`step`	RunSystem object算法
`release`	Release resources and allow changes toSystem objectproperty values and input characteristics
`reset`	Reset internal states ofSystem object

Examples

collapse all

CPM Modulate and Demodulate Signal with Gray Mapping and Bit Inputs

Open Live Script

Create CPM modulator, and CPM demodulator System objects.

cpmmodulator = comm.CPMModulator(8,...'BitInput',true,...'SymbolMapping','Gray'）;cpmdemodulator = comm.CPMDemodulator (8,...'BitOutput',true,...'SymbolMapping','Gray'）;

Create an error rate calculator System object™, that accounts for the delay caused by the Viterbi algorithm.

delay = log2(cpmdemodulator.ModulationOrder)...* cpmdemodulator.TracebackDepth; errorRate = comm.ErrorRate('ReceiveDelay',delay);

Transmit 100 3-bit words and print the error rate results.

forcounter = 1:100 data = randi([0 1],300,1); modSignal = cpmmodulator(data); noisySignal = awgn(modSignal,0); receivedData = cpmdemodulator(noisySignal); errorStats = errorRate(data,receivedData);endfprintf('Error rate = %f\nNumber of errors = %d\n',...errorStats(1),errorStats(2))

Error rate = 0.004006 Number of errors = 120

Apply GFSK Modulation and Demodulation

Open Live Script

Using thecomm.CPMModulatorandcomm.CPMDemodulatorSystem objects, apply Gaussian frequency-shift keying (GFSK) modulation and demodulation to random bit data.

Create a GFSK modulator and demodulator pair.

gfskMod = comm.CPMModulator(...'ModulationOrder',2,...'FrequencyPulse','Gaussian',...'BandwidthTimeProduct',0.5,...'ModulationIndex',1,...'BitInput',true); gfskDemod = comm.CPMDemodulator(...'ModulationOrder',2,...'FrequencyPulse','Gaussian',...'BandwidthTimeProduct',0.5,...'ModulationIndex',1,...'BitOutput',true);

Generate random bit data and apply GFSK modulation. Use a scatter plot to view the constellation.

numSym = 100; x = randi([0 1],numSym*gfskMod.SamplesPerSymbol,1); y = gfskMod(x); eyediagram(y,16)

数字Eye Diagram contains 2 axes objects. Axes object 1 with title Eye Diagram for In-Phase Signal contains an object of type line. This object represents In-phase. Axes object 2 with title Eye Diagram for Quadrature Signal contains an object of type line. This object represents Quadrature.

Demodulate the GFSK-modulated data. To verify that the demodulated signal data is equal to the original data, account for the delay introduced by the Gaussian filtering in the GFSK modulation and demodulation processes.

z = gfskDemod(y); delay = finddelay(x,z); isequal(x(1:end-delay),z(delay+1:end))

ans =logical1

More About

expand all

CPM Demodulation Method

The CPM demodulation method process consists of a correlator followed by a maximum-likelihood sequence detector (MLSD) that searches the paths through the state trellis for the minimum Euclidean distance path. When the modulation index is rational (h=m/p), a finite number of phase states exist in the symbol. The implementation uses the Viterbi algorithm to perform MLSD.

{h_i} is a sequence of modulation indices that moves cyclically through a set of indices {h₀,h₁,h₂、……h_H-1}.

h_i=m_i/p_iis the modulation index in proper rational form.
m_iis the numerator of the modulation index.
p_i的分母是the modulation index.
m_iandp_iare relatively prime positive numbers.
The least common multiple (LCM) of{p₀,p₁,p₂、……p_H-1}is denoted asp.
h_i=m'_i/p.

{h_i} determines the number of phase states,

$n u m P h a s e S t a t e s = {\begin{cases} p, for all even m'_{i} \\ 2 p, for any odd m'_{i} \end{cases}},$

and affects the number of trellis states,

numStates=numPhaseStates×M^(L-1),

Lis the pulse length.
Mis the modulation order.

CPM Method

The input to the demodulator is a baseband representation of the modulated signal:

$\begin{array}{l} s (t) = \exp [j 2 π \sum_{i = 0}^{n} α_{i} h_{i} q (t - i T)], and \\ n T < t < (n + 1) T . \end{array}$

{α_i} is a sequence ofM-ary data symbols selected from the alphabet ±1, ±3, ±(M–1).
Mmust have the form 2^kfor some positive integerk, whereMis the modulation order and specifies the size of the symbol alphabet.
{h_i} is a sequence of modulation indices, andh_imoves cyclically through a set of indices {h₀,h_i,h₂, ...,h_H-1}. WhenH=1, only one modulation index exists,h₀, which is denoted ash.

Pulse Shape Filtering

The CPM method uses pulse shaping to smooth the phase transitions of the modulated signal. The functionq(t) is the phase response obtained from the frequency pulse,g(t), through this relation: $q (t) = \int_{- \infty}^{t} g (t) d t$ .

The specified frequency pulse shape corresponds to these pulse shape expressions forg(t）.

Pulse Shape	Expression
Rectangular	$g (t) = {\begin{matrix} \frac{1}{2 L T}, & 0 \leq t \leq L T \\ 0 & otherwise \end{matrix}$
Raised Cosine	$g (t) = {\begin{matrix} \frac{1}{2 L T} [1 - \cos (\frac{2 π t}{L T})], & 0 \leq t \leq L T \\ 0 & otherwise \end{matrix}$
Spectral Raised Cosine	$g (t) = \frac{1}{L_{main} T} \frac{\sin (\frac{2 π t}{L_{main} T})}{\frac{2 π t}{L_{main} T}} \frac{\cos (β \frac{2 π t}{L_{main} T})}{1 - {(\frac{4 β}{L_{main} T} t)}^{2}}, 0 \leq β \leq 1$
Gaussian	$\begin{matrix} g (t) = \frac{1}{2 T} {Q [2 π B_{b} \frac{t - \frac{T}{2}}{\sqrt{\ln 2}}] - Q [2 π B_{b} \frac{t + \frac{T}{2}}{\sqrt{\ln 2}}]}, where \\ Q (t) = \int_{t}^{\infty} \frac{1}{\sqrt{2 π}} e^{- τ^{2} / 2} d τ \end{matrix}$
Tamed FM (tamed frequency modulation)	$\begin{array}{l} g (t) = \frac{1}{8} [g_{0} (t - T) + 2 g_{0} (t) + g_{0} (t + T)], where \\ g_{0} (t) \approx \frac{1}{T} [\frac{\sin (\frac{π t}{T})}{\frac{π t}{T}} - \frac{π^{2}}{24} \frac{2 \sin (\frac{π t}{T}) - \frac{2 π t}{T} \cos (\frac{π t}{T}) - {(\frac{π t}{T})}^{2} \sin (\frac{π t}{T})}{{(\frac{π t}{T})}^{3}}] \end{array}$

L_mainis the main lobe pulse duration in symbol intervals.
βis the roll-off factor of the spectral raised cosine.
B_bis the product of the bandwidth and the Gaussian pulse.
The duration of the pulse,LT, is the pulse length in symbol intervals. As defined by the expressions, the spectral raised cosine, Gaussian, and tamed FM pulse shapes have infinite length. For all practical purposes,LTspecifies the truncated finite length.

For more information on pulse shape filtering, see[1].

Integer-Valued and Binary-Valued Output Signals

When you set theBitOutputproperty tofalse:

The object outputs an integer column vector of length equal toN/SamplesPerSymbol, whereNis the length of the input signal and indicates the number of input baseband modulated symbols. The output values are odd integers in the range [–(ModulationOrder–1), (ModulationOrder–1)].
You cannot set theOutputDataTypeproperty to'logical'.

When you set theBitOutputproperty totrue:

The object outputs a binary column vector of length equal tok×(N/SamplesPerSymbol), wherek= log2(ModulationOrder)andNis the number of input baseband modulated symbols (specifically, the length of the input signal).
TheSymbolMappingproperty determines how the object maps integers in the range [0,ModulationOrder– 1] tok-length bit word. The binary word mapping options are natural binary-coded ordering or Gray-coded ordering.
You can set theOutputDataTypeproperty to only'double'or'logical'.
The object follows this process.
1. Map each demodulated symbol to an odd integerL在the range [–(ModulationOrder–1), (ModulationOrder–1)].
2. MapLto the nonnegative integer (L+ModulationOrder–1)/2.
3. Map each nonnegative integer to ak-length binary word. The binary word mapping options are natural binary-coded ordering or Gray-coded ordering, as specified by theSymbolMappingproperty.

Traceback Depth and Output Delays

The traceback depth is the number of trellis branches used to construct each traceback path. Traceback depth influences the output delay, which is the number of zero symbols that precede the first meaningful demodulated value in the output.

The optimal traceback depth setting depends on minimum squared Euclidean distance calculations. Alternatively, you can choose a typical value, dependent on the number of states, using thefive-times-the-constraint-lengthrule, which corresponds to5log₂(numStates).

For a binary raised cosine pulse shape with a pulse length of 3 andh=2/3, applying this rule(5log₂(3×2²) = 18gives a result that is close to the optimum value of 20.

参考文献

[1]Anderson, John B., Tor Aulin, and Carl-Erik Sundberg.Digital Phase Modulation. New York: Plenum Press, 1986.

Extended Capabilities

C/C++ Code Generation
Generate C and C++ code using MATLAB® Coder™.

Usage notes and limitations:

SeeSystem Objects in MATLAB Code Generation(MATLAB Coder).

Version History

Introduced in R2012a

comm.CPMDemodulator

Description

Creation

Syntax

Description

Properties

ModulationOrder—Modulation order4(default) |power of two scalar

BitOutput—Option to output data as bits0orfalse(default) |1ortrue

SymbolMapping—Symbol encoding'Binary'(default) |'Gray'

Dependencies

ModulationIndex—Modulation index0.5(default) |nonnegative scalar|column vector

FrequencyPulse—Type of frequency pulse shaping'Rectangular'(default) |'Raised Cosine'|'Spectral Raised Cosine'|'Gaussian'|'Tamed FM'

MainLobeDuration—Main lobe duration1(default) |positive integer

Dependencies

RolloffFactor—Roll-off factor0.2(default) |scalar in the range [0, 1]

Dependencies

BandwidthTimeProduct—Product of bandwidth and symbol time of Gaussian pulse shape0.3(default) |positive scalar

Dependencies

PulseLength—Length of frequency pulse shape1(default) |positive integer

SymbolPrehistory—Symbol prehistory1(default) |scalar|vector

InitialPhaseOffset—Initial phase offset0(default) |scalar

SamplesPerSymbol—Number of samples per input symbol8(default) |positive integer

TracebackDepth—Traceback depth for Viterbi algorithm16(default) |positive integer

OutputDataType—Data type of output'double'(default) |'logical'|'int8'|'int16'|'int32'

Usage

Syntax

Description

Input Arguments

x—CPM-modulated signalcolumn vector

Output Arguments

y— Output signalcolumn vector | matrix

Object Functions

Common to All System Objects

Examples

CPM Modulate and Demodulate Signal with Gray Mapping and Bit Inputs

Apply GFSK Modulation and Demodulation

More About

CPM Demodulation Method

CPM Method

Pulse Shape Filtering

Integer-Valued and Binary-Valued Output Signals

Traceback Depth and Output Delays

参考文献

Extended Capabilities

C/C++ Code GenerationGenerate C and C++ code using MATLAB® Coder™.

Version History

See Also

Objects

Blocks

`ModulationOrder`—Modulation order
4(default) |power of two scalar

`BitOutput`—Option to output data as bits
`0`or`false`(default) |`1`or`true`

`SymbolMapping`—Symbol encoding
`'Binary'`(default) |`'Gray'`

`ModulationIndex`—Modulation index
`0.5`(default) |nonnegative scalar|column vector

`FrequencyPulse`—Type of frequency pulse shaping
`'Rectangular'`(default) |`'Raised Cosine'`|`'Spectral Raised Cosine'`|`'Gaussian'`|`'Tamed FM'`

`MainLobeDuration`—Main lobe duration
`1`(default) |positive integer

`RolloffFactor`—Roll-off factor
`0.2`(default) |scalar in the range [0, 1]

`BandwidthTimeProduct`—Product of bandwidth and symbol time of Gaussian pulse shape
`0.3`(default) |positive scalar

`PulseLength`—Length of frequency pulse shape
`1`(default) |positive integer

`SymbolPrehistory`—Symbol prehistory
`1`(default) |scalar|vector

`InitialPhaseOffset`—Initial phase offset
`0`(default) |scalar

`SamplesPerSymbol`—Number of samples per input symbol
`8`(default) |positive integer

`TracebackDepth`—Traceback depth for Viterbi algorithm
`16`(default) |positive integer

`OutputDataType`—Data type of output
`'double'`(default) |`'logical'`|`'int8'`|`'int16'`|`'int32'`

`x`—CPM-modulated signal
column vector

`y`— Output signal
column vector | matrix

C/C++ Code Generation
Generate C and C++ code using MATLAB® Coder™.