Basic Operations On Matrices: Ij MXN Ij MXN
Basic Operations On Matrices: Ij MXN Ij MXN
THEORY :
1. ADDITION OF MATRICES :
If A[aij]mxn and B[bij]mxn are two matrices of the same order, then
their sum A + B is a matrix, and each element of that matrix is the sum
of the corresponding elements, i.e. A + B = [a ij + bij]mxn
𝑎1 𝑏1 𝑎2 𝑏2 𝑎1 + 𝑎2 𝑏1 + 𝑏2
[ ] + [ ] = [ ]
𝑐1 𝑑1 𝑐2 𝑑2 𝑐1 + 𝑐2 𝑑1 + 𝑑2
Properties of Matrix Addition:
2. SUBTRACTION OF MATRICES :
A – B = A + (-B)
Consider the two matrices, A and B, of order 2 x 2. Then, the difference
is given by:
𝑎1 𝑏1 𝑎2 𝑏2 𝑎1 − 𝑎2 𝑏1 − 𝑏2
[ ] - [ ] = [ ]
𝑐1 𝑑1 𝑐2 𝑑2 𝑐1 − 𝑐2 𝑑1 − 𝑑2
3. MATRIX MULTIPLICATION :
(AB)ij = Cij = ∑𝑛
𝑟 = 1(𝑎)𝑖𝑟. (𝑏)𝑟𝑗
AB ≠ BA.
(b) Matrix multiplication is associative, i.e. (AB)C = A(BC).
A.O = O.A = O
EXP NO. 21011A0253 PAGE NO.
i.e. the product of the matrix with a null matrix is always a null
matrix.
(f) If AB = 0 (It does not mean that A = 0 or B = 0, again, the product
of two non-zero matrices may be a zero matrix).
(g) If AB = AC , B ≠ C (Cancellation law is not applicable).
(h) There exists a multiplicative identity for every square matrix, such
as
AI = IA = A
PROCEDURE :
1. Open MATLAB => New script
2. Type the program & save the file with .m extension.
3. Click Run & view the output in command window.
PROGRAM :
clc;
close all;
A = [5 2 -9; -1 9 4; 3 -4 7];
B = [4 6 -2; -9 -1 -8; 2 -5 0];
disp('Matrix A :');
A
disp('Matrix B :');
B
C = A + B;
disp('Sum of A & B is');
C
D = A - B;
disp('Difference of A & B is');
D
E = A*B;
disp('Product of A & B is');
E
EXP NO. 21011A0253 PAGE NO.
OUTPUT :
Matrix A :
A =
5 2 -9
-1 9 4
3 -4 7
Matrix B :
B =
4 6 -2
-9 -1 -8
2 -5 0
Sum of A & B is
C =
9 8 -11
-10 8 -4
5 -9 7
Difference of A & B is
D =
1 -4 -7
8 10 12
1 1 7
Product of A & B is
E =
-16 73 -26
-77 -35 -70
62 -13 26
EXP NO. 21011A0253 PAGE NO.
CALCULATIONS :
5 2 −9 4 6 −2
A = −1 9 4 , B = −9 −1 −8
3 −4 7 2 −5 0
5 − 4 2 − 6 −9 + 2 1 −4 −7
D = A - B = −1 + 9 9+1 4+8 = 8 10 12
3 − 2 −4 + 5 7 − 0 1 1 7
E = A x B =
−16 73 −26
= −77 −35 −70
62 −13 26
RESULT
THEORY :
∞
F(s) = L[f(t)] = ∫0 𝑓(𝑡)𝑒^(−𝑠𝑡). 𝑑𝑡
An integral formula for the inverse Laplace transform, called the
Mellin's inverse formula, the Bromwich integral, or the Fourier–Mellin
integral, is given by the line integral:
PROCEDURE :
PROGRAM :
clc;
close all;
syms f t w s;
f1 = 4*t^3;
z = laplace(f1);
disp('Laplace transform of f1 ');
z
f2 = 3*sin(2*w*t);
v = laplace(f2);
disp('Laplace transform of f2 ');
v
f3 = 5*exp(-3*t);
m = laplace(f3);
disp('Laplace transform of f3 ');
m
y1 = ilaplace(z);
disp('Inverse Laplace transform of z ');
y1
figure(1);
ezplot(y1);
y2 = ilaplace(v);
disp('Inverse Laplace transform of v ');
y2
figure(2);
ezplot(y2);
y3 = ilaplace(m);
disp('Inverse Laplace transform of m ');
y3
figure(3);
ezplot(y3);
EXP NO. 21011A0253 PAGE NO.
OUTPUT :
Laplace transform of f1
z =
24/s^4
Laplace transform of f2
v =
(6*w)/(s^2 + 4*w^2)
Laplace transform of f3
m =
5/(s + 3)
OUTPUT GRAPHS :
EXP NO. 21011A0253 PAGE NO.
CALCULATIONS :
3
1. f(t) = 4t
3! 24
L{f(t)} = L{4t3} = 4. =
s^4 𝑠^4
2. f(t) = 3*sin(2wt)
2w 6w
L{f(t)} = L{3*sin(2wt)} = 3. =
s^2 + (2w)^2 s^2 + 4.w^2
3. f(t) = 5.𝒆−𝟑𝐭
1 5
L{f(t)} = L{5.𝒆−𝟑𝐭 } = 5. =
s−(−3) s+3
RESULT :
Hence the Laplace transform & inverse Laplace transform of given signals
are determined using MATLAB.
EXP NO. 21011A0253 PAGE NO.
THEORY :
Poles and Zeros of a transfer function are the frequencies for which
the value of the denominator and numerator of transfer function becomes
infinite and zero respectively. The values of the poles and the zeros of
a system determine whether the system is stable, and how well the system
performs. Control systems, in the most simple sense, can be designed
simply by assigning specific values to the poles and zeros of the system.
𝐶(𝑠)
G(s) =
𝑅(𝑠)
Z-Transform :
The bilateral or two-sided Z-transform of a discrete-time signal
x[n] for n ≥ 0 is the formal power series defined as:
, where z = r.𝐞−𝐣𝛚
PROCEDURE :
EXP NO. 21011A0253 PAGE NO.
PROGRAM :
clc;
close all;
num = input("Enter the numerator coefficients : ");
den = input("Enter the denominator coefficients : ");
G = tf(num,den);
p = roots(den);
z = roots(num);
G
figure(1);
sgrid;
pzplot(G);
figure(2);
zplane(z,p);
OUTPUT :
OUTPUT GRAPHS:
EXP NO. 21011A0253 PAGE NO.
CALCULATIONS :
1. s - Plane :
For zeros,
s2 + 7s + 18 = 0
For poles,
s3 - 4s2 - 16s + 30 = 0
2. z – Plane :
For zeros,
z2 + 7z + 18 = 0
For poles,
z3 – 4z2 – 16z + 30 = 0
RESULT :
Hence the poles & zeros of given transfer function are determined &
located on s-plane & z-plane using MATLAB.
EXP NO. 21011A0253 PAGE NO.
CIRCUIT PARAMETERS :
THEORY :
During the positive half cycle of the input AC supply (0 to π), the
polarity of the secondary voltage across terminal A is positive with
respect to terminal B. This causes the diodes D1 and D2 to be forward-
biased and diodes D3 and D4 to be reversed-biased.
EXP NO. 21011A0253 PAGE NO.
During the negative half cycle of the input AC supply (π to 2π), the
polarity of the secondary voltage across terminal B is positive with
respect to terminal A.
In both cycles, the load current flows through the positive half
cycles of the input AC voltage and in the same direction through the load
resistor RL. The polarity of the voltage across RL is the same as the
direction of load current through conducting states of diodes D1, D2, and
D3, D4. The output pulse can either be completely positive or negative.
1 𝜋 2𝑉𝑚
i. Average load voltage VL = ∫ 𝑉𝑚
𝜋 0
𝑠𝑖𝑛(𝑤𝑡) 𝑑(𝑤𝑡) =
𝜋
,
𝑉𝐿
& load current IL = .
𝑅
𝐼𝐿
ii. Current through each diode Id =
2
iii. Average Output Power Po = VL*IL
8
iv. Efficiency η = = 0.8105 or 81.05 %
𝜋^2
𝑉(𝑅𝑀𝑆)
v. Ripple Factor γ = √( )^2 − 1
𝑉(𝐷𝐶)
𝑉𝐿 ∗ 𝐼𝐿
vi. Power Factor cosɸ =
𝑉(𝑅𝑀𝑆) ∗ 𝐼(𝑅𝑀𝑆)
PROCEDURE :
1. Open MATLAB => SIMULINK => Blank Model & create model as per circuit
diagram.
3. Add Mean block with supply fundamental frequency to output voltage &
current & display in Display block.
4. Add Powergui block, then save the model & run with desired stop time
to view output in scopes.
EXP NO. 21011A0253 PAGE NO.
OUTPUT GRAPHS :
EXP NO. 21011A0253 PAGE NO.
OUTPUT GRAPHS :
EXP NO. 21011A0253 PAGE NO.
OUTPUT GRAPHS :
EXP NO. 21011A0253 PAGE NO.
OUTPUT GRAPHS :
EXP NO. 21011A0253 PAGE NO.
CALCULATIONS :
For 1ɸ rectifier with R-load without filter,
2𝑉𝑚 2 𝑥 325.269
Load Voltage, VL = = = 207.07 V
𝜋 𝜋
𝑉𝐿 207.07
Load Current, IL = = = 20.707 A
𝑅 10
𝐼𝐿 20.707
Current through each diode, Id = = = 10.354 A
2 2
𝑉(𝑅𝑀𝑆) 230
Ripple Factor γ = √( )^2 − 1 = √( )^2 − 1 = 1.211
𝑉(𝐷𝐶) 207.07
𝑉𝐿 ∗ 𝐼𝐿 207.07 𝑥 20.707
Power Factor cosɸ = = = 0.81
𝑉(𝑅𝑀𝑆) ∗ 𝐼(𝑅𝑀𝑆) 230 𝑥 23
RESULT :
Hence the operation of 1ɸ diode bridge rectifier is studied using MATLAB
Simulink & the output waveforms are obtained.
EXP NO. 21011A0253 PAGE NO.
CIRCUIT PARAMETERS :
THEORY :
Smaller Filter Size: Due to the reduced output ripple, the size of the
output filter components can be smaller, which in turn leads to a more
compact overall design.
PROCEDURE :
1. Open MATLAB => SIMULINK => Blank Model & create model as per circuit
diagram.
3. Add Mean block with supply fundamental frequency to output voltage &
current & display in Display block.
4. Add Powergui block, then save the model & run with desired stop time
to view output in scopes.
EXP NO. 21011A0253 PAGE NO.
OUTPUT GRAPHS :
EXP NO. 21011A0253 PAGE NO.
OUTPUT GRAPHS :
EXP NO. 21011A0253 PAGE NO.
OUTPUT GRAPHS :
EXP NO. 21011A0253 PAGE NO.
OUTPUT GRAPHS :
EXP NO. 21011A0253 PAGE NO.
CALCULATIONS :
For 3ɸ rectifier with R-load without filter,
RESULT :
Hence the operation of 3ɸ diode bridge rectifier is studied using MATLAB
Simulink & the output waveforms are obtained.
EXP NO. 21011A0253 PAGE NO.
CIRCUIT PARAMETERS :
THEORY :
The total instantaneous power absorbed by the three loads Z1, Z2 and
Z3, is equal to the sum of the powers measured by the two wattmeters, W1
and W2.
EXP NO. 21011A0253 PAGE NO.
Therefore, the total power measured by the two wattmeters W 1 and W2 will
be obtained by adding the equations (1) & (2)
PROCEDURE :
1. Open MATLAB => SIMULINK => Blank Model & create model as per circuit
diagram.
3. Add Powergui block, then save the model & run with desired stop time
to view output in scopes.
EXP NO. 21011A0253 PAGE NO.
I. 3ɸ BALANCED LOAD :
MODEL DIAGRAM :
EXP NO. 21011A0253 PAGE NO.
CALCULATIONS :
𝑋𝐿 100
=> L = = = 0.318 H = 318 mH
2𝜋𝑓 2𝜋 ∗ 100
𝑉𝑝ℎ 230
IL = IPh = = = 1.627 A
𝑍 √100^2 + 100^2
= 793.82 W
= 793.82 W
CALCULATIONS :
Phase A : R = XL = 100 Ω
𝑋𝐿 100
=> L = = = 0.318 H = 318 mH
2𝜋𝑓 2𝜋 ∗ 50
Phase B : R = XL = 50 Ω
𝑋𝐿 50
=> L = = = 0.159 H = 159 mH
2𝜋𝑓 2𝜋 ∗ 50
Phase C : R = XL = 200 Ω
𝑋𝐿 200
=> L = = = 0.637 H = 637 mH
2𝜋𝑓 2𝜋 ∗ 50
𝑉𝑎 230 ∠0°
IAB = IA = = = 1.627 ∠-45° A
𝑍𝑎 100√2 ∠45°
𝑉𝑐 230 ∠−240°
ICB = IC = = = 0.813 ∠75° A
𝑍𝑐 200√2 ∠45°
= (480.6 + j 709.9) W
RESULT :
Hence the measurement of active & reactive power of 3ɸ circuit for
balanced & unbalanced loads using two-wattmeter method is demonstrated
using MATLAB Simulink & the results are verified successfully.
EXP NO. 21011A0253 PAGE NO.
CIRCUIT PARAMETERS :
1. Va = 160 V
2. Ra = 0.5 Ω
3. La = 0.1 H
4. J = 5 kg-m2
5. B = 0.01 N-m/(rad/s)
𝐿𝑎
6. Electrical time constant, Ta =
𝑅𝑎
7. Ke.phi (ke.ɸ) = 1.6
𝐽
8. Mechanical time constant, Tm =
𝐵
9. Repeating Sequence Stair :
Time = [0 1 2]
Output Value = [1 2 3]
Sample time = 2
THEORY :
Tm ∝ ɸ.Ia
Tm = ke.ɸ.Ia
The motor back EMF is proportional to speed, i.e.
𝑑𝛳
Eb ∝
𝑑𝑡
EXP NO. 21011A0253 PAGE NO.
𝑑𝛳
=> Eb = kb.
𝑑𝑡
𝑑 2𝛳 𝑑𝛳
Tm = ke.ɸ.Ia = J. 2 + B. + TL
𝑑𝑡 𝑑𝑡
where
J = Moment of Inertia of motor & load referred to motor shaft
B = Coefficient of Friction of motor & load referred to motor shaft
ϴ = Angular displacement of motor shaft
PROCEDURE :
1. Open MATLAB => SIMULINK => Blank Model & create model as per circuit
diagram.
3. Save the model & run with desired stop time to view output in scopes.
EXP NO. 21011A0253 PAGE NO.
MODEL DIAGRAM :
EXP NO. 21011A0253 PAGE NO.
OUTPUT GRAPH :
RESULT :
Hence the characteristics of feedback control systems is studied upon
simulating closed loop control of DC motor by using MATLAB Simulink.
EXP NO. 21011A0253 PAGE NO.
DESIGN OF 1ɸ AC - DC CONVERTERS
AIM : To design a 1ɸ AC-DC converter & study its operation for R & R-L
loads.
CIRCUIT PARAMETERS :
THEORY :
𝑉𝑚
Average load voltage, Vo = (1 + 𝑐𝑜𝑠 𝛼)
𝜋
𝑉𝑚 𝑠𝑖𝑛 2𝛼
RMS load voltage, VRMS = √𝜋−𝛼 +
√2𝜋 2
2. 1ɸ Full-wave Full-Converter :
1ɸ full converters utilize four SCRs to deliver a constant DC output
voltage. Like dual converters, full converters offer an uninterrupted
output voltage and support bidirectional power flow. With the
implementation of four SCRs, the full converter is capable of providing a
continuous output voltage for both positive and negative half-cycles of
the AC input voltage.
A 1ɸ full converter can operate in two quadrants: the first quadrant
(Q1) for forward motion and the fourth quadrant (Q4) for reverse braking.
During the first quadrant operation, the input voltage is supplied to the
load in the same polarity as the output voltage, while during the fourth
quadrant operation, the input voltage is applied to the load in the
opposite polarity as the output voltage. This allows the load to either
dissipate energy as heat or return energy to the power source, making
full converters ideal for bidirectional power flow applications like DC
motor drives, where both motoring and braking operations are needed.
Full converters can generate higher levels of harmonic distortion in
the input current, which may lead to issues such as increased power
losses and electromagnetic interference. Therefore, it is crucial to
carefully evaluate the trade-offs among different types of rectifiers
when designing power electronics systems to achieve optimal efficiency
and performance.
EXP NO. 21011A0253 PAGE NO.
2 𝑉𝑚
Average load voltage, Vo = 𝑐𝑜𝑠 𝛼
𝜋
PROCEDURE :
1. Open MATLAB => SIMULINK => Blank Model & create model as per circuit
diagram.
3. Add Mean & RMS blocks with supply fundamental frequency to output
voltage & display in Display block.
4. Add Powergui block, then save the model & run with desired stop time
to view output in scopes.
EXP NO. 21011A0253 PAGE NO.
MODEL DIAGRAM :
OUTPUT GRAPHS :
1. Voltage Waveforms :
EXP NO. 21011A0253 PAGE NO.
2. Current Waveforms :
CALCULATIONS :
𝑉𝑚 325.269
Vo,Avg = (1 + 𝑐𝑜𝑠 𝛼) = (1 + 𝑐𝑜𝑠 60°) = 155.30 V
𝜋 𝜋
= 206.296 V
EXP NO. 21011A0253 PAGE NO.
MODEL DIAGRAM :
OUTPUT GRAPHS :
1. Voltage Waveforms :
EXP NO. 21011A0253 PAGE NO.
2. Current Waveforms :
CALCULATIONS :
𝑉𝑚 325.269
Vo,Avg = (1 + 𝑐𝑜𝑠 𝛼) = (1 + 𝑐𝑜𝑠 60°) = 155.30 V
𝜋 𝜋
= 206.296 V
EXP NO. 21011A0253 PAGE NO.
MODEL DIAGRAM :
OUTPUT GRAPHS :
1. Voltage Waveforms :
EXP NO. 21011A0253 PAGE NO.
2. Current Waveforms :
CALCULATIONS :
2 𝑉𝑚 2 ∗ 325.269
Vo,Avg = 𝑐𝑜𝑠 𝛼 = 𝑐𝑜𝑠 60° = 103.54 V
𝜋 𝜋
Vo,RMS = Vs = 230 V
EXP NO. 21011A0253 PAGE NO.
MODEL DIAGRAM :
OUTPUT GRAPHS :
1. Voltage Waveforms :
EXP NO. 21011A0253 PAGE NO.
2. Current Waveforms :
CALCULATIONS :
2 𝑉𝑚 2 𝑥 325.269
Vo,Avg = 𝑐𝑜𝑠 𝛼 = 𝑐𝑜𝑠 60° = 103.53 V
𝜋 𝜋
Vo,RMS = Vs = 230 V
RESULT :
Hence the operation of 1ɸ AC-DC converters are studied by using MATLAB
Simulink & the output waveforms are obtained.
EXP NO. 21011A0253 PAGE NO.
DESIGN OF 1ɸ INVERTERS
AIM : To design a 1ɸ multiple sinusoidal PWM inverter & study its
performance characteristics.
CIRCUIT PARAMETERS :
THEORY :
FULL-BRIDGE INVERTER :
The circuit comprises four diodes & four controlled switches, which
can be BJT, IGBT, MOSFET or Thyristors, play a crucial role in the
inversion process. The power circuit of a 1ɸ full bridge inverter is
constructed with precision, featuring four thyristors labelled T1 to T4,
four diodes D1 to D4 & a two wire DC input power source denoted as Vs.
The four diodes, also known as freewheeling or feedback diodes,
facilitate the redirection of stored energy in the load back to the DC
source, particularly beneficial for non-pure resistive loads. This
feedback mechanism enhances efficiency and stability.
EXP NO. 21011A0253 PAGE NO.
In this mode, T1 & T2 are triggered, while T3 & T4 are OFF. The current
starts flowing through T1 & T2, hence voltage across A & B V AB = +VS.
2. Mode 2 :
EXP NO. 21011A0253 PAGE NO.
In this mode, T3 & T4 are triggered, while T1 & T2 are OFF. The current
starts flowing through T3 & T4, hence voltage across A & B V AB = -VS. The
whole process is repeated for the next cycles.
PROCEDURE :
1. Open MATLAB => SIMULINK => Blank Model & create model as per circuit
diagram.
3. Add Powergui block, then save the model & run with desired stop time
to view output in scopes.
EXP NO. 21011A0253 PAGE NO.
MODEL DIAGRAM :
EXP NO. 21011A0253 PAGE NO.
OUTPUT GRAPHS :
EXP NO. 21011A0253 PAGE NO.
CALCULATIONS :
𝐴𝑟 0.9
Modulation index, m = = = 0.9
𝐴𝑐 1
RESULT :
Hence the 1ɸ multiple sinusoidal PWM inverter is designed using MATLAB
Simulink & the performance characteristics are studied.