njit-etd2003-081 - New Jersey Institute of Technology

265
end
waitbar(n/(length(signal)-4*hlf-1))
end
kernel(fftlen-hlf+1:fftlen,:)=flipud(conj(kernel(2:hlf+1,:)));
TFDs=2/fftlen*abs(fft(kernel,fftlen));
close(h)
else
distribution=['Wigner'];
for n=0:skip:skip*(Timeslice-1);
V 1=(signal(I+n:I+hlf+n). *s ignal_conjugate(I+n:- 1 :I-hlf+n));
kernel(1 :11lf+ 1,(n/skip)+ 1 )—(V 1 +(sqrt(- 1 )*V 1))';
kernel(fftlen-
hlf+1:fftlen,(n/skip)+1)=flipud((conj(V1(2:hlf+1))+[imag(V1(2:hlf+1))+(sqrt(-
1)*real(V1(2:hlf+1)))]));
waitbar(n/(length(signal)-4*hlf-1))
end
TFDs=2/fftlen*abs(fft(flipud(kernel)));
close(h)
end
% (Algorithm to calculate instantaneous frequency)
if Question_3=='y';
[r,C]=size(TFDs);
for i=1:C,
W=TFDs(:,i);
Y=(1:r)';
M=W. *Y;
S=sum(M);
F=sum(W);
E(i)=S/F;
end
elseif Question_3=='n';
end
%end % correct mismatched end (DAN)
(Algorithm to extract certain frequency)
if Question_6=='y';
symvag=sum(TFDs(LFC,1:Timeslice));
vagal=sum(TFDs(HFC, 1 : Time slice));
symtopar=symvag./vagal;
end
elseif Question_6='n';
end
(Three Dimensional Graphics)