Computational Astronomy 16 12 Fitting interpolation and Signal

Computational Astronomy 제 16 장 자료 분석 1/2 (Fitting, interpolation and Signal processing) 전산천문학 2 2007. 가을

Mission Computational Astronomy • Fitting Ø linear regression Ø curve fitting • Interpolation Ø 1 D interpolation Ø Bilinear interpoloation • 시간에 따라 나오는 신호 혹은 물리량을 signal processing 을 통하여 분석할 수 있다.

무작위수 발생시키기 • 무작위수(random number) Ø randomu : 균일분포 Ø randomn : 정규분포 Ø x=randomu(seed, 30) • seed : 무작위수의 시작점 • 시작점이 같으면 같은 분포의 무 작위수를 지정함 • • • Seed 1=1 L X=randomu(seed 1, 32) Print, seed 1 Y= randomu(seed 1, 32) Plot, x, y Computational Astronomy

Linear Regression • X 와 Y가 선형 관계에 있으면, Ø Yi = a + b * xi + εi Ø 이 때, a 와 b 는 선형 회귀 방법 (linear regression)으로 구할 수 있 으며, 상관계수 (Correlation coefficient) r 도 구할 수 있다. Ø 기울기 = regress(xarr, yarr, const=a, correation=r, [/double]) Ø 상관계수 r 구하기 • correlate(xarr, yarr, /double) • LIN_REGRESS. PRO • 선형 관계 외에도 곡선 및 다항식으로 fitting해주는 함수도 있다. Ø curvefit, svdfit 등 참고 Computational Astronomy

Computational Astronomy Interpolation (내삽) • 유한한 개수의 목록화 된 자료들의 연속적인 값을 구하고자 할 때 • 내삽 (Interpolation) Ø 주어진 값들의 사이에 있는 값을 구하기 • 외삽 (extrapolation) Ø 주어진 범위 바깥의 값 을 구하기 • 1 -D interpolation • INTERPOL function Ø Ø Linear (no keyword) Linear square fit quadratic 3 pt average quadratic spline

Interpolation • 2 -D interpolation Ø bilinear (2 darr, xarr, yarr) Ø xarr, yarr는 내삽할 자료 의 인덱스 값이며, 정수 로 주어진 인덱스 값의 사이값인 실수값이 들어 간다. Computational Astronomy

삼각점 내삽 • TRIGRID 함수 Ø Ø KEYWORD NX, NY /EXTRAPOLATE / QUINTIC Computational Astronomy