Пример выполнения задания

Пусть задано отображение Эно: ; .    

Программа для построения отображения и фазового портрета:

a=1.4;

b=.3;

x(1)=0;

n=500;

y(1)=0;

for i=2:n

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

y(i)=b*x(i-1);

 end

plot(x),grid;

title('x of Henon descret map');

xlabel('i');ylabel('x');

 figure;

plot(y),grid;

title('Y of Henon descret map');

xlabel('i');ylabel('Y');

figure;

i=2:n;

plot(x(i-1),x(i)),grid;

title('phase portrait for Henon descret map');

xlabel('x(i)');ylabel('x(i+1)');

 

 

Программа для построения бифуркационной диаграммы при b=0.3:

x0=0;

y0=0;

b=.3;

amin=0;

amax=1.4;

astep=0.001;

n=1000; nplot=500;

nstep=(amax-amin)/astep+1;

t=zeros(nstep+1,1);

z=zeros(nstep+1,nplot);

for j=1:nstep+1

 x=zeros(nstep+1,1);

 x(1)=x0;

 t(j)=(j-1)*astep+amin;

 a=t(j);y(1)=y0;

 for i=1:n

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

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

if (i>n-nplot)

    z(j,i-(n-nplot))=y(i+1);

end

 end

end

plot(t,z,'b.','MarkerSize',2), grid

title('Bifurcation diagram for Henon map')

xlabel('a','Fontsize',12), ylabel('Population','Fontsize',12)

 

Результаты работы программы показаны на рисунках 1-3.

 

Рисунок 1.- Отображение Эно               Рисунок 2.- Бифуркационное дерево

 

 

 Рисунок 3.-  Фазовый портрет      

 

 

Лабораторная работа № 4. Вейвлет-анализ динамики системы

Цель работы: научиться выполнять вейвлет-анализ в среде Matlab.

Задание. Выполнить дискретный и непрерывный вейвлет-анализ динамики системы.  

 

Таблица 1. - Варианты заданий

№	Файл данных
1	1.txt
2	2.txt
3	3.txt
4	4.txt
5	5.txt
6	6.txt