distubplotpi.m

过程系统建模和PID整定ggg gjkjkjjjhgg g 、

%  ====================================
%    互偶水槽PID整定比较和响应指标计算
%  ====================================
%    特别注意：
%           将干扰的检测点设为2，切记。
%   清屏，清理工作空间，命令窗口
clear all;
close all;
clc

%  PI整定效果数据加载
wwa=load('responsedist');    %  加载响应曲线法整定PID的数据                                 
wwb=load('SMCdist') ;        %  加载SMC法整定PID的数据
wwc=load('boundPIdist');

wwa(:,2)=wwa(:,5);    %  T2高度传感器检测数据放到wwa（：，2）                         
wwb(:,2)=wwb(:,5);    %  T2高度传感器检测数据放到wwa（：，2）
wwc(:,2)=wwc(:,5);    %  T2高度传感器检测数据放到wwa（：，2）

%   确定wwa加干扰时间为t10
for i=1:length(wwa(:,3))-1
    if wwa(i,3)~=wwa(i+1,3)
        t10=wwa(i,1);
    end
end

%   确定wwb加干扰时间为t20
for i=1:length(wwb(:,3))-1
    if wwb(i,3)~=wwb(i+1,3)
        t20=wwb(i,1);
    end
end

%   确定wwc加干扰时间为t30
for i=1:length(wwc(:,3))-1
    if wwc(i,3)~=wwc(i+1,3)
        t30=wwc(i,1);
    end
end


x1=wwa(:,1)-t10;   %  修正的X轴
x2=wwb(:,1)-t20;   %  修正的X轴
x3=wwc(:,1)-t30;   %  修正的X轴

%  两种PID响应曲线比较（数据未滤波）
figure;
plot(x1,wwa(:,2),'b',x2,wwb(:,2),'r',x3,wwc(:,2),'g')    
grid;
legend('响应曲线法整定PI干扰响应','SMC整定PI干扰响应','稳定边界法PI干扰响应',2);
title('PI整定效果响应曲线（未滤波）')
xlabel ('Time(s)');
ylabel('Respones(V)')