njit-etd2003-081 - New Jersey Institute of Technology

njit-etd2003-081 - New Jersey Institute of Technology njit-etd2003-081 - New Jersey Institute of Technology

archives.njit.edu
from archives.njit.edu More from this publisher
20.01.2015 Views

261 6. The program will run and output the appropriate time-frequency plots and the activity plots (LF, HF and LF/HF: non-normalized and normalized) of the channel that was analyzed using the selected time-frequency distribution.

262 This program provides the STFT and the bilinear time-frequency analysis (Wigner, Smoothed Wigner, Choi-William, Born-Jordan-Cohen and Rihaczek-Margenau) with 3D power spectra versus time and frequency, contour plots and activity plots using the selected time-frequency distribution as outputs. % This program is a modified version of the original work of Adib Mansour (1994) % by Rindala Saliba (1996) % Modified by Douglas Newandee 7/14/2001 for better 3D Graphs and to work with % MATLAB 6.0 and later. clear clc choose=menu('SELECT YOUR DISTRIBUTION', 'Short Time Fourier Transform',... `Wigner','Smoothed Wigner',... `Choi Williams',... `Born_Jordan_Cohen (RID)','Rihaczek_Margenau'); clc disp(") disp([blanks(30) 'WELLCOME TO']) disp([blanks(21) 'TIME FREQUENCY DISTRIBUTION']) disp([blanks(33) 'WORLD']) disp(") disp(") SIGNAL=input('Please enter name of datafile with no extension --> ','s'); filename=[SIGNAL,'.asc']; eval(['load ' filename]); original_rawdata=eval(SIGNAL)';clear SIGNAL; Question 1=input('Which column do you want to analyze [1 2 3] --> ','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,:); end Signa=rawdata; Lwin=input('Please enter the length of analysis window --> '); fftlen=input('Please enter the length of FFT analysis --> '); skip=input('Please enter the number of the skip points --> '); sample=input('Please enter the sample rate of the data --> '); Question_2=input('Do you want to detrend the input signal [y/n] --> ','s');

261<br />

6. The program will run and output the appropriate time-frequency plots and the<br />

activity plots (LF, HF and LF/HF: non-normalized and normalized) <strong>of</strong> the channel that<br />

was analyzed using the selected time-frequency distribution.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!