Scribd
Upload
4 views

Tutorial 2

Uploaded by

jayakailas34
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
4 views

Tutorial 2

Uploaded by

jayakailas34
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 7

%example 1

syms I1 I2 I3 I4 I5 I6
eqns = [I1 - I2 + I3 == 0, I1 - I2 + I4 == 0, I4 - I5 + I6 == 0, I3 - I5 +
I6 == 0, 3*I1 + 2*I2 == 14, 2*I2 + 4*I3 + 2*I4 + I5 == 25, I5 + I6 == 8]

eqns =

[A,B] = equationsToMatrix(eqns)

A =

B =

P = [A,B]

P =

k1 = rank(A)

k1 =
6

k2 = rank(P)

k2 =
6

1
n = size(A,2)

n =
6

if k1==k2 && k1==n


disp("System has unique solution")
sol=linsolve(A,B)
elseif ((k1==k2)&& k1<n)
disp("System has infinitely many solutions")
sol=solve(eqns,[I1,I2,I3,I4,I5])
else k1~=k2
disp("System has no solution")
end

System has unique solution


sol =

%example 2
syms i1 i2 i3
eqns = [3*i1 + 2*i2 == 7, 2*i2 + 4*i3 == 8, i1 - i2 + i3 ==0]

eqns =

[A,B] = equationsToMatrix(eqns)

A =

B =

P = [A,B]

P =

2
k1 = rank(A)

k1 =
3

k2 = rank(P)

k2 =
3

n = size(A,2)

n =
3

if k1==k2 && k1==n


disp("System has unique solution")
sol=linsolve(A,B)
elseif ((k1==k2)&& k1<n)
disp("System has infinitely many solutions")
sol=solve(eqns,[x,y])
else k1~=k2
disp("System has no solution")
end

System has unique solution


sol =

%example 3
syms i1 i2 i3
eqns = [4*i1 + i2 == 16, i2 + 4*i3 == 8, i1 - i2 + i3 == 0]

eqns =

[A,B] = equationsToMatrix(eqns)

A =

B =

3
P = [A,B]

P =

k1 = rank(A)

k1 =
3

k2 = rank(P)

k2 =
3

n = size(A,2)

n =
3

if k1==k2 && k1==n


disp("System has unique solution")
sol=linsolve(A,B)
elseif ((k1==k2)&& k1<n)
disp("System has infinitely many solutions")
sol=solve(eqns,[x,y])
else k1~=k2
disp("System has no solution")
end

System has unique solution


sol =

%example 4
syms x1 x2 x3 x4 x5
eqns = [x1 + x2 == 20, x3 - x4 == -20, x5 - x4 == -10, x2 + x3 == 20]

eqns =

[A,B] = equationsToMatrix(eqns)

4
A =

B =

P = [A,B]

P =

k1 = rank(A)

k1 =
4

k2 = rank(P)

k2 =
4

n = size(A,2)

n =
5

if k1==k2 && k1==n


disp("System has unique solution")
sol=linsolve(A,B)
elseif ((k1==k2)&& k1<n)
disp("System has infinitely many solutions")
sol=solve(eqns,[x1,x2,x3,x4])
else k1~=k2
disp("System has no solution")
end

System has infinitely many solutions


sol = struct with fields:
x1: x5 - 10
x2: 30 - x5
x3: x5 - 10
x4: x5 + 10

5
%example 5
syms x1 x2 x3 x4
eqns = [x1 - x2 == 200, x4 -x2 == 100, x3 - x4 == 200, x3 - x1 == 100]

eqns =

[A,B] = equationsToMatrix(eqns)

A =

B =

P = [A,B]

P =

k1 = rank(A)

k1 =
4

k2 = rank(P)

k2 =
4

n = size(A,2)

n =
4

if k1==k2 && k1==n


disp("System has unique solution")
sol=linsolve(A,B)
elseif ((k1==k2)&& k1<n)

6
disp("System has infinitely many solutions")
sol=solve(eqns,[x1,x2,x3,x4])
else k1~=k2
disp("System has no solution")
end

System has unique solution


sol =

You might also like