_________________________________________________________________________

DEPARTMENT OF AVIONICS ENGINEERING

SUBJECT : CONTROL SYSTEMS LAB

SUBJECT CODE : 305312
LAB NO : 01

TITLE : INTRODUCTION TO MATLAB

SUBMITTED TO : Ms. Anila Ali Ahmed

SEMESTER : 5th
SECTION : B

Marks Obtained
Group Member Group Member Group Member
1 2 3
NAME DUA NASEEM
REGISTRATION # 220701014
LAB REPORT MARKS
PERFORMANCE MARKS
TOTAL MARKS

DEADLINE: 1100Hrs, 25th October 2024

DATE OF SUBMISSION: 24th October 2024

________________________________________________________________________
Experiment # 01 Page 1 of 9
Control Systems Lab
_________________________________________________________________________

Questions For Lab Report

Exercise 1:

Use Matlab command to obtain the following

a) Extract the fourth row of the matrix generated by magic(6)
b) Show the results of ‘x’ multiply by ‘y’ and ‘y’ divides by ‘x’. Given x =
[0:0.1:1.1] and y = [10:21]
c) Generate random matrix ‘r’ of size 4 by 5 with number varying between -8 and 9

CODE:

a) m=magic(6)
m4=m(4,:)
disp(m4)

OUTPUT:

CODE:

b) x = [0:0.1:1.1]
y = [10:21]
mul=x.*y
div=y./x
disp(mul)
disp(div)

OUTPUT:

________________________________________________________________________
Experiment # 01 Page 2 of 9
Control Systems Lab
_________________________________________________________________________

CODE:
c) r=randi([-8,9],4,5)
disp(r)
OUTPUT:

Exercise 2:

Use MATLAB commands to get exactly as the figure shown below

x=pi/2:pi/10:2*pi;
y=sin(x);
z=cos(x);

________________________________________________________________________
Experiment # 01 Page 3 of 9
Control Systems Lab
_________________________________________________________________________

CODE:

x=pi/2:pi/10:2*pi ; y=sin(x); z=cos(x)

subplot(2,1,1)
plot(x,y,'--+b' % Plot the x and z values with red '*'
%markers and a dashed linxlim([1.5,6.5]) % Set the x-axis

limits between 1.5 and 6.5

title("sine curve")
xlabel("angle")
ylabel("sinx")
ax = gca % Get the current axes object
grid on
ax.GridLineStyle = '--' % Set the grid line style to dashed
subplot(2,1,2)
plot(x,z,'--*r') % Plot the x and z values with red '*'
%markers and a dashed line

title("cos curve")
xlabel("angle")
ylabel("cosx")
ax = gca
grid on
ax.GridLineStyle = '--'

OUTPUT:

________________________________________________________________________
Experiment # 01 Page 4 of 9
Control Systems Lab
_________________________________________________________________________

sine curve
1

0.5
sinx
0

-0.5

-1
1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5
angle
cos curve
1

0.5
cosx

-0.5

-1
1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5
angle

Exercise 3:

Define Following Matrices in MATLAB

[ ]
1.2 3.04 7.3 3.2
−1 −12 1.2 5
A= 2.7 1 −2 4 ,
6.6 4.51 0 1
−3 −2.2 1 6

[]
1.3
−1
B=
2.5
1.6

a. Compute the size of A.

b. Compute AB.
c. Extract elements of A that are in 2nd to 4th row and 1st & 3rd coloumn.
d. Solve Ax=B for

CODE:

A=[1.2,3.04,7.3,3.2;
-1,-12,1.2,5;
2.7,1,-2,4;
________________________________________________________________________
Experiment # 01 Page 5 of 9
Control Systems Lab
_________________________________________________________________________

6.6,4.51,0,1;
-3,-2.2,1,6]

B=[1.3;-1;2.5;1.6]

size(A)
mult=A*B
disp(mult)
A(2:4,[1,3])

% to solve the linear equation AX=B matrix dimensions must

agree so ;
A = [1.2, 3.04, 7.3, 3.2;
-1, -12, 1.2, 5;
2.7, 1, -2, 4;
6.6, 4.51, 0, 1]
B = [1.3; -1; 2.5; 1.6]
x = A \ B
disp('Solution x for Ax = B:')
disp(x)

OUTPUT:

________________________________________________________________________
Experiment # 01 Page 6 of 9
Control Systems Lab
_________________________________________________________________________

Exercise 4:

a. Draw the graph of following function:

𝑧 = 𝑒−2𝑡 cos(3𝑡 + 1) − 𝜋 ≤ 𝑡 ≤ 5𝜋
________________________________________________________________________
Experiment # 01 Page 7 of 9
Control Systems Lab
_________________________________________________________________________

CODE:

t = linspace(-pi, 5*pi, 1000)

z = exp(-2*t) .* cos(3*t + 1)
plot(t, z, 'b-', 'LineWidth', 1)
title('Graph of z = e^{-2t} * cos(3t + 1)')
xlabel('t')
ylabel('z')
grid on

OUTPUT:
-2t
Graph of z = e * cos(3t + 1)
200

150

100

50

-50
z

-100

-150

-200

-250

-300
-5 0 5 10 15 20
t

b. The deflection of a cantilever beam is the distance its end moves in response to a
force applied at the end. The following table gives the deflection 𝑥 that was
produced in a particular beam by the given applied force 𝑓. Plot the graph:

Force 𝑓 (pound) 0 100 200 300 400 500 600 700 800
Deflection 𝑥 0 0.09 0.18 0.28 0.37 0.46 0.55 0.65 0.74
(inches)

CODE:

force = [0, 100, 200, 300, 400, 500, 600, 700,

800] % in pounds

________________________________________________________________________
Experiment # 01 Page 8 of 9
Control Systems Lab
_________________________________________________________________________

deflection = [0, 0.09, 0.18, 0.28, 0.37, 0.46,

0.55, 0.65, 0.74] % in inches
plot(force, deflection, 'm-o', 'LineWidth', 2,
'MarkerSize', 6)
title('Deflection of Cantilever Beam')
xlabel('Force (pounds)')
ylabel('Deflection (inches)')
grid on

OUTPUT:

Deflection of Cantilever Beam

0.8

0.7

0.6
Deflection (inches)

0.5

0.4

0.3

0.2

0.1

0
0 100 200 300 400 500 600 700 800
Force (pounds)

________________________________________________________________________
Experiment # 01 Page 9 of 9
Control Systems Lab