Scribd
Upload
11 views

Tutorial 1

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
11 views

Tutorial 1

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/ 9

%example 1

A = [3 3 4; 2 -3 4; 0 -1 1]

A = 3×3
3 3 4
2 -3 4
0 -1 1

rank(A)

ans =
3

%example 2
B = [1 2 5 3; 4 -3 1 0; -2 3 6 7; 0 2 4 6]

B = 4×4
1 2 5 3
4 -3 1 0
-2 3 6 7
0 2 4 6

rank(B)

ans =
4

%example 3
C = [5 -2 3; 0 1 5]

C = 2×3
5 -2 3
0 1 5

rank(C)

ans =
2

%example 4
D = [6 -2 3; 2 1 -2; 3 0 -1]

D = 3×3
6 -2 3
2 1 -2
3 0 -1

rank(D)

ans =
3

%example 5
E = [1 6 4 0; 3 2 2 -1; 1 3 4 -1]

1
E = 3×4
1 6 4 0
3 2 2 -1
1 3 4 -1

rank(E)

ans =
3

%Question 2

% example 1

syms x y z
eqns = [3*x + 4*y - z == 2, 2*x - y + 3*z == 1, x + 2*y + z == 3]

eqns =

[A,B] = equationsToMatrix(eqns)

A =

B =

P = [A,B]

P =

k1 = rank(A)

k1 =
3

k2 = rank(P)

k2 =
3

n = size(A,2)

n =

2
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 2

syms x y z
eqns = [x+y-z==2,2*x+2*y-2*z==4,3*x+3*y-3*z==6]

eqns =

[A,B]= equationsToMatrix(eqns)

A =

B =

P=[A,B]

P =

k1=rank(A)

3
k1 =
1

k2=rank(P)

k2 =
1

n=size(A,3)

n =
1

if (k1==k2 && k1==3)


disp("System has a unique solutions")
solution=linsolve(A,B)
elseif (k1==k2 && k1<3)
disp("System have infinitely many solutions")
solution=solve(eqns,[x,y])
else(k1~=k2)
disp("System does not have solution")
end

System have infinitely many solutions


solution = struct with fields:
x: z + 2
y: 0

%example 3

syms x y z w
eqns = [x+2*y-z+w==4,2*x+4*y-2*z+2*w==8,x-y+z==1]

eqns =

[A,B] = equationsToMatrix(eqns)

A =

B =

P = [A,B]

P =

4
k1 = rank(A)

k1 =
2

k2 = rank(P)

k2 =
2

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)
disp("System has infinitely many solutions")
sol=solve(eqns,[x,y,z])
else k1~=k2
disp("System has no solution")
end

System has infinitely many solutions


sol = struct with fields:
x: 2 - w/3
y: 1 - w/3
z: 0

%example 4

syms x y z
eqns = [5*x + 2*y - z ==3, 2*x + 4*y -2*y == 7, 3*x + 6*y - 3*z == 10]

eqns =

[A,B] = equationsToMatrix(eqns)

A =

B =

5
P = [A,B]

P =

k1 = rank(A)

k1 =
3

k2 = rank(P)

k2 =
4

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

ans = logical
1
System has no solution

%example 5

syms x y z
eqns = [x + 2*y + 2*z == 1, 2*x + 2*y + 3*z == 3, x - y + 3*z == 5]

eqns =

[A,B] = equationsToMatrix(eqns)

6
A =

B =

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 6

7
syms x y z
eqns = [2*x - 3*y + 5*z == 1, 3*x + y - z == 2, x + 4*y - 6*z == 1]

eqns =

[A,B] = equationsToMatrix(eqns)

A =

B =

P = [A,B]

P =

k1 = rank(A)

k1 =
2

k2 = rank(P)

k2 =
2

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 infinitely many solutions


sol = struct with fields:
x: 7/11 - (2*z)/11
y: (17*z)/11 + 1/11

8
9

You might also like