Skip to content

Fractional order proportional derivative controller tuner

License

Notifications You must be signed in to change notification settings

cnpcshangbo/FOPD-tuner

Repository files navigation

FOPD-tuner

Fractional order proportional derivative controller tuner

Related Publications:

Bo Shang, Jianxin Liu, Tiebiao Zhao and YangQuan Chen, "Fractional order robust visual servoing control of a quadrotor UAV with larger sampling period," 2016 International Conference on Unmanned Aircraft Systems (ICUAS), Arlington, VA, 2016, pp. 1228-1234. PDF

1. License

FOPD-tunner is released under a MIT license.

If you use FOPD-tunner in an academic work, please cite:

@inproceedings{shang2016fractional,
  title={Fractional order robust visual servoing control of a quadrotor UAV with larger sampling period},
  author={Shang, Bo and Liu, Jianxin and Zhao, Tiebiao and Chen, YangQuan},
  booktitle={Unmanned Aircraft Systems (ICUAS), 2016 International Conference on},
  pages={1228--1234},
  year={2016},
  organization={IEEE}
}

2. Prerequisites

We have tested the library in MATLAB 2018a and Windows 7, but it should be easy to run in other platforms.

Control model

screenshot 2017-08-15 17 51 47

MATLAB System Identification Toolbox

Convert the identified model to this format:

screenshot 2017-08-15 20 25 37

Then we can get the parameter k and tau.

FOTF Toolbox

https://www.mathworks.com/matlabcentral/fileexchange/60874-fotf-toolbox

Oustaloup

3. Downloading FOPD-tunner library and run examples

Clone the repository:

git clone https://github.com/cnpcshangbo/FOPD-tunner.git

4. Run examples

Tune IOPID controller

Enter the folder named "1. iopid_tune_pid", open "iopid_tune.m", update the plant parameters identified with MATLAB System Identification Toolbox and set cross-over frequency omega_c and phase margin phi.

%% Setting plant parameters
global K T1 wc phi %inputs
K = -0.8592; %plant gain
T1 = 1.0710; %plant time constant
wc = 1.2; %crossoverfrequency
phi = 83.9; %phase margin

then run it to tune IOPID controller and get related parameters (integer order PID parameters kp, ki, kd, ).

Tune FOPD controller

Setting parameters

Open file "cal_lambda_ki_kp_bode.m". Set fractional order lambda to -1~0.

Note: When lambda is negative, a FOPI controller becomes a FOPD controller.

Plotting ki against lambda

Run the first cell of "cal_lambda_ki_kp_bode.m". You will see the three curves. Zoom the curves to find crosspoints. The coordinates of the crosspoints are lambda and ki.

Verify the crosspoints

Run the "Solution 1" cell of "cal_lambda_ki_kp_bode.m" to double check if the crosspoint is the REAL solution of the equations. If ki is complex number, there's no real solution. It means the phase curve at cross-over frequency can not be flat under current cross-over frequency omega_c and phase margin phi. Either change them or use optimization method to make the derivate of phase the smallest at cross-over frequency.

Verify the tuned FOPD controller

Run the "Verify the controller 1" cell of "cal_lambda_ki_kp_bode.m" to verify if the phase curve is flat at crossover frequency. At the same time kp will be displayed. The we can get the transfer function of the tuned FOPD controller.

Controller analysis with Simulink

Open the file "run_simulink.m", update the parameters and run the script. image image

Co-Simulation with Gazebo

We use the Robotics System Toolbox to cooperate with Gazebo. The Gazebo part follows Erlerobotics simulation.

Follow this document to take off the drone:

http://docs.erlerobotics.com/simulation/vehicles/erle_copter/tutorial_5 drone taken off

Launch MATLAB

Install Robotics System Toolbox. Initial ROS in MATLAB.

ROS with MATLAB node

Use the ROS Publish block in Simulink to publish topic for velocity setpoint. ROS Publish block

Use Simulink to move a drone in Gazebo

About

Fractional order proportional derivative controller tuner

Topics

Resources

License

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published