main_danchunxin.m

基于单纯形无约束方法的PID参数整定。在myfun.m中输入目标函数和控制对象的传递函数

%main_danchunxin.m
clear all
close all
clc
N = 50;
%input
r_input = ones(1,N);
y_out = zeros(1,N);
u_ctrl = zeros(1,N);
e_error = r_input - y_out;
B0 = [2,1,1]';
B1 = [1,-1,7]';
B2 = [1,9,-1]';
B3 = [-10,1,1]';
Init_B = [B0,B1,B2,B3];
figure;
[fx1,fy1,iter,out1,u_con1] = DanchunXin(Init_B,3,N,u_ctrl,y_out,r_input,e_error,1);
plot(out1);grid on
hold on
[fx2,fy2,iter2,out2,u_con2] = DanchunXin(Init_B,3,N,u_ctrl,y_out,r_input,e_error,0);
plot(out2,'r')
hold on
plot(ones(1,N),'k');
tr1 = find(abs(out1-0.5)==min(abs(out1-0.5)))
tr2 = find(abs(out2-0.5)==min(abs(out2-0.5)))
figure;
plot(u_con1);grid on
hold on
plot(u_con2,'r')