How to remove the AWGN noise from data?

34 views (last 30 days)
Sadiq Akbar
Sadiq Akbar on 29 Sep 2022
Commented: Sadiq Akbar on 4 Oct 2022
Let's say we have a vector data Uo=[20 30 40 50];
If we add AWGN noise to it, it becomes data U i.e.,
U=awgn(Uo,30);
Now we see that U and Uo are different. We want to get back our original data Uo. I tried as:
Uo=U-awgn(30);
But it gives error as:
Error using awgn (line 63)
Not enough input arguments.
What commands should I use to get back my vector Uo?

Answers (2)

Chunru
Chunru on 30 Sep 2022
Uo=[20 30 40 50]; % signal
U=awgn(Uo,30) % add noise to signal
U = 1×4
19.9837 29.9773 40.0475 50.0504
% Uo=U-awgn(30);
% if you don't know the noise (which is U-Uo), you are not able to fully
% recover the signal.
% However, it is possible to "filter" out noise if signal and noise are
% different in certain characteristics.
% For your case, signal is a low-pass signal and lowpass filter can help to
% remove some noise.
Uo = [20:1:100]/30;
U = awgn(Uo, 10);
b = fir1(20, 0.02); % low pass filter
Ufiltered = filtfilt(b, 1, U);
t = 0:length(Uo)-1;
plot(t, Uo, 'r', t, U, 'b', t, Ufiltered, 'k')
legend("Signal", "Signal+noise", "Filtered")
  2 Comments
Chunru
Chunru on 30 Sep 2022
The plot is for u to see how good the filtered result as an approximation. You can always comment it out when you are satisfied the results.
For filter to work. You need sufficient data. If you have only 4 points, you may not be able to do much with the added noise.

Sign in to comment.


Star Strider
Star Strider on 30 Sep 2022
I would simply use a moving average filter, for example movmean, or choose one of the options in smoothdata, to eliminate normally-distributed additive noise.
  15 Comments
Sadiq Akbar
Sadiq Akbar on 4 Oct 2022
Thanks a lot for your kind response dear Star Strider. Can you demonstrate it here as I am not too much expert in Matlab?

Sign in to comment.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!