function [wrap,t]=my_wrap(x,N)
% Figure the wrap of the input signal
% s, the input signal which type is float
% Fs, the sampling rate
% N, one data out of N data

% Author: Crystal Shi
% ShenZhen HYT Science&Technology CO.,LTD.
% Date: 2007.05.31

% If the input signal is fftpower, convert it into dB first
% dBs = 20*log10(s);
% [val,pos] = max(dBs);
% x = dBs - val;
x = x';
L = length(x); % the length of the input signal 
len = L / N; % each part contains len datas 

% store the input signal in parts
data = zeros(N,len); % the matrix that store the input signal in parts 
max_data=zeros(1,len);
for i=1:1:len
    data(:,i)=x(1,i*N-N+1:i*N);
%     dt=data(:,i)';
    max_data(1,i) = data((find(data(:,i)==max(data(:,i)))),i);
    % Find out the maxminum values in each column
     
end

wrap=max_data;

int = Fs/(length(max_data)-1);
t = -Fs/2:int:Fs/2;
% % figure(4)
plot(t/1e3,max_data,'r')
% grid on