Gutierrez Experiment 1
Gutierrez Experiment 1
Experiment No. 1
A Laboratory Report
ECE 514 – Feedback and Control System
Submitted to:
Submitted by:
December 9, 2021
Control System
Experiment 1
Title:
Objectives:
• To be solve the linear constant differential equation by Laplace Transformation Method with the
use of MATLAB.
Discussion:
A system represented by a differential equation is difficult to model as a block diagram. Thus, we now lay
the groundwork for the Laplace transform, with which we can represent the input, output, and system as
separate entities. Further, their interrelationship will be simply algebraic. Let us first define the Laplace
transform and the show how it simplifies the representation of physical system (Nilsson, 1996).
L [f(t)] = F(s) = 𝑠𝑡 𝑑𝑡
Where s = σ + jω, a complex variable. Thus, knowing f(t) and that the integral in the equation above exist,
we can find a function,F(s), that is called the Laplace transform of f(t).
The notation for the lower limit means that even if f(t) is discontinuous at t = 0, we can start the integration
prior to the discontinuous as long as the integral converges. Thus, we can find the Laplace transform of
impulse functions. This property has distinct advantages when applying the Laplace transform to the
solution of differential equations where the initial conditions before the discontinuity. Using the Laplace
transform we need only know the initial condition before the discontinuity.
The inverse Laplace transform, which allows us to find f(t) given F(s), is
L -1 [F(s)] = f(t)u(t)
Where u(t) =
1 t>0 =
0 t < 0 Matlab
command reviews:
POLY Convert roots to polynomial.
POLY(A), when A is an N by N matrix, is a row vector with N+1 elements which are the coefficients of
the characteristic polynomial, DET(lambda*EYE(SIZE(A))-A).
POLY(V), when V is a vector, is a vector whose elements are the coefficients of the polynomial whose
roots are the elements of V. For vectors ROOTS and POLY are inverse function of each other, up to
ordering, scaling, and round off error.
ROOTS© computes the roots of the polynomial whose coefficients are the elements of the vectors C. If C
has N+ components, the polynomial is C(1)*X^N+…+C(N)*X+C(N+1).
PZMAP(SYS) computes the poles and (transmission) zero of the LTI models SYS and plots them
in the complex plane. The poles are plotted as x’s and the zeros are plotted as o’s.
PZMAP (SYS1, SYS2 ,….) shows the poles and zeros of multiple LTI models SYS1, SYS2,… on a single
plot. You can specify distinctive colors for each model, as in pzmap(sys1,sys2,’y’,sys3, ,’g’)
[P,Z] = PZMAP(SYS) returns the poles and zeros of the system in two column vectors P and Z. No plot is
drawn on the screen.
S = SYM(A) constructs an object S, of class ‘sym’, from A. If the input argument is a string, the result is a
symbolic number or variable. If the input argument is a numeric scalar or matrix, the result is a symbolic
representation of the given numeric values.
X = sym(‘x’) creates the symbolic variable with name ‘x’ and stores the result in x. x = sym(‘x’,’real’) also
assumes that x is real, so that conj(x) is equal to x. alpha = sym(‘alpha’) and r = sym (‘Rho’, ‘real’) are
other examples. Similarly, k = sym(‘k’, ‘positve’) makes k a positive (real) variable. X = sym(‘x’, ‘unreal’)
makes x a purely formal variable with no additional properties.
L = LAPLACE(F) is the Laplace transform of the scalar sym F with default independent variable t. The
default return is a function of s. If F = F(s), then LAPLACE returns a function of t: L = L(t). By definition
L(s) = int (F(t)* exp(-s*t),0,int), where integration occurs with respect to t.
L=LAPLACE(F,w,z) makes L a function of z instead of the default s (integration with respect to w).
Laplace transform:
MATLAB Command:
Roots:
Laplace transform:
MATLAB Command:
Laplace transform:
MATLAB Command:
𝑑2𝑐(𝑡) 𝑑𝑐(𝑡)
a.) 𝑑𝑡2 + 6 𝑑𝑡 + 13 c(t) =r(t)
Laplace Transform:
Inverse Laplace:
MATLAB Command:
𝑑2𝑐(𝑡) 𝑑𝑐(𝑡)
b.) 𝑑𝑡2 + 10 𝑑𝑡 + 25 c(t) =r(t)
Laplace Transform:
Inverse Laplace:
MATLAB Command:
𝑑3𝑐(𝑡) 𝑑2𝑐(𝑡) 𝑑𝑐(𝑡) 𝑑2𝑐(𝑡) 𝑑𝑐(𝑡)
c.) 𝑑𝑡3 + 3 𝑑𝑡2 +7 𝑑𝑡 +7 c(t) = 𝑑𝑡2 +4 𝑑𝑡 +3r(t)
Laplace Transform:
Inverse Laplace:
MATLAB Command: