Hi Pako,
As per my understanding, you would like to calculate signal power, noise power, and signal-to-noise ratio(SNR) of the AWGN waveform within the specified bandwidth.
I would suggest you to design a 'bandpass filter' with the desired bandwidth, sampling frequency, and filter order using the MATLAB function 'designfilt'.
You can filter the waveform from the filter created before using the MATLAB function 'filter'. Singal power and noise power can be determined by calculating the 'rms' value of the respective waveforms using the MATLAB function 'rms'.
SNR can be calculated from the below mentioned formula by calculating the signal and noise powers:
snr = 10 * log10(signal_power / noise_power);
Here is an example showing how you can do this:
% Generating the waveforms. This is a sample waveform.
t = (0:0.1:60)';
waveform = sawtooth(t);
waveform_with_noise = awgn(waveform, 0.8,'measured'); %adding awgn to the waveform
% Design of bandpass filter
fs = 1000; % Sampling Frequency
fpass = [1 400]; % Lower and upper cutoff frequencies
filterOrder = 10; % Filter Order
bpFilter = designfilt('bandpassfir', 'FilterOrder', filterOrder, 'CutoffFrequency1', fpass(1), 'CutoffFrequency2', fpass(2), 'SampleRate', fs);
% Filtering waveforms with bandpass filter
filteredwaveform = filter(bpFilter , waveform);
filteredwaveform_with_noise = filter(bpFilter , waveform_with_noise);
% Calculating Signal Powers
signal_power = rms(filteredwaveform)
noise_power = rms(filteredwaveform_with_noise - filteredwaveform)
% Determine the SNR from signal and noise powers
snr = 10 * log10(signal_power / noise_power)
Please refer to the following documentation links for more information on:
- 'designfilt' function: https://www.mathworks.com/help/signal/ref/designfilt.html
- 'filter' function: https://www.mathworks.com/help/matlab/ref/filter.html
- 'rms' function : https://www.mathworks.com/help/matlab/ref/rms.html
