warning_id = 'Simulink:blocks:TDelayDirectThroughAutoSet';
warning('off',warning_id)

freq_linolg = 10.^(-2:0.01:6);
io(1) = linio([noiseModel '/In1'],1,'input');
io(2) = linio([noiseModel '/From1'],1,'output');
opt = linearizeOptions;
opt.UseExactDelayModel = 'on';

SW = [0 1 1 1 0 0 0 0];
disp(' ')
disp(' ')
disp(' ')
disp('linearlizing FSS total path OLG')
%linearize
fss = linearize(noiseModel,io,opt);

SW = [0 1 0 0 0 0 0 0];
disp('linearlizing FSS PZT path OLG')
%linearize
pzt = linearize(noiseModel,io,opt);

SW = [0 0 1 0 0 0 0 0];
disp('linearlizing FSS EOM path OLG')
%linearize
eom = linearize(noiseModel,io,opt);
lin_model_c_fss = mnbode_list(fss,freq_linolg,'c');
lin_model_c_pzt = mnbode_list(pzt,freq_linolg,'c');
lin_model_c_eom = mnbode_list(eom,freq_linolg,'c');

%plot
figure()
mnbode(freq_linolg,abs(lin_model_c_fss),mnangle(lin_model_c_fss),...
    freq_linolg,abs(lin_model_c_pzt),mnangle(lin_model_c_pzt),...
    freq_linolg,abs(lin_model_c_eom),mnangle(lin_model_c_eom),...
   'title','FSS OLG','linestyle',{'-','--','--'},'legend',...
   {'total','PZT','EOM'},'ylim',[1e-1 1e6]);


io(1) = linio([noiseModel '/In5'],1,'input');
io(2) = linio([noiseModel '/From5'],1,'output');
opt = linearizeOptions;
opt.UseExactDelayModel = 'on';
SW = [1 1 1 1 0 1 1 0];
disp('linearlizing IMC total path OLG')
%linearize
imc = linearize(noiseModel,io,opt);

SW = [1 1 1 1 0 1 0 0];
disp('linearlizing IMC AOM path OLG')
%linearize
aom = linearize(noiseModel,io,opt);

SW = [1 1 1 1 0 0 1 0];
disp('linearlizing IMC suspension path OLG')
%linearize
sus = linearize(noiseModel,io,opt);

warning('on',warning_id)

lin_model_c_imc = mnbode_list(imc,freq_linolg,'c');
lin_model_c_aom = mnbode_list(aom,freq_linolg,'c');
lin_model_c_sus = mnbode_list(sus,freq_linolg,'c');




% close all
%plot
figure()
mnbode(freq_linolg,abs(lin_model_c_imc),mnangle(lin_model_c_imc),...
    freq_linolg,abs(lin_model_c_aom),mnangle(lin_model_c_aom),...
    freq_linolg,abs(lin_model_c_sus),mnangle(lin_model_c_sus),...
   'title','IMC OLG','linestyle',{'-','--','--'},'legend',...
   {'total','AOM','suspension'},'ylim',[1e-3 1e6],...
   'linewidth',{2, 1, 1});