A rather large number of SDF demonstrations have been developed. These can serve as valuable illustrations of the possibilities. Almost every star is illustrated in the demos. Because of the large number, the demos are organized into a set of palettes. Certain demos may appear in more than one palette. A top-level palette, shown in figure

`butterfly`

Use sines and cosines to compute a curve known as the butterfly curve, invented by T. Fay. The curve is plotted in polar form.`chaoticNoise`

Chaotic Markov map example with a nonlinear feedback loop.`comparison`

Compare two sinusoidal signals using the`Test`

star.`complexExponential`

Generate and plot a complex exponential.`delayTest`

Illustrates the use of initializable delays.`lmsFreqDetect`

Illustrate the use of the LMS algorithm to estimate the dominant sinusoidal frequency in the input signal.`freqPhaseOffset`

Impose frequency jitter and phase offset on a sinusoid using the`freqPhase`

SDF block.`gaussian`

Generate a Gaussian white noise signal, and plot its histogram and estimated autocorrelation.`integrator`

Demonstrate the features of the integrator star, such as limiting, leakage, and resetting.`Modulo`

Demonstrate modulus computation for float and integer data types.`muxDeMux`

Demonstrate the`Mux`

and`DeMux`

stars, which perform multiplexing and demultiplexing. Contrast with the`scramble`

demo below.`quantize`

Demonstrate the use of the`Quantizer`

star.`scramble`

This system rearranges the order of samples of signal using the`Commutator`

and`Distributor`

stars. Note that because these are multirate stars, one iteration involves more than one sample. Contrast with the`muxDeMux`

demo above.`sinMod`

Modulate a sinusoid by multiplying by another sinusoid.`tbus`

Illustrate the bus facility in Ptolemy, in which multiple signals are combined onto a single graphical connection.

`analytic`

Use a`FIRCx`

star filter to reduce the sample rate of a sinusoid by a factor of 8/5, and at the same time produce a complex approximately analytic signal (one that has no negative frequency components).`broken`

Give an example of an inconsistent SDF system. It fails to run, generating an error message instead.`downSample`

Convert from the digital audio tape sampling rate (48 kHz) to the compact disc sampling rate (44.1 kHz). The conversion is performed in multiple stages for better performance.`filterBank`

Implement an eight-level perfect reconstruction one-dimensional filter bank based on the biorthogonal wavelet decomposition.`filterBank-NonUniform`

Implement a simple split of the frequency domain into two non-uniform frequency bands.`interp`

Use an FIR filter to upsample by a factor of 8 and linearly interpolate between samples.`multirate`

Upsample a sinusoidal signal by a ratio of 5/2 using a polyphase lowpass interpolating FIR filter.`upSample`

Convert from the compact disc sampling rate (44.1 kHz) to the digital audio tape sampling rate (48 kHz). The conversion is performed in multiple stages for better performance.

`constellation`

A 16-QAM signal is sent through a baseband equivalent channel that simulates the following impairments: frequency offset, phase jitter and white Gaussian noise.`DTMFCodec`

Dual-Tone Modulated Frequency Demo. Generate touch tones and decode the based on the Goertzel Algorithm.`eye`

Plot an eye diagram for a binary antipodal signal with a raised-cosine pulse shape and user controlled noise.`lmsDTMFCodec`

Dual-Tone Modulated Frequency Demo. Generate touch tones

and decode them based on the LMS Algorithm.`lossySpeech`

Illustrate the effect on speech of a zero-substitution policy in a network (such as ATM) with 48 byte packets and a variable loss probability. Note that this demo requires audio capability and will probably only work on Sun workstations.`lossySpeechPrevCell`

Illustrate the effect on speech of a previous cell substitution policy in a network (such as ATM) with 48 byte packets and a variable loss probability. Note that this demo requires audio capability and will probably only work on Sun workstations.`modem`

Baseband model of a 16-QAM modem.`pseudoRandom`

Generate a pseudo-random sequence of zeros and ones using a maximal-length shift register and test its randomness by estimating it autocorrelation.`pulses`

Generate raised cosine and square-root raised cosine pulses and demonstrate matched filtering with the square-root raised cosine pulse.`xmitber`

Bit Error determination through simulation at various noise levels.`xmit2rec`

Simple 2-level PAM communication system (matched filtering at the receiver).`xmit4rec`

Simple 4-level PAM communication system (no filtering at the receiver).

`qam`

Produce a 16-point quadrature amplitude modulated (QAM) signal and displays the eye diagram for the in-phase part, the constellation, and the modulated transmited signal.`QAM4withDFE`

This is a model of a digital communication system that uses quadrature amplitude modulation (QAM) and a fractionally spaced decision feedback equalizer.`codeDecode`

Encode and decode a 16-QAM signal using differential encoding for the quadrant and Gray coding for the point within the quadrant.`plldemo`

Simulate a fourth-power optical phase-locked loop with laser phase noise and additive Gaussian white noise operating on a complex baseband envelope model of the signal.`telephoneChannelTest`

Assuming a sampling rate of 8 kHz, a sinusoid at 500 Hz is transmitted through a simulation of a telephone channel with additive Gaussian noise, nonlinear distortion, and phase jitter.

`adaptFilter`

An LMS adaptive filter converges so that its transfer function matches that of a fixed FIR filter.`allPole`

Two realizations of an all-pole filter are shown to be equivalent. One uses an FIR filter in a feedback path, the other uses the`BlockAllPole`

star.`animatedLMS`

An LMS adaptive filter is configured as in the`adaptFilter`

demo, but this time the filter taps are displayed as they adapt.`animatedLMSCx`

A complex LMS adaptive filter is configured as in the`adaptFilter`

demo, but in addition, user-controlled noise is added to the feedback loop using an on-screen slider to control the amount of noise. The filter taps are displayed as they adapt.`cep`

Given the coefficients of any polynomial, this demo uses the cepstrum to find a minimum-phase polynomial. Thus, given the coefficients of the denominator polynomial of an unstable filter, this demo will compute the coefficients of a stable denominator polynomial that has the same magnitude frequency response.`chaos`

This is a simple demonstration of chaos, in which the phase-space plot of the famous Henon map is given.`convolve`

Convolve two rectangular pulses in order to demonstrate the`Convolve`

star.`dft`

Compute a discrete Fourier transform of a finite signal using the`FFT`

star. The magnitude and phase (unwrapped) are plotted.`doppler`

A sine wave is subjected to four successive amounts of doppler shift. The doppler shift is accomplished by the`phaseShift`

galaxy, which forms an analytic signal (using a Hilbert transform) that modulates a complex exponential.`dtft`

Demonstrate the`DTFT`

star, showing how it is different from the`FFTCx`

star. Specifically, the range, number, and spacing of frequency samples is arbitrary.`freqsample`

This system designs FIR filters using the frequency sampling method. Samples of the frequency response are converted into FIR filter coefficients.`iirDemo`

Two equivalent implementations of IIR filtering.`lattice`

Demonstrate the use of lattice filters to synthesize an auto-regressive (AR) random process.`latticeDesign`

Use of Levinson-Durbin algorithm to design a lattice filter with a specified transfer function.`levinsonDurbin`

Use the Levinson-Durbin algorithm to estimate the parameters of an AR process.`linearPrediction`

Perform linear prediction on a test signal consisting of three sinusoids in colored, Gaussian noise. Two mechanisms (Burg's algorithm and an LMS adaptive filter) for linear prediction are compared.`overlapAddFFT`

Convolution is implemented in the frequency domain using overlap and add.`phasedArray`

Simulate a plane wave approaching a phased array with four sensors. The plane wave approaches from angles starting from head on and slowly rotating 360 degrees. The response of the antenna is plotted as a function of direction of arrival in polar form.`powerSpectrum`

Compare three methods for estimating a power spectrum of a signal with three sinusoids plus colored noise. The three methods are the periodogram method, the autocorrelation method, and Burg's method.`timeVarSpec`

A time-varying spectrum is computed using the autocorrelation method and displayed using a waterfall plot.`window`

Generate and display four window functions and the magnitude of their Fourier transforms. The windows displayed are the Hanning, Hamming, Blackman, and steep Blackman.

`ptplay`

is in your path, and that it accepts data of an appropriate format and will play it over a workstation speaker at an 8 kHz sample rate. If you are using a Sun SPARCStation, these conditions will most likely be satisfied, if your path is correct. The `ptplay`

program has also been used on SGI Indigos and HP 700s and 800s. If you are on an HP, you may need other publicly available software.The samples are written into a file before they are played. Since a large number of samples must be generated, these demos can take some time to run. By contrast, the CGC domain has some audio demos that generate sounds in real time at 44.1kHz, assuming a reasonably fast workstation. For further information about playing audio files, see
"Sounds" on page 2-38.

`chirpplay`

Chirp generator that plays on the workstation speaker.`fmplay`

Sound generator using FM modulation that plays on the workstation speaker.`speech`

Read a speech signal from a file, and encode it at two bits per sample using adaptive differential pulse code modulation with a feedback-around-quantizer structure. The signal is then reconstructed from the quantized data. The original and reconstructed speech are played over the workstation speaker.`KSchord`

Simulation of plucked string sounds using the Karplus-Strong algorithm.`vox`

Coarticulation with an Adaptive Vocoder. The resulting FM synthesized sound is played over the workstation speaker.`blockVox`

A block processed version of the vox demo.`lossySpeech`

Illustrate the effect on speech of a zero-substitution policy in a network (such as ATM) with 48 byte packets and a variable loss probability. This demo also appears in the basic demos palette`lossySpeechPrevCell`

Illustrate the effect on speech of a previous cell substitution policy in a network (such as ATM) with 48 byte packets and a variable loss probability. This demo also appears in the basic demos palette.`perfectReconstuction`

Eight-channel perfect reconstruction one-dimensional analysis/synthesis filterbank. The incoming speech signal is split into eight adjacent frequency bins and then reconstructed. The original and reconstructed speech are played over the workstation speaker.`subbandcoding`

Four channel subband speech coding with APCM at 16kps.

`BlendImage`

Combine two images and display the result.`bwDither`

Demonstrate four different forms of black and white dithering: error diffusion, clustered dither, dispersed dither, and use custom mask.`cntrastEnhance`

Contrast enhancement by histogram modification.`ColorImage`

Convert an RGB (red-green-blue) format color image to YUV (luminance-hue-saturation) format and back, and then display it on the workstation screen.`CompareMedian`

Median filter an image to reduce artifacts due to interleaved scanning in video sequences.`DctImage`

Perform discrete cosine transform (DCT) coding of an image sequence.`DpcmImage`

Perform differential pulse code modulation (DPCM) on an image sequence.`EdgeDetect`

Demonstrate four different forms of edge detection: Sobel, Roberts, Prewitt, and Frei-Chen.`MC_DCT`

Perform motion compensation and DCT encoding of video.`MotionComp`

Perform motion compensation video coding.

- The Vec_Quan icon in the image processing palette brings up a sub-palette that has several vector quantization demonstrations in figure
5-28
:

`fullVQCodebk`

Generate a codebook for full search vector quantization.`fullVQ`

Full search vector quantization using codebook generated by`fullVQCodebk`

.`SGVQCodebk`

Generate codebooks for shape-gain vector quantization.`SGVQ`

Shape-gain vector quantization using codebook from`SGVQCodebk`

.`MRVQCodeBk`

Generate codebooks for mean-removed vector quantization using independent quantizer structure.`MRVQmeanCB`

Generate codebook for mean-removed vector quantization.`MRVQshapeCB`

Generate the shape codebook for mean-removed quantization using alternate structure. This universe uses the codebook generated by`MRVQmeanCB`

.`MRVQ`

Mean-removed vector quantization.

`fixConversion`

Illustration of the different masking options available.`fixFIR`

Effect of filter tap precision on the frequency response.`fixIIRdf`

Comparison of a fourth-order direct-form IIR filter implemented with floating-point arithmetic and a similar filter implemented with fixed-point arithmetic.`fixMpyTest`

Testing of fixed-point multiplication over a range of numbers by comparison against floating-point multiplication. The results should be the same.

`animatedLMS`

See "Digital signal processing demos" on page 5-55.`animatedLMSCx`

See "Digital signal processing demos" on page 5-55.`buttons`

Demonstrate`TkButtons`

.`phased_Array`

Demonstrate`TkSlider`

by creating a vertical array of radar sensors that can be move in the horizontal plane. Note that small movements of the sensors radically change the polar gain plot. This simulation demonstrates the importance of sensor calibration to performance of the sensor array.`sinWaves`

Demonstrate`TkBarGraph`

by generating and displaying a complex exponential.`tclScript`

Demonstrate`TclScript`

by generating two interactive X window follies that consist of circles that move in the same playing field.`tkMeter`

Demonstrate`TkMeter`

by creating three bar meters. The first oscillates sinusoidally. The second displays a random number between zero and one. The third displays a random walk.`tkShowValues`

Demonstrate`TkShowValues`

and`TkText`

by displaying the ASCII form of two ramp sequences.`xyplot`

Demonstrate the dynamic plotting capabilities of the`xyplot`

star.

`MatrixTest1`

Demonstrate the use of the Matrix stars that have one input. These include the operations inverse, transpose, and multiply by a scalar gain for all matrix types. Also conjugate and Hermitian transpose are available for the complex matrix type.`MatrixTest2`

Demonstrate the use of some simple Matrix stars with two inputs. These include multiply, add, and subtract.`MatrixTest3`

Demonstrate the use of the Matrix conversion stars. These convert between the scalar particles and the matrix particles as well as between the various matrix types.`initDelays`

Illustrate the use of initializable delays with the`matrix`

class.`Kalman_M`

Compare the convergence properties of a Kalman filter to those of an LMS filter when addressing the problem of adaptive equalization of a process in noise.`SVD_MUSIC_1`

Show the use of singular-value decomposition (SVD) and the Multiple-Signal Characterization (MUSIC) algorithm to identify the frequency of a single sinusoid in a signal that has two different signal to noise ratios.`SVD_MUSIC_2`

Demonstrate the use of the Multiple-Signal Characterization (MUSIC) algorithm to identify three sinusoids in noise that have frequencies very close to each other.

`matlab_hilb`

This demo uses MATLAB as a signal source to produce a Hilbert matrix. The Hilbert matrix is an ill-conditioned matrix used to test the robustness of numerical linear algebra routines. The matrix element (i,j) has the value of 1 / (i + j - 1). The matrix values appear similar to the coefficients of a discrete Hilbert transformer.`matlab_eig`

This demo shows the use of MATLAB to perform eigendecomposition of a 2 x 2 Hermitian symmetric complex matrix. A matrix of eigenvectors and a matrix of eigenvalues are produced. The eigenvalues are real because the input matrix is Hermitian symmetric`sombrero`

This demo is an entire universe composed of a cascade of four MATLAB stars. The MATLAB stars are used a signal source and a signal sink. The overall system generates and plots a mathematical model of a two-dimensional sinc function that resembles a sombrero.`filterPrototype`

This system uses a halfband lowpass filter prototype for the lowpass and highpass filters. All parameters are computed using MATLAB.

`demoscript`

An interactive tutorial that leads a user through a session that runs a simple universe.`sinescript`

This demo runs the same sine wave modulation universe three times, each time with a different frequency.`xmitber`

This demo runs a bit error determination universe at various noise levels and then plots the output.