제목: 가우스 함수 [] | 작성일: 2010. 5. 7. 17:05 |
주로 통계학, 천문학등지에서 자주 사용됨
ae^-{(x-b)^2/2c^2)} - 표준정규분포도
항상 값은 1 - 그 까닭은 확율의 총 합은 1이기때문
가우스 함수 그래프
; x-y gaussian function
; a= Graph's maxium value
; b= Maxium value of x condition
; c= Graph's width
; exp => e=2.71828182.....
; by wiinemo.kor@gmail.com
x=16*findgen(101)/100-8
a=1 & b=1 & c=1
y=a*exp(-1*((x-b)^2/2*c^2))
plot, x, y, title='X - Y Gaussian Function', color='ffffff'x
a=0.5 & b=0 & c=3
y=a*exp(-1*((x-b)^2/2*c^2))
oplot, x, y, color='ff00ff'x, linestyle=1
a=1.1 & b=3 & c=2
y=a*exp(-1*((x-b)^2/2*c^2))
oplot, x, y, color='343fad'x, linestyle=2
end
2차원 가우스 함수
pro gauss_2d
gaussz=fltarr(101,101)
;fltarry는 2차원 100,100선언 (값은 없음)
;만약 fltarry 선언을 하지 않으면
;아래의 gaussz[i,j]에서 에러가 남
a=1.
b1=2.
b2=3.
c1=0.8
c2=0.5
x=16*findgen(101)/100-8
y=16*findgen(101)/100-8
for i=0,100 do begin
for j=0,100 do begin
gaussz[i,j]=a*exp(-1*((x[i]-b1)^2/(2*c1)+(y[j]-b2)^2/(2*c2^2)))
endfor
endfor
device, decomposed=0
loadct, 0
;shade_surf,gaussz,x,y,shade=(bytscl(gaussz)),background=128
shade_surf, gaussz,x,y
;surface,gaussz,x,y
;save,gaussz,x,y,filename='gauss_2d.sav'
;stop
end