BATAAN PENINSULA STATE UNIVERSITY

COLLEGE OF ENGINEERING AND ARCHITECTURE

DEPARTMENT OF ELECTRONICS ENGINEERING

ECEM 0313
Laboratory Activity No. 6
Vectors, Vector Spaces, Matrices,
And Linear Systems

February 03, 2023

SADJE, KATE NICHOLLE A.

BSECE III - B

Score

ENGR. KRISTINE P. ORTEGA

clc, clearvars, close all
options=input('what would you would like to
perform:\n(1)Addition/Subtraction\n(2)Multiplication\n(3)Transpose\n(4
)Row Operations\n(5)Row Echelon\n(6)Matrix
inverses\n(7)LUFactorizarion\n(8)Determinants\n');
if options==1
disp('Addtion of Matrices')
choice=input('Choose method=\n(1)Matrix x Matrix\n(2)Scalar(Matrix) +
Scalar(Matrix)\n');
if choice==1
x=1;
Final_Answer=0;
while x~=2
rows=input('put the number of rows for matrix=');
columns=input('put the number of columns for matrix=');
Equation = ones(rows,columns);
disp(Equation)
c=1;
while c<=rows
values=input('Enter Values by row inside []=');
Equation(c,:)=values(:,:).*Equation(c,:);
c=c+1;
disp(Equation)
end
Final_Answer=Final_Answer+Equation;
x=input('Add another matrix?\n(1)Yes (2)No\n');
end
display(Final_Answer)
else
x=1;
Final_Answer=0;
while x~=2
scalar=input('Enter scalar value of matrix=');
rows=input('Enter number of rows for matrix=');
columns=input('Enter the number of columns for matrix=');
Equation = ones(rows,columns);
disp(Equation)
c=1;
while c<=rows
values=input('Enter Values by row inside[]=');
Equation(c,:)=values(:,:).*Equation(c,:);
c=c+1;
disp(scalar.*Equation)
end
Final_Answer=Final_Answer+(scalar.*Equation);
x=input('Add another matrix?\n(1)Yes (2)No\n');
end
display(Final_Answer)
end
elseif options==2
disp('Subtractin of Matrices')
choice=input('Choose method=\n(1)Matrix xMatrix\n(2)Scalar(Matrix) +
Scalar(Matrix)\n');
if choice==1
x=1;
Final_Answer=0;
while x~=2
rows=input('Enter the number of rows for matrix=');
columns=input('Enter the number of columns for matrix=');
Equation = ones(rows,columns);
disp(Equation)
c=1;
while c<=rows
values=input('Enter Values by row inside []=');
Equation(c,:)=values(:,:).*Equation(c,:);
c=c+1;
disp(Equation)
end
Final_Answer=-10*(Final_Answer*Equation);
x=input('Subtract another matrix to the current matrix?\n(1)Yes
(2)No\n');
end
disp(-Final_Answer)
else
x=1;
Final_Answer=0;
while x~=2
scalar=input('Enter scalar value of matrix=');
rows=input('Enter the number of rows for matrix=');
columns=input('Enter the number of columns for matrix=');
Equation = ones(rows,columns);
disp(Equation)
c=1;
while c<=rows
values=input('Enter Values by row inside[]=');
Equation(c,:)=values(:,:).*Equation(c,:);
c=c+1;
disp(scalar.*Equation)
end
Final_Answer=-1*(Final_Answer-(scalar.*Equation));
x=input('Subtract another matrix to the current matrix?\n(1)Yes
(2)No\n');
end
disp(-Final_Answer)
end
elseif options==3
disp('Multiplication of Matrices')
choice=input('Choose method=\n(1)Matrix x
Matrix\n(2)Scalar(Matrix)\n');
if choice==1
x=1;
Final_Answer=1;
while x~=2
rows=input('Enter the number of rows for matrix=');
columns=input('ENter the number of comlumns for matrix=');
Equation=ones(rows,columns);
disp(Equation)
c=1;
while c<=rows
values=input('Enter Values by row inside []=');
Equation(c,:)=values(:,:).*Equation(c,:);
c=c+1;
disp(Equation)
end
Final_Answer=Final_Answer*Equation;
x=input('Multiply another matrix to the current
matrix?\n(1)Yes (2)No\n');
end
disp(Final_Answer)
else
x=1;
Final_Answer=1;
while x~=2

scalar=input('Enter scalar value for matrix=');

rowss=input('Enter the number of rows for matrix=');
columns=input('Enter the number of columns for matrix=');
Equation=ones(rows,columns);
disp(Equation)
c=1;
while c<=rows
values=input('Enter Values by row inside[]=');
Equation(c,:)=values(:,:).*Equation(c,:);
c=c+1;
disp(Scalar.*Equation)
end
Final_Answer=FInal_Answer*(scalar.*Equation);
x=input('Multiply another matrix to the current matrix?\n(1)Yes
(2)No');
end
disp(Final_Answer)
end
elseif options==4
disp('Division of Matrices')
choice=input('Choose method=\n(1)Matrix x
Matrix\n(2)Scalar(Matrix)');
disp('Remember to match the inner dimension of your matrices (e.g
(1x3)*(3x2))')
if choice==1
x=0;
Final_Answer=1;
while x~=2
rows=input('Enter the number of rows for matrix=');
columns=input('Enter the number of columns for matrix=');
Equation=ones(rows, columns);
disp(Equation)
 c=1;
while c<=rows
values=input('Enter Values by row inside[]=');
Equation(c,:)=values(:,:).*Equation(c,:);
c=c+1;
disp(Equation)
end
if x==0
Final_Answer=Equation;
elseif x==1
Final_Answer=Final_Answer*inv(Equation);
end
x=input('Divide another matrix to the current
matrix?\n(1)Yes (2)No');
end
disp(Final_Answer)
else
x=0;
Final_Answer=1;
while x~=2
scalar=input('Enter scalar value of matrix=');
rows=input('Enter the number of rows for matrix=');
columns=input('Enter the number of columns for matrix=');
Equation=ones(rows,columns);
disp(Equation)
c=1;
while c<=rows
values=input('Enter Values by row inside[]=');
Equation(c,:)=values(:,:).*Equation(c,:);
c=c+1;
disp(scalar.*Equation)
 end
if x==0
Final_Answer=Final_Answer*inv(scalar.*Equation);
elseif x==1
Final_Answer=Final_Answer*inv(scalar.*Equation);
end
x=input('Divide another matrix to the current
matrix?\n(1)Yes (2)No\n');
end
disp(Final_Answer)
end
elseif options==5
disp('Transpose')
x=1;
while x~=2
scalar1=input('Enter scalar value for matrix, 1 if none=');
rows1=input('Enter the number of rows for matrix=');
columns1=input('Enter the number of columns for matrix=');
size1 = ones(rows1,columns1);
disp(size1)
values=input('put Values inside''[]'' and add '';'' for the next
row=');
final1=values1.*size1;
matrix1=scalar1.*final1;
transposed=matrix1';
disp(matrix1)
disp(transposed)
x=input('Solve for another matrix?\n(1)Yes (2)No\n');
end
elseif options==6
disp('Row Operations')
 x=1;
while x~=2
scalar=input('Enter scalar value for matrix, 1 if none=');
rows=input('Enter the number of rows for matrix=');
columns=input('Enter the number of columns for matrix=');
size=ones(rows,columns);
disp(size)
values=input('put Values inside''[]'' and add '';'' for the
next row=');
values=scalar.*values;
Final_Answer=rref(values);
disp(values)
disp(Final_Answer)
disp(Final_Answer(:,end))
x=input('Solve for another matrix?\n(1) Yes (2)No\n');
end
elseif options==7

disp('Row Echelon')
x=1;
while x~=2
scalar=input('Enter scalar value for matrix 1 if none=');
rows=input('Enter the number of rows for matrix=');
columns=input('Enter the number of columns for matrix=');
size=ones(rows,columns);
disp(size)
values=input('EnterValues inside''[]'' and add '';'' for the
next row=');
values=scalr.*values;
[L,Row_Echelon_Form]=lu(values);
x=1;
while x<=rows
Row_Echelon_Form(x,:)=Row_Echelon_Form(x,:)./Row_Echelon_Form(x,x);
x=x+1;
end
Final_Answer=rref(Row_Echelon_Form);
disp(values)
disp(Row_Echelon_Form)
disp(Final_Answer(:,end))
x=input('Solve for another matrix?\n(1)Yes (2)No\n');
end
elseif options==8
disp('Matrix Inverses')
x=1;
while x~=2
scalar=input('Enter scalar value for matrix, 1 if none=');
rows=input('Enter the number of rows for 1st matrix=');
columns=input('Enter the number of columns for 1st matrix=');
size = ones(rows,columns);
disp(size)
values=input('put values inside''[]'' and add '';'' for the next
row=');
values=scalar.*values;
Inverse=inv(values);
disp(Inverse)
x=input('Solve for another matrix?\n(1)Yes (2)No\n');
end
elseif options==9
disp('LU Factorization of matrixes')
scalar=input('Enter scalar value for matrix, 1 if none=');
rows=input('Enter the number of rows for 1st matrix=');
columns=input('Enter the number of columns for 1st matrix=');
size = ones(rows,columns);
disp(size)
values=input('Enter values inside''[]'' and add '';'' for the next
row=');
values=scalar.*values;
[L,U]=lu(values);
option=input('Is there a value after the equal sign of the linear
equation or matrix?\n(1)Yes (2)No');
if option==1
B=input('Enter values after the equal sign of the equation or other
matrix inside''[]'' and add '';'' for the next row=');
L=[L B];
rl=rref(L);
U=[U rl(:,end)];
ru=rref(U);
L(:,end)=[];
U(:,end)=[];
disp(L([4 1 3 2],:))
disp(U)
disp(ru)
else
disp(L([4 1 3 2],:))
disp
end
elseif options==10
disp('Determinants of a matrix: (Must be a square Matrix)')
x=1;
while x~=2
scalar=input('Enter scalar value for matrix, 1 if none=');
rows=input('Enter the number of rows for 1st matrix=');
columns=input('Enter the number of columns for 1st matrix=');
 size = ones(rows,columns);
disp(size)
values=input('Enter values inside''[]'' and add '';'' for the
next row=');
values=scalar.*values;
Determinant=det(values);
disp(Determinant)
x=input('Solve for another matrix?\n(1)Yes (2)No\n');
end
else
disp('None of the choices')
end