Control Tutorials for MATLAB and Simulink (2024)

Once appropriate mathematical models of a system have been obtained, either in state-space or transfer function form, we may then analyze these models to predict how the system will respond in both the time and frequency domains. To put this in context, control systems are often designed to improve stability, speed of response, steady-state error, or prevent oscillations. In this section, we will show how to determine these dynamic properties from the system models.

Key MATLAB commands used in this tutorial are: tf , ssdata , pole , eig , step , pzmap , bode , linearSystemAnalyzer

Run Live Script Version in MATLAB Online

Related Tutorial Links

  • Time Resp Activity
  • Freq Resp Activity

Related External Links

Contents

  • Time Response Overview
  • Frequency Response Overview
  • Stability
  • System Order
  • First-Order Systems
  • Second-Order Systems

Time Response Overview

The time response represents how the state of a dynamic system changes in time when subjected to a particular input. Since the models we have derived consist of differential equations, some integration must be performed in order to determine the time response of the system. For some simple systems, a closed-form analytical solution may be available. However, for most systems, especially nonlinear systems or those subject to complicated inputs, this integration must be carried out numerically. Fortunately, MATLAB provides many useful resources for calculating time responses for many types of inputs, as we shall see in the following sections.

The time response of a linear dynamic system consists of the sum of the transient response which depends on the initial conditions and the steady-state response which depends on the system input. These correspond to the homogenous (free or zero input) and the particular solutions of the governing differential equations, respectively.

Frequency Response Overview

All the examples presented in this tutorial are modeled by linear constant coefficient differential equations and are thus linear time-invariant (LTI). LTI systems have the extremely important property that if the input to the system is sinusoidal, then the steady-state output will also be sinusoidal at the same frequency, but, in general, with different magnitude and phase. These magnitude and phase differences are a function of the frequency and comprise the frequency response of the system.

The frequency response of a system can be found from its transfer function in the following way: create a vector of frequencies (varying between zero or "DC" to infinity) and compute the value of the plant transfer function at those frequencies. If Control Tutorials for MATLAB and Simulink (1) is the open-loop transfer function of a system and Control Tutorials for MATLAB and Simulink (2) is the frequency vector, we then plot Control Tutorials for MATLAB and Simulink (3) versus Control Tutorials for MATLAB and Simulink (4). Since Control Tutorials for MATLAB and Simulink (5) is a complex number, we can plot both its magnitude and phase (the Bode Plot) or its position in the complex plane (the Nyquist Diagram). Both methods display the same information, but in different ways.

Stability

For our purposes, we will use the Bounded Input Bounded Output (BIBO) definition of stability which states that a system is stable if the output remains bounded for all bounded (finite) inputs. Practically, this means that the system will not "blow up" while in operation.

The transfer function representation is especially useful when analyzing system stability. If all poles of the transfer function (values of Control Tutorials for MATLAB and Simulink (6) for which the denominator equals zero) have negative real parts, then the system is stable. If any pole has a positive real part, then the system is unstable. If we view the poles on the complex s-plane, then all poles must be in the left-half plane (LHP) to ensure stability. If any pair of poles is on the imaginary axis, then the system is marginally stable and the system will tend to oscillate. A system with purely imaginary poles is not considered BIBO stable. For such a system, there will exist finite inputs that lead to an unbounded response. The poles of an LTI system model can easily be found in MATLAB using the pole command, an example of which is shown below:

s = tf('s');G = 1/(s^2+2*s+5)pole(G)
G = 1 ------------- s^2 + 2 s + 5 Continuous-time transfer function.ans = -1.0000 + 2.0000i -1.0000 - 2.0000i

Thus this system is stable since the real parts of the poles are both negative. The stability of a system may also be found from the state-space representation. In fact, the poles of the transfer function are the eigenvalues of the system matrix Control Tutorials for MATLAB and Simulink (7). We can use the eig command to calculate the eigenvalues using either the LTI system model directly, eig(G), or the system matrix as shown below.

[A,B,C,D] = ssdata(G);eig(A)
ans = -1.0000 + 2.0000i -1.0000 - 2.0000i

System Order

The order of a dynamic system is the order of the highest derivative of its governing differential equation. Equivalently, it is the highest power of Control Tutorials for MATLAB and Simulink (8) in the denominator of its transfer function. The important properties of first-, second-, and higher-order systems will be reviewed in this section.

First-Order Systems

First-order systems are the simplest dynamic systems to analyze. Some common examples include mass-damper systems and RC circuits.

The general form of the first-order differential equation is as follows

(1)Control Tutorials for MATLAB and Simulink (9)

The form of a first-order transfer function is

(2)Control Tutorials for MATLAB and Simulink (10)

where the parameters Control Tutorials for MATLAB and Simulink (11) and Control Tutorials for MATLAB and Simulink (12) completely define the character of the first-order system.

DC Gain

The DC gain, Control Tutorials for MATLAB and Simulink (13), is the ratio of the magnitude of the steady-state step response to the magnitude of the step input. For stable transfer functions, the Final Value Theorem demonstrates that the DC gain is the value of the transfer function evaluated at Control Tutorials for MATLAB and Simulink (14) = 0. For first-order systems of the forms shown, the DC gain is Control Tutorials for MATLAB and Simulink (15).

Time Constant

The time constant of a first-order system is Control Tutorials for MATLAB and Simulink (16) which is equal to the time it takes for the system's response to reach 63% of its steady-state value for a step input (from zero initial conditions) or to decrease to 37% of the initial value for a system's free response. More generally, it represents the time scale for which the dynamics of the system are significant.

Poles/Zeros

First-order systems have a single real pole, in this case at Control Tutorials for MATLAB and Simulink (17). Therefore, the system is stable if Control Tutorials for MATLAB and Simulink (18) is positive and unstable if Control Tutorials for MATLAB and Simulink (19) is negative. Standard first-order system have no zeros.

Step Response

We can calculate the system time response to a step input of magnitude Control Tutorials for MATLAB and Simulink (20) using the following MATLAB commands:

k_dc = 5;Tc = 10;u = 2;s = tf('s');G = k_dc/(Tc*s+1)step(u*G)
G = 5 -------- 10 s + 1 Continuous-time transfer function.

Control Tutorials for MATLAB and Simulink (21)

Note: MATLAB also provides a powerful graphical user interface for analyzing LTI systems which can be accessed using the syntax linearSystemAnalyzer('step',G).

If you right-click on the step response graph and select Characteristics, you can choose to have several system metrics overlaid on the response: peak response, settling time, rise time, and steady-state.

Settling Time

The settling time, Control Tutorials for MATLAB and Simulink (22), is the time required for the system output to fall within a certain percentage (i.e. 2%) of the steady-state value for a step input. The settling times for a first-order system for the most common tolerances are provided in the table below. Note that the tighter the tolerance, the longer the system response takes to settle to within this band, as expected.

10%5%2%1%
Ts=2.3/a=2.3TcTs=3/a=3TcTs=3.9/a=3.9TcTs=4.6/a=4.6Tc

Rise Time

The rise time, Control Tutorials for MATLAB and Simulink (23), is the time required for the system output to rise from some lower level x% to some higher level y% of the final steady-state value. For first-order systems, the typical range is 10% - 90%.

Bode Plots

Bode diagrams show the magnitude and phase of a system's frequency response, Control Tutorials for MATLAB and Simulink (24), plotted with respect to frequency Control Tutorials for MATLAB and Simulink (25). We can generate the Bode plot of a system Control Tutorials for MATLAB and Simulink (26) in MATLAB using the syntax bode(G) as shown below.

bode(G)

Control Tutorials for MATLAB and Simulink (27)

Again the same results could be obtained using the Linear System Analyzer GUI, linearSystemAnalyzer('bode',G).

Bode plots employ a logarithmic frequency scale so that a larger range of frequencies are visible. Also, the magnitude is represented using the logarithmic decibel unit (dB) defined as:

(3)Control Tutorials for MATLAB and Simulink (28)

As with the frequency axis, the decibel scale allows us to view a much larger range of magnitudes on a single plot. Also, as we shall see in subsequent tutorials, when components and controllers are placed in series, the transfer function of the overall system is the product of the individual transfer functions. Using the dB scale, the magnitude plot of the overall system is simply the sum of the magnitude plots of the individual transfer functions. The phase plot of the overall system is also just the sum of the individual phase plots.

The low frequency magnitude of the first-order Bode plot is Control Tutorials for MATLAB and Simulink (29). The magnitude plot has a bend at the frequency equal to the absolute value of the pole (ie. Control Tutorials for MATLAB and Simulink (30)), and then decreases 20 dB for every factor of ten increase in frequency (slope = -20 dB/decade). The phase plot is asymptotic to 0 degrees at low frequencies, and asymptotic to -90 degrees at high frequencies. Between frequency 0.1a and 10a, the phase changes by approximately -45 degrees for every factor of ten increase in frequency (-45 degrees/decade).

We will see in the Frequency Methods for Controller Design Section how to use Bode plots to calculate closed-loop stability and performance of feedback systems.

Second-Order Systems

Second-order systems are commonly encountered in practice, and are the simplest type of dynamic system to exhibit oscillations. Examples include mass-spring-damper systems and RLC circuits. In fact, many true higher-order systems may be approximated as second-order in order to facilitate analysis.

The canonical form of the second-order differential equation is as follows

(4)Control Tutorials for MATLAB and Simulink (31)

The canonical second-order transfer function has the following form, in which it has two poles and no zeros.

(5)Control Tutorials for MATLAB and Simulink (32)

The parameters Control Tutorials for MATLAB and Simulink (33), Control Tutorials for MATLAB and Simulink (34), and Control Tutorials for MATLAB and Simulink (35) characterize the behavior of a canonical second-order system.

DC Gain

The DC gain, Control Tutorials for MATLAB and Simulink (36), again is the ratio of the magnitude of the steady-state step response to the magnitude of the step input, and for stable systems it is the value of the transfer function when Control Tutorials for MATLAB and Simulink (37). For the forms given,

(6)Control Tutorials for MATLAB and Simulink (38)

Damping Ratio

The damping ratio Control Tutorials for MATLAB and Simulink (39) is a dimensionless quantity charaterizing the rate at which an oscillation in the system's response decays due to effects such as viscous friction or electrical resistance. From the above definitions,

(7)Control Tutorials for MATLAB and Simulink (40)

Natural Frequency

The natural frequency Control Tutorials for MATLAB and Simulink (41) is the frequency (in rad/s) that the system will oscillate at when there is no damping, Control Tutorials for MATLAB and Simulink (42).

(8)Control Tutorials for MATLAB and Simulink (43)

Poles/Zeros

The canonical second-order transfer function has two poles at:

(9)Control Tutorials for MATLAB and Simulink (44)

Underdamped Systems

If Control Tutorials for MATLAB and Simulink (45), then the system is underdamped. In this case, both poles are complex-valued with negative real parts; therefore, the system is stable but oscillates while approaching the steady-state value. Specifically, the natural response oscillates with the damped natural frequency, Control Tutorials for MATLAB and Simulink (46) (in rad/sec).

k_dc = 1;w_n = 10;zeta = 0.2;s = tf('s');G1 = k_dc*w_n^2/(s^2 + 2*zeta*w_n*s + w_n^2);pzmap(G1)axis([-3 1 -15 15])

Control Tutorials for MATLAB and Simulink (47)

step(G1)axis([0 3 0 2])

Control Tutorials for MATLAB and Simulink (48)

Settling Time

The settling time, Control Tutorials for MATLAB and Simulink (49), is the time required for the system ouput to fall within a certain percentage of the steady-state value for a step input. For a canonical second-order, underdamped system, the settling time can be approximated by the following equation:

(10)Control Tutorials for MATLAB and Simulink (50)

The settling times for the most common tolerances are presented in the following table:

10%5%2%1%
Ts=2.3/(zeta*w_n)Ts=3/(zeta*w_n)Ts=3.9/(zeta*w_n)Ts=4.6/(zeta*w_n)

Percent Overshoot

The percent overshoot is the percent by which a system's step response exceeds its final steady-state value. For a second-order underdamped system, the percent overshoot Control Tutorials for MATLAB and Simulink (51) is directly related to the damping ratio by the following equation. Here, Control Tutorials for MATLAB and Simulink (52) is a decimal number where 1 corresponds to 100% overshoot.

(11)Control Tutorials for MATLAB and Simulink (53)

For second-order underdamped systems, the 1% settling time, Control Tutorials for MATLAB and Simulink (54), 10-90% rise time, Control Tutorials for MATLAB and Simulink (55), and percent overshoot, Control Tutorials for MATLAB and Simulink (56), are related to the damping ratio and natural frequency as shown below.

(12)Control Tutorials for MATLAB and Simulink (57)

(13)Control Tutorials for MATLAB and Simulink (58)

(14)Control Tutorials for MATLAB and Simulink (59)

Overdamped Systems

If Control Tutorials for MATLAB and Simulink (60), then the system is overdamped. Both poles are real and negative; therefore, the system is stable and does not oscillate. The step response and a pole-zero map of an overdamped system are calculated below:

zeta = 1.2;G2 = k_dc*w_n^2/(s^2 + 2*zeta*w_n*s + w_n^2);pzmap(G2)axis([-20 1 -1 1])

Control Tutorials for MATLAB and Simulink (61)

step(G2)axis([0 1.5 0 1.5])

Control Tutorials for MATLAB and Simulink (62)

Critically-Damped Systems

If Control Tutorials for MATLAB and Simulink (63), then the system is critically damped. Both poles are real and have the same magnitude, Control Tutorials for MATLAB and Simulink (64). For a canonical second-order system, the quickest settling time is achieved when the system is critically damped. Now change the value of the damping ratio to 1, and re-plot the step response and pole-zero map.

zeta = 1;G3 = k_dc*w_n^2/(s^2 + 2*zeta*w_n*s + w_n^2);pzmap(G3)axis([-11 1 -1 1])

Control Tutorials for MATLAB and Simulink (65)

step(G3)axis([0 1.5 0 1.5])

Control Tutorials for MATLAB and Simulink (66)

Undamped Systems

If Control Tutorials for MATLAB and Simulink (67), then the system is undamped. In this case, the poles are purely imaginary; therefore, the system is marginally stable and the step response oscillates indefinitely.

zeta = 0;G4 = k_dc*w_n^2/(s^2 + 2*zeta*w_n*s + w_n^2);pzmap(G4)axis([-1 1 -15 15])

Control Tutorials for MATLAB and Simulink (68)

step(G4)axis([0 5 -0.5 2.5])

Control Tutorials for MATLAB and Simulink (69)

Bode Plot

We show the Bode magnitude and phase plots for all damping conditions of a second-order system below:

bode(G1,G2,G3,G4)legend('underdamped: zeta < 1','overdamped: zeta > 1','critically-damped: zeta = 1','undamped: zeta = 0')

Control Tutorials for MATLAB and Simulink (70)

The magnitude of the bode plot of a second-order system drops off at -40 dB per decade in the limit, while the relative phase changes from 0 to -180 degrees. For underdamped systems, we also see a resonant peak near the natural frequency, Control Tutorials for MATLAB and Simulink (71) = 10 rad/s. The size and sharpness of the peak depends on the damping in the system, and is charaterized by the quality factor, or Q-Factor, defined below. The Q-factor is an important property in signal processing.

(15)Control Tutorials for MATLAB and Simulink (72)


Published with MATLAB® 9.2

Control Tutorials for MATLAB and Simulink (2024)

FAQs

How do I find answers in MATLAB? ›

To view all of your solutions, go to a Problem page and click View my solutions. You can view your solutions in a list or in the Solution Map. If using the list view, you can review the display by selecting a Sort by option.

Is MATLAB Simulink hard to learn? ›

Is MATLAB Hard to Learn? MATLAB is designed for the way you think and the work you do, so learning is accessible whether you are a novice or an expert. The Help Center is always available to guide you with robust documentation, community answers, and how-to videos.

How much time is required to learn MATLAB? ›

If you're a novice programmer, you can expect it to take a little longer than if you were a more seasoned programmer. Someone who can afford to devote all their time to MATLAB can finish learning the language in two weeks. If you have a lot of other responsibilities, however, it will take you longer to complete.

How do you get a long answer in MATLAB? ›

To format the way numbers display, do one of the following:
  1. On the Home tab, in the Environment section, click Preferences. Select MATLAB > Command Window, and then choose a Numeric format option.
  2. Use the format function, for example: format short format short e format long.

How do I find Solver in MATLAB? ›

Select solvers in the Solver pane of model configuration parameters. All solvers provided by MATLAB® and Simulink follow a similar naming convention: ode , followed by two or three numerals indicating the orders of the solver.

Is MATLAB harder than Python? ›

The Difference in Technical Computing:

They are both used for the same type of work, such as numerical analysis, data visualization, and scientific computation. When it comes to syntax and readability, Python is often easier to read and understand than MATLAB.

What is the salary of MATLAB Simulink? ›

Matlab Developer salary in India ranges between ₹ 2.0 Lakhs to ₹ 9.4 Lakhs with an average annual salary of ₹ 5.5 Lakhs. Salary estimates are based on 343 latest salaries received from Matlab Developers. 1 - 5 years exp. 1 - 5 years exp.

Is MATLAB enough for a job? ›

Conclusion. The industry has some familiar buzz that learning MATLAB will not be a good opportunity for a better career. But this is not fully true. Yes, it is an acceptable reason that salary or company structure will never be able to touch available popular jobs on other programming technologies.

What coding language is MATLAB most like? ›

Language Comparison

The language of Python and MATLAB can be used interactively (a single command at a time) or to develop large-scale applications.

Can I learn MATLAB on my own? ›

Get Started with Introductory Videos

See common applications and workflows, and discover new product capabilities. Get started with MATLAB by walking through an example. This video shows you the basics, and it gives you an idea of what working in MATLAB is like.

What should I learn first MATLAB or Python? ›

The OOP in MATLAB is more advanced and complex, which to some can be more confusing. That being said, MATLAB is generally a more advanced language while Python is more of a beginner's language. Therefore, just because MATLAB may be more complex and confusing at first, with practice, it will become easier to grasp.

Is MATLAB in high demand? ›

Matlab careers are actually on the rise today. It's a very popular programming language. It can be used by a developer, engineer, programmer, scientist, etc. to collect and sort out data, and develop apps, software, and sites.

Can I learn MATLAB in 1 month? ›

If I generalize my answer highly, It may take at least 3 months to learn matlab and may take maximum 3 years to become an expert.

What should I learn before learning MATLAB? ›

Also, it would be beneficial to have some basic knowledge on linear algebra and calculus. However, it is not essential to have any prior programming experience or a strong mathematical background as this course is designed with beginners in mind.

How to check results in MATLAB? ›

View Results in Command Window

The Summary Report link provides access to the Model Advisor Command-Line Summary report. You can review additional results in the Command Window by calling the DisplayResults parameter when you run the Model Advisor.

How do you find the step response in MATLAB? ›

[ y , tOut ] = step( sys , tFinal ) computes the step response from t = 0 to the end time t = tFinal . [ y , tOut ] = step( sys , t ) returns the step response of a dynamic system model sys at the times specified in the vector t .

How do I find something in MATLAB code? ›

Search Using Find Dialog Box

The Find dialog box opens. The search begins at the current cursor position. MATLAB finds the text you specified and highlights it. MATLAB beeps when a search for Find Next reaches the end of the Command Window, or when a search for Find Previous reaches the top of the Command Window.

How do you search text in MATLAB? ›

Search for Text

To locate where the text occurs, use the strfind function, which returns starting indices. Find and extract text using extraction functions, such as extract , extractBetween , extractBefore , or extractAfter . Optionally, include the boundary text.

Top Articles
Latest Posts
Recommended Articles
Article information

Author: Margart Wisoky

Last Updated:

Views: 6391

Rating: 4.8 / 5 (78 voted)

Reviews: 85% of readers found this page helpful

Author information

Name: Margart Wisoky

Birthday: 1993-05-13

Address: 2113 Abernathy Knoll, New Tamerafurt, CT 66893-2169

Phone: +25815234346805

Job: Central Developer

Hobby: Machining, Pottery, Rafting, Cosplaying, Jogging, Taekwondo, Scouting

Introduction: My name is Margart Wisoky, I am a gorgeous, shiny, successful, beautiful, adventurous, excited, pleasant person who loves writing and wants to share my knowledge and understanding with you.