ttt400

1. Configure G=1/3*((5*Pp-2*P1)/(P1-2*Pp)-P1*c1^2/Pp*cp^2); p> p>

Cambiar a G=1/3*((5*Pp-2*P1)/(P1+2*Pp)-P1*c1^2/(Pp*cp^2));

El error en la expresión provoca una oscilación de alta frecuencia y el tamaño del paso de integración es muy pequeño (aproximadamente 10 ^ -8), por lo que no se puede calcular durante mucho tiempo.

2. A juzgar por los resultados, 50 segundos de tiempo de simulación son suficientes (ha entrado en un estado estable).

3. Cuando ocurre ocupado, puede considerar dejar que la función ode* no devuelva parámetros. En este caso, la curva de respuesta se dibujará automáticamente y la operación se puede detener en cualquier momento a través de la GUI. Facilitar la observación del proceso de solución.