Circuit Simulation Laboratory
Assignment 3
Using MATLAB solve the following circuit analysis problems:
a) Use the mesh analysis to find the current flowing through the resistor RB of fig.3(a) . In
addition, find the power supplied by the 10-volt voltage source.
b) In fig. 3(b), as RL varies from 0 to 50KΩ, plot the power dissipated by the load. Verify
that the maximum power dissipation by the load occurs when RL is 10 KΩ
c) For the circuit shown in fig.3(c), find the nodal voltages V 1, V 2 , and V3 .
Fig.3(a)
1
Circuit Simulation Laboratory
Fig.3(b)
Fig.3(c)
Part A).
Fig.3(a)
Applying mesh analysis to loop1,
10 = 10 * ( I1 - I2 ) + 30 * ( I2- I3 )
2
Circuit Simulation Laboratory
40*I1 - 10*I2 - 30*I3 = 10
4*I1 - 1*I2 - 3*I3 = 10 ................... Eqn(1)
Applying mesh analysis to loop2,
0 = 15*I2 + 5 * ( I2 - I3 ) + 10 * ( I2 – I1 )
-10*I1 + 30*I2 - 5*I3 = 0
-2*I1 + 6*I2 - 1*I3 = 0 ................... Eqn(2)
Applying mesh analysis to loop3,
0 = 30*I3 + 5*( I3 – I2 ) + 30*( I3 – I1 )
-30*I1 – 5*I2 + 65*I3 = 0
-6*I1 – 1*I2 + 13*I3 = 0 ................... Eqn(3)
From eqn (1), (2), (3)
Representing in matrix form:
4 -1 -3 I1 1
-2 6 -1 I2 = 0
-6 -1 13 I3 0
Hence, V = Z*I
Where , V= Voltage matrix
I= Current matrix
Z= Impedance matrix
Therefore we get , I = Z-1 * V
3
Circuit Simulation Laboratory
%Use Mesh analysis to find current in resistor Rb. Also find power
supplied by %voltage source.
Z=[4 -1 -3 ; -2 6 -1 ; -6 -1 13 ]; %Impedance matrix
V=[1 ; 0 ; 0]; %Voltage matrix
I=inv(Z)*V; %Current matrix
Ir=I(3)-I(2); %Current in resistor Rb
P=10*I(1); %Power supplied by voltage source
Part B).
Maximum Power Transfer theorem:-
The maximum power transfer theorem states that, to obtain maximum external power from a
source with a finite internal resistance, the resistance of the load must equal the resistance of the
source as viewed from its output terminals.
Fig.3(b)
%Maximum Power Transfer Theorem
R=0:1000:50000; %Rl varies from 0K to 50kohm
Vs=10; %Voltage source
Rs=10000; %Source resistance
k=length(R);
for i=1:k
P(i)=(Vs/(Rs+R(i)))^2*R(i);
end
4
Circuit Simulation Laboratory
plot(R,P,'*');
title('Maximum Power Transfer Theorem');
xlabel('Load in Ohms');
ylabel('Power in Watts');
[val,ind]=max(P);
Maxpower= val %To display maximum power dissipated by load
Matchresistor= R(ind) %To display value of Rl for which maximum
%power is dissipated i.e. Rs=Rl
Output:
5
Circuit Simulation Laboratory
Fig 3.1 Maximum power Transfer theorem
Part C).
6
Circuit Simulation Laboratory
Fig.3(c)
Applying nodal analysis at node 1,
(V1 -V2) + ( V1 - V3) = 5
10 20
0.15 V1 - 0.1 V2 + 0.05 V3 = 5 ................... Eqn(1)
Applying nodal analysis at node 2,
(V2 -V1) + ( V2 - V3) + V2 = 0
10 40 50
-0.1 V1 + 0.145 V2 - 0.025 V3 = 0 ................... Eqn(2)
Applying nodal analysis at node 3,
(V3 -V1) + ( V3 - V2) = 2
20 40
-0.05 V1 - 0.025 V2 + 0.075 V3 = 2 ................... Eqn(3)
From eqn (1), (2), (3)
Representing in matrix form:
7
Circuit Simulation Laboratory
0.15 -0.1 0.55 V1 5
-0.1 0.145 -0.025 V2 = 0
-0.05 -0.025 0.075 V3 2
Hence, I = Y*V
Where , V= Voltage matrix
I= Current matrix
Y= Admittance matrix
Therefore we get , V = Y-1 * I
%Circuit analysis node voltages
y=[0.15 -0.1 -0.05 ; -0.1 0.145 -0.025 ; -0.05 -0.025 0.075];
i=[5 0 2 ]';
v=inv(y)*i
Output:
8
Circuit Simulation Laboratory
9
Circuit Simulation Laboratory
INFORMATION:
Length:
INV:
10
Circuit Simulation Laboratory
Max:
Conclusion:-
Using MATLAB we have solved the circuit analysis problems. We have calculated the value of
voltages and current by using Mesh and nodal analysis respectively. For given circuit the load
resistance for which the maximum power will be delivered to load is found out and the graph of
maximum power dissipated by the load is plotted.
11