function Sys_ID(Ws,Vs,Fs,n)
%
% Author: K.Srikanthan
% Hardware and Systems Group - Automotive MBD
% Tata Elxsi Limited
%
% Inputs:
% Ws - Output structure of the motor or the plant (Recorded)
% Vs - Input structure of the motor or the plant (Recorded)
% Fs - Sampling frequency
% n - order of the system - (must be entered as 2 for motor)
W = Ws.signals.values;
V = Vs.signals.values;
[Num, Den] = Find_Parameters(W,V,Fs,n);
Num_coeff = mat2str(Num);
Den_coeff = mat2str(Den);
simulink;
sfnew('System_Identification')
delete_block('System_Identification/Chart')
add_block('simulink/Discrete/Discrete Transfer Fcn', 'System_Identification/Transfer_Function')
set_param('System_Identification/Transfer_Function','Numerator',Num_coeff,'Denominator',Den_coeff, 'SampleTime','1/Fs')
set_param('System_Identification/Transfer_Function','Position',[250 150 550 250])
add_block('simulink/Sources/From Workspace','System_Identification/Input','variablename','Vs','Sampletime','1/Fs','interpolate','off','outputafterfinalvalue','holdingfinalvalue','position',[50 185 100 215]);
add_block('simulink/Sinks/To Workspace','System_Identification/Output','variablename','Wcap','Sampletime','1/Fs','saveformat','structurewithtime','Position',[700 185 750 215]);
add_block('simulink/Sources/From Workspace','System_Identification/Recorded_Output','variablename','Ws','Sampletime','1/Fs','interpolate','off','outputafterfinalvalue','holdingfinalvalue','position',[350 95 400 125]);
add_block('simulink/Signal Routing/Mux','System_Identification/Mux','Position',[650 100 655 135])
add_block('simulink/Sinks/Scope','System_Identification/Scope_output','Position',[700 100 730 135],'limitdatapoints','off')
add_line('System_Identification','Input/1','Transfer_Function/1')
add_line('System_Identification','Transfer_Function/1','Output/1')
add_line('System_Identification','Recorded_Output/1','Mux/1')
add_line('System_Identification','Transfer_Function/1','Mux/2','autorouting','on')
add_line('System_Identification','Mux/1','Scope_output/1','autorouting','on')
sim('System_Identification',Vs.time(end))