Programación iterativa de problemas

Gracias a la amabilidad de mi amigo de arriba, copié palabra por palabra la respuesta a una pregunta que respondí hace dos años y medio (Nº 6876467041069324). De hecho, ese código fue escrito por mí para observar el proceso iterativo de ciertos tipos de líneas de gravedad de Hénon y no cumplía con los requisitos del problema en ese momento.

El trasfondo de esta pregunta es un modelo de sistema dinámico de tiempo discreto llamado mapeo de Hénon sobre el caos.

Entre ellos, el clásico diagrama de Hénon lleva una? =?1.4,b? =? 0.3, estos son los datos dados en la pregunta.

(1) Escriba funciones según sea necesario:

¿Función? 【incógnita,? y]? =?henon(x0,?¿y tú?n)

x? =?[0;? ceros(N,1)*NaN];

y? =?x;

¿Para qué? ¿I? =?1?:?Normal

x(i+1)? =?1?+?y(yo)? -?1.4*x(i)^2;

y(i+1)? =?0.3?*?x(一);

Los parámetros de entrada de end son el valor inicial x0, y0 y el número de iteraciones, y la salida es el resultado de la iteración.

(2) Función de llamada

x0? =?0.1;? ¿y0? =?0.2;

[x,? y]? =?henon(x0,?¿Y tú?100);

disp([x(end)?Y(end)]) ¿Cuál es el resultado? xk = -1,2493, yk = 0,3805.

(3) Función de llamada

x0? =?0;? ¿y0? =?0;

[x,? y]? =?henon(x0,?¿y tú?30000);

plot(x,y,?'.')

xlabel? incógnita,? ¿Ilabell? y