njit-etd2003-081 - New Jersey Institute of Technology

273
Z=zeros([p,k]);
for i=1:k
g=x(n+l).*conj(x(n-l));
g(2 *L)=0;
y=w.*g; % Apply window to g or kernel.
end
Y=2/m*abs(fft(y,m)); % evidently because it's analytic
% we only need 2/N
Z(:,i)=Y(1 :p);
n=n+skip;
LFC=1:4;
symvag=sum(Z(LFC,1:k));
HFC=4:20;
vagal=sum(Z(HFC,1:k));
symtopar=symvag./vagal;
% This section is for determining the instantaneous frequency
% from the Wigner distribution.
[r,c]=size(Z);
for i=1:c;
W=Z(:,i);
Y=(1:r)';
M=W. *Y;
S=sum(M);
F=sum(W);
E(i)=S/F;
end
%Determines the size of the matrix created by wgjka4.m.
%Repetition based on number of columns.
%W is assigned the values of the column.
%Y represents the frequencies.
%Sum of the Wigner Distribution values.
% Plotting commands
J=[ 12.8+(1:k)*
subplot(211);
plot(A,dtrendata);
title('IBI detrended');
xlabel('time');
ylabel('amplitude');
subplot(212);
plot(J,E);
xlabel('time');