njit-etd2003-081 - New Jersey Institute of Technology
njit-etd2003-081 - New Jersey Institute of Technology njit-etd2003-081 - New Jersey Institute of Technology
271 B.2.3.2 Frequency to Indices Conversions The equation for converting the frequency ranges to the corresponding indices used in the analysis program is shown below: Freq (Hz) Index Freq (Hz) _ Index Freq (Hz) Index 0.039 1 0.351 9 0.663 17 0.078 2 0.390 10 0.702 18 0.117 3 0.429 11 0.741 19 0.156 4 _ 0.468 12 0.780 20 0.195 5 _ 0.507 13 0.819 21 0.234 6 0.546 14 0.858 22 0.273 7 0.585 15 0.897 23 0.312 8 0.624 16 0.936 24 Frequency LF Range: HF Range: Indices 0.04-0.10 Hz 0.04-0.15 Hz 0.15-0.4 Hz 0.15-0.8 Hz 1:3 1:4 4:10 4:20 B.2.3.3 MATLAB Time-Frequency Analysis Matrix Row represents Frequency Column represents Time
272 B.2.3.4 Program to Generate SympathoVagal Activity Plots Using Wigner Distribution This .m MATLAB program calculates the vagal tone and sympatho-vagal ratio. It is an implementation of time-frequency analysis using the Wigner distribution. The MATLAB code for the program is given below: function symparl(rawdata,top) % SYMPAR1(iibi,gtitle) % iibi is the interpolated interbeat interval calculated % by pslwsu, either in Matlab or in S-plus. Gtitle is the % title of the output graphs. % Written by Sanjay Fernando. % last revision: 7/16/01 by Douglas Newandee (MATLAB 6.) rawdata=rawdata(:); order=5; freq=0.03; sample=20; nfreq=freq/sample; dtrendata=rawdata; [row,col]=size(dtrendata); I=1:row; I=I(:); A=(I/sample)/60; % This part is for testing Janse Kaiser Wigner calculation % algorithm with no window fs=1000; m=512; % The size of the fft we will be computing. skip=25; % Number of points we skip to get the next segment. p=60; % The number of freq vals we will be plotting k=fix((row-m)/skip); % the number of spectra we compute w=ones(size(1:m)); % window specification. Can be changed. w=w(:); x=hilbert(dtrendata); % Forms the analytic function of x L=m/2; l=-(L-1):(L-1); n=L;
- 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 and 290: 260 B.2.2 Time — Frequency Analys
- Page 291 and 292: 262 This program provides the STFT
- 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: 270 4. The program creates five out
- 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
271<br />
B.2.3.2 Frequency to Indices Conversions The equation for converting the<br />
frequency ranges to the corresponding indices used in the analysis program is shown<br />
below:<br />
Freq (Hz) Index Freq (Hz) _ Index Freq (Hz) Index<br />
0.039 1 0.351 9 0.663 17<br />
0.078 2 0.390 10 0.702 18<br />
0.117 3 0.429 11 0.741 19<br />
0.156 4 _ 0.468 12 0.780 20<br />
0.195 5 _ 0.507<br />
13 0.819 21<br />
0.234 6 0.546 14 0.858 22<br />
0.273 7 0.585 15 0.897 23<br />
0.312 8 0.624 16 0.936 24<br />
Frequency<br />
LF Range:<br />
HF Range:<br />
Indices<br />
0.04-0.10 Hz<br />
0.04-0.15 Hz<br />
0.15-0.4 Hz<br />
0.15-0.8 Hz<br />
1:3<br />
1:4<br />
4:10<br />
4:20<br />
B.2.3.3 MATLAB Time-Frequency Analysis Matrix<br />
Row represents Frequency<br />
Column represents Time