cstr_rk4.m

本程序用来实现CSTR模型的模糊PID控制仿真

% 本程序用四阶龙格－库塔方法实现CSTR模型的仿真的递推计算
% 
%   完成日期：2007年4月22日

function [Ca, T, Tc, H] = CSTR_RK4(Q, Qc, Ca, T, Tc, H)
    %   仿真步长  
    h = 0.05;
    
    [K1,L1,M1,N1] = CSTR_MD(Q, Qc, Ca, T, Tc, H);
    [K2,L2,M2,N2] = CSTR_MD(Q, Qc, Ca+0.5*h*K1, T+0.5*h*L1, Tc+0.5*h*M1, H+0.5*h*N1);
    [K3,L3,M3,N3] = CSTR_MD(Q, Qc, Ca+0.5*h*K2, T+0.5*h*L2, Tc+0.5*h*M2, H+0.5*h*N2);
    [K4,L4,M4,N4] = CSTR_MD(Q, Qc, Ca+h*K3, T+h*L3, Tc+h*M3, H+h*N3);

    Ca = Ca + h/6*(K1+2*K2+2*K3+K4);
    T  = T + h/6*(L1+2*L2+2*L3+L4);
    Tc = Tc + h/6*(M1+2*M2+2*M3+M4);
    H  = H + h/6*(N1+2*N2+2*N3+N4);
end