Example: Signal Generation

Functions > Reading and Writing Files > Example: Signal Generation

Signal Synthesis

Consider a dual-tone multi-frequency (DTMF) signal. These signals are commonly used for touch tone telephone dialing and remote control functions. DTMF signals are composed of one high frequency and one low frequency sinusoidal signal. A table of common tones is shown below.

Digital Phone Keypad		High Tone Group
Digital Phone Keypad		1209 Hz	1336 Hz	1477 Hz	1633 Hz
Low	697 Hz	1	2	3	A
Tone	770 Hz	4	5	6	B
Tone	852 Hz	7	8	9	C
Group	941 Hz	*	0	#	D

For example, to generate the tone for the number 5 button, mix the 770 Hz signal from the Low Tone Group with the 1336 Hz signal from the High Tone Group.

In this example, use PTC Mathcad to generate the # tone.

1. Define the number of channels, sample rate and resolution.

2. Set the time duration and the number of data points.

3. Calculate the time interval and generate two arrays, one for indexing and one for time.

4. Set the amplitude of the signal and choose the two frequencies required to generate the # tone.

5. Use the sin function to generate the signal, and plot it.

<region id="ID0EUT41" actualWidth="389.62000000000006" actualHeight="31.791200000000003" top="1132.8000000000002" left="38.400000000000006" xmlns="http://schemas.mathsoft.com/worksheet50">
 <math resultRef="22">
 <define xmlns="http://schemas.mathsoft.com/math50">
 <id labels="VARIABLE" xml:space="preserve">Signal</id>
 <apply>
 <vectorize />
 <parens>
 <apply>
 <mult />
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">Amplitude</id>
 <parens>
 <apply>
 <plus />
 <apply>
 <id labels="FUNCTION" xml:space="preserve" label-is-contextual="true">sin</id>
 <apply>
 <mult />
 <apply>
 <mult />
 <apply>
 <mult />
 <real>2</real>
 <id labels="CONSTANT" xml:space="preserve" label-is-contextual="true">π</id>
 </apply>
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">
 f<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">1</Subscript>
 </id>
 </apply>
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">τ</id>
 </apply>
 </apply>
 <apply>
 <id labels="FUNCTION" xml:space="preserve" label-is-contextual="true">sin</id>
 <apply>
 <mult />
 <apply>
 <mult />
 <apply>
 <mult />
 <real>2</real>
 <id labels="CONSTANT" xml:space="preserve" label-is-contextual="true">π</id>
 </apply>
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">
 f<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">2</Subscript>
 </id>
 </apply>
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">τ</id>
 </apply>
 </apply>
 </apply>
 </parens>
 </apply>
 </parens>
 </apply>
 </define>
 <resultFormat>
 <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="false" />
 </resultFormat>
 </math>
</region>

<region id="ID0ECU41" actualWidth="576.621459277374" actualHeight="249.600005722046" top="1171.2000000000003" left="38.400000000000006" height="249.600005722046" width="576.621459277374" xmlns="http://schemas.mathsoft.com/worksheet50">
 <plot background-type="white" origin-positioning="true">
 <xyPlot>
 <title class="- topic/title " wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
 <legend />
 <traces>
 <trace resultRef="24">
 <traceStyle color="#FF068149" symbol="none" line-weight="1" line-style="Solid">lines</traceStyle>
 </trace>
 </traces>
 <graph-size width="438.221459277374" height="192.000005722046" />
 <axes>
 <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="0.01">
 <axisLine position="origin" positionticmark="5" legendWidth="60.6233333333333" />
 <axisGrid>
 <gridFrequency>11</gridFrequency>
 <gridLabels display="true" />
 <gridLines />
 <tickMarks display="true" />
 </axisGrid>
 <axisLabel />
 <markers />
 <numberFormat>
 <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" />
 </numberFormat>
 <plotEquations>
 <plotEquation>
 <math resultRef="27">
 <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">τ</id>
 </math>
 <math resultRef="28">
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </math>
 </plotEquation>
 </plotEquations>
 <xyDomain scale-type="linear" auto-scale="true">
 <startValue>
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </startValue>
 <endValue resultRef="26">
 <real xmlns="http://schemas.mathsoft.com/math50">.01</real>
 </endValue>
 </xyDomain>
 </xAxis>
 <yAxis rank="1" legend-position="PlotBoundaryLeft" start="-20000" end="20000">
 <axisLine position="origin" positionticmark="0" legendWidth="96.66" />
 <axisGrid>
 <gridFrequency>11</gridFrequency>
 <gridLabels display="true" />
 <gridLines />
 <tickMarks display="true" />
 </axisGrid>
 <axisLabel />
 <markers />
 <numberFormat>
 <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" />
 </numberFormat>
 <plotEquations>
 <plotEquation>
 <math resultRef="29">
 <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">Signal</id>
 </math>
 <math resultRef="30">
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </math>
 </plotEquation>
 </plotEquations>
 <xyDomain scale-type="linear" auto-scale="true">
 <startValue>
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </startValue>
 <endValue>
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </endValue>
 </xyDomain>
 </yAxis>
 </axes>
 </xyPlot>
 </plot>
</region>

6. Use the WRITEWAV function to write the signal to a WAV file using the specified parameters.

If the specified bit resolution is 1-8, the data is written to the file as unsigned byte data. The limits of unsigned byte data are 0-256 (28). If the bit resolution is 9-16, word data (two bytes) is written to the file. The limits on word data are –32768 - +32767 (215).

7. Try changing the frequencies for the DTMF signal and then play the sound file. The sound is similar to a touch-tone telephone.

Exponential Ramp in Frequency

Another useful signal in audio testing is the exponential ramp in frequency.

1. Define an exponential ramp in frequency signal.

2. Define the sample rate, the resolution, and the number of channels.

3. Define the following parameters for the new signal.

<region id="ID0E4441" actualWidth="281.95000000000005" actualHeight="29.811200000000007" top="2131.2000000000003" left="38.400000000000006" xmlns="http://schemas.mathsoft.com/worksheet50">
 <math resultRef="48">
 <define xmlns="http://schemas.mathsoft.com/math50">
 <id labels="VARIABLE" xml:space="preserve">NumPts</id>
 <eval>
 <apply>
 <mult />
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">SampleRate</id>
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">
 τ<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">f</Subscript>
 </id>
 </apply>
 <unitOverride>
 <placeholder />
 </unitOverride>
 </eval>
 </define>
 <resultFormat>
 <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="false" />
 </resultFormat>
 </math>
</region>

4. Use the sin function to generate the signal, then use the length function to find its length

<region id="ID0EGC51" actualWidth="282.08000000000004" actualHeight="29.6" top="2400.0000000000005" left="38.400000000000006" xmlns="http://schemas.mathsoft.com/worksheet50">
 <math resultRef="58">
 <define xmlns="http://schemas.mathsoft.com/math50">
 <id labels="VARIABLE" xml:space="preserve">Signal</id>
 <apply>
 <vectorize />
 <parens>
 <apply>
 <mult />
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">Amplitude</id>
 <parens>
 <apply>
 <id labels="FUNCTION" xml:space="preserve" label-is-contextual="true">sin</id>
 <apply>
 <mult />
 <apply>
 <mult />
 <real>2</real>
 <id labels="CONSTANT" xml:space="preserve" label-is-contextual="true">π</id>
 </apply>
 <apply>
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">er</id>
 <id labels="VARIABLE" xml:space="preserve" label-is-contextual="true">τ</id>
 </apply>
 </apply>
 </apply>
 </parens>
 </apply>
 </parens>
 </apply>
 </define>
 <resultFormat>
 <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="false" />
 </resultFormat>
 </math>
</region>

5. Use WRITEWAV to write the signal to a new file.

6. Plot the signal.

<region id="ID0EZF51" actualWidth="576.621459277374" actualHeight="345.600008583069" top="2649.6000000000004" left="38.400000000000006" height="345.600008583069" width="576.621459277374" xmlns="http://schemas.mathsoft.com/worksheet50">
 <plot background-type="white" origin-positioning="true">
 <xyPlot>
 <title class="- topic/title " wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
 <legend />
 <traces>
 <trace resultRef="64">
 <traceStyle color="#FF2E3192" symbol="none" line-weight="1" line-style="Solid">lines</traceStyle>
 </trace>
 </traces>
 <graph-size width="438.221459277374" height="288.000008583069" />
 <axes>
 <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="1">
 <axisLine position="origin" positionticmark="2" legendWidth="60.6233333333333" />
 <axisGrid>
 <gridFrequency>6</gridFrequency>
 <gridLabels display="true" />
 <gridLines />
 <tickMarks display="true" />
 </axisGrid>
 <axisLabel />
 <markers />
 <numberFormat>
 <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" />
 </numberFormat>
 <plotEquations>
 <plotEquation>
 <math resultRef="69">
 <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">τ</id>
 </math>
 <math resultRef="70">
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </math>
 </plotEquation>
 </plotEquations>
 <xyDomain scale-type="linear" auto-scale="true">
 <startValue>
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </startValue>
 <secondTickValue resultRef="66">
 <real xmlns="http://schemas.mathsoft.com/math50">0.2</real>
 </secondTickValue>
 <endValue resultRef="68">
 <real xmlns="http://schemas.mathsoft.com/math50">1</real>
 </endValue>
 </xyDomain>
 </xAxis>
 <yAxis rank="1" legend-position="PlotBoundaryLeft" start="-10000" end="10000">
 <axisLine position="origin" positionticmark="0" legendWidth="96.66" />
 <axisGrid>
 <gridFrequency>5</gridFrequency>
 <gridLabels display="true" />
 <gridLines />
 <tickMarks display="true" />
 </axisGrid>
 <axisLabel />
 <markers />
 <numberFormat>
 <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" />
 </numberFormat>
 <plotEquations>
 <plotEquation>
 <math resultRef="77">
 <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">Signal</id>
 </math>
 <math resultRef="78">
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </math>
 </plotEquation>
 </plotEquations>
 <xyDomain scale-type="linear" auto-scale="true">
 <startValue resultRef="72">
 <apply xmlns="http://schemas.mathsoft.com/math50">
 <mult />
 <real>-1</real>
 <apply>
 <pow />
 <real>10</real>
 <real>4</real>
 </apply>
 </apply>
 </startValue>
 <secondTickValue resultRef="74">
 <apply xmlns="http://schemas.mathsoft.com/math50">
 <mult />
 <real>-5</real>
 <apply>
 <pow />
 <real>10</real>
 <real>3</real>
 </apply>
 </apply>
 </secondTickValue>
 <endValue resultRef="76">
 <apply xmlns="http://schemas.mathsoft.com/math50">
 <mult />
 <real>1</real>
 <apply>
 <pow />
 <real>10</real>
 <real>4</real>
 </apply>
 </apply>
 </endValue>
 </xyDomain>
 </yAxis>
 </axes>
 </xyPlot>
 </plot>
</region>