end

hold off;

axis([0 5.5 0 0.6]);

ПРИЛОЖЕНИЕ Б

Программная реализация построения фазовых портретов волн пластической деформации

clear all;

alpha=1; beta=1; gamma=1;

f=@(t,y)[-y(1)+alpha*y(2);beta*(1-y(2)-gamma*y(1)*y(2))];

for i=0:1/5:1,

for j=0:1/5:1,

[T,Y]=ode45(f,[0 100],[i j]);

plot(Y(:,1),Y(:,2));

hold on;

end

end

hold off;

axis([0 1 0 1]);

pause;

clear all;

alpha=20.48;

beta=0.043;

gamma=25;

f=@(t,y)[-y(1)+alpha*y(2);beta*(1-y(2)-gamma*y(1)*y(2))];

for i=0:1/5:1,

for j=0:1/5:1,

[T,Y]=ode45(f,[0 100],[i j]);

plot(Y(:,1),Y(:,2));

hold on;

end

end