Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- N = 600;
- x = 1:N;
- mean = 61;
- scale = 3;
- noisy = random('Rayleigh', scale, 1, N) + mean;
- drop = 0.5/100; % packet drop chance
- for i = x
- if drop > rand()
- noisy(i) = 3*noisy(i);
- end
- end
- alpha = 0.5;
- b = alpha;
- a = [1, -(1-alpha)];
- filtered5= filter(b, a, noisy);
- alpha = 0.1;
- b = alpha;
- a = [1, -(1-alpha)];
- filtered1= filter(b, a, noisy);
- plot(x, noisy, '-b', x, filtered5, '-r', x, filtered1, '-g');
- legend('input', 'filter a=0.5', 'filter a=0.1')
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement