Filter input signal through multipath Rayleigh fading channel
The comm.RayleighChannel
System object™ filters an input signal through a multipath Rayleigh fading channel. For more
information on fading model processing, see Methodology for Simulating
Multipath Fading Channels.
To filter an input signal using a multipath Rayleigh fading channel:
Create the comm.RayleighChannel
object and set its properties.
Call the object with arguments, as if it were a function.
To learn more about how System objects work, see What Are System Objects?.
creates a
frequencyselective or frequencyflat multipath Rayleigh fading channel System
object. This object filters a real or complex input signal through the multipath
channel to obtain the channelimpaired signal.rayleighchan
= comm.RayleighChannel
sets properties using one or more namevalue pair arguments. Enclose each property
name in quotes. For example,
rayleighchan
= comm.RayleighChannel(Name,Value)comm.RayleighChannel('
sets the input signal sample rate to 2.SampleRate
',2)
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.
SampleRate
— Input signal sample rate1
(default)  positive scalarInput signal sample rate in Hz, specified as a positive scalar.
Data Types: double
PathDelays
— Discrete path delay0
(default)  scalar  row vectorDiscrete path delay in seconds, specified as a scalar or row vector.
When PathDelays
is a scalar, the channel is
frequency flat.
When PathDelays
is a vector, the channel is
frequency selective.
Data Types: double
AveragePathGains
— Average of the discrete paths0
(default)  scalar  row vectorAverage gains of the discrete paths in decibels, specified as a scalar or row
vector. AveragePathGains
must be of the same
size as the PathDelays
property.
Data Types: double
NormalizePathGains
— Normalize average path gains to 0 dBtrue
or 1
(default)  false
or 0
Normalize average path gains to 0 dB, specified as a logical
1
(true
) or 0
(false
).
When NormalizePathGains
is
true
, the fading processes are normalized so that
the total power of the path gains, averaged over time, is 0 dB.
When NormalizePathGains
is
false
, the total power of the path gains is not
normalized.
Data Types: logical
MaximumDopplerShift
— Maximum Doppler shift for all channel paths0.001
(default)  nonnegative scalarMaximum Doppler shift for all channel paths, specified as a nonnegative scalar. Units are in hertz.
The maximum Doppler shift limit applies to each channel path. When you set this
property to 0
, the channel remains static for the entire input.
You can use the reset
object function to
generate a new channel realization. The MaximumDopplerShift
property value must be smaller than SampleRate
/10/f_{c} for each
path. f_{c} represents the cutoff frequency
factor of the path. For most Doppler spectrum types, the value of
f_{c} is 1. For Gaussian and biGaussian
Doppler spectrum types, f_{c} is dependent
on the Doppler spectrum structure fields. For more details about how
f_{c} is defined, see Cutoff Frequency Factor.
Data Types: double
DopplerSpectrum
— Doppler spectrum shape for all channel pathsdoppler('Jakes')
(default)  Doppler spectrum structure  1byN_{P} cell array of Doppler
spectrum structuresDoppler spectrum shape for all channel paths, specified as a Doppler spectrum
structure or a 1byN_{P} cell array of
Doppler spectrum structures. These Doppler spectrum structures must be outputs of
the form returned from the doppler
function.
N_{P} is the number of discrete delay
paths, that is, the length of the PathDelays
property.
When DopplerSpectrum
is defined by a single
Doppler spectrum structure, all paths have the same specified Doppler
spectrum.
When DopplerSpectrum
is defined by a cell array
of Doppler spectrum structures, each path has the Doppler spectrum
specified by the corresponding structure in the cell array.
Options for the spectrum type are specified by the
specType
input to the doppler
function. If you set the
FadingTechnique
property to 'Sum of
sinusoids'
, you must set DopplerSpectrum
to
doppler('Jakes')
.
To enable this property, set the MaximumDopplerShift
property to a positive value. The MaximumDopplerShift
property defines the maximum
Doppler shift value permitted when specifying the Doppler spectrum.
Data Types: struct
 cell
FadingTechnique
— Channel model fading technique'Filtered Gaussian noise'
(default)  'Sum of sinusoids'
Channel model fading technique, specified as 'Filtered Gaussian
noise'
or 'Sum of sinusoids'
.
Data Types: char
 string
NumSinusoids
— Number of sinusoids used48
(default)  positive integerNumber of sinusoids used to model the fading process, specified as a positive integer.
To enable this property, set the FadingTechnique
property to 'Sum of sinusoids'
.
Data Types: double
InitialTimeSource
— Source to control start time of fading process'Property'
(default)  'Input port'
Source to control the start time of the fading process, specified as
'Property'
or 'Input port'
.
When InitialTimeSource
is set to
'Property'
, use the
InitialTime
property to set the initial time
offset.
When InitialTimeSource
is set to 'Input
port'
, specify the start time of the fading process by using
the itime
input to the System object. The input value can change in consecutive calls to the
System object.
To enable this property, set the FadingTechnique
property to 'Sum of sinusoids'
.
Data Types: string
 char
InitialTime
— Initial time offset0
(default)  nonnegative scalarInitial time offset, in seconds, for the fading model, specified as a
nonnegative scalar. InitialTime
must be greater than end time
of the last frame. InitialTime
is rounded up to the nearest
sample position when it is not a multiple of
1/SampleRate
.
To enable this property, set the FadingTechnique
property to 'Sum of sinusoids'
and the
InitialTimeSource
property to
'Property'
.
Data Types: double
RandomStream
— Source of random number stream'Global stream'
(default)  'mt19937ar with seed'
Source of random number stream, specified as 'Global stream'
or 'mt19937ar with seed'
.
'Global stream'
— The current global random number
stream is used for normally distributed random number generation. In this
case, the reset
object function
resets the channel filters only.
'mt19937ar with seed'
— The mt19937ar algorithm is
used for normally distributed random number generation. In this case, the
reset
object function resets the channel filters
and reinitializes the random number stream to the value of the
Seed
property.
Data Types: string
 char
Seed
— Initial seed of mt19937ar random number stream73
(default)  nonnegative integerInitial seed of mt19937ar random number stream generator algorithm, specified
as a nonnegative integer. When the reset
object function is
called, the mt19937ar random number stream is reinitialized to the
Seed
property value.
To enable this property, set the RandomStream
property
to 'mt19937ar with seed'
.
Data Types: double
PathGainsOutputPort
— Option to output path gainsfalse
or 0
(default)  true
or 1
Option to output path gains, specified as 0
(false
) or 1
(true
).
Set this property to 1
(true
) to output the
channel path gains of the underlying fading process.
Data Types: logical
Visualization
— Channel visualization'Off'
(default)  'Impulse response'
 'Frequency response'
 'Impulse and frequency responses'
 'Doppler spectrum'
Channel visualization, specified as 'Off'
, 'Impulse
response'
, 'Frequency response'
, 'Impulse
and frequency responses'
, or 'Doppler spectrum'
.
For more information, see Channel Visualization.
Data Types: string
 char
SamplesToDisplay
— Percentage of samples to display'25%'
(default)  '10%'
 '50%'
 '100%'
Percentage of samples to display, specified as '25%'
,
'10%'
, '50%'
, or
'100%'
. Displaying fewer samples improves (decreases) the
display update rate at the expense of decreasing the visualized precision.
To enable this property, set the Visualization
property
to 'Impulse response'
, 'Frequency
response'
, or 'Impulse and frequency
responses'
.
Data Types: string
 char
PathsForDopplerDisplay
— Path for Doppler display1
(default)  positive integer Path for Doppler display, specified as a positive integer. Use this property to
select the discrete path used in constructing a Doppler spectrum plot. The
specified path must be an element of {1, 2, ...,
N_{P}}. In this set, N_{P} is
the number of discrete delay paths, that is, the length of the PathDelays
property value.
To enable this property, set the Visualization
property
to 'Doppler spectrum'
.
Data Types: double
x
— Input signalInput signal, specified as an N_{S}by1 vector, where N_{S} is the number of samples.
Data Types: single
 double
Complex Number Support: Yes
itime
— Initial time0
 nonnegative scalarInitial time in seconds, specified as a nonnegative scalar. The
itime
input must be greater than the end time of the
last frame. When itime
is not a multiple of
1/SampleRate
, it is rounded up to the nearest sample
position.
Data Types: single
 double
y
— Output signalOutput signal, returned as an
N_{S}by1 vector of complex values with
the same data precision as the input signal x
.
N_{S} is the number of
samples.
pathgains
— Path gainsPath gains, returned as an
N_{S}byN_{P}
array. N_{S} is the number of samples.
N_{P} is the number of discrete delay
paths, that is, the length of the PathDelays
property value. pathgains
contains complex values with the same precision as the input signal
x
.
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)
This example shows how to produce the same multipath Rayleigh fading channel response by using two different methods for random number generation. The multipath Rayleigh fading channel System object™ includes two methods for random number generation. You can use the current global stream or the mt19937ar algorithm with a specified seed. By interacting with the global stream, the System object can produce the same outputs from these two methods.
Create a PSK modulator System object to modulate randomly generated data.
pskModulator = comm.PSKModulator; channelInput = pskModulator(randi([0 pskModulator.ModulationOrder1],1024,1));
Create a multipath Rayleigh fading channel System object, specifying the random number generation method as the my19937ar algorithm and the random number seed as 22.
rayleighchan = comm.RayleighChannel(... 'SampleRate',10e3, ... 'PathDelays',[0 1.5e4], ... 'AveragePathGains',[2 3], ... 'NormalizePathGains',true, ... 'MaximumDopplerShift',30, ... 'DopplerSpectrum',{doppler('Gaussian',0.6),doppler('Flat')}, ... 'RandomStream','mt19937ar with seed', ... 'Seed',22, ... 'PathGainsOutputPort',true);
Filter the modulated data by using the multipath Rayleigh fading channel System object.
[chanOut1,pathGains1] = rayleighchan(channelInput);
Set the System object to use the global stream for random number generation.
release(rayleighchan);
rayleighchan.RandomStream = 'Global stream';
Set the global stream to have the same seed that was specified when creating the multipath Rayleigh fading channel System object.
rng(22)
Filter the modulated data by using the multipath Rayleigh fading channel System object again.
[chanOut2,pathGains2] = rayleighchan(channelInput);
Verify that the channel and path gain outputs are the same for each of the two methods.
isequal(chanOut1,chanOut2)
ans = logical
1
isequal(pathGains1,pathGains2)
ans = logical
1
This example shows how to create a frequencyselective multipath Rayleigh fading channel System object and display its impulse and frequency responses.
Set the sample rate to 3.84 MHz. Specify path delays and gains by using the ITU pedestrian B channel configuration. Set the maximum Doppler shift to 50 Hz.
fs = 3.84e6; % Hz pathDelays = [0 200 800 1200 2300 3700]*1e9; % sec avgPathGains = [0 0.9 4.9 8 7.8 23.9]; % dB fD = 50; % Hz
Create a multipath Rayleigh fading channel System object by using the previously defined properties and to visualize the impulse response and frequency response plots.
rayleighchan = comm.RayleighChannel('SampleRate',fs, ... 'PathDelays',pathDelays, ... 'AveragePathGains',avgPathGains, ... 'MaximumDopplerShift',fD, ... 'Visualization','Impulse and frequency responses');
Generate random binary data and pass it through the multipath Rayleigh fading channel. The impulse response plot enables you to identify the individual paths and their corresponding filter coefficients. The frequency response plot shows the frequencyselective nature of the ITU pedestrian B channel.
x = randi([0 1],1000,1); y = rayleighchan(x);
The example shows how the channel state is maintained in cases in which data is discontinuously transmitted. Create a multipath Rayleigh fading channel System object and pass data through it using the sumofsinusoids technique.
Set the channel properties.
fs = 1000; % Sample rate (Hz) pathDelays = [0 2.5e3]; % Path delays (s) pathPower = [0 6]; % Path power (dB) fD = 5; % Maximum Doppler shift (Hz) ns = 1000; % Number of samples nsdel = 100; % Number of samples for delayed paths
Define the overall simulation time and three time segments for which data is to be transmitted. In this case, the channel is simulated for 1s with a 1000 Hz sampling rate. One 1000sample continuous data sequence is transmitted at time 0. Three 100sample data packets are transmitted at times 0.1 s, 0.4 s, and 0.7 s, respectively.
to0 = 0.0; to1 = 0.1; to2 = 0.4; to3 = 0.7; t0 = (to0:ns1)/fs; % Transmission 0 t1 = to1+(0:nsdel1)/fs; % Transmission 1 t2 = to2+(0:nsdel1)/fs; % Transmission 2 t3 = to3+(0:nsdel1)/fs; % Transmission 3
Generate random binary data corresponding to the previously defined time intervals.
d0 = randi([0 1],ns,1); d1 = randi([0 1],nsdel,1); d2 = randi([0 1],nsdel,1); d3 = randi([0 1],nsdel,1);
Create a frequencyflat multipath Rayleigh fading System object, specifying the sumofsinusoids fading technique. So that results can be repeated, specify a seed value. Use the default InitialTime
property setting so that the fading channel will be simulated from time 0. Enable the path gains output.
rayleighchan1 = comm.RayleighChannel('SampleRate',fs, ... 'MaximumDopplerShift',5, ... 'RandomStream','mt19937ar with seed', ... 'Seed',17, ... 'FadingTechnique','Sum of sinusoids', ... 'PathGainsOutputPort',true);
Create a clone of the multipath Rayleigh fading channel System object. Set the source for initial time so that the fading channel offset time can be specified as an input argument when using the System object.
rayleighchan2 = clone(rayleighchan1);
rayleighchan2.InitialTimeSource = 'Input port';
Pass random binary data through the first multipath Rayleigh fading channel System object, rayleighchan1
. Data is transmitted over all 1000 time samples. For this example, only the complex path gain is needed.
[~,pg0] = rayleighchan1(d0);
Pass random data through the second multipath Rayleigh fading channel System object, rayleighchan2
, where the initial time offsets are provided as input arguments.
[~,pg1] = rayleighchan2(d1,to1); [~,pg2] = rayleighchan2(d2,to2); [~,pg3] = rayleighchan2(d3,to3);
Compare the number of samples processed by the two channels by using the info
object function. The rayleighchan1
object processed 1000 samples, while the rayleighchan2
object only processed 300 samples.
G = info(rayleighchan1); H = info(rayleighchan2); [G.NumSamplesProcessed H.NumSamplesProcessed]
ans = 1×2
1000 300
Convert the path gains into decibels.
pathGain0 = 20*log10(abs(pg0)); pathGain1 = 20*log10(abs(pg1)); pathGain2 = 20*log10(abs(pg2)); pathGain3 = 20*log10(abs(pg3));
Plot the path gains for the continuous and discontinuous cases. The gains for the three segments match the gain for the continuous case. The alignment of the two highlights that the sumofsinusoids technique is suited to the simulation of packetized data, as the channel characteristics are maintained even when data is not transmitted.
plot(t0,pathGain0,'r') hold on plot(t1,pathGain1,'b') plot(t2,pathGain2,'b') plot(t3,pathGain3,'b') grid xlabel('Time (sec)') ylabel('Path Gain (dB)') legend('Continuous','Discontinuous','location','nw') title('Continuous and Discontinuous Transmission Path Gains')
This example shows how to use the ChannelFilterCoefficients
property returned by the info
object function of the comm.RayleighChannel
System object to reproduce the multipath Rayleigh fading channel output across multiple frames.
Create a multipath Rayleigh fading channel System object, defining two paths and specifying data to pass through the channel.
rayleighchan = comm.RayleighChannel('SampleRate',1000,'PathDelays',[0 1.5e3], ... 'AveragePathGains',[0 3],'PathGainsOutputPort',true)
rayleighchan = comm.RayleighChannel with properties: SampleRate: 1000 PathDelays: [0 0.0015] AveragePathGains: [0 3] NormalizePathGains: true MaximumDopplerShift: 1.0000e03 DopplerSpectrum: [1x1 struct] Show all properties
data = randi([0 1],600,1);
Pass data through the channel. Assign the ChannelFilterCoefficients
property value to the variable coeff
. Within a for loop, calculate the fractional delayed input signal at the path delay locations stored in coeff
, apply the path gains, and sum the results for all paths. Compare the output of the multipath Rayleigh fading channel System object (chanout1
) to the output reproduced using the path gains and the ChannelFilterCoefficients
property of the multipath Rayleigh fading channel System object (chanout2
).
chaninfo = info(rayleighchan); coeff = chaninfo.ChannelFilterCoefficients; Np = length(rayleighchan.PathDelays); state = zeros(size(coeff,2)1, size(coeff, 1)); nFrames = 10; chkChan = zeros(nFrames,1); for jj = 1 : nFrames data = randi([0 1],600,1); [chanout1,pg] = rayleighchan(data); fracdelaydata = zeros(size(data,1),Np); % Calculate the fractional delayed input signal. for ii = 1:Np [fracdelaydata(:,ii),state(:,ii)] = ... filter(coeff(ii,:),1,data,state(:,ii)); end % Apply the path gains and sum the results for all paths. Compare the % channel outputs, chanout1 and chanout2. chanout2 = sum(pg .* fracdelaydata,2); chkChan(jj) = isequal(chanout1,chanout2); end chkChan'
ans = 1×10
1 1 1 1 1 1 1 1 1 1
Verify that the autocorrelation of the path gain output from the Rayleigh channel System object is a Bessel function. Based on the results in [1] and Appendix A of [2], we know that when the autocorrelation of the path gain outputs is a Bessel function, and the Doppler spectrum is Jakesshaped.
Initialize simulation parameters.
Rsym = 9600; % Input signal symbol rate (symbols/sec) sps = 10; % Number of samples per input symbol Fs = sps*Rsym; % Input signal sampling frequency (samples/sec) Ts = 1/Fs; % Input signal sampling period (sec) numsym = 1e6; % Total number of input symbols to simulate numsamp = numsym*sps; % Total number of channel samples to simulate fd = 100; % Maximum Doppler frequency shift (Hz) num_acsamp = 5000; % Number of samples of autocovariance of complex fading process calculated numtx = 1; % Number of transmit antennas numrx = 1; % Number of receive antennas numsin = 48; % Number of sinusoids frameLen = 10000; numFrames = numsamp/frameLen;
Configure a Rayleigh channel System object.
chan = comm.RayleighChannel('FadingTechnique','Sum of sinusoids', ... 'NumSinusoids',numsin,'RandomStream','mt19937ar with seed', ... 'PathDelays',0,'AveragePathGains',0,'SampleRate',Fs, ... 'MaximumDopplerShift',fd,'PathGainsOutputPort',true);
Apply DPSK modulation to a random bit stream.
tx = randi([0 1],numsamp,numtx); % Random bit stream dpskSig = dpskmod(tx,2); % DPSK signal
Pass modulated signal through the channel.
outsig = zeros(numsamp,numrx); pg_rx = zeros(numsamp,numrx,numtx); for frmNum = 1:numFrames [outsig((1:frameLen)+(frmNum1)*frameLen,:),pathGains] = chan(dpskSig((1:frameLen)+(frmNum1)*frameLen,:)); for i = 1:numrx pg_rx((1:frameLen)+(frmNum1)*frameLen,i,:) = pathGains(:,:,:,i); end end
Using the channel path gains received per antenna, compute the autocovariance of the fading process for each transmitreceive path.
autocov = zeros(frameLen+1,numrx,numtx); autocov_normalized_real = zeros(num_acsamp+1,numrx,numtx); autocov_normalized_imag = zeros(num_acsamp+1,numrx,numtx); for i = 1:numrx % Compute autocovariance of simulated complex fading process for j = 1:numtx autocov(:,i,j) = xcov(pg_rx(:,i,j),num_acsamp); % Real part of normalized autocovariance autocov_normalized_real(:,i,j) = real(autocov(num_acsamp+1:end,i,j)/autocov(num_acsamp+1,i,j)); % Imaginary part of normalized autocovariance autocov_normalized_imag(:,i,j) = imag(autocov(num_acsamp+1:end,i,j)/autocov(num_acsamp+1,i,j)); end end
Using the besselj
function compute the theoretical autocovariance of complex fading process.
Rrr = zeros(1,num_acsamp+1); for n = 1:1:num_acsamp+1 Rrr(n) = besselj(0,2*pi*fd*(n1)*Ts); end Rrr_normalized = Rrr/Rrr(1);
Display the autocovariance comparing the results from the Rayleigh channel System object and the besselj function.
subplot(2,1,1) plot(autocov_normalized_real,'b') hold on plot(Rrr_normalized,'r') hold off legend('Using comm.RayleighChannel','Using Bessel function of the first kind') title('Autocovariance of real part of Rayleigh process') subplot(2,1,2) plot(autocov_normalized_imag) legend('Using comm.RayleighChannel') title('Autocovariance of imaginary part of Rayleigh process')
Compute mean square error.
act_mse_real = sum((autocov_normalized_realrepmat(Rrr_normalized.',1,numrx,numtx)).^2,1)/size(autocov_normalized_real,1)
act_mse_real = 7.0043e08
act_mse_imag = sum((autocov_normalized_imag0).^2,1)/size(autocov_normalized_imag,1)
act_mse_imag = 4.1064e07
References
1. Dent, P., G.E. Bottomley, and T. Croft. “Jakes Fading Model Revisited.” Electronics Letters 29, no. 13 (1993): 1162. https://doi.org/10.1049/el:19930777.
2. Pätzold, Matthias. Mobile Fading Channels. Chichester, UK: John Wiley & Sons, Ltd, 2002. https://doi.org/10.1002/0470847808.
The cutoff frequency factor, f_{c}, is determined for different Doppler spectrum types.
For any Doppler spectrum type other than Gaussian and biGaussian, f_{c} equals 1.
For a doppler
('Gaussian')
spectrum type,
f_{c} equals NormalizedStandardDeviation
$$\text{\hspace{0.17em}}\times \text{\hspace{0.17em}}\sqrt{2\mathrm{log}2}$$.
For a doppler
('BiGaussian')
spectrum type:
If the PowerGains
(1)
and
NormalizedCenterFrequencies
(2)
field values are both 0
, then
f_{c} equals NormalizedStandardDeviation
(1)
$$\text{\hspace{0.17em}}\times \text{\hspace{0.17em}}\sqrt{2\mathrm{log}2}$$.
If the PowerGains
(2)
and
NormalizedCenterFrequencies
(1)
field values are both 0
, then
f_{c} equals NormalizedStandardDeviation
(2)
$$\text{\hspace{0.17em}}\times \text{\hspace{0.17em}}\sqrt{2\mathrm{log}2}$$.
If the NormalizedCenterFrequencies
field value is
[0,0]
and the
NormalizedStandardDeviation
field has two identical
elements, then f_{c} equals NormalizedStandardDeviation
(1)
$$\text{\hspace{0.17em}}\times \text{\hspace{0.17em}}\sqrt{2\mathrm{log}2}$$.
In all other cases, f_{c} equals 1.
The comm.RayleighChannel
System object enables visualization of the channel impulse response, frequency response,
and Doppler spectrum.
Note
The displayed and specified path gain locations can differ by as much as 5% of the input sample time.
The visualization display speed is controlled by the combination of
the SamplesToDisplay
property and the
Playback > Reduce Updates to Improve
Performance scope menu item. Reducing the percentage
of samples to display and enabling reduced updates can speed up the
rendering of the visualization scope.
After the visualization scopes are manually closed, calls to the System object are executed at its normal speed.
Code generation is available only when the
Visualization
property is set to
'Off'
.
Impulse Response
The impulse response scope displays the path gains, channel filter coefficients, and
interpolated path gains. The path gains occur at time instances which correspond to the
specified PathDelays
property and might not be aligned with the
input sampling time. The channel filter coefficients are used to model the channel. They
are interpolated from the actual path gains and are aligned with the input sampling
time. In cases in which the path gains are aligned with the sampling time, the path
gains overlap the filter coefficients. Sinc interpolation is used to connect the channel
filter coefficients and is shown as the interpolated path gains. These points are used
solely for display purposes and not used in subsequent channel filtering. For a flat
fading channel (one path), the sinc interpolation curve is not displayed. For all
impulse response plots, the frame and sample numbers are shown in the upperleft corner
of the scope.
The impulse response plot shares the same toolbar and menus as the dsp.ArrayPlot
System object.
When the specified path gains align with the sample rate, the path gains and the channel filter coefficients overlap. In this figure the impulse response shows the path gains overlap the filter coefficients.
If the specified path gains are not aligned with the sample rate, the path gains and the channel filter coefficients do not overlap. In this figure, the filter coefficients are equally distributed. The path gains do not overlap with the channel filter coefficients because the path gains are not aligned with the sample rate.
This figure shows the impulse response for a frequencyflat channel. In this case, the interpolated path gains are not displayed.
Frequency response
The frequency response scope displays the multipath Rayleigh fading channel spectrum
by taking a discrete Fourier transform of the channel filter coefficients. The frequency
response plot shares the same toolbar and menus as the dsp.SpectrumAnalyzer
System object.
The yaxis limits of the plot are computed based on the
NormalizePathGains
and AveragePathGains
properties of the comm.RayleighChannel
System object.
This table shows other selected default spectrum settings. These settings can be changed from their default values by using the View > Spectrum Settings menu.
Spectrum Settings  Value 

Main options>Type 

Main options>Window length  Channel filter length 
Main options>NFFT 

Window options>Window 

Trace options>Units 

This figure shows the frequency response plot for a frequency selective channel.
Doppler Spectrum
The Doppler spectrum plot displays the theoretical Doppler spectrum and the
empirically determined data points. The theoretical data is displayed as a line for the
case of nonstatic channels and as a point for static channels. The empirical data is
shown as * symbols. Before the empirical plot is updated, the internal buffer must be
completely filled with filtered Gaussian samples. The empirical plot is the running mean
of the spectrum calculated from each full buffer. For nonstatic channels, the number of
input samples needed before the next update is displayed in the upperleft corner of the
plot. The number of samples needed is a function of the sample rate and the maximum
Doppler shift. For static channels, the text "Reset fading channel for next
update
" is displayed.
[1] Oestges, Claude, and Bruno Clerckx. MIMO Wireless Communications: From RealWorld Propagation to SpaceTime Code Design. 1st ed. Boston, MA: Elsevier, 2007.
[2] Correia, Luis M., and European Cooperation in the Field of Scientific and Technical Research (Organization), eds. Mobile Broadband Multimedia Networks: Techniques, Models and Tools for 4G. 1st ed. Amsterdam ; Boston: Elsevier/Academic Press, 2006.
[3] Kermoal, J.P., L. Schumacher, K.I. Pedersen, P.E. Mogensen, and F. Frederiksen. “A Stochastic MIMO Radio Channel Model with Experimental Validation.” IEEE Journal on Selected Areas in Communications 20, no. 6 (August 2002): 1211–26. https://doi.org/10.1109/JSAC.2002.801223.
[4] Jeruchim, Michel C., Philip Balaban, and K. Sam Shanmugan. Simulation of Communication Systems. Second edition. Boston, MA: Springer US, 2000.
[5] Patzold, M., ChengXiang Wang, and B. Hogstad. “Two New SumofSinusoidsBased Methods for the Efficient Generation of Multiple Uncorrelated Rayleigh Fading Waveforms.” IEEE Transactions on Wireless Communications 8, no. 6 (June 2009): 3122–31. https://doi.org/10.1109/TWC.2009.080769.
Usage notes and limitations:
To generate C code, set the DopplerSpectrum
property to a
single Doppler spectrum structure.
Code generation is available only when the Visualization
property is 'Off'
.
See System Objects in MATLAB Code Generation (MATLAB Coder).
You have a modified version of this example. Do you want to open this example with your edits?
You clicked a link that corresponds to this MATLAB command:
Run the command by entering it in the MATLAB Command Window. Web browsers do not support MATLAB commands.
Choose a web site to get translated content where available and see local events and offers. Based on your location, we recommend that you select: .
Select web siteYou can also select a web site from the following list:
Select the China site (in Chinese or English) for best site performance. Other MathWorks country sites are not optimized for visits from your location.