작성일:
 
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