Листинг 2. %Решение уравненея Ван-дер-Поля с внешней силой

%Решение уравненея Ван-дер-Поля с внешней силой

global w b L;

w=6;

b=5;

L=1;

x0=[2 0];

[T,Y]=ode23('vdp1',[0 100],x0);

subplot(211);

plot(T,Y(:,1),'k-',T,Y(:,2),'r');

xlabel('Время Т');

ylabel('');

legend('у1','у2');

grid on;

subplot(212);

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

41.Найти (написать программу) фазовую траекторию отображения Икеды:

X_n+1=1-ax²-by_n, y_n+1=x_n.

x=0; % задаём начальные значения переменных

y=0;

b=0.3; % задаём параметрические константы

for a=0:0.01:1.5; % открываем первый цикл для

%расчёта состояния системы при разных значениях

%параметрической константы

for i=1:100; %открываем второй цикл для расчёта состояний

%системы при каждом значении константы

y(i+1)=x(i);

x(i+1)=1-a.*x(i).^2-b.*y(i);

plot(a,y(i+1),'k.');

hold on

end;

end;

grid on

42.Покажите хаотические режимы (написать программу вывода колебаний) генератора хаоса:

dx/dt=-y-z

dy/dt=x+ay

dz/dt=b+(x-r)z

a=2.3; b=1; r=0.3