njit-etd2003-081 - New Jersey Institute of Technology

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');