Unit 2.3: Elementary Signals#

The preparatory reading for this section is Chapter 1 of [Karris, 2012] which

  • begins with a discussion of the elementary signals that may be applied to electrical circuits

  • introduces the unit step, unit ramp and dirac delta functions

  • presents the sampling and sifting properties of the delta function and

  • concludes with examples of how other useful signals can be synthesised from these elementary signals.

Additional information has been adapted from Section 1.4 of [Hsu, 2020].

Follow along at cpjobling.github.io/eg-150-textbook/signals_and_systems/elementary_signals

QR Code

Introduction#

Consider the network shown in below where the switch is closed at time t=T and all components are ideal.

Network with a switch which is closed at t = T.

Express the output voltage Vout as a function of the unit step function, and sketch the appropriate waveform.

Slido poll 1#

Solution

Before the switch is closed at t<T:

Vout=0.

After the switch is closed for t>T:

Vout=Vs.

We imagine that the voltage jumps instantaneously from 0 to Vs volts at t=T seconds as shown below.

The step function

We call this type of signal a step function.

What is happening at t=0?

Mathematically, the value of the signal at t=0 is undefined and the time derivative is infinite. In practice, we could imagine that the signal will be somewhere on the dotted line. But, we don’t know where.

The Unit Step Function#

u0(t)={0t<01t>0

The unit step function

In MATLAB#

In MATLAB, we use the heaviside function (named after Oliver Heaviside) to model the unit step signal.

cd('/Users/eechris/code/src/github.com/cpjobling/eg-150-textbook/signals_and_systems/elementary_signals')
edit plot_heaviside
plot_heaviside
../../_images/496844206370e5109044ed21a2b18987bc8b7e089b094179d081d079e4ccb12a.png
ans = 1

Note that, so that it can be plotted, MATLAB defines the Heaviside function slightly differently from the mathematically ideal unit step:

heaviside(t)={0t<01/2t=01t>0

Simple Signal Operations#

Amplitude Scaling#

Sketch Au0(t) and Au0(t)

syms t;
u0(t) = heaviside(t); % rename heaviside function for ease of use
A = 2; % so signal can be plotted
fplot(A*u0(t),[-1,1],'LineWidth',2),ylim([-0.2,2.2]),grid,title('Amplitude scaling $$Au_0(t)$$','interpreter','latex')
../../_images/b40feb40e948b27996cbf55afe8322ba8b8eafb8fdebf46ad4aeaa9dcad45fbd.png

Note that the signal is scaled in the y direction.

fplot(-A*u0(t),[-1,1],'LineWidth',2),ylim([-2.2,0.2]),grid,...
title('Amplitude scaling and mirroring $$-Au_0(t)$$','interpreter','latex'),...
xlabel('t')
../../_images/891a2f5712f0bd1c78f39f884754787bebb77a7117fe28ba645aae6b7eb7f537.png

Note that, because of the sign, the signal is mirrored about the x axis as well as being scaled by 2.

Time Reversal#

Sketch u0(t)

fplot(A*u0(-t),[-1,1],'LineWidth',2),ylim([-0.2,2.2]),grid,title('Time reversal $$Au_0(-t)$$','interpreter','latex'),xlabel('t')
../../_images/9746a3208c510a2e78a094ed7e2e95ac8e120de98c88bb57ce626a5a00ba50f5.png

The sign on the function argument t causes the whole signal to be reversed in time. Note that another way of looking at this is that the signal is mirrored about the y axis.

Time Delay and Advance#

Sketch u0(tT) and u0(t+T)

T = 1; % again to make the signal plottable.
fplot(u0(t - T),[-1,2],'LineWidth',2),ylim([-0.2,1.2]),grid,title('Time delay $$u_0(t - T)$$','interpreter','latex'),xlabel('t')
../../_images/c12a12cf12d2eb816892aae038f03eb7cf71501cb48c86ed020363f84f1f0f27.png

This is a time delay … note for u0(tT) the step change occurs T seconds later than it does for uo(t).

fplot(u0(t + T),[-2,1],'LineWidth',2),ylim([-0.2,1.2]),grid,title('Time advance $$u_0(t + T)$$','interpreter','latex'),xlabel('t')
../../_images/99ef4b4b57059fe15bbe4026b3b2048030edfa9aad5e8daf2c6e61e6e24f4ef3.png

This is a time advance … note for u0(t+T) the step change occurs T seconds earlier than it does for uo(t).

Examples#

We will work through some examples in class. See Examples 3.

Synthesis of Signals from the Unit Step#

Unit step functions can be used to represent other time-varying functions such as rectangular pulses, square waves and triangular pulses. See Examples 3 for the examples that we will look at in class.

The Ramp Function#

Integrator (RC) circuit

In the circuit shown above is is a constant current source and the switch is closed at time t=0.

When the current through the capacitor ic(t)=is is a constant and the voltage across the capacitor is

vc(t)=1Ctic(τ)dτ

where τ is a dummy variable.

Since the switch closes at t=0, we can express the current ic(t) as

ic(t)=isu0(t)

and if vc(t)=0 for t<0 we have

vc(t)=isCtu0(τ)dτ=isC00dτ0+isC0t1dτ

So, the voltage across the capacitor can be represented as

vC(t)=isCtu0(t)

Note that in this as in other examples throughout these notes, and in published tables of transforms, the inclusion of u0(t) in vc(t) acts as a “gating function” that limits the definition of the signal to the causal range 0t<.

To sketch the wave form, let’s arbitrarily let C and is be one and then plot with MATLAB.

C = 1; is = 1;
vc(t)=(is/C)*t*u0(t);
fplot(vc(t),[-1,4],'LineWidth',2),grid,title('A ramp function'),xlabel('t')
../../_images/407bcc21d64bc7277ca9e622cfb65fba99fd4bd4f85c9b6a48258cbe94555689.png

This type of signal is called a ramp function. Note that it is the integral of the step function (the resistor-capacitor circuit implements a simple integrator circuit).

The unit ramp function is defined as

u1(t)=tu0(τ)dτ

so

u1(t)={0t<0tt0

and

u0(t)=ddtu1(t)

Note

Higher order functions of t can be generated by the repeated integration of the unit step function.

For future reference, you should determine u2(t), u3(t) and un(t) for yourself and make a note of the general rule:

un1=1nddtun(t)

Details are given in equations 1.26—1.29 in [Karris, 2012].

The Dirac Delta Function#

A differentiating function (RL network)

In the circuit shown above, the switch is closed at time t=0 and iL(t)=0 for t<0. Express the inductor current iL(t) in terms of the unit step function and hence derive an expression for vL(t).

Solution

vL(t)=LdiLdt

Because the switch closes instantaneously at t=0

iL(t)=isu0(t)

Thus

vL(t)=isLddtu0(t).

The unit Impulse Function#

The unit impulse function δ(t), is the derivative of the unit step.

δ(t)=ddtu0(t)

which is tricky to compute because u0(t) is discontinuous at t=0 but it must have the properties

tδ(τ)dτ=u0(t)

To solve this problem we need to invent a function that represents the derivative of the unit step function. This function is called the unit impulse function 𝛿(𝑡), also known as the Dirac delta function (named after Paul Dirac).

Traditionally, δ(t) is often defined as the limit of a suitably chosen conventional function having unity area over an infinitesimal time interval as shown in Fig. Fig. 21.

Visualisation of the Dirac delta function as the limit of a conventional function with unit area.

Fig. 21 Visualisation of the Dirac delta function as the limit of a conventional function with unit area.#

The Dirac delta posseses the following properties

δ(t)={0t0t=0
ϵϵδ(t)dt=1

The unit impulse function plays a fundamental role in systems analysis.

Sketch of the delta function#

Continuing the example, and replacing the derivative of the unit step u0(t) with the unit impulse δ(t)

Vout(t)=VL(t)=isLδ(t)

The delta function

Note when we draw the unit impulse we show the height of δ(t) as one so the height of the impulse in the figure is isL.

MATLAB Confirmation#

syms is L;
vL(t) = is * L * diff(u0(t))
vL(t)=Lisδ(t)

Note that we can’t plot dirac(t) in MATLAB with fplot. The best we can do is a stem plot.

L = 1; is = 2;
stem(0,L*is),ylim([-0.2,2.2]),title('Impulse $$v_L(t) = L*i_s*\delta(t)$$','interpreter','latex'),grid,xlabel('t')
../../_images/a08518d124a137a1c130e1228e1e5a27669e02ce383f73f92e7a1508d76c74b1.png

Important properties of the delta function#

Sampling Property#

The sampling property of the delta function states that

f(t)δ(ta)=f(a)δ(ta)

or, when a=0,

f(t)δ(t)=f(0)δ(t)

Multiplication of any function f(t) by the delta function δ(t) results in sampling the function at the time instants for which the delta function is not zero.

The study of descrete-time (sampled) systems is based on this property.

You should work through the proof for youself.

Sifting Property#

The sifting property of the delta function states that

f(t)δ(tα)dt=f(α)

That is, if multiply any function f(t) by δ(tα), and integrate from to +, we will get the value of f(t) evaluated at t=α.

You should also work through the proof for yourself.

Higher Order Delta Fuctions#

the nth-order delta function is defined as the nth derivative of u0(t), that is

δn(t)=dndtn[u0(t)]

The function δ(t) is called the doublet, δ(t) is called the triplet and so on.

By a procedure similar to the derivation of the sampling property we can show that

f(t)δ(ta)=f(a)δ(ta)f(t)δ(ta)

Also, derivation of the sifting property can be extended to show that

f(t)δn(tα)dt=(1)ndndtn[f(t)]|t=α

Quiz and take aways#

Summary#

In this chapter we have looked at some elementary signals and the theoretical circuits that can be used to generate them.

Unit 2.3: Take aways#

  • You should note that the unit step is the heaviside function u0(t).

  • Many useful signals can be synthesized by use of the unit step as a “gating function” in combination with other signals

  • That unit ramp function u1(t) is the integral of the step function.

  • The Dirac delta function δ(t) is the derivative of the unit step function. We sometimes refer to it as the unit impulse function.

  • The delta function has sampling and sifting properties that will be useful in the development of time convolution and sampling theory.

Exercises#

We will do some of these in class. See Examples 3.

References#

[Hsu20]

Hwei P. Hsu. Schaums outlines signals and systems. McGraw-Hill, New York, NY, 2020. ISBN 9780071634724. Available as an eBook. URL: https://www.accessengineeringlibrary.com/content/book/9781260454246.

[Kar12] (1,2)

Steven T. Karris. Signals and systems with MATLAB computing and Simulink modeling. Orchard Publishing, Fremont, CA., 2012. ISBN 9781934404232. Library call number: TK5102.9 K37 2012. URL: https://ebookcentral.proquest.com/lib/swansea-ebooks/reader.action?docID=3384197.

Next Time#

Systems and Classification of Systems