>> % Simulate data of two loops
>> sys=tf([0 0.5],[1 -0.5],1); % generate process dynamics
>> t=0:100;
>> u1=sin(0.1*t);
>> y1=lsim(sys,u1); % simulate loop1 data
>> t=0:500;
>> u2=cos(0.2*t);
>> y2=lsim(sys,u2); % simulate loop2 data
>>%--------------------------------------------------------

>> gen_cmpct_Data_stiction    % run data conversion code
 
Consider a control-loop with control valve
 
1. Input data (OP): N x 1 matrix
2. Output data (PV): N x 1 matrix
3. Sampling time (sec)
4. Comment on the control loop
 
How many control loops? 2
[Loop No 1] Input data (OP) vector: u1
[Loop No 1] Output data (PV) vector: y1
[Loop No 1] Sampling time (sec) (Just enter if unknown): 
[Loop No 1] Comment on the control loop: loop1
 
[Loop No 2] Input data (OP) vector: u2
[Loop No 2] Output data (PV) vector: y2
[Loop No 2] Sampling time (sec) (Just enter if unknown): 
[Loop No 2] Comment on the control loop: loop2
 
A file name to save compact data (e.g. test_cmpct_data_stiction): test_data

% Now we are ready to run stiction detection/quantification code main_UA_stiction.m and load the data generated