DETERMINATION of the APPROPRIATE CUTOFF FREQUENCY in the DIGITAL FILTER DATA SMOOTHING PROCEDURE By
Total Page:16
File Type:pdf, Size:1020Kb
'DETERMINATION OF THE APPROPRIATE CUTOFF FREQUENCY IN THE DIGITAL FILTER DATA SMOOTHING PROCEDURE by BING YU B.S., Peking Institute of Physical Education, 1982 A MASTER'S THESIS Submitted in partial fulfillment of the requirements for the degree MASTER OF SCIENCE Department of Physical Education and Leisure Studies KANSAS STATE UNIVERSITY 1988 Approved by: Major Professo 3# AllSDfl 5327b7 ']cl ACKNOWLEDGEMENTS The author wishes to acknowledge the assistance and support of the entire graduate faculty of Kansas State University's Department of Physical Education and Leisure Studies. Special thanks go to committee members Dr. Stephan Konz and Dr. Kathleen Williams for their unique perspectives and editorial assistance. Most of all, I would like to thank my major professor, Dr. Larry Noble, for his integrity, his enthusiasm for knowledge, and the tremendous amount of time and assistance he has given me over the past two years. 11 DEDICATION This thesis is dedicated to my parents, Dr. Gou-Rei Yu and Ming-Hua Lu, to my wife Wei Li, to her parents, Dr. Ping Li and Dr. Xiu-Zhang Yu, and to all of the other folks of my family and her family for their understanding of my absence when my son Charlse Alan Yu was born. Their constant support and encouragement are deeply appreciated. 111 . CONTENTS ACKNOWLEDGMENTS ii DEDICATION iii LIST OF FIGURES vi LIST OF TABLES ix Chapter 1 INTRODUCTION 1 Statement Of The Problem 2 Definitions 3 2 REVIEW OF RELATED LITERATURE 7 The Nature Of Errors 7 Sources Of Errors 10 Data Smoothing Techniques Used In Sport Biomechanics 15 Finite difference technique 15 Least square polynomial approxination. ... 17 Spline function 18 Fourier series 19 Digital filter process 21 Determination Of Cutoff Frequency For Digital Filter 22 Summary 2 6 3. METHOD 27 Standard Data 27 Raw Data 29 Digital Filter Process 31 IV . Evaluation Of Filtered Results 33 Fourier Fitness 33 Evaluation Of Fourier Fitness Results 36 Harmonic Analysis 37 Statistical Analysis 37 Application of Regression equation 38 4 RESULTS AND DISCUSSION 39 Harmonic Analysis Of Measured Signal 39 Selection Of Variables 41 Dependent variable 41 Independent variable 41 Regression Analysis 61 Application Of Regression Equations 82 5. SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 88 Summary 88 Conclusions 91 Recommendations 91 REFERENCES 93 APPENDIX A. A Computer Program For Digital Filter With Different Methods Of Determination Of Cutoff Frequency 9 5 . LIST OF FIGURES Figure 1. Measured and actual sinusoidal signals 9 2. The first derivative of the measured and actual signals 11 3. The second derivative of the measured and actual signals 12 4. Frequency spectrum of measured signal before and after filtered and filter response 24 5. Flow chart of the procedure 28 6. Harmonic spectrum of measured signal 40 7 A symmetry plot for the best cutoff frequency (Fc) 42 8 A symmetry plot for the square root of the best cutoff frequency (PFc) 43 9 The relationship between the square root of the best cutoff frequency (PFc) and sampling frequency (Fs) 44 10. The relationship between the square root of the best cutoff frequency (PFc) and the square root of sampling frequency (PFs) 46 11. The relationship between the normalized harmonic amplitude before overlap harmonic (NHA1) and overlap harmonic (Ho) .... 51 12. A symmetry plot for the normalized harmonic amlitude before overlap harmonic (NHA1) 52 13. A symmetry plot for the -0.125th power transformation of the normalized harmonic amplitude before overlap harmonic (PNHA1) .... 53 14. The relationship between the 0.12 5th power transformation of the normalized harmonic VI amplitude before overlap harmonic (PNHA1) and the difference between the square root of the best cutoff frequency and the mean of the square root of the best cutoff frequency at the same sampling frequency (DPFc) 54 15. The relationship between the normalized overlap harmonic amplitude (NHAo) and overlap harmonic (Ho) 55 16. A symmetry plot for the normalized overlap harmonic (NHAo) 57 17. A symmetry plot for the -0.125th power transformation of the normalized overlap harmonic amplitude (PNHAo) 58 18. The relationship between the -0.125th power transformation of the normalized overlap harmonic (PNHAo) and the difference between the square root of the best cutoff frequency and the mean of the square root of the best cutoff frequency at the same sampling frequency (DPFc) 59 19. The relationship between mean error and the standard deviation of the five normalized harmonic amplitudes after the overlap harmonic (ME) 60 20. A symmetry plot for the standard deviation of the five normalized harmonic amplitudes after overlap harmonic (ME) 62 21. A symmetry plot for the 0.2 5th power transformation of the standard deviation of the five normalized harmonic amplitudes after overlap harmonic (PME) 63 22. The relationship between the 0.25th power transformation of the standard deviation of the five normalized harmonic amplitudes after overlap harmonic (PME) and the difference between the best cutoff frequency and the mean of the best cutoff frequency at the same sampling frequency (DPFc) 64 23. Residuals vs the actual value of the square root of the best cutoff frequency for the first best model 70 vii 24. Residuals vs the predicted value of the square root of the best cutoff frequency for the first best model 71 25. Residuals vs the square root of sampling frequency for the first best model 72 26. Residuals vs the -0.125th power transformation of normalized harmonic amplitude before the overlap harmonic for the first best model 73 27. The relationship between the normalized harmonic amplitude before the overlap harmonic (NHA1) and mean error 74 28. The relationship between normalized harmonic amplitude before overlap harmonic (NHA1) and the standard deviation of the five normalized harmonic amplitude after the overlap harmonic (ME) 75 29. Residuals vs the actual value of the square root of the best cutoff frequency for the second best model 79 30. Residuals vs the predicted value of the square root of the best cutoff frequency for the second best model 80 31. Residuals vs the square root of sampling frequency for the second best model 81 32. A comparison of analog and filtered angular displacement data for the first best model 8 3 33. A comparison of analog and filtered angular displacement data for the second best model 84 34. A comparison of analog and filtered finite angular acceleration data for the first best model 85 35. A comparison of analog and filtered finite angular acceleration data for the second best model 8 6 vm . LIST OF TABLES Table 1 Sources of errors 13 2 C values and corresponding mean random errors 3 3. Sampling frequencies 30 4. Statistical results for Equation 35 65 5. Statistical results for Equation 36 67 6. Statistical results for Equation 37 68 7. Statistical results for Equation 39 77 IX CHAPTER INTRODUCTION Cinematography is the most widely used technique in sport biomechanics research to provide a digital record and/or an image record of overt human body movements. These film records are analyzed quantitatively to obtain linear and angular coordinate-time data for total body or segmental movements. Typically, the basic coordinate-time function of a motion does not provide sufficient information to fully describe the activity. Therefore, these data have to be further treated mathematically to determine the corresponding velocity and acceleration functions. Knowledge of the velocity and acceleration patterns provides more descriptive information about the movement under investigation. It is therefore critical in most studies that precise estimates of velocity and acceleration be obtained. For precise estimates of velocity and acceleration, various mathematical techniques have been employed to smooth raw data to eliminate or reduce the effect of errors which arise in filming and digitizing procedures. Winter (1974) advocated a second-order recursive digital filter for smoothing film analysis data. This digital filter has been suggested as the best among the smoothing techniques used in biomechanics studies (Li & Yu, 1983; Winter, 1979; Wood, 1982) , and is widely used in biomechanics studies at present. In the digital filter smoothing procedure, there is a set of coefficients which control the "degree" of smoothing or the sharpness of cutoff. This set of coefficients is determined according to the cutoff frequency of the raw data. At the present time, the cutoff frequency is subjectively determined by the researcher. There is no method to objectively determine the cutoff frequency. The determination of cutoff frequency has a great influence on the smoothing quality of digital filter. Therefore, if a precise estimate of velocity and acceleration is expected from the data smoothed by the digital filter, it is essential that the problem of determination of cutoff frequency be solved. Statement Of The Problem The purpose of this study was to develop a method for objective determination of the proper cutoff frequency of the digital filter used in data smoothing procedures in sport biomechanics. A set of coordinate-time data of freely falling movement was used as the standard data in which no error was involved. Different sets of computer-generated random numbers were used as random errors, and mixed into the standard data. These sets of artificial raw data with random errors were smoothed by the digital filter procedure using varying cutoff frequencies. The filtered outputs were compared with the standard data to determine the best 2 cutoff harmonic. Multiple regression procedures were used to analyze the relationship between the best cutoff frequency and the characteristic variables of raw data. It was expected that a multiple regression equation could be established to determine the best cutoff harmonic from the characteristic variables of raw data. The independent variables selected to reflect the characteristics of errors involved in raw data were sampling frequency, normalized harmonic amplitudes of the signal, and corresponding harmonics.