%fft for X signal
fftx=fft(X,NFFT);
% absolute value of the signal
Xn = abs(fftx);
A = mean(Xn);  % I took the mean to get single value.
phase = unwrap(angle(fftx));
phase_mean = mean(alpha);

[b,a] = butter(2, [ 2440/(fs/2),2441/(fs/2)],'bandpass');
X_filtered = filter(b,a,X);
% Applying FFT to filtered signal
fft_filtered = fft( X_filtered,NFFT);
Xn_filtered = abs(fft_filtered);
A_filtered = mean(Xn_filtered);
phase_filtered = unwrap(angle(fft_filtered));
phase_filtered_mean = mean(phase_filtered);

L = length(X); % 50050 is length of signal
NFFT = 2^nextpow2(length(X)); % Zero padding to nearest N power 2
%Define frequency axis
fs = 1e7; % Sampling frequency
df = fs/NFFT; % frequency resolution
dt = 1/df; % time resolution
X_detrend = detrend(X,0); % Removing DC Offset