Introduction to signal processing: sine wave
and complex signals
E. Tikhonov and M. Sneps-Sneppe
Abstract — The article describes one of the powerful tools for the signals analysis while transmission - complex mathematics. Here the basic concepts (complex numbers, complex signal as expansion of a real signal to a complex-valued function) and some operations (exponentiation, absolute value and angular phase representation) are introduced.
As an example complex mathematical methods are used for the case of classical electrical circuits with resistors, capacitors, coils. Actions can be performed both analytically and on a computer, using, for example, MatLab or any other computing software.
Illustrations of the signal passage through electrical circuits are supplied with a commented Matlab code, which is used to their producing.
The article is supposed to be used as a teaching aid, for self-learning the basics of the signal processing. It is considered as a first introduction part of a series on modern signal processing technologies.
Three main approaches to the exploring of the case are also demonstrated: ideological (experiment and its principal explanation), mathematical and computer modeling.
Keywords — complex numbers, complex signals, AC circuits, impedance, MatLab.
I. INTRODUCTION
This article is an introduction to a series related to the signal processing tasks, actual for radio astronomy measurements and satellite data collection. The works are focused on applications of the "Ventspils International Radio Astronomy Center" of the Ventspils University of Applied Sciences, for example, for signal processing at the Reconfigurable Communication Subsystem, obtained from a 16-meter diameter parabolic antenna [1]. Digital section consists: in-flight reconfigurable FPGA baseband processor; subsystem managing microcontroller; radiation resistant FRAM memory. Analog section: 802.11a WLAN ADC/DAC frontend; 802.11a WLAN RF transceiver; RF PA and LNA.
Oa-Bcard Computet
FRAU Mamory
Fig. 1. Reconfigurable Communication Subsystem block diagram [1]
Manuscript received February 13, 2019.
Eugene Tikhonov is with Ventspils University of Applied Sciences,
101a Inzenieru Street, LV-3601, Ventspils, Latvia (e-mail: [email protected]).
Manfred Sneps-Sneppe is with Ventspils University of Applied
Sciences, 101a Inzenieru Street, LV-3601, Ventspils, Latvia (e-mail: [email protected]).
The rest of the paper is the following. In Section 2, electrical response to alternative voltage of the basic elements like resistance, capacitor and coil is discussed. Section 3 suggests mathematical description of alternative current and Section 4 apply it to a previously explored elements behavior. Sections 5 and 6 introduce complex numbers and signals (like complex sine wave) correspondently. Section 6 applies this mathematical tool for each element. And in Section 7 this method is used to analyze RC and RL-circuits analytically and using MatLab computing.
II. Resistance, capacitor and coil (R, C, L)
ALTERNATING CURRENT
We borrow the description of schemes from the textbook [2].
If we apply the voltage of a simple alternative (harmonic) current source to the basic electrical elements and then explore the current flowing through them using oscilloscope, we get different results, shown on Figure 2.
a) Voltage and current on the resistor I :: I
AX U U(t)
V| M
0.02 0.04 0.06 0.08 0.1 Time, sec. 0.12 0.14
g 1 b) Voltage and current on the capacitor
; (V); l / xv/ NY/ NY,
D 0.02 0.04 0.06 0.08 0.1 Time, sec. 0.12 0.14
S 1 c) Voltage and current on the coil
U
S 0 VV \Y/ \V/
U(t)
^C
U(t)
Ml
Fig 2. Time diagrams of alternating current on basic elements for applied harmonic voltage with amplitude 1 V, frequency 50 Hz simulation: resistor with resistance 2 kOm (a), capacitor with capacitance 5 mF (b), coil with inductance 10 mH (c) comparing with the voltage.
Matlab code:
% Time parameters
t0 = 0; % start time
t1 = 8/50; % finish time
steps = 1000; % time accuracy step
% Cosine generator signal parameters
U0 = 1; % amplitude
w0 = 50*pi; % angular frequency
p0 = pi; initial phase
% RCL parameters
R = 2e3; % 2 kOm
C = 5e-3; % 5 mF
L = 10e-3; % 10 mH
% Processing
SI R
t = linspace(t0, t1, steps); % timeline
% voltage
U = U0*cos(w0*t + p0);
% resistance current
I_R = (U0/R)*cos(w0*t + p0);
% capacitor current
I_C = (U0/w0/C)*cos(w0*t + p0 + pi/2); % coil current
I_L = (U0*L*w0)*cos(w0*t + p0 - pi/2); % Plot a: resistor
figure();subplot(3,1,1);hold on; grid on; plot(t, U); % plot voltage
plot(t, 1e3*I_R); % plot current title('a) Voltage and current on the resistor'); xlabel( 'Time, sec.');ylabel('U_R (V); I_R (mA)'); legend('U_R', 'I_R');
% Plot b: capacitor
subplot(3,1,2);hold on; grid on;
plot(t, U); % plot voltage
plot(t, I_C); % plot current
title('b) Voltage and current on the capacitor');
xlabel('Time, sec.') ;
ylabel('U_C (V); I_C (A)');
legend('U_C', 'I_C');
% Plot c: coil
subplot(3,1,3);hold on; grid on; plot(t, U); % plot voltage
plot(t, I_L); % plot current
title('c) Voltage and current on the coil'); xlabel( 'Time, sec.');ylabel('U_L (V) / I_L (A)'); legend('U_L', 'I_L');
In all the cases the shape of the current is the same as for initial voltage. But it is scaled in all of them and shifted in time for the capacitor and the coil.
If we do some more tests then we explore:
— time shift doesn't depend on the parameters of the element - just on the type;
— a increase in frequency leads to a proportional change in the scale of the current for the same coil (increases) and inductive (reduces) and doesn't affect the resistor;
— increasing the nominal of the element (resistance, capacitance, inductance) proportionally changes the scale: reduces (resistance, capacitance) or increases (coil).
Let's try to explain these properties and offer the way to predict the currents and voltages of the complex schemes with these types of elements.
III. Sine wave (alternative voltage)
The voltage of a harmonic current source (that we could interpret as a signal) has a shape of one of the basic analogous (i.e. continuous, without breaks) signals called "sine wave". It's a mathematical curve that describes a smooth periodic oscillation (most commonly used in communications). One of its common form is:
s(t) = A cos (oit + ç0) = A cos q>(t)
where:
— "A" is called amplitude (corresponds to maximum value of s)
— "ai" is called angular frequency (corresponds to changes measured in radians per second". To use
common frequency in Hz we should divide to 2 n: f = -
' 2 n
— " is called initial phase (specified initial value of s for t=0, in radians). For example, if ç0 = n/2 then s(t) = A cos (œt + ^/2) = A sin (œt)
— "ç = p(t)" is called full phase or current phase.
If all but t parameters always stay the same, we say that signal is stationary. In practice some of these parameters changes to add some more information - that is called "modulation".
As we have just one changing parameter (t or <p(t) derived from it) we could make differens 2D visualizations of it:
b) Signal phase 'pit), 'f e [—n -i- it]
a) Signal s(t)
ip e [0 - 2jr]
c) Signal a(tp), <p f [—7r jr]
d) Signal s{tp), ip € [—jr jr]
lA=l ! /o) 1 / 1 / 1 / 1 N 1 \i = 0, ^1,... _ v „ „ l\ 1 1 \ 1 1 \ 1 ..
1 / 1 r 1 / 1 / 1 / 1/ 1 ^ \ 1 \ 1 \ 1 \ 1 \J 1
Phase, ip = ip(i) = wt+^o e [—tt : n]
Fig. 3. Sine wave different visualizations: s(t) as a time-function (a), phase <p(t) as a time-function (b), s(tp) as a phase-function (c) on Cartesian planes, and s(tp) as a phase-function on the polar plane (d).
Because the sine/cosine takes on the same values for any x = x0 + 2n x i (where i is integer), polar description (signal as a function of a phase) is not definitely reversible. For each function value there are multiple possible phases giving it. Several interpretations of that are possible:
— suggest that phase abruptly decreases to 0 when reaches 2n (see b, c on Figure 3);
— imagine sine function as an infinite number of phases (indistinguishable), that are "active" when their value is in range [0 - 2n] (see b on Figure 3);
— consider only representation in polar coordinates (see d on Figure 3) - with obvious round cycle;
The interesting and important fact is that rate of sine wave change is also sine wave (scaled and time shifted):
ds(t) d cos(œt + ç0)
dt
dt
= —Am sin(o)t + q>0)
= Am cos (2 + + Vo) = œs (f +--)
2u>'
IV. Sine wave passing throught resistor,
CAPACITOR, COIL IN CIRCUITS
When sine wave signal (sinusoidal alternating voltage) is applied to active (ohmic) resistance, value current is proportional to fall of potential on this element and has the zero time shift, i.e., they coincide in phase (see Figure 2a and 4a).
b) c)
a)
Ik Uk
<P,=o°
Ic
h Tt=-90°
__Uc^
Ul
<P=90°
Ic (t) = cOr:
dt
= C-
dt
= œC x U(t + —)
starting to move throught the coil under the force of the applied external voltage U(t), it generates a magnetic field around it according to Maxwell's laws. This magnetic field creates a self-induction vsoltage E(t) in the coil, directed against the external voltage and equals to that to that slowing down the current rise on the coil. When the external field becomes zero, the current arising stops (and so the self-induction does), the current has the maximum value. Then the current decreasing to zero begins. After that the same cycle in the negative values takes the place.
In mathematical form:
E(t) = —U(t) = —L
dl(t) dt
= — cos(œt + ç)
So,
Fig. 4. Vector diagram for of alternating voltage and current value on basic elements: resistor (a), capacitor (b), coil (c)
This is displayed in physics as Ohm's law as a definition for resistor resistance:
UR (t) = R x I(t)
The situation is different when an external sine voltage is applied to a capacitor (two metal plates separated by a dielectric). Through an ideal capacitor the flow of electrons (that is an electric current in fact) does not pass. But when external voltage is changed, the first capacitor's plane immidiately accumulates additional electrons and the second plate emits them compensating the new voltage level. The greater the change in voltage (over time), the greater is the number of accumulated and emmited electons (over that time). In oher words, the current is proportional to a voltage change rate.
Therefore the "charging" current flows in the circuit before and after the capacitor during the voltage growth to the peak. When external voltage begins to decrease (from the maximum value to 0), a similar reversed "discharge" apears. Then the same process occurs in the opposite direction (to a negative maximum and again back to zero).
Mathematically, this current dependence on changes in external voltage is written as the definition of inductance (assuming the current is linearly dependent on voltage change rate):
tdUc(t) d cos(œt + n
If n \
U(t)= -CJc(t- 20Ù
Thus, the sinusoidal form of the signal is preserved. This fact has fundamental importance for the linearity of systems and for using of complex voltages that can be represented as sum of sinusoids due to the Fourier transform. But charging/discharging of a capacitor causes a shift between the voltage and the current (phase shift (p = 90° corresponds to — — time shift), which is shown in the timeline (Figure 2b) and vector diagrams (Figure 4b). In addition, the amplitude of current and voltage does not depend only on the capacitance C, but also on the frequency of change of the external voltage.
Similarly, when a flow of electrons (electical current)
I(t) =--sin(wt + q>) = — cos fwt + m--)
lm lm V 2>
= 1 U (t n ^ Lw V 2u>'
So the current is ahead of the voltage in the coil in 90
degrees (or — seconds), or we could say the voltage is behind of the current. Sinusoidal waveform is preserved, the proportionality coefficient of the amplitudes of voltage and current also depends not only on the characteristics of the coil (inductance), but also on the frequency of the signal, that is, the rate of its change:
uloo = lm x IL(t + 2^)
Composite schemes with these elements can be analyzed using vector form (analytically or graphically), or the complex representation of numbers and signals could be uased to keep the convenient simple form of Ohm's law: Uc (t) = Xc x I(t) UL (t) = XL x I(t)
V. Complex numbers
Complex number is the mathematical abstraction, that has the form: x = a + jb, where a and b are real numbers (called "real part" and "imaginary part" respectively) and j is a solution of the equation j2 = 1. No real number correspond a j, so it is called "imaginary". Imaginary and real part of complex number are completely independent (one cannot be received from the other with any mathematical operations).
Complex numbers are the powerful tool to operate with different physical processes. Basic operations with these type of numbers are obvious from the definition (addition, multiplication and so on).
Real part of x (a) usually noted Re(x) or Re x x and imaginary part (b) noted Im(x)or Im x. The number x* = a — jb = Re(x) — jlm(x) is called "complex conjugate". Because of x is a function of 2 independent parameters (a and b), it could be shown on the "complex" plane:
Complex number z
a) Complex z(t) in 2d
Fig. 5. Complex number z = Re(z) + jlm(z)on complex plane (polarproperties: length\z\ = ^Re(z)2 + Im(z)2, phase (p = arctg(Jm 2)) and complex conjugate z* = Re(z) -jlm(z).
Also it could be written as a vector:
= (re(xa
W = ( Im.(x) )
^b^ \Im(x) J for processing as common Cartesian vectors. Alternative we could use the polar coordinates - length and angle: |x|\ _ / Va2 + b2 <p ) \arctg(b/a)) In this representation we could use all the methods of analytical geometry.
Exponentiation: Euler's formula states that, for any real number x:
eix = cos X + j sin X where e is the base of the natural logarithm. Complex conjugate for this number is cos x — j sin x = cos(—x) + j sin(— x) = e-ix.
This type of representation of trigonometric functions make many operations with them easier.
VI. Complex sine wave
Sine wave could be written as
s(t) = A cos(œt + po) = Re(A x ej(Mt+^o))
= Re(A x ej<po x ejMt) = Re(Z x ejMt) = Re{Z(t)),
or:
Re Z(t)+ Im Z(t) s(t)=-j-'
where Z = A x called "complex amplitude" and
Z(t) = Z x ejMt = A x ej<po x ejMt = A(cos(œt + y0) + j sin(o)t + ç0)) called "complex sine wave".
1
s 0 °= -1
Time, t b) Complex z(t) in polar
120 60 lm(z) 1.5
. 1 0.5 \ t=0 1 ~
z(t) ^ ^^^
270
c) Complex Re z / Im z in 3d
Fig. 6. Complex sine wave z(t) = X e^03t (with A = 1,
& = 2n, = u/q): Cartesian (a), polar (b), three-dimensions Re-Im (c) and Abs-Angle (d) representations. Real part corresponds to real sine wave s(t) = cos(&t + with the same parameters (colored red)
lm(z)
Re(z) = s(t)
0
Time, t
< -
0
2
90
150
30
330
240
300
Angle z
Time, t
Matlab code (some design features are skipped):
% Complex sine wave figure('units', 'normalized' 0 0.8 1]);
'outerposition', [0
%% Constants definition A = 1; w = 2*pi; ph0 = pi/8;
scale = 2; N_counts = 100; t_max = 3; t min = 0;
amplitude angle velocity initial phase
accuracy
max time to plot min time to plot
%% Processing
t = linspace(0, t_max, N_counts); % timeline t_step = (t_max-t_min) / N_counts; % time step
sine_c = A*exp(j*(w*t+ph0)); Re = real(sine_c); Im = imag(sine_c);
o = @(t) zeros(size(t)); e = @(t) ones(size(t)); ph_full = @(t) w*t+ph0;
complex sine wave real part of z(t) imaginary of z(t)
zero-vector
one-vector
full-phase
%% Figure a1: Re/Im in 2d subplot(4, 2, 1); hold on; grid on;
plot(t, Re, 'r'); plot(t, Im, 'b');
plot Re plot Im
%% Figure a2: Abs / angle in 2d
subplot(4, 2, 3); hold on;
plot(t, abs(sine_c), 'm'); % plot Abs
% plot angle
for i = ceil((w*t_min+ph0)/2/pi)-1:floor((w*t_max+ph0)/2/pi) % before and after pi t_start = min(max((pi t_max);
t_end = max(min((pi* t_min);
t_temp = t_min:t_step:t_max; plot(t_temp,ph_full(t_temp)-2*pi*i, [0 0.7 0], 'LineStyle', ':'); % -pi - pi
t_temp = t_start:t_step:t_end; plot(t_temp,ph_full(t_temp)-2*pi*i, [0 0.7 0], 'LineStyle', '-'); end
*i-1)-ph0)/w, t_min) i+1)-ph0)/w, t_max)
'Color' 'Color'
%% Figure b: Abs/angle in 2d subplot(4, 2, [5,7]); scale_trick=polar([0 0], [0 A*2]); % polar scale trick
set(scale_trick, 'Visible', 'off'); hold on;
polar(angle(sine_c), abs(sine_c), 'k'); % Polar
%% Figure c: Re/Im in 3d
subplot(4, 2, [2,4]); hold on; grid on;
view(30, 26);
% Complex sine curve
plot3(t, Re, Im, 'k', 'LineWidth', 3); plot3(t, e(t)*scale, Im, 'b'); % Im projection plot3(t, Re, -e(t)*scale, 'r'); % Re projection plot3(-1*e(t), Re, Im, 'g'); % Polar projection
% polar vectors
quiver3(t, o(t), o(t), o(t), Re, Im, 0, 'Color', [0.5 0.5 0.5]);
%% Figure d: abs / angle
subplot(4, 2, [6,8]); hold on; grid on; view(30, 26);
% Compex sine wave
plot3(t, angle(sine_c), abs(sine_c), 'k',
'LineWidth', 3); % Abs projection
plot3(t, e(t)*1.5*pi, abs(sine_c), 'b'); % Angle projection
plot3(t, angle(sine_c), -e(t)*scale, 'g');
Matlab code for the Figutre 7 (some design features are skipped):
% Complex sine forming
figure('units', 'normalized', 'outerposition', [0 0 0.5 1])
%% Constants definition w = 2*pi; % angle velociry
ph0 = pi/8; % initial phase
A = 1; % amplitude
t_max = 2; % time to plot
N_counts = 60; % accuracy
%% Signal
t = linspace(0, t_max, N_counts);% timeline o = zeros(size(t)); % zero-vector
sine_c = A*exp(j*(w*t+ph0)); sine_r = real(sine_c); sine_i = imag(sine_c);
comlex sine wave
%% Figure a: real
subplot(2, 2, 1); hold on; grid on; xlim([0 t_max]); ylim([-1 1]); zlim([-1 1]); view(-40, 40);
plot3(t, sine_r, o, 'r'); % sine curve
quiver3(t, o, o, o, sine_r, o, 0, 'r'); % ar-field
%% Figure b: imaginary part subplot(2, 2, 2); hold on; grid on; xlim([0 t_max]); ylim([-1 1]); zlim([-1 1]); view(-40, 40);
plot3(t, o, sine_i, 'b'); % sine curve
quiver3(t, o, o, o, o, sine_i, 0, 'b'); % ar-field
%% Figure c: real + imaginary subplot(2, 2, 3); hold on; grid on; xlim([0 t_max]); ylim([-1 1]); zlim([-1 1]); view(-3,20);
plot3(t, sine_r, o, 'r'); % sine curves
plot3(t, sine_r, sine_i, 'b');
quiver3(t, o, o, o, sine_r, o, 0, 'r'); % ar-field quiver3(t, sine_r, o, o, o, sine_i, 0, 'b')
%% Figure d: result polar subplot(2, 2, 4); hold on; grid on; xlim([0 t_max]); ylim([-1 1]); zlim([-1 1]); view(-30,15);
quiver3(t, o, o, o, sine_r, sine_i, 0, 'k'); plot3(t, sine_r, sine_i, 'k'); % arrow field
a) Real part of z(t) = s(t)
Re z(t)
-1 0
Time, t
b) Imaginary part of z(t) = s(t)
-0.5
Re z(t)
-1 0
1 c) Summarize tlie vectors Re z(t) + j Im z(t)
d) Resulting complex sine wave z(t)
Re z(t) ° Time, t
Fig. 7. Complex sine wave forming z(t) = A Re z(t) + j A Im z(t) = A cos(a>t + <p0) + j A sin(a>t + <p0) forming process (with M = 2n, tp0 = u/q)
So, the real and imaginary parts in complex sine wave are not independent. They are bind so that imaginary part exactly repeats the real part with a phase shift n/2 corresponding to time delay n/2^
The real part of complex sine wave corresponds to a real sine signal. And imaginary part is synchronous addition to that, and is needed to provide this in a complex space.
VII. COMPLEX IMPEDANCE OF THE RESISTOR, CAPACITOR, COIL
The phase difference between the voltage and the current of the coil or capacitor can be taken into account, while maintaining the convenient form of Ohm's law (voltage equal to the current multiplied by the resistance). Of couse, we should imagine complex sine wave generator U(t) to provide that it's real part will correspond to a real sine signal.
We assume that the resistance has a phase: Rc = Xc = \Xc\ei<P
The phase delay of the voltage on the capacitor by -90 degrees from the current ((p = -n/2, that is, 1/4 of the full period) means:
Uc = Ic X \XC\e* = Ic X \XC= Ic X — e-%
ti>L
According to Euler's formula:
eJ(-2) = cos 72) + j sin 72) = -j =-
Therefore, Ohm's law for a capacitor can be written in the form of capacitance:
Uc (t) = Xc X Ic (t) = x ic (t) = X Ic (t)
Similarly, due to the fact that the amplitudes of the voltage and current on the coil are related as \XL \ = mL, and the current shift from the voltage is +90 degrees: UL(t)= XL X IL(t)= juL X IL(t)
These complex resistance analogs are called "reactance" -resistance that a capacitor or inductance exerts on alternating current. Reactive resistance in contrast to active resistance does not consume energy, but accumulates it and returns it to the circuit.
Any real elements contain both active and reactive (inductive and capacitive) resistance and can be represented as their sum, the so-called impedance [3]:
Z = (Re Z + j Im Z) = R + jX UA(t) = Z X IA(t) = (Re Z + j Im Z) X IA(t) = \Z\eVz X IA(t)
The impedance at X = 0 corresponds to the active resistance R, at X> 0 has a capacitive and when X<0 -inductive type.
Calculations of parallel and serial connections of active elements (coils and capacitors), as well as passive (resistors), can be used in this complex form, as well as in the vector form:
-'parallel
=!Zi
i=1 n
- = Y1
4-iZi
Matlab code: complex method gives the same result as on Figure 2.
% Time parameters
t0 = 0; % start time
t1 = 8/50; % finish time
steps = 1000; % time accuracy step
% Cosine generator signal parameters U0 = 1; % amplitude
w0 = 50*pi; % angular frequency
p0 = pi; % initial phase
% complex sine wave c = @(t) U0*exp(j*(w0*t + p0));
% RCL parameters
R = 2e3; % 2 kOm
C = 5e-3; % 5 mF
L = 10e-3; % 10 mH
% Impedance X_R = R;
X_C = 1/(j*w0*C); X_L = j*w0*L;
% Processing
t = linspace(t0, t1, steps);
U = c(t); I_R = U/X_R I_C = U/X_C IL = U/X L
voltage
resistance current capacitor current coil current
D" f^l 1
" rr «* \ i
Fig. 8. RC circuit and its time and vector diagrams
The time diagram shows the four curves I, U, UC and UR. It is seen that the current I is ahead of the voltage UC by 900, but the phase shift between the common voltage U and the current I is in the range between 0 and 90 degrees. Summary impedance:
Z = Zfl + Zc = R—j± Though voltage in the circuit is U, the current is
u u
Ir = ¡c = ~=~
Voltage drop on resistor is: R
UR = R x IR = U-^ = U-
= U-
= u(l+-
1 -j
wRC
+J:
œRC œRC - j
œRC
% Plot a: resistor
figure();subplot(3,1,1);hold on; grid on; plot(t, real(U)); % plot voltage
plot(t, 1e3*real(I_R)); % plot current title('a) Voltage and current on the resistor'); xlabel('Time, sec.'); ylabel('U_R (V); I_R (mA)'); legend('U_R', 'I_R');
% Plot b: capacitor
subplot(3,1,2);hold on; grid on;
plot(t, real(U)); % plot voltage
plot(t, real(I_C)); % plot current
title('b) Voltage and current on the capacitor');
xlabel('Time, sec.') ;
ylabel('U_C (V); I_C (A)');
legend('U_C', 'I_C');
% Plot c: coil
subplot(3,1,3);hold on; grid on;
plot(t, real(U)); % plot voltage
plot(t, real(I_L)); % plot current
title('c) Voltage and current on the coil');
xlabel('Time, sec.');
ylabel('U_L (V) / I_L (A)');
legend('U_L', 'I_L');
VIII. Example. Analysis of RC and RL circuits
A. RC circuit Figure 7 shows the series connection of resistance R and capacitor C. The task is to determine the particular voltages uR and uC. The circuit has a resistance R and a reactance capacitor Xc, which introduces a phase shift between current I and voltage U.
(œRC)2 — 1 J (œRC)2
Voltage drop on capacity is:
i 1 a>RC
Uc = U — UR = U[-^-r- — j
c R \1 — (uRC)2 J
So the real voltages are:
n )
(coRC )2
i)
Re U,
= + ) 1
Re Ur = U-
(a>RC )2
B. RL-circuit
Figure 8 shows the series connection of resistance R and inductance L. The task is to determine the particular voltages UR and UL. The circuit has a resistance R and a reactance inductance XL, which also introduces a phase shift between current I and voltage U, but it has a different character.
Figure 9. RL-circuit and its time and vector diagrams
Figure 8 shows that the reactive voltage UL is 90° ahead of the current I, the active voltage UR is in the same phase as the current I, the phase shift between the common voltage U and the current I is in the range between 0 and 90 degrees.
Summary impedance:
Z = ZR + ZL = R + juL
Though voltage in the circuit is U, the current is _ -U - U
Ir = II = Z = R + juL Voltage drop on resistor is:
R 1
UR = R x IR = U „ . . , = U
R + ja>L
= U
1+ i
,u>L T
(ml/r)2 J 1
n)
Voltage drop on inductance is: ju>L
u>L
u, = UJ . , = u[—v 'RJ 2 + j- 2
■ ■ (<ol/r) 1 -(œL,)
R + ju>L y i So the real voltages are: Re(UR) =
'R) U
Re(UL ) = U
- (œl/r) r(ml/r)2 1 - (ml/r)2
a) Voltages and current in RC-circuit
0.02
0.05
0.03 0.04
Time, sec. b) Voltages and current in RL-circuit
0.03 0.04 0.05 0.06
Time, sec.
Figure 10. Time diagrams simulation for RC (a) and RL (b) circuits for applied harmonic voltage with amplitude 1 V, frequency 50 Hz: resistor with resistance 2 Om, capacitor with capacitance 5 mF, coil with inductance 10 mH
Matlab code:
% Processing
t = linspace(t0, t1, steps); U = c(t);
I_RC = U/(X_R+X_C); U_RC_R = X_R * I_RC; U_RC_C = X_C * I_RC;
I_RL = U/(X_R+X_L); U_RL_R = X_R * I_RL; URLL = XC * I RL;
voltage
RC-current
U_R in RC circuit
U_C in RC circuit
RC-current
U_R in RL circuit
U_C in RL circuit
% Plot a: RC figure();subplot(2,1,1); plot(t, real(U)); plot(t, real(I_RC)); plot(t, real(U_RC_R)); plot(t, real(U_RC_C));
hold on; grid on; % plot voltage % plot current % plot U_R % plot U_C
% Time parameters
t0 = 0; % start time
t1 = 4/50; % finish time
steps = 1000; % time accuracy step
% Cosine generator signal parameters
U0 = 1; % amplitude
w0 = 50*pi; % angular frequency
p0 = pi; % initial phase
% complex sine wave
c = @(t) U0*exp(j* (w0*t + p0));
% RCL parameters
R = 2; % 2 Om
C = 5e-3; % 5 mF
L = 10e-3; % 10 mH
% Impedance
XR = R;
XC = 1/(j*w0*C);
XL = j*w0*L;
title('a) Voltages and current in RC-circuit'); xlabel( 'Time, sec.');ylabel('U(V); I (A)'); legend('U', 'I', 'U_R', 'U_C');
% Plot b: RL
subplot(2,1,2); hold on; grid on; plot(t, real(U)); % plot voltage
plot(t, real(I_RL)); % plot current
plot(t, real(U_RL_R)); % plot U_R plot(t, real(U_RL_L)); % plot U_L title('b) Voltages and current in RL-circuit'); xlabel( 'Time, sec.');ylabel('U(V); I (A)'); legend('U', 'I', 'U_R', 'U_L');
IX. Conclusion
The basic concept tool of signal processing analysis was discussed: the complex mathematical method basics, like complex sine wave signal and related operations. The in following parts some technologies based on this idea are planned to discuss: spectrum, analytical (complex) representation of a signal, information (digital) signal, modulations, spectrum modifications (like Direct Sequence Spread Spectrum, Frequency Hopping, Time Hopping, Chirp Spread), etc.
References
[1] J. Sateet all. Concept of spectrally efficient communication subsystem// Space Research Review, Vol.4, January 2016, 52-61.
[2] Werner Dzieia, Josef Kammerer, Wolfgang Oberthur, Hans-Jobst Siedler, Peter Zastrow Elektotechnische Grundlagen der Mirkoelektronik. Lehrbuch. 2001, ISBN 3-7905-0707-5.
[3] Horowitz, Paul; Hill, Winfield. The Art of Electronics, 7th printing 2016 with corrections. Cambridge University Press. ISBN 978-0-52180926-9.
[4] Charles K. Alexander; Matthew N.O. Sadiku. Fundamentals of Electric Circuits - 5th ed. p.cm. The McGraw-Hill Companies, Inc. NewYork., 2013, ISBN 978-0-07-338057-5.
0
0