%% ABOUT THIS FILE
% ----------------------------------------------------------------------
% Type-A with Dummy Mass for KAGRA
% Coded by K. Okutomi on 2016.12.14
% ----------------------------------------------------------------------

%% PRELIMINARY
clear all;  % Clear workspace
close all;  % Close plot windows
addpath('../../utility');  % Add path to utilities
addpath('controlmodel');
addpath('susmodel');
g = 9.81;

%% IMPORT SUSPENSION MODEL
matfile='TypeA161114mdl';
load([matfile,'.mat']);

%%
%%%%% NO-CONTROL MODEL %%%%%
%% IMPORT SERVO FILTERS
addpath('parameter');         % Add path to servo
TypeA_paramNoCtrl161214;      % NO CONTROL MODE
rmpath ('parameter');         % Remove path to servo

%% IMPORT SIMULINK MODEL
mdlfile='TypeA_ctrl161209';
st     =linmod(mdlfile);
invl   =strrep(st.InputName, [mdlfile,'/'],'');
outvl  =strrep(st.OutputName,[mdlfile,'/'],'');
sysc0  =ss(st.a,st.b,st.c,st.d,'inputname',invl,'outputname',outvl);

%% FREQUENCY
freq =logspace(-2,1,1001);
freq2=logspace(-3,1,1001);