>질문에 깊은 감사를 드립니다. 재미있게 잘 읽었습니다. 앞으로도 재미있는 것들로 얘길 나눌 수 있다면 좋겠습니다.

제가 궁금한 점들은,
1. 대체로 올라온 코드들을 보면 시간을 위에서 t를 정의한 것처럼 따로 정의해주고 있었는데요, 저희 데이터는 시간에 대한 정보가 있습니다. 이 상황에서도 굳이 sampling frequency를 이용해서 시간을 다시 정의해주어야 할지요.

> 계측된 자료의 시간 간격을 dt로 하면, 샘플링 주파수 fs = 1/dt를 구할 수 있습니다.
> dt = time(2,1) - time(1,1);
> fs = 1/dt;


2. [P, F] = spectrum(wxf,4096,0,hanning(4096),fs);
이 부분은 구글에서 검색된 것을 그대로 적은 것입니다. 저는 4096과 0이 무엇을 의미하는지 모르겠는데요, 의미를 알면 저희 처지에 맞는 값을 입력할 수 있을 것 같습니다.


> spectrum은 Power Spectral Density Function을 생성하는 매트랩 함수입니다. 요즘 버전에서는 spectrum대신에 spectrum.welch를 사용합니다. 4096은 FFT(fast Fourier transform)을 사용할 때 전체 데이터를 4096개씩(즉, 2^12개씩) 잘라서 사용하는 것을 의미하고, 0은 잘라 쓴 데이터 사이에 overlapping(중복)되는 데이터의 개수가 0개임을 뜻합니다.

> 혹시라도 자세한 내용이 필요할까 싶어, 구조동역학(김두기, 2009)의 1-4.4절 푸리에변환과 3-4.2절 신호처리를 첨부하였습니다.
> 재미있고 보람 있는 결과가 있다면 좋겠습니다. 주말이 겹쳐 답변이 늦었습니다.

================== 원본글 내용 ==================

저는 강남세브란스 병원 신경과에서 강사로 일하고 있는 의사입니다. 이번에 저희는 손가락의 움직임을 자이로 센서로 분석하는 연구를 준비하고 있는데요, 이런 종류의 일을 해본 적이 전혀 없어서 여쭤보게 되었습니다. 아는 것이 없어서, 좀 황당한 질문이라도 이해해주시면 감사하겠습니다.

저희 연구에서는 사람이 일부러 움직이는 동작을 분석하려고 하고 있고, 손떨림은 배제하려고 하고 있습니다.기존 연구에서는 손가락의 자발적인 움직임과 불수의적인 떨림을 주파수 스펙트럼을 통해 구분한 곳이 몇 군데 있었습니다. 저희도 그런 방법을 시도해보려고 합니다.

제가 매트랩이나 이런 종류의 프로그램에 대해서는 아는 것이 없어서, 구글로 검색되는 것 중 그럴 듯한 것을 써보고 있습니다만, 이것이 맞는 내용인지 몰라 여기까지 질문을 적게 되었습니다.

저희 데이터는 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이 무엇을 의미하는지 모르겠는데요, 의미를 알면 저희 처지에 맞는 값을 입력할 수 있을 것 같습니다.