PCA SPE T2 Contribution Simulink matlab

matlab不可多得的代码 利用simulink仿真数据 涉及到PCA，SPE及其控制限，T2及其控制限，以及在超限时原始变量的贡献图。不管对算法或者GUI的制作都有很大的帮助意义。本人使用的时matlab2017b，不保证其他版本能否打开，谨慎下载！！

function xprime=colamod（tXU） 
%
% colamod - This is a nonlinear model of a distillation column with
%           NT-1 theoretical stages including a reboiler （stage 1） plus a
%           total condenser （“stage“ NT）. The liquid flow dynamics are
%           modelled by a simple linear relationship.
%           Model assumptions: Two components （binary separation）; constant
%           relative volatility; no vapor holdup; one feed and two products;
%           constant molar flows （same vapor flow on all stages）; 
%           total condenser
%
%           The model is based on column A in Skogestad and Postlethwaite
%           （1996）. The model has 82 states.
%
% Inputs:    t    - time in [min].
%            X    - State the first 41 states are compositions of light
%                   component A with reboiler/bottom stage as X（1） and 
%                   condenser as X（41）. State X（42）is holdup in reboiler/
%                   bottom stage and X（82） is hold-up in condenser. 
%            U（1） - reflux L
%            U（2） - boilup V
%            U（3） - top or distillate product flow D
%            U（4） - bottom product flow B
%            U（5） - feed rate F
%            U（6） - feed composition zF.
%            U（7） - feed liquid fraction qF.
%
% Outputs:   xprime - vector with time derivative of all the states 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%------------------------------------------------------------
% The following data need to be changed for a new column.
% These data are for “colmn A“.
% Number of stages （including reboiler and total condenser: 
    NT=41; 
% Location of feed stage （stages are counted from the bottom）:
    NF=21;
% Relative volatility
    alpha=1.5;
% Nominal liquid holdups
    M0（1）=0.5;      	% Nominal reboiler holdup （kmol）
    i=2:NT-1; M0（i）=0.5*ones（1NT-2）;% Nominal stage （tray） holdups （kmol）
    M0（NT）=0.5;      	% Nominal condenser holdup （kmol）
% Data for linearized liquid flow dynamics （does not apply to reboiler and condenser）:
    taul=0.063;     	% time constant for liquid dynamics （min）
    F0=1;	 	% Nominal feed rate （kmol/min） 
	load cola_init
    qF0 = Uinit（7）;%1; 		% Nominal fraction of liquid in feed 
    L0=Uinit（1）;%2.70629;     	% Nominal reflux flow （from steady-state data）
    L0b=L0 + qF0*F0;	% Nominal liquid flow below feed （kmol/min）
    lambda=0;		% Effect of vapor flow on liquid flow （“K2-effect“）
    V0=Uinit（2）;%3.20629;
    V0t=V0+（1-qF0）*F0;% Nominal vapor flows - only needed if lambda is nonzero 
% End data which need to be changed
%------------------------------------------------------------

% Splitting the states
x=X（1:NT）‘;                          % Liquid composition from btm to top
M=X（NT+1:2*NT）‘;                     % Liquid hold up from btm to top

% Inputs and disturbances
LT = U（1）;                            % Reflux
VB = U（2）;                            % Boilup
D  = U（3）;                            % Distillate
B  = U（4）;             

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        860  2018-01-17 14:50  A benchmark software for multivariate statistical process control\.addOnmetadata\metadata.xml

     文件    8684315  2018-01-17 14:51  A benchmark software for multivariate statistical process control\.addOnmetadata\MSPC.1.6.zip

     文件      14693  2018-01-17 14:51  A benchmark software for multivariate statistical process control\.addOnmetadata\previewImage.png

     文件       4511  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\colamod.m

     文件       1746  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\colas.m

     文件      12281  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\colas_PI.m

     文件     121099  2018-01-17 15:33  A benchmark software for multivariate statistical process control\code\colas_PI.mdl

     文件     114873  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\colas_PI.mdl.r2009b

     文件        774  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\cola_init.mat

     文件        774  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\cola_init_PO1.mat

     文件        705  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\cola_init_PO2.mat

     文件        700  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\cola_init_PO3.mat

     文件        699  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\cola_init_PO4.mat

     文件       1426  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\cola_lv.m

     文件       2115  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\CreateFeedSignal.m

     文件    8539087  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\example.cdb

     文件       2802  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\InterpTemp.m

     文件       1304  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\license.txt

     文件     126218  2018-01-17 16:26  A benchmark software for multivariate statistical process control\code\MSPC_main.fig

     文件     127020  2018-01-17 16:26  A benchmark software for multivariate statistical process control\code\MSPC_main.m

     文件       7441  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\MSPC_scores.fig

     文件      12614  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\MSPC_scores.m

     文件       5884  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\MSPC_SPE.fig

     文件       7814  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\MSPC_SPE.m

     文件       8554  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\MSPC_T2.fig

     文件      10897  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\MSPC_T2.m

     文件       4116  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\MSPC_variance.fig

     文件       4356  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\MSPC_variance.m

     文件       4258  2018-01-17 14:51  A benchmark software for multivariate statistical process control\code\qFactor.m

     文件     453850  2018-01-17 15:41  A benchmark software for multivariate statistical process control\code\Sim_20180117145154.cdb

............此处省略10个文件信息