njit-etd2003-081 - New Jersey Institute of Technology
njit-etd2003-081 - New Jersey Institute of Technology njit-etd2003-081 - New Jersey Institute of Technology
281 %* orient landscape; %* print -dps wvletf.eps; subplot(2,1,2), contour(t, freq, abs(coef),20); axis([0 100 0 2]); grid on; xlabel('Time (sec)'); ylabel('Frequency (Hz)'); title(['Contour of ',eval('filename'),' using ',eval('wave'),' wavelet']); % print -dps wvletf.eps; figure; subplot(2,1,1); %plot(A,x); plot(t,x); title(['Detrended IIBI Plot of ',eval('filename'),' using ',eval('wave'),' wavelet']); xlabel('Time (sec)'); ylabel('Amplitude'); subplot(2,1,2); plot(t,E); %plot(A,E); axis([0 50 0 2]); xlabel('Time (sec)'); ylabel('Frequency (Hz)'); title(['Instantaneous frequency Plot of',eval('filename')]); figure; subplot(3,1,1); plot(J,symvag); %plot(A,symvag); %gtext(top); title(['Symp. and Parasymp Mixture (LF) Plot of ',eval('filename'),' using ',eval('wave'),' wavelet']); %xlabel('time (sec)'); ylabel('Amplitude'); subplot(3,1,2); plot(J,vagal); %plot(A,vagal); title('Parasympathetic range (HF)'); %xlabel('time (sec)'); ylabel('Amplitude'); subplot(3,1,3)
282 plot(J,symtopar); %plot(A,symtopar); %gtext(top); title('Ratio of LF to HF) xlabel('Time (sec)'); ylabel('Amplitude'); %order=input('Please enter the order of the lowpass filter. '); %freq=input('Please enter the cuttoff frequency for LPF. '); %sample=input('Please enter the sample rate of the data. '); %order=12; %freq=0.03; %nfreq=freq/sf; %[poles,zeros]=butter(order,nfreq); %dtrnd_nsymvag=filtfilt(poles,zeros,n_symvag); %dtrnd_nvagal=filtfilt(poles,zeros,n_vagal); %dtrnd_nsymtopar=filtfilt(poles,zeros,n_symtopar); figure subplot(3,1,1); plot(J,n_symvag); %plot(A,n_symvag); %gtext(top); title(['Normalized Symp and Parasymp Mixture (NLH) Plot of ',eval('filename'),' using ',eval('wave'),' wavelet']); %xlabel('time (sec)'); ylabel('Amplitude'); subplot(3,1,2); plot(J,n_vagal); %plot(A,n_vagal); title('Normalized Parasympathetic range (NHF)'); %xlabel('time (sec)'); ylabel('Amplitude'); %figure subplot(3 13) plot(J,n_symtopar); %plot(A,symtopar); %gtext(top); title('Ratio of NLF to NHF') xlabel('Time (sec)'); ylabel('Amplitude');
- 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 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: 280 vagal=sum(TFDs(HFC,1:k)); 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
281<br />
%* orient landscape;<br />
%* print -dps wvletf.eps;<br />
subplot(2,1,2), contour(t, freq, abs(coef),20);<br />
axis([0 100 0 2]);<br />
grid on;<br />
xlabel('Time (sec)');<br />
ylabel('Frequency (Hz)');<br />
title(['Contour <strong>of</strong> ',eval('filename'),' using ',eval('wave'),' wavelet']);<br />
% print -dps wvletf.eps;<br />
figure;<br />
subplot(2,1,1);<br />
%plot(A,x);<br />
plot(t,x);<br />
title(['Detrended IIBI Plot <strong>of</strong> ',eval('filename'),' using ',eval('wave'),' wavelet']);<br />
xlabel('Time (sec)');<br />
ylabel('Amplitude');<br />
subplot(2,1,2);<br />
plot(t,E);<br />
%plot(A,E);<br />
axis([0 50 0 2]);<br />
xlabel('Time (sec)');<br />
ylabel('Frequency (Hz)');<br />
title(['Instantaneous frequency Plot <strong>of</strong>',eval('filename')]);<br />
figure;<br />
subplot(3,1,1);<br />
plot(J,symvag);<br />
%plot(A,symvag);<br />
%gtext(top);<br />
title(['Symp. and Parasymp Mixture (LF) Plot <strong>of</strong> ',eval('filename'),' using ',eval('wave'),'<br />
wavelet']);<br />
%xlabel('time (sec)');<br />
ylabel('Amplitude');<br />
subplot(3,1,2);<br />
plot(J,vagal);<br />
%plot(A,vagal);<br />
title('Parasympathetic range (HF)');<br />
%xlabel('time (sec)');<br />
ylabel('Amplitude');<br />
subplot(3,1,3)
Change language