[A,B,C,D]=power_analyze('study2');
eig(A)


 
freq=0:1500;
w=2*pi*freq;
% [mag1,phase1]=bode(A,B,C,D,2,w);
% semilogy(freq,mag1(:,2));

[A,B,C,D]=power_analyze('study2');
[mag10,phase10]=bode(A,B,C,D,2,w);
semilogy(freq,mag1(:,2),freq,mag10(:,2));