function measured = sinusoidal_with_noise ( time, A, f, phase, alpha, seed )
% �³Œ·”g + ƒmƒCƒY�iˆê—l—��”�j ‚ðŒvŽZ
% A sin (2 pi f t + delta) + U(-alpha,alpha)
    rng(seed, 'twister');
    measured = A*sin(2*pi*f*time - phase) + (2*alpha*rand(size(time)) - alpha);
end
