twostate.m

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

%  ====================================
%     计算两个稳态值即：Y（0）和Y（inf）
%  ====================================
%  输入：y1——滤波后的曲线Y值
%        wwa——响应曲线全部数据
%  输出：mm1——Y(0)
%        mm3——Y(inf)

function [mm1,mm3]=twostate(y1,wwa)

%  求两稳态输出
%  Y(inf)=mm3
%  Y(0)=mm1

%  将Y*1000，然后将其统计出来在temp1
L1=length(y1);
temp2=round(y1*1000);
temp1=zeros(max(temp2+1000),1);
L2=length(temp1);
for i1=1:L1
    if(temp2(i1,1)>0)
       temp1(temp2(i1,1)+1000,1)=temp1(temp2(i1,1)+1000,1)+1;
    end
end

%  对temp1数据提取稳态值所对应的Y*1000
i=0;
for i2=1:L2    
    if temp1(i2,1)>=100 
        i=i+1;
        sum1(i)=i2-1000;
     end
end


%  求稳态出平均值作为Y（0）和Y（inf）

mm=0;
for i3=1:length(sum1)-1
      mm=sum1(1,i3)+mm;
      if (abs(sum1(i3+1)-sum1(i3))>=20)       % 20为稳态区分阀值
        mm1=mm/i3;
        mm=0;
        mm2=0;
         for i4=(i3+1):length(sum1)
             mm2=sum1(i4)+mm2;
         end
         mm3=mm2/(i4-i3);
      end 
end
mm1=mm1/1000;
mm3=mm3/1000;