Cálculo de programación de primavera

Oye, soy viejo y Matlab está oxidado.

Hermano, trabajé duro toda la tarde para ti por capricho, ¡por favor dame puntos extra!

El código del programa es el siguiente:

β= 0.1

w0 = 1

delta_t = 25/1000;

t(1)= 0;

Porque i=1:1000

t(I 1)= delta _ t * I;

End

Umbral = 0.9;

x(1)= 100; Para resolver el problema de manera más conveniente, finalmente cambie la escala.

temp_v 1 = 0;

n = 1; esquina de x

Porque i=1:1000

Para j =- 100:0.001:100

temp_x = j;

temp_v2 =(temp_x-x(n))/delta_t;

temp_a =(temp_v2-temp_v1) /delta_t;

Si abs(temp_a 2*beta*temp_v2 w0^2*temp_x)lt; umbral

n = n 1;

x(n )= temp_x;

temp_v 1 = temp_v2;

Pausa;

Fin

Para i=1:1001 vuelve a la escala real.

x(I)= x(I)/100;

Fin

plot(t, x);

xlabel( ' Tiempo ');

ylabel('Distancia X');

Título ('Curva de vibración amortiguada del oscilador de resorte'); oscilador La curva se puede obtener mediante cálculo directo.