nr. 477 - 2011 - Institut for Natur, Systemer og Modeller (NSM)

166 matlab Code
w_7,V_34,’.b’,w_7,V_35,’.b’,k,0,’black’)
axis([0 3.5e-5 -2 2 ])
Title(’Plot of eigenvalues for the Jacobian evaluated in the stable healthy rest state ’)
xlabel(’f_1’)
ylabel(’real(eig)’)
%------------------
figure
subplot(2,1,1); plot(w,real(V_1),’.b’,w,real(V_2),’.b’,
w,real(V_3),’.b’,w,real(V_4),’.b’,w,real(V_5),’.b’,
w_2,real(V_6),’.b’,w_2,real(V_7),’.b’,w_2,real(V_8),’.b’,
w_2,real(V_9),’.b’,w_2,real(V_10),’.b’,
w_3,real(V_11),’.r’,w_3,real(V_12),’.r’,w_3,real(V_13),’.r’,
w_3,real(V_14),’.r’,w_3,real(V_15),’.r’,w_4,
real(V_16),’.r’,w_4,real(V_17),’.r’,w_4,real(V_18),’.r’,
w_4,real(V_19),’.r’,w_4,real(V_20),’.r’,
w_8,real(V_36),’.r’,w_8,real(V_37),’.r’,w_8,real(V_38),’.r’,
w_8,V_39,’.r’,w_8,real(V_40),’.r’,k,0,’black’)
axis([0 3.5e-5 -4 4 ])
Title(’Plot of the eigenvalues for the Jacobian evaluated in...
the upper inflamed region’)
xlabel(’f_1’)
ylabel(’real(eig)’)
subplot(2,1,2); plot(w,imag(V_1),’.b’,w,imag(V_2),’.b’,
w,imag(V_3),’.b’,w,imag(V_4),’.b’,w,imag(V_5),’.b’,
w_2,imag(V_6),’.b’,w_2,imag(V_7),’.b’,w_2,imag(V_8),’.b’,
w_2,imag(V_9),’.b’,w_2,imag(V_10),’.b’
,w_3,imag(V_11),’.r’,w_3,imag(V_12),’.r’,w_3,imag(V_13),’.r’,
w_3,imag(V_14),’.r’,w_3,imag(V_15),’.r’,
w_4,imag(V_16),’.r’,w_4,imag(V_17),’.r’,w_4,imag(V_18),’.r’,
w_4,imag(V_19),’.r’,w_4,imag(V_20),’.r’,
w_8,imag(V_36),’.r’,w_8,imag(V_37),’.r’,w_8,imag(V_38),’.r’,
w_8,V_39,’.r’,w_8,imag(V_40),’.r’,k,0,’black’)%
%axis([0 3.5e-5 -4 4 ])
axis([2e-5 3e-5 -1 1 ])
xlabel(’f_1’)
ylabel(’img(eig)’)
figure
%
plot(real(V_4),imag(V_4),’*b’,real(V_5),imag(V_5),’*b’,
real(V_9),imag(V_9),’*b’, real(V_10),imag(V_10),’*b’,