DC Motor Speed: PID Controller Design

Contents

  • Proportional control
  • PID control
  • Tuning the gains

From the main problem, the dynamic equations in the Laplace domain and the open-loop transfer function of the DC Motor are the following.

(1)$$ s(Js + b)\Theta(s) = KI(s) $$

(2)$$ (Ls + R)I(s) = V(s) - Ks\Theta(s) $$

(3)$$ P(s) = \frac{\dot{\Theta}(s)}{V(s)} = \frac{K}{(Js + b)(Ls + R) + K^2}  \qquad [\frac{rad/sec}{V}] $$

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

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

For a 1-rad/sec step reference, the design criteria are the following.

  • Settling time less than 2 seconds
  • Overshoot less than 5%
  • Steady-state error less than 1%

Now let's design a controller using the methods introduced in the Introduction: PID Controller Design page. Create a new m-file and type in the following commands.

J = 0.01;
b = 0.1;
K = 0.01;
R = 1;
L = 0.5;
s = tf('s');
P_motor = K/((J*s+b)*(L*s+R)+K^2);

Recall that the transfer function for a PID controller is:

(4)$$ C(s) = K_{p} + \frac {K_{i}} {s} + K_{d}s = \frac{K_{d}s^2 + K_{p}s + K_{i}} {s} $$

Proportional control

Let's first try employing a proportional controller with a gain of 100, that is, C(s) = 100. To determine the closed-loop transfer function, we use the feedback command. Add the following code to the end of your m-file.

Kp = 100;
C = pid(Kp);
sys_cl = feedback(C*P_motor,1);

Now let's examine the closed-loop step response. Add the following commands to the end of your m-file and run it in the command window. You should generate the plot shown below. You can view some of the system's characteristics by right-clicking on the figure and choosing Characteristics from the resulting menu. In the figure below, annotations have specifically been added for Settling TimePeak Response, and Steady State.

 

         t = 0:0.01:5;
         step(sys_cl,t)
         grid
         title('Step Response with Proportional Control')
         

 

From the plot above we see that both the steady-state error and the overshoot are too large. Recall from the Introduction: PID Controller Design page that increasing the proportional gain $K_p$ will reduce the steady-state error. However, also recall that increasing $K_p$ often results in increased overshoot, therefore, it appears that not all of the design requirements can be met with a simple proportional controller.

This fact can be verified by experimenting with different values of $K_p$. Specifically, you can employ the Control System Designer by entering the command controlSystemDesigner(P_motor) or by going to the APPS tab and clicking on the app icon under Control System Design and Analysis and then opening a closed-loop step response plot from the New Plot tab of the Control System Designer window as shown below.

After that you can right-click on the plot and select Edit Compensator. You can then vary the control gain in the Compensator Editor window and see the resulting effect on the closed-loop step response as shown below.

A little experimentation verifies what we anticipated, a proportional controller is insufficient for meeting the given design requirements; derivative and/or integral terms must be added to the controller.

PID control

Recall from the Introduction: PID Controller Design page, adding an integral term will eliminate the steady-state error to a step reference and a derivative term will often reduce the overshoot. Let's try a PID controller with small $K_i$ and $K_d$. Modify your m-file so that the lines defining your control are as follows. Running this new m-file gives you the plot shown below.

Kp = 75;
Ki = 1;
Kd = 1;
C = pid(Kp,Ki,Kd);
sys_cl = feedback(C*P_motor,1);
step(sys_cl,[0:1:200])
title('PID Control with Small Ki and Small Kd')

Inspection of the above indicates that the steady-state error does indeed go to zero for a step input. However, the time it takes to reach steady-state is far larger than the required settling time of 2 seconds.

Tuning the gains

In this case, the long tail on the step response graph is due to the fact that the integral gain is small and, therefore, it takes a long time for the integral action to build up and eliminate the steady-state error. This process can be sped up by increasing the value of $K_i$. Go back to your m-file and change $K_i$ to 200 as in the following. Rerun the file and you should get the plot shown below. Again the annotations are added by right-clicking on the figure and choosing Characteristics from the resulting menu.

 

         Kp = 100;
         Ki = 200;
         Kd = 1;
         C = pid(Kp,Ki,Kd);
         sys_cl = feedback(C*P_motor,1);
         step(sys_cl, 0:0.01:4)
         grid
         title('PID Control with Large Ki and Small Kd')
         

 

As expected, the steady-state error is now eliminated much more quickly than before. However, the large $K_i$ has greatly increased the overshoot. Let's increase $K_d$ in an attempt to reduce the overshoot. Go back to the m-file and change $K_d$ to 10 as shown in the following. Rerun your m-file and the plot shown below should be generated.

 

         Kp = 100;
         Ki = 200;
         Kd = 10;
         C = pid(Kp,Ki,Kd);
         sys_cl = feedback(C*P_motor,1);
         step(sys_cl, 0:0.01:4)
         grid
         title('PID Control with Large Ki and Large Kd')
         

 

As we had hoped, the increased $K_d$ reduced the resulting overshoot. Now we know that if we use a PID controller with $K_p$ = 100, $K_i$ = 200, and $K_d$ = 10, all of our design requirements will be satisfied.

Matlabsolutions.com provides guaranteed satisfaction with a commitment to complete the work within time. Combined with our meticulous work ethics and extensive domain experience, We are the ideal partner for all your homework/assignment needs. We pledge to provide 24*7 support to dissolve all your academic doubts. We are composed of 300+ esteemed Matlab and other experts who have been empanelled after extensive research and quality check.

Matlabsolutions.com provides undivided attention to each Matlab assignment order with a methodical approach to solution. Our network span is not restricted to US, UK and Australia rather extends to countries like Singapore, Canada and UAE. Our Matlab assignment help services include Image Processing Assignments, Electrical Engineering Assignments, Matlab homework help, Matlab Research Paper help, Matlab Simulink help. Get your work done at the best price in industry.

Machine Learning in MATLAB

Train Classification Models in Classification Learner App

Train Regression Models in Regression Learner App

Distribution Plots

Explore the Random Number Generation UI

Design of Experiments

Machine Learning Models

Logistic regression

Logistic regression create generalized linear regression model - MATLAB fitglm 2

Support Vector Machines for Binary Classification

Support Vector Machines for Binary Classification 2

Support Vector Machines for Binary Classification 3

Support Vector Machines for Binary Classification 4

Support Vector Machines for Binary Classification 5

Assess Neural Network Classifier Performance

Naive Bayes Classification

ClassificationTree class

Discriminant Analysis Classification

Ensemble classifier

ClassificationTree class 2

Train Generalized Additive Model for Binary Classification

Train Generalized Additive Model for Binary Classification 2

Classification Using Nearest Neighbors

Classification Using Nearest Neighbors 2

Classification Using Nearest Neighbors 3

Classification Using Nearest Neighbors 4

Classification Using Nearest Neighbors 5

Linear Regression

Linear Regression 2

Linear Regression 3

Linear Regression 4

Nonlinear Regression

Nonlinear Regression 2

Visualizing Multivariate Data

Generalized Linear Models

Generalized Linear Models 2

RegressionTree class

RegressionTree class 2

Neural networks

Gaussian Process Regression Models

Gaussian Process Regression Models 2

Understanding Support Vector Machine Regression

Understanding Support Vector Machine Regression 2

RegressionEnsemble



matlab assignment help


matlab assignment help