Hola amigos, estoy trabajando con una ecuación de Matthieu no lineal, que simula un una viga con un motor desequilibrado. La he pasado a frecuencia para una comprobación de datos empiricos y la querría representar. He generado este código:
for t=1:1000;
Y1(1)=0;
x=(t/(4*122.13)); %frecuencia natural 122.13
X1(t)=x;
y = fzero(@(y) ecuale(x,y),Y1(1))
Y1(t)=y;
end
plot(X1,abs(Y1))el problema es que hay una zona de la ecuación con múltiples raices reales, en concreto 3, y la función fzero solo me calcula una (¿o estoy equivocado?). Podrían aconsejarme para suavizar la gráfica y obtener el resto de raices.
Muchas gracias, la función con la que trabajo es :
function xdot=ecuale(x,y)
global F a b
F=[(0.4*0.04)/(10.16*122.13)];
a=(3*1.76e9)/(4*122.13^2);
b=2*0.007;
xdot=[((1-x^2)*y)+(a*y^3)]^2+[b*x*y]^2-F;