njit-etd2003-081 - New Jersey Institute of Technology
njit-etd2003-081 - New Jersey Institute of Technology njit-etd2003-081 - New Jersey Institute of Technology
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');
- Page 239 and 240: Figure 5.50 HF coherence of COPD (1
- Page 241 and 242: 212 For better presentation of the
- Page 243 and 244: Figure 5.53 Coherence and partial c
- Page 245 and 246: 216 2. Interpretation of the transf
- Page 247 and 248: 218 covariances of the parameters,
- Page 249 and 250: 220 deviations are interpreted as A
- Page 251 and 252: 222 Figure 5.58 Bode plot of the HR
- Page 253 and 254: 224 In this section a simple model
- Page 255 and 256: 226 The data for all 47 COPD subjec
- Page 257 and 258: 228 Figure 5.60 Normal and COPD cla
- Page 259 and 260: 230 Figure 5.61 Normal and COPD cla
- Page 261 and 262: 232 Figure 5.62 Normal classificati
- Page 263 and 264: 234 5.7 Cluster Analysis The purpos
- Page 265 and 266: 236 Figure 5.64 Severity classifica
- Page 267 and 268: 238 both the normal and COPD subjec
- Page 269 and 270: 240 In summary, COPD subjects had h
- Page 271 and 272: APPENDIX A EXERCISE PHYSIOLOGY A.1
- Page 273 and 274: 244 A.3 Figure Out Your Target Hear
- Page 275 and 276: APPENDIX B ANALYSIS PROGRAM LISTING
- Page 277 and 278: 248 4) Click on file, close to exit
- Page 279 and 280: 250 • TN 11
- Page 281 and 282: 252 B.1.2 Partial Coherence Between
- Page 283 and 284: 254
- Page 285 and 286: 256 Block Diagram !rime of record K
- Page 287 and 288: 258
- Page 289: 260 B.2.2 Time — Frequency Analys
- Page 293 and 294: 264 G(:j+1)=G(:,j+1)/(2*sum(G(:j+1)
- Page 295 and 296: 266 T=(length(Signa)/sample)/(Times
- Page 297 and 298: 268 subplot(3, 1,3), plot(T,E); xla
- Page 299 and 300: 270 4. The program creates five out
- Page 301 and 302: 272 B.2.3.4 Program to Generate Sym
- Page 303 and 304: 274 ylabel('frequency'); title('Ins
- Page 305 and 306: 276 The program will run and output
- Page 307 and 308: 278 axis([0 1 0 2]); grid on; xlabe
- Page 309 and 310: 280 vagal=sum(TFDs(HFC,1:k)); symto
- Page 311 and 312: 282 plot(J,symtopar); %plot(A,symto
- Page 313 and 314: 284 4. Remove the constant levels a
- Page 315 and 316: 286 Make sure the agreement is quit
- Page 317 and 318: 288 B.2.6 Principal Components Anal
- Page 319 and 320: 290 Columns 12 through 15 'LF_pcoh_
- Page 321 and 322: 292 I= 1.0000 -0.0000 -0.0000 -0.00
- Page 323 and 324: 294 variances = 3.4083 1.2140 1.141
- Page 325 and 326: 296 B.2.7 Cluster Analysis Program
- Page 327 and 328: 298 end [R,C]=size(Data); if length
- Page 329 and 330: 300 B.2.8 Cross-correlation Program
- Page 331 and 332: 302 C.3 Partial coherence of HR and
- Page 333 and 334: 304 [13] Madwed, J., and R. Cohen.
- Page 335 and 336: 306 [41] Mallat, S. G., "A Theory f
- Page 337: [70] Tazebay, M.V., R.T. Saliba and
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.