clc; clear all; close all;
format short

%tiempo de simulación (segundos)
ti=0; h=0.0025; tf = 5; 
ts=ti:h:tf;  %vector tiempo
opciones=odeset('RelTol',1e-3,'InitialStep',2.5e-3,'MaxStep',2.5e-3);
[t,x]=ode45('filtro',ts,0,opciones);
u=5*sin(t)+random('Normal',t,t);
subplot(2,2,1); plot(t,u)
subplot(2,2,2); plot(t,x)