임시 블로그 이름

MATLAB의 연산 결과를 자동으로 저장해 보자. 매우 오래 걸리는 연산을 하거나, 매번 수행하는 연산 결과를 알아서 자동으로 따로따로 저장하게 만들고 싶을때 유용하다. 날짜별로 폴더를 만들고, 그 밑에 파일 이름을 시간으로 해서 저장을 해보자. 예시:clear all;close all; s = 0; for j= 1: 1000 s = s + j; % 1에서 1000까지 합을 계산하는 연산end DIRNAME = datestr(now,'yyyymmdd');FILENAME = sprintf('./%s/SAVE_%s.mat',DIRNAME,datestr(now,'HHMMSS')); mkdir(DIRNAME); save(FILENAME); datestr 함수는 날짜와 시간을 문자열로 반환해 주는 함수고, now..

(이미지는 위키피디아 optical flow 설명 페이지에서 퍼왔다. http://en.wikipedia.org/wiki/File:Opticfloweg.png) 영상에서 Motion Estimation을 하는 방법 중 아마 가장 유명한것이 Block-Matching-Algorithm (BMA) 와 이 Optical Flow가 아닐까 한다. BMA는 뭐 그냥 블록을 나누고 두 영상에서 가장 비슷한 블록을 찾아내는 거니까 이해하는데 별 큰 어려움이 없고, 근데 이 Optical Flow는 좀 아리송 하다. 설명을 해야하는데, 귀찮기도 하고 위키피디아에 설명이 잘 되어있으니까 그걸 가지고 하겠다. 아래는 그냥 위키피디아에 있는걸 요점 정리만 해서 쓴거다. 원본은 여기 있다. http://en.wikipedia..

Matlab에서 주파수 분석을 할때 일일이 fft를 하고, 다시 fftshift를 하고, 거기다 abs를 취해서 plot을 하기가 귀찮아서 만들어 보았다. 사용법은 간단하다. >> freqshow(signal) >> freqshow(signal, Fs, 'option', 'scale') signal은 주파수 분석을 할 신호이고 Fs는 sampling frequency, 'option'은, 'amp'와 'phase', 'both' 중에 하나를 선택하면 되는데, 각각 amplitude만 보여줄건지, phase만 보여줄건지, 둘 다 보여줄건지 선택하는 것이다. 기본값은 'amp'이다. 'scale'은 'normal'과 'log'가 있다. amplitude spectrum의 스케일을 조정하는 것이다. 첫 줄 처..

블로그 유입 키워드를 보면 matlab lowpass filter가 꽤 있다. 아마도 컴퓨터나 전기전자 공학 학생들인것 같고, 처음에 이런 필터를 어떻게 설계해야 하는지 찾다 찾다 여기 까지 왔으리라... 사실 이런 필터를 설계하는 방법은 여러가지가 있고, 직접 필터 설계하는 방법 및 공식을 코드로 작성할 수 도 있다. 그렇지만 매번 이렇게 코드를 작성하는 것은 너무 귀찮은 일이고 magnitude response, phase response, pole-zero plot등을 보기 위해서는 또 코드를 작성해야 할 것이다., 매틀랩에서는 다행히도 좋은 툴을 제공하고 있다. 커맨드 라인에서 fdatool을 입력해 보자. (단 signal processing toolbox가 설치되어 있어야 한다.) 왼쪽 아래 ..