Control Tutorials for MATLAB and Simulink (2024)

Key MATLAB commands used in this tutorial are: tf , controlSystemDesigner

Related Tutorial Links

  • Freq Resp Intro
  • Freq Resp Activity
  • Lead/Lag Freq Resp

Related External Links

Contents

  • Drawing the uncompensated system's Bode plot
  • Adding an integrator
  • Gain and phase margin specifications and controller design

From the main problem, the open-loop transfer function of the DC Motor is given as follows.

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

The structure of the control system has the form shown in the figure below.

Control Tutorials for MATLAB and Simulink (2)

For the original problem setup and the derivation of the above equations, please refer to the DC Motor Position: System Modeling page.

With a 1-radian step reference, the design criteria are the following.

  • Settling time less than 0.040 seconds
  • Overshoot less than 16%
  • No steady-state error, even in the presence of a step disturbance input

First create a new m-file and type in the following commands (refer to main problem for the details of getting these commands).

J = 3.2284E-6;b = 3.5077E-6;K = 0.0274;R = 4;L = 2.75E-6;s = tf('s');P_motor = K/(s*((J*s+b)*(L*s+R)+K^2));

Drawing the uncompensated system's Bode plot

The main idea of frequency-based design is to use the Bode plot of the open-loop transfer function to estimate the system's closed-loop response. Adding a controller to the system changes the open-loop Bode plot, thereby changing the closed-loop response. Let's first draw the Bode plot for the original open-loop transfer function. Throughout this section, we will use the interactive Control System Designer to design our controller. This tool can be opened for the uncompensated plant by adding the command controlSystemDesigner('bode',P_motor) to your m-file or by going to the APPS tab of the MATLAB toolstrip and clicking on the app icon under Control System Design and Analysis.

One window titled Control System Designer will open initially having the form shown in the figure below. In the window, you will be able to see both the open-loop Bode plot and the closed-loop step response of the transfer function passed via the controlSystemDesigner function. If the string 'bode' is omitted from the function call, the default initial window includes the root locus plot, in addition to the open-loop Bode plot and closed-loop step response plot.

Control Tutorials for MATLAB and Simulink (3)

From previous pages we know that the open-loop plant in unity feedback does not satisfy the given system requirements. We will, therefore, add different compensators to our system in order to reshape the open-loop Bode plot to give the desired closed-loop performance.

Adding an integrator

Now let's add a pure integrator compensator to our system. The architecture of the feedback system can be modified within the Control System Designer window. For the default format, the compensator C is placed in series with the plant G in the forward path of the block diagram. This corresponds to the architecture we are assuming, as shown in the schematic at the top of this page. Placing an integral compensator in this position will achieve the goal of zero steady-state error in the presence of a step disturbance input.

The form of our controller is then modified under the Compensator Editor window that we can open by right-clicking on the Bode plot and then selecting Edit Compensator. Specifically, an integrator is added to our controller by right-clicking in the Dynamics section of the window and selecting Add Pole/Zero > Integrator from the resulting menu. The figure below illustrates what the Compensator Editor window should look like.

Control Tutorials for MATLAB and Simulink (4)

The Bode plot that is already open should then change as shown below to reflect the addition of the integrator.

Control Tutorials for MATLAB and Simulink (5)

Gain and phase margin specifications and controller design

Recall that we additionally desire an overshoot of less than 16%. Overshoot of a step response correlates to a system's phase margin. We will apply the relationship that phase margin approximately equals 100 times the damping ratio (Control Tutorials for MATLAB and Simulink (6)) to estimate the desired phase margin of our system. We will additionally use our 0.040 second settling time requirement to estimate the desired bandwidth frequency (Control Tutorials for MATLAB and Simulink (7)) of our closed-loop system. Add the following code to your m-file to estimate the desired phase margin and bandwidth. Note that these relationships are intended for a second-order canonical system, which we do not have, but the results will serve as a good starting point for our design.

Mp = 0.16;ts = 0.04;zeta = -log(Mp) / sqrt( pi^2 + (log(Mp))^2 );PM = 100*zetaWbw = (4/(ts*zeta))*sqrt((1 - 2*zeta^2) + sqrt(4*zeta^4 - 4*zeta^2+2))
PM = 50.3868Wbw = 251.5743

The above results demonstrate that we desire to have at least 50 degrees of phase margin in an attempt to meet the overshoot requirement. Also, a closed-loop bandwidth of 250 rad/sec is desired to additionally achieve the settle time requirement. We will estimate the closed-loop bandwidth as the point on the open-loop Bode plot where the magnitude is between -6 and -7.5 dB. From the Bode plot we see that we must add about 110 degrees of phase and 80 dB of gain at a frequency of 250 rad/sec in order to move the gain crossover frequency to 250 rad/sec and provide 50 degrees of phase margin. The magnitude plot will then also lie between -6 and -7.5 dB at a frequency greater than 250 rad/sec.

From the Bode phase plot we can see that there is a pole near -60 as indicated by the blue x on the Bode plot near the frequency 60 rad/sec (as well as the change in the magnitude and phase plots). We will, therefore, begin to modify our compensator by adding a zero at s = -60 in order to flatten out the phase curve. A real zero can be added again by right-clicking in the Dynamics section of the Compensator Editor window and selecting Real Zero. The location of the zero can then be set by typing -60 into the Location cell of the Edit Selected Dynamics portion of the window and pressing Enter. Your compensator now has the form of a PI controller and the compensated open-loop Bode plot should have the modified form shown below where the red o shows the location of the new zero. To achieve the limits on the axes shown below, right-click on the Bode plot and select Properties from the resulting menu.

Control Tutorials for MATLAB and Simulink (8)

From the Bode plot we can see that we need at least 50 additional degrees of phase if we move our crossover frequency to 250 rad/sec. Let's then add a lead compensator to add exactly 50 degrees of phase at this frequency. The following equations provide the pole and zero locations of the lead compensator.

a = (1 - sin(PM*pi/180))/(1 + sin(PM*pi/180));T = 1/(Wbw*sqrt(a));zero = -1/Tpole = -1/(a*T)
zero = -90.6050pole = -698.5222

The pole and zero of the lead compensator can then be added to our system in the same manner that the integrator and real zero were added above. Alternatively, we can add a lead compensator by specifiying the desired maximum phase and frequency location within the Control System Designer window. This is accomplished also by right-clicking in the Dynamics section of the Compensator Editor window. Selecting the Lead compensator in the Dynamics section of the window then opens four cells in the Edit Selected Dynamics portion of the window. These four cells allow you to specify the lead portion of the compensator either in terms of the pole and zero locations, or in terms of the phase and frequency location. If we enter 50 into the Max Delta Phase (deg) cell and 250 into the at Frequency cell and hit Enter, we will get a pole and zero location that closely (though not exactly) matches those we calculated above. The modified Bode plot should then have a similar appearance to the one shown below.

Control Tutorials for MATLAB and Simulink (9)

This new Bode plot now shows that the phase margin is about right at 250 rad/sec, but the gain is too small by about 55 dB if we wish for the gain crossover to occur at 250 rad/sec. Note that each pole and zero we have added to the compensator has the form

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

rather than (s+z) or (s+p). This default form has a DC gain of 1. Therefore, your magnitude plot will be shifted up or down compared to the figure shown above if the terms of your compensator have the form (s+z) and (s+p). Converting from decibels determines that a gain of 600 will shift the magnitude plot up approximately 55 dBs. This change in gain can be accomplished by typing 600 into the cell multiplying the compensator at the top of the Compensator Editor window and pressing Enter. The resulting Bode plot is shown below. As indicated on the plot, the phase margin is 49.8 degrees and the gain crossover frequency is 236 rad/sec which are close to our goals of 50 degrees and 250 rad/sec.

Control Tutorials for MATLAB and Simulink (11)

Let's now check the resulting step response of the closed-loop system by moving to the IOTransfer_r2y:step tab. If you have accidentally closed this tab, you can re-open it from the Control System Designer window by clicking on the New Plot menu and selecting New Step. In response, a new window titled New Step to plot will appear. From the Select Responses to Plot menu, then choose IOTransfer_r2y and click the button Plot. The response of the output y of the closed-loop system for a step reference r will then appear in the Control System Designer window. See the figure below for an illustration of this process.

Control Tutorials for MATLAB and Simulink (12)

The resulting step response plot is shown below, where the details of the peak response and the settling time are added from the right-click menu under Characteristics.

Control Tutorials for MATLAB and Simulink (13)

From the above figure, it can be seen that the overshoot is too large and the settling time is better than expected. This is attributable to the fact that our system does not actually have the form of a canonical second-order system and because we did not exactly achieve a phase margin of 50 degrees and a bandwidth of 250 rad/sec.

In order to reduce the overshoot, we will adjust the lead portion of our compensator to provide more phase margin. Specifically, within the Compensator Editor window you can change the value of the Max Delta Phase (deg). By increasing this value, you should see the bump in the phase plot increase and the overshoot in the step response plot decrease. This change will of course also affect the gain crossover frequency. Therefore, we need to change the value of the compensator gain to shift the magnitude plot such that the gain crossover frequency is still near 250 rad/sec. This can be accomplished by manually typing in different values for the gain in the Compensator Editor window as we have done previously. Another option is to use the Bode plot since plots in the Control System Designer allow for graphical tuning. Specifically, if you place your cursor over the magnitude plot it will turn into a hand and you can "grab" the plot and move it up or down. The gain in the compensator will automatically change to reflect the movement of the plot. In the end, a choice of 70 degrees of additional phase at 250 rad/sec and a gain of 315 give us a phase margin of 70 degrees and an overshoot less than 16% as shown below.

Control Tutorials for MATLAB and Simulink (14)

Control Tutorials for MATLAB and Simulink (15)

From the above step response we see that even though the overshoot is fine, the settling time is too long. We can speed up the system response by trying a slightly higher bandwidth. This can again be accomplished from the Compensator Editor window by increasing the frequency where the peak phase of the lead compensator is positioned. Through trial and error we choose an at Frequency of 350 rad/sec. We then use the graphical tuning technique to "grab" the magnitude plot of the Bode diagram to shift the crossover frequency to a location that provides approximately 70 degrees of phase margin. A gain of 670 seems to work well.

Throughout the tuning, the Bode diagram and step response plot changed according to the changes being made to the compensator. Your final plots should be similar to the two shown below.

Control Tutorials for MATLAB and Simulink (16)

Control Tutorials for MATLAB and Simulink (17)

Now we have the performance we desire. We have less than 16% overshoot and a settling time less than 40 milliseconds. Furthermore, the integral action of the compensator provides zero steady-state error to a step reference, even in the presence of a step disturbance input.


Published with MATLAB® 9.2

Control Tutorials for MATLAB and Simulink (2024)

References

Top Articles
Latest Posts
Article information

Author: Terrell Hackett

Last Updated:

Views: 5736

Rating: 4.1 / 5 (72 voted)

Reviews: 87% of readers found this page helpful

Author information

Name: Terrell Hackett

Birthday: 1992-03-17

Address: Suite 453 459 Gibson Squares, East Adriane, AK 71925-5692

Phone: +21811810803470

Job: Chief Representative

Hobby: Board games, Rock climbing, Ghost hunting, Origami, Kabaddi, Mushroom hunting, Gaming

Introduction: My name is Terrell Hackett, I am a gleaming, brainy, courageous, helpful, healthy, cooperative, graceful person who loves writing and wants to share my knowledge and understanding with you.