njit-etd2003-081 - New Jersey Institute of Technology

276
The program will run and output the 9 appropriate time-frequency and activity plots
of the channel that was analyzed.
The MATLAB program given below calculates the wavelet coefficients and plots
them in 3D as power intensity versus time and frequency as well as the regular and
normalized activity plots.
% d_wvlet.m
% Matlab Wavelet Toolbox 2.0.
% Modified by D. Newandee
% 7-13-2001: using time - freq spectrogram and change to run with MATLAB 6.0
% 10-23-2001: add calculation of normalized sympathovagal ratio using chosen
% wavelet.
% This program will accept any wavelet from the Matlab Wavelet Toolbox 2.0 by
% entering the wavelet name in " for the 'wave' variable below.
clear all;
close all;
SIGNAL=input('Please enter name of datafile with no extension --> ','s');
filename=[SIGNAL,'.asc'];
eval(['load ' filename]);
original_rawdata=eval(SIGNAL)';clear SIGNAL;
Question_l=input('Which channel do you want to analyze [1=02, 2=BP, 3=resp,
4=ECG] --> ','s');
if Question_1=='1'
rawdata=original_rawdata(1,:);
elseif Question_1=='2'
rawdata=original_rawdata(2,:);
elseif Question_1=='3'
rawdata=original_rawdata(3,:);
elseif Question 1=='4'
rawdata=original_rawdata(4,:);
end
sf=input('Please enter the sampling rate of the data --> ');
T=input('Please enter the window size in seconds --> ');
nfft=input('Please enter the fft size --> '); % must be larger than window size and in ^^2
wave = input('Please enter name of wavelet used --> Vs');