This example shows how to use the Sensitivity Analysis tool to explore the behavior of a PI controller for a DC motor. The controller is susceptible to variations caused by component tolerances, and the impact on controller reliability is explored.
You explore the controller reliability by characterizing the components using probability distributions. You use the distributions to generate random samples and perform Monte-Carlo evaluation of the controller design at these sample points. You evaluate the impact of the component tolerances on the controller behavior, and use statistical analysis to determine which components have the most influence on whether the controller meets its requirements. This analysis guides the selection of component tolerances.
This example requires Statistics and Machine Learning Toolbox™.
The controller enables the DC motor's angular position to match a desired reference value. The load on the motor is subject to disturbances, and the controller needs to reject these disturbances. The Simulink model can be used to probe how well the controller rejects a step disturbance at 1 second.
The gains of the PI controller,
Ki, are set using resistors in the circuit below:
R4 are 47 kOhm, 180 kOhm, 10 kOhm and 10 kOhm respectively. These were chosen to set
Ki to values that enable the controller to meet the requirements for disturbance rejection. However in practice, the actual resistor values will differ from the nominal ones, within a tolerance. This raises concern about whether the actual controller will still satisfy the requirements. To explore the effect of different resistance values, use the Sensitivity Analysis tool. In the Simulink Analysis menu, click Sensitivity Analysis... to open the tool.
The controller needs to maintain the motor at a reference position in the presence of disturbances. If a step disturbance occurs, the motor needs to deviate no more than 20 degrees, and needs to settle back to within 5 degrees of the reference position by 4 seconds after the disturbance.
Load previously specified design requirements for disturbance rejection. In the tool, click Open Session and select Open from model workspace in the drop-down menu.
You can plot the requirements and verify that they are met when the resistances have the nominal values. In the Results area in the data browser, right-click on the
LowerBound requirement, and select Plot and Simulate. Do the same for the
The motor position satisfies the disturbance rejection requirements, when the resistances are at their nominal values. However in practice, the actual resistor values will differ from the nominal ones, and we need to determine whether the controller will still meet the requirements. Click Select Parameters and make a new parameter set. This creates
ParamSet in the Parameter Sets area of the tool. Specify that
R4 are in the parameter set, and click OK.
Click Generate Values and generate random values. For repeatable results, reset the state of the random number generator in MATLAB®.
In the Generate Random Parameters dialog box, specify 500 samples to generate.
Specify the probability distribution for each parameter. Standard precision resistors match their nominal component value within a tolerance of 5%. This could be modeled using a uniform probability distribution. However, because resistors that measure within 1% of the nominal value are separated out and sold as higher-priced precision resistors, the 5% resistors can be more accurately modeled by a probability distribution with a well that excludes values within 1% of nominal. This can be modeled using a piecewise linear probability distribution if Statistics and Machine Learning Toolbox™ is available.
Specify the distribution of
R1 as piecewise linear with 4 points. Specify the x values as [0.95 0.99 1.01 1.05] times 47 (the nominal value of the resistor). Specify the Fx values as [0 0.5 0.5 1]; these are the values of the cumulative distribution function corresponding to each x value. Similarly, set the distributions of
R4 to piecewise linear with 4 points, the x values as [0.95 0.99 1.01 1.05] times the nominal values (180, 10, and 10, respectively), and the Fx values as [0 0.5 0.5 1].
Click OK to generate parameter values. The generated values are stored in the
ParamSet variable in the Parameter Set area of the tool. (Note that due to the random number generator, the specific values in the table below may differ from what you get when running the example.)
To plot the parameter set click
ParamSet in the Parameter Sets area of the tool browser. In the Plots tab, select Scatter Plot in the plot gallery. The plot shows the histogram of the generated parameters on the diagonal and pair-wise parameter scatter plots off the diagonal. Each marker on the plot represents one row of the
ParamSet table, with each row being simultaneously displayed on all the scatter plots. You can use the View tab to arrange the layout of the table and plot so they are both visible.
Evaluate the requirements for each row of parameter values in the table to see if the requirements are satisfied. In the Sensitivity Analysis tab, click Select for Evaluation. By default, all requirements are selected to be evaluated. Click Evaluate Model to evaluate the
LowerBound requirements for each row of parameter values in
ParamSet. Note you can speed up evaluation by using parallel computing if you have the Parallel Computing Toolbox™, or by using fast restart. For more information, see "Use Parallel Computing for Sensitivity Analysis" and "Use Fast Restart Mode During Sensitivity Analysis" in the Simulink Design Optimization™ documentation.
A results scatter plot showing each requirement vs. each parameter is updated during model evaluation. At the end of evaluation a table with the evaluation results is shown. Each row in the evaluation result table contains values for
R4 and the resulting requirement values
LowerBound. The evaluation results are stored in the
EvalResult variable in the Results area of the tool. You can use the View tab to arrange the layout of the table and plot so they are both visible.
You can sort the evaluation results table by clicking on the column headers in the table. The
LowerBound requirement is still met, as indicated by the fact that all evaluation results for the signal bound requirement are negative. That is not the case for the
UpperBound requirement, which has several positive values. By selecting the rows of the table with these positive values, you can also see the corresponding points highlighted in the scatter plot.
Using 5% tolerance components resulted in violation of the
UpperBound requirement. Precision components with 1% tolerance would satisfy the design requirements, but they are more costly, so it is desirable to use only as many as necessary. You can use statistical analysis to identify the components that most influence the design requirements.
In the Statistics tab, select a variety of analyses to be done, including Correlation and Standardized Regression methods, and Linear and Ranked types of processing. Click Compute Statistics. The analysis result is stored in
StatsResult in the Results area of the tool, and a tornado plot shows the analysis results. For each requirement, the tornado plot shows the most influential parameters at the top, and the others in decreasing order of the magnitude of their influence on the requirement. For the
R4 have the most influence, so we will try replacing these by higher precision 1% components.
Explore the use of 1% component tolerances only for resistors
R4. In the Sensitivity Analysis tab, click Generate Values and generate gridded values. For
R2, specify that the nominal value is to be perturbed by plus-and-minus 5%. For
R4, specify that the nominal value is to be perterbed by plus-and-minus 1%.
Click Overwrite to generate the new parameter values. To plot the parameter set click
ParamSet in the Parameter Sets area of the tool browser. In the Plots tab, select Scatter Plot in the plot gallery.
In the Sensitivity Analysis tab, click Evaluate Model. The requirements are evaluated for each row in the table of parameter values, and results are stored in
EvalResults_1 shown in the Results area of the tool. The evaluation results scatter plot and the evaluation results table show that both requirements are met for all combinations of component values.
The Sensitivity Analysis tool was used to explore the effect of standard precision components on the design requirements of a PI controller. With standard precision components, some requirements were found to be violated. Statistical analysis was used to identify which parameters most influence the requirements. The analysis resulted in replacement of only two of the four components with most costly high-precision components.
Close the model.