2012-06-09 주파수 스펙트럼에 대해서 여쭤보고 싶습니다
작성자 : 이명준
(2022-10-01)
조회수 : 2083
저는 강남세브란스 병원 신경과에서 강사로 일하고 있는 의사입니다. 이번에 저희는 손가락의 움직임을 자이로 센서로 분석하는 연구를 준비하고 있는데요, 이런 종류의 일을 해본 적이 전혀 없어서 여쭤보게 되었습니다. 아는 것이 없어서, 좀 황당한 질문이라도 이해해주시면 감사하겠습니다.
저희 연구에서는 사람이 일부러 움직이는 동작을 분석하려고 하고 있고, 손떨림은 배제하려고 하고 있습니다.기존 연구에서는 손가락의 자발적인 움직임과 불수의적인 떨림을 주파수 스펙트럼을 통해 구분한 곳이 몇 군데 있었습니다. 저희도 그런 방법을 시도해보려고 합니다.
제가 매트랩이나 이런 종류의 프로그램에 대해서는 아는 것이 없어서, 구글로 검색되는 것 중 그럴 듯한 것을 써보고 있습니다만, 이것이 맞는 내용인지 몰라 여기까지 질문을 적게 되었습니다.
저희 데이터는 2열짜리 표로 되어 있는데 첫 줄은 시간, 둘째 줄은 각속도입니다. 다음은 제가 매트랩에 입력한 코드입니다.
data=load('Lt_FT.txt');
time=data(:,1);
wxf=data(:,2);
z=length(wxf);
t=0: 1/fs : (z-1)/fs;
[P,F]=spectrum(wxf,4096,0,hanning(4096),fs);
plot(F,P);
제가 궁금한 점들은,
1. 대체로 올라온 코드들을 보면 시간을 위에서 t를 정의한 것처럼 따로 정의해주고 있었는데요, 저희 데이터는 시간에 대한 정보가 있습니다. 이 상황에서도 굳이 sampling frequency를 이용해서 시간을 다시 정의해주어야 할지요.
2. [P, F] = spectrum(wxf,4096,0,hanning(4096),fs);
이 부분은 구글에서 검색된 것을 그대로 적은 것입니다. 저는 4096과 0이 무엇을 의미하는지 모르겠는데요, 의미를 알면 저희 처지에 맞는 값을 입력할 수 있을 것 같습니다.