function symmetric_noise = create_symmetric_noise(signal,width) m=mean(signal); distance_from_mean=(signal-m).^2; max_distance=max(distance_from_mean); noise_spread=signal*NaN; for k=1:length(signal) noise_spread(k)=randn*width*exp(-distance_from_mean(k)/(3*max_distance)) end symmetric_noise=noise_spread;