Ei 7211-Circuit Simulation Lab List of Experiments
Ei 7211-Circuit Simulation Lab List of Experiments
Ei 7211-Circuit Simulation Lab List of Experiments
LIST OF EXPERIMENTS
1. Generation of Continuous Time(CT) and Discrete Time (DT) signals
(i) Standard signals: – impulse, step, ramp, exponential
(ii) Periodic and a-periodic signals
(iii) Deterministic and random signals
2. CT and DT system characterization:
(i) Linearity
(ii) Time invariance
(iii) Causality
(iv) Stability
3. Time response & Frequency response of DT systems
4. Discretization and Reconstruction of signals
(i) Sampling and aliasing effects
(ii) A/D conversion
(iii) D/A conversion
5. Statistical analysis of random signals
6. Verification Kirchhoff’s laws, Thevenin’s and Norton’s theorems.
7. Verification of Superposition, Maximum Power transfer and Reciprocity
theorems.
8. Response of RL, RC and RLC circuits for step input.
9. Frequency response of Series and Parallel resonance circuits.
10.Determination of self and mutual inductances and coupling coefficient of
coupled coils.
11. Power and power factor measurement in three phase circuits by two wattmeter
method.
12. Determination of Z, Y and H parameters of a two port network.
Exp No:1 Generation of Continuous Time(CT) and Discrete Time signals
Aim:
To generate i) standard signals (Impulse, step, ramp and exponential) ii) periodic and
aperiodic signals and iii) deterministic and random signals using Matlab software.
Software Required
Matlab software
Theory:
Program:-
y7=sin(2*pi*50*t);
figure;
subplot(2,2,1);
plot(t,y7);
axis([0 0.1 -2 2]);
xlabel('time');
ylabel('amplitude');
title(' sinsoidal wave signal');
subplot(2,2,2);
stem(y7);
xlabel('n');
ylabel('amplitude');
title('sin wave sequence');
% Random signal:
t=0:1/fs:0.1;
x=rand(1000,1)
subplot(2,2,1);
plot(t,x)
subplot(2,2,2);
stem(x);
xlabel('n');
ylabel('amplitude');
title('sin wave sequence');
xlabel('time');
ylabel('amplitude');
title(' sinsoidal wave signal');
Result:
Various signals & sequences generated using Matlab software.
Exp. No:2 DT system Characterization
Aim:
To test the following properties of CT and DT systems:-
i) Linearity ii) Time invariance iii) Causality and iv)Stability.
Software Required:
Mat lab software 7.0
Theory:
Linearity property: Any system is said to be linear if it satisfies the Superposition principal.
superposition principal state that Response to a weighted sum of input signal equal to the
corresponding weighted sum of the outputs of the system to each of the individual input signals.
Two input sequences x[n] and x[n - D], are generated and corresponding
output sequences y1[n], y2[n] are plotted.
if (d==0)
fprintf('The system is time invariant\n');
else
fprintf('The system is time variant\n');
% stability
%y(n)-y(n-1)+.9y(n-2)=x(n);
b=[1];
a=[1,-1,.9];
figure;
zplane(b,a);
% Causality
N=-10:1:10
b=[1];
a=[1,-1,.9];
x=[zeros(1,10) 1 zeros(1,10)];
y1=filter(a,b,x)
subplot(2,1,1)
stem (n1,y1)
xlabel (‘sampling’)
ylabel(‘Amplitude’)
Result:
The Linearity,time varying property, causality and stability of a given
Discrete System are verified.
Exp. No. 3 UNIT SAMPLE, UNIT STEP, SINUSOIDAL RESPONSE and FREQUENCY
RESPONSE OF THE GIVEN LTI SYSTEM
AIM:
To compute the Unit sample, unit step, sinusoidal response and frequency response of the
given LTI system.
Software Required:
Theory:
A discrete time system performs an operation on an input signal based on predefined
criteria to produce a modified output sign al. The input signal x(n) is the system excitation, and
y(n) is the system response. The transform operation is shown as, If the input to the system is unit
impulse i.e. x(n) = δ(n) then the output of the system is known as impulse response denoted by
h(n) where, h(n) = T[δ(n)]we know that any arbitrary sequence x(n) can be represented as a
weighted sum of discrete impulses.
Program:
%given difference equation y(n)-y(n-1)+.9y(n-2)=x(n);
b=[1];
a=[1,-1,.9];
n =0:3:100;
%generating impulse signal
x1=(n==0);
%impulse response
h1=filter(b,a,x1);
subplot(3,1,1);
stem(n,h1);
xlabel('n');
ylabel('h(n)');
title('impulse response');
%generating step signal
x2=(n>0);
% step response
s=filter(b,a,x2);
subplot(3,1,2);
stem(n,s);
xlabel('n');
ylabel('s(n)')
title('step response');
%Frequency response:
num=[1 0]
den=[1 -0.5]
ww=-pi:0.01:pi
[H]=freqz(num,den,ww)
Figure;
Plot(ww,abs(H))
Result
The Unit sample, unit step and sinusoidal response of the given LTI system is computed verified.
Output:
Aim:
To study the sampling effect, A/D conversion and/A conversion.
Software Required:
Theory:
Sampling Theorem:
A bandlimited signal can be reconstructed exactly if it is sampled at a rate atleast twice the
maximum frequency component in it.
Under sampling :-
If the sampling frequency fs < 2fm, then the sampling is called under sampling.
Effect of under sampling is called aliasing.
Over sampling:-
If the sampling frequency fs ˃ 2fm,,then the sampling is called over sampling. No aliasing
problem occurs.
Critical sampling:-
If the sampling frequency fs = 2fm,,then the sampling is called critical sampling
Program:
Clc;
Clear all;
Close all;
t=-10:.01:10;
T=4;
fm=1/T;
x=cos(2*pi*fm*t);
subplot(2,2,1);
plot(t,x);
xlabel('time');
ylabel('x(t)');
title('continous time signal');
grid;
n1=-4:1:4;
fs1=1.6*fm;
fs2=2*fm;
fs3=8*fm;
x1=cos(2*pi*fm/fs1*n1);
subplot(2,2,2);
stem(n1,x1);
xlabel('time');
ylabel('x(n)');
title('discrete time si
gnal with fs<2fm');
hold on;
subplot(2,2,2);
plot(n1,x1);
grid;
n2=-5:1:5;
x2=cos(2*pi*fm/fs2*n2);
subplot(2,2,3);
stem(n2,x2);
xlabel('time');
ylabel('x(n)');
title('discrete time si
gnal with fs=2fm');
hold on;
subplot(2,2,3);
plot(n2,x2)
grid;
n3=-20:1:20;
x3=cos(2*pi*fm/fs3*n3);
subplot(2,2,4);
stem(n3,x3);
xlabel('time');
ylabel('x(n)');
title('discrete time signal with fs>2fm')
hold on;
subplot(2,2,4);
plot(n3,x3)
grid;
Software Required:
Matlab software 7.0 and above.
Theory:
Program:
clc;
clear all;
close all;
%generates first set of 2000 sa
mples of Gaussian distributed
random numbers
x1=randn(1,2000);
%generates second set
of 2000 samples of Gaussian distributed
random numbers
x2=randn(1,2000);
%plot the joint distribution of
both the sets using dot.
plot(x1,x2,'.');
title('scatter plot of gaussian di
stributed random numbers');
%generates two sets of 2000 samples of uniform distributed
random numbers
x3=rand(1,2000);
x4=rand(1,2000);
figure;
plot(x3,x4,'.');
title('scatter plot of uniform
distributed random numbers');
figure;
subplot(2,1,1);
%plot the histogram graph of x2
Department of ECE, MRCET
BS
Lab Manual
75
hist(x2);
title(' gaussian distribution');
subplot(2,1,2);
%plot the histogram graph of x4
hist(x4);
title(' uniform distribution');
ymu=mean(x2)
ymsq=sum(x2.^2)/length(x2)
ysigma=std(x2)
yvar=var(x2)
yskew=skewness(x2)
ykurt=kurtosis(x2)
Output:
ymu =
0.0172
ymsq =
0.9528
ysigma =
0.9762
yvar =
0.9530
yskew =
-0.0041
ykurt =
2.9381
Output
Exp.No:6a Verification of Kirchhoff’s law
AIM:
To verify the Kirchhoff’s current law and voltage law, both theoretically and experimentally in the
given circuit.
APPARATUS REQUIRED:
CIRCUIT DIAGRAM:
(i) KCL
+ - + -
A B A
A C
100 220
+
(0-50mA) + (0-50mA)
E1=10V A (0-50mA) E1=5V
-
I1
100 I2
F D
100 E 220
(ii)KVL
A 100 B 220 C
+ V - - V +
E1=10V +n + E1=5V
(0-10V) (0-10V)
(0-10V) V 100
-
100 220
F D
The algebraic sum of all the currents at junction is equal to zero (or) the sum of the current flowing
towards the junction is equal to the sum flowing away from junction.
PROCEDURE:
Compare the measured values with the theoretical values obtained through calculations
LOOP2:
220I2 + 200I2 + 100(I2-I1) = -5
-100I1 + 540I2 = -5
From the two equations
We have,
I1 = 32.34mA
I2 = -3.289mA
I3 = I1 - I2
I3 = 35.529mA
TABULATION: KCL
The algebraic sum of product of the sum of current and resistance and the algebraic sum of the voltage
sources in a closed loop or mesh is zero.
PROCEDURE:
LOOP1:
LOOP2:
LOOP1:
LOOP2:
RESULT:
Thus the Kirchhoff’s current and voltage laws are verified. The theoretical and experimental
values are approximately equal.
To verify the Thevenin’s and Norton’s theorem both experimentally and theoretically
APPARATUS REQUIRED:
(i)THEVENINS THEOREM:
Any two terminals of a network composed of linear passive and active components may be
replaced by an equivalent voltage source in series with an equivalent resistance Rth the voltage source Voc
is equal to potential difference between two terminals caused by active network with no external
resistance. The Rth is an equivalent resistance looking back into network terminals with all sources in
network active.
PROCEDURE:
1. TO FIND IL:
1k
Ω 220Ω
220Ω
+ - 10V
A
B D
0-50mA
1KΩ
1k
Ω
C
2. To FIND Rth:
A
1k
Ω 220Ω
B D
1KΩ
1k
Ω
C
3. To find Voc.
A
1k
Ω 220Ω
10V
B D
1KΩ
1k
Ω
C
4. EQUIVALENT CIRCUIT:
+ -
A
Rth=680.3Ω
(0-50mA)
RL=220Ω
VOC=-3.2V
TABULATION-THEVENINS THEOREM:
1000∗220 1000∗1000
Rth = +
1000+220 1000+1000
Rth = 680.32Ω
By solving we get,
I1 = 3.22mA
I2 = 6.775mA
I3 = 14.33mA
𝑉𝑜𝑐
IL = 𝑅𝑡ℎ+𝑅𝑙
3.2
= 226+680.3
IL = 3.55mA
VOC = VB - VD
𝑅4 𝑅3
= 𝑅1+𝑅4
* 10 – 𝑅2+𝑅3 ∗ 10
VOC = -3.2 V
(ii)NORTONS THEOREM:
Any two terminals of a network composed of linear passive and active elements may be replaced
by an equivalent constant current source.
The current source is equal to short circuit current developed when terminals of original network
and short circuited the parallel back into original network terminal with all sources in network mode
inactive.
PROCEDURE:
1. TO FIND IL:
A
1k
Ω 220Ω
220Ω
+ - 10V
A
B D
0-50mA
1KΩ
1k
Ω
C
2. TO FIND Rth:
A
1k
Ω 220Ω
B D
1KΩ
1k
Ω
C
3. TO FIND ISC:
A
1k
Ω 220Ω
+ - 10V
A
B D
0-50mA
1KΩ
1k
Ω
C
Rth=680.3Ω 220Ω
ISC=4.6A
+
A (0-50mA)
-
TABULATION-NORTONS THEOREM
1000∗220 1000∗1000
Rth = +
1000+220 1000+1000
Rth = 680.32Ω
V = IR
Isc= V/R = 3.2 / 680.3 = 4.7mA
In ABDA, 1000I1 + 220(I1 – I3) + 220(I1 – I2) = 0
By solving we get,
I1 = 3.22mA
I2 = 6.775mA
I3 = 4.33mA
I1 + IL = I2
IL = 3.5mA
RESULT:
Thus Thevenin’s and Norton’s theorem are verified theoretically and experimentally.
To verify the Superposition theorem, both theoretically and experimentally in the given circuit.
APPARATUS REQUIRED
If any multisource complex networks consisting of bilateral linear elements, the voltage across or
current through any given circuit of element is equal to sum of individual voltage or current produced
independently through that current by each source acting independently when all remaining sources are
replaced by their internal resistance.
PROCEDURE:
CIRCUIT DIAGRAM:
A + - B + - C D
A A
120Ω +
220Ω 120Ω
(0-50mA)
A (0-50mA)
(0-50mA)
E1 -
I1 I2 I3
E2
10V
H G F E
A’ + - B’ + - C’ D’
A A
120Ω 220Ω 120Ω
+
(0-50mA)
A (0-50mA)
(0-50mA)
I1 I2 I3
E2
H’ G’ F’ 10V E’
(0-50mA)
A (0-50mA)
(0-50mA)
E1 -
10V
220Ω 220Ω 220Ω
I1 I2 I3
TABULATION:
1. READINGS FOR CIRCUIT 1:
CIRCUIT 1:
340I1 – 220I2 = 10
-220I1 + 660I2 – 220I3 = 0
-220I1 + 560I3 = 10
By solving we get,
I1 = 45mA
I2 = 24.1mA
I3 = 27.3mA
CIRCUIT 2:
340I1’– 220I2’ = 0
-220I1’ + 660I2’ – 220I3’ = 0
-220I2’ + 560I3’ = 10
By solving we get,
I1 = 5.89mA
I2 = 9.1mA
I3 = 21.4mA
CIRCUIT 3:
340I1’’– 220I2’’ = 10
-220I1’’ + 660I2’’ – 220I3’’ = 0
-220I2’’ + 560I3’’ = 0
By solving we get,
I1 = 39.2mA
I2 = 15mA
I3 = 5.87mA
RESULT:
APPARATUS REQUIRED:
A load resistance being connected to a DC network receives maximum power when the load
resistance is equal to the internal resistance of a source network as sum of load terminals.
CIRCUIT DIAGRAM:
1KΩ 1KΩ
Load Resistance, RL
5V
1KΩ +
F -
D
B) E
TO FIND RTH (CD):
A B C
1KΩ 1KΩ
1KΩ
F D
C) E
TO FIND VOC:
A B C
1KΩ 1KΩ
5V
1KΩ
F
E D
D)
EQUIVALENT CIRCUIT:
RTh=1.5K
+
2.5V A (0-50mA)
-
PROCEDURE:
PL = I2L * RL
4. Draw the graph by taking PL on y-axis and RL on x-axis where maximum power transfer takes
place.
THEORETICAL:
Remove the resistor RL across CD and calculate the resistance Rth by looking back from the open
terminals with all voltage levels removed and replaced by internal resistance.
Draw the Thevenin’s equivalent circuit as shown in (D) and calculate the values of IL from
Thevenins equivalent circuit.
Find the value of maximum power transfer theoretically from Thevenins equivalent circuit using
formula IL2 * RL * Voc / Rth
THEORETICAL CALCULATION:
1∗1
Rth = 1 + 1+1
Rth = 1.5KΩ
1∗5
Vth = 2.5
= 2.5V
IL = Vth/ 3*103
IL = 0.833mA
TABULATION:
Maximum power:
REFERENCE:
(ii)RECIPROCITY THEOREM:
The ratio of excitation voltage to response current with a single excitation applied at on one
terminal pair and the response observed at the other is the same when the pair of terminals is interchanged
for excitation and response
PROCEDURE
B
A C
1KΩ 1KΩ
1KΩ
1KΩ
+
5V
A (0-10mA)
-
F D
E
A + - B C
A
1KΩ 1KΩ
(0-10mA)
(0-30V)
1KΩ
F E D
THEORETICAL CALCULATION
LOOP1:
2I1 – I2 = 5
-I1 + 3I2 = 0
I1 = 3mA
I2 = 1mA
LOOP2:
2I1’– I2’ = 0
-I1’ + 3I2’ = 5
I1’= 1mA
I2’ = 2mA
TABULATION
RESULT:
Thus the Maximum power transfer and Reciprocity theorem are verified both experimentally and
theoretically.
10. Simulation of DC Circuit
Aim:-To study the DC transient response of RLC series circuit.
Apparatus Required:-
iii. P-SPICE software
iv. Computer kit
Theory:-
When a network is switched from one condition to another by change in applied voltage
or by change in anyone of the circuit elements during a period of time .Branch current and
voltage change from the initial value to new values. The time interval is called “Transient
Period”. The response or output of the network during transient period is called transient
response of the network.
Circuit Diagram:-
R=2Ω L=50uh C=10uf
1 2
3
=
vinDC
Program:-
RLC series circuit for step input
Vin 1 0 pwl (0 0 1ns 1v 1ms 1v)
R 1 2 2 ohm
L 2 3 50μH
C 3 0 10μF
Tran 1μs 400μs
64
Probe
end
Model Graph:-
Result:-
Name Current
PROCEDURE:-
1. Make connections as per the circuit diagram.
PROCEDURE:-
1. Make connections as per the circuit diagram
2. Give 2V Peak to peak square wave supply through function generator with suitable
frequency.
3. Take out put across inductor in RL Circuit, across capacitor in RC Circuits.
4. Calculate the time constant from CRO.
5. For deferent values of T and V Calculate corresponding (L/R) Values.
6. Compare the time constant theoretically and practically.
OBSERVATIONS:-
To experimentally find the resonant frequency of series and parallel RLC circuit and to compare
with the value predicted using the inductance and capacitance in the circuit.
APPARATUS REQUIRED:
CRO
1KΩ
+
-
100mH
Signal Generator
+
0.1µF
PARALLEL RESONANCE CIRCUIT:
0.1µF
1KΩ 100mH
+ +
CRO
- 27Ω
THEORY:
An RLC circuit is an electrical circuit consisting of a resistor (R) and an inductor (L) and a capacitor
(C) connected in series or parallel.
A resonant circuit also called tuning circuit consists of an inductor and capacitor together with a
voltage or current source. It is one of the best and most important circuits used in electronics. For
example a resonant circuit in one of its many forms allows us to select a desired radii or television
signals flow from the vast member of signals that are member of signal around us any time. The
resonant frequency can be formed from the condition
XL = X C
2πf0L = 1/2πf0C
f0 = 1/2π LC
FREQUENCY RESPONSE:
It is the plot of the magnitude of the output voltage of a resonant circuit as function of frequency. The
response of the curve starts at zero, reaches a maximum value in the vicinity of a natural frequency and
that drops again to zero as frequency becomes infinite.
PROCEDURE:
THEORITICAL:
At resonance f0 = 1/2π LC
=1592.3 Hz
TABULATION:
V0
0.707 V0
Current
V0
Vmin
RESULT:
Thus the characteristics of series and parallel resonance circuit are studied.
Exp.No:10 Determination of self and mutual inductances and
coupling coefficient of coupled coils.
Aim:-
To determine self, mutual inductance of coils a, b, c and coupling co-efficient between
„a‟ and „b‟ and „a‟ and „c‟.
1 Voltmeter MI 0-300V 1
2 Voltmeter MI 0-150V 1
3 Voltmeter MI 0-75V 1
Ammeter MI 0-1A 1
Apparatus Required:-
Circuit Diagram:-
48
Procedure:-
1. Connections are made as per the circuit diagram
2. Ensure that the Auto Transformer at 0V position before giving supply to the circuit.
3. Switch ON the supply and apply Va =90V and note down Vb, Vcand the ammeter
readings.
4. Similarly apply Va=195 and Va=200V and note down Vb, Vc and ammeter readings at
„V‟ excitation of coil „A‟.
5. Above steps are repeated by exciting coil B (95,100,105) and C (45, 50, 55).
Precautions:-
1. Before switching a coil which should be kept at „0‟ position.
2. While exciting a coil should be ensure that the applied voltage does not exceed its
reading.
Observation and Calculation:-
Average of Laa=
Average of Mba=
Average of Mca=
Average of Lbb=
Average of Mab=
Average of Mcb=
3) When coil C excited
Average of Lcc=
Average of Mac=
Average of Mbc=
Coupling Coefficient:-
Kab= (Mab+Mba)/2 .
Kca= (Mca+Mac)/2 .
Kbc= (Mbc+Mcb)/2 .
Result:-
The self inductance, mutual inductance of a, b, c and coefficient of coupling between „a‟ and „b‟
and „a‟ and „c‟ is determined.
Coupling Coefficient between „a‟ and „b‟ is Kab=
Coupling Coefficient between „b‟ and „c‟ is Kbc=
Coupling Coefficient between „c‟ and „a‟ is Kca=
Exp. 11. Power and power factor measurement in three phase circuits by
two wattmeter
method.
1. Power drawn by a 3-phase inductive load using 2-wattmeter method and 2. To calculate the
power factor of the load from the readings of two wattmeters.
APPRATUS REQUIRED:
1. Wattmeter Dynamometer 2
5/10A,
150/300/600 V
THEORY:
Power consumed by a 3-phase balanced or unbalanced load (star connected) can be measured by
using 2-wattmeters properly connected in the circuit.The current coil of the wattmeter are
connected in series with the load in any two line. Where as the pressure coils are connected
between these two lines and the third line. The phasor diagram of this circuit assuming balanced
lagging load has been shown in the figure. Under running conditions the power consumed by the
three phase system is the sum o the two individual wattmeters.
Where,
When the load power factor is less then 0.5 then wattmeter 2 will show the correct deflection and
first wattmeter will show the reverse deflection. In the first wattmeter the current coil or voltage
coil connection is reversed. Thus the wattmeter pointer direction is corrected. The net power is
obtained by adding the two wattmeter readings.
(𝑊1 − 𝑊2)
Ø = 𝑡𝑎𝑛−1 √3 ×
(𝑊1 + 𝑊2)
(𝑾𝟏 − 𝑾𝟐)
𝒄𝒐𝒔Ø = 𝐜𝐨𝐬[𝒕𝒂𝒏−𝟏 √𝟑 × ]
(𝑾𝟏 + 𝑾𝟐)
CIRCUIT DIAGRAM:
PROCEDURE:
OBSERVATION TABLE:
PRECAUTIONS:
Result:
1. Net power consumed, P =___________Watts.
2. Power Factor = _____________________
AIM: To obtain experimentally Z parameters and Y parameters of a given two port network.
APPARATUS:
S.No Name of the Equipment Type Range Quantity
1 Bread board 1
2 DMM MC/Digital
Regulated power supply
unit(RPS)
3 DC 30v/2A 2
PROCEDURE:
1. Open Circuiting Output Terminals (I2 = 0):
Connections are made as per the circuit diagram shown in fig (2). Output terminals are kept
Open via a voltmeter. Supply is given to input port. Note the readings of ammeter as I1 and
Voltmeter as V2.
2. Short circuiting output terminals (V2 = 0):
Connections are made as per the circuit diagram shown in fig (4). Output terminals are short
circuited via an ammeter. Supply is given to input port. Note the readings of ammeters as I1 and
I2.
When I1=0
V1 I2 V2
S.No.
When I2=0
V1 I2 V2
S.No.
When V1=0
I2 I1 V2
S.No.
When V2=0
I2 I1 V2
S.No.
RESULT TABLE:
Z Parameters Y Parameters
Theoretical
Practical
H-Parameters:-
Theoretical Procedure: - In this V2 and I1 are taken as independent variables, while V1 and I2
are taken as dependent variables.
V1=h11I1+h12V2
I2=h21I1+h22V2
Analytical Procedure:-
For determining h11 and h21
Thus port 2 is shorted. Therefore
V2=0 I1=V1/Req
h12=V1/V2 h22=I2/V2
h11
h12
h21
h22