1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
|
// To calculate current in each branch using loop analysis.
clc;
clear;
// MESH Equations for the given network.
//3*i1-i2+0*i3=11
//-i1+10*i2-2*i3=0
//0*i1+-2*i2+5*i3=13
//Voltage supplies are 11V and 13V
R=[3 -1 0;-1 10 -2; 0 -2 5];
E=[11;0;13];
// Loop Currents
I=inv(R)*E;
i1=I(1);
i2=I(2);
i3=I(3);
ia=i1-i2; // Assumed direction from Mesh 1
ib=i2-i3; // Assumed direction from Mesh 2
disp('A',ib,'ib (through 2 resistor between 7 ohm and 3 ohm resistor) =','A',ia,'ia(through 1 ohm resistor) =','A',i3,'i3 =','A',i2,'i2 =','A',i1,'i1 =','The Calculated Loop Currents are')
disp('The Negative sign indicates that the assumed direction of flow of current should be reveresed')
// To obtain the magnitude of direction.
if(i1<0)
i1=abs(i1);
end
if(i2<0)
i2=abs(i2);
end
if(i3<0)
i3=abs(i3);
end
if(ia<0)
ia=abs(ia);
end
if(ib<0)
ib=abs(ib);
end
disp('A',i1,'The Current through 2 ohm resistor on the 11V side =')
disp('A',i2,'The Current through 7 ohm resistor =')
disp('A',i3,'The Current through 3 ohm resistor on the 13V side =')
disp('A',ia,'The Current through 1 ohm resistor =')
disp('A',ib,'The Current through 2 ohm resistor between the 7 and 3 ohm resistors =')
|
