- #1

joshmccraney

Gold Member

- 1,957

- 102

Code:

```
global delta alpha beta gamma OMEG
delta = 0; % DAMPING
alpha = 50; % STIFFNESS
beta = 0; % RESTORATION NONLINEARITY
OMEG = 1.2; % DRIVING FORCE ANGULAR FREQ
gamma = 0; % FORCING
Fs = 100; % Sampling frequency
T = 1/Fs; % Sampling period
L = 10000; % Length of signal
[t, x] = ode45(@duffing,(0:L-1)*T,[0 1]);
plot(t,x(:,2))
```

Code:

```
function xdot = duffing(t,x)
global delta alpha beta gamma OMEG
xdot(1) = -delta*x(1) - alpha*x(2) - beta*x(2)^3 + gamma*cos(OMEG*t);
xdot(2) = x(1);
xdot = xdot';
end
```