게시판 답변이 생성됐습니다
-
글쓴이글
-
inmyta27회원
ㄹㄹ
inmyta27회원답변 너무 감사합니다. 제가 올린 코드에서 데이터는
두가지 축에는 서로 관계가 없기 때문에 말씀해주신것처럼 current 함수를 사용하는게 좋을꺼 같은데
current에 대해 잘 몰라서
axis 축을 따르고 싶은 데이터는 그럼yaxis=axis(…………..) 이후에
d=plot(……………../current) 을 하면 d는 yaxis 축을 따르게 되는 건가요?inmyta27회원감사합니다 ~ 덕분에 해결하였습니다.
inmyta27회원idl 8.3버전을 사용하고 있으며, 위에서 쉽게 설명하기위해 M2, M4,M6값을 지정해서 계산하는 방식으로 설명하였지만 실제로는 아래와 같이 M2, M4, M6가 계산됩니다.
delta_d=0.2 ; [mm]
M_2=fltarr(orb, ncount)
M2=fltarr(orb)
M_4=fltarr(orb, ncount)
M4=fltarr(orb)
M_6=fltarr(orb, ncount)
M6=fltarr(orb)for aa=0, orb-1 do begin
for i=0, ncount-1 do begin
M_2[aa,i]=(d[aa,i]^2)*nd[aa,i]*delta_d
M2[aa]=M2[aa]+M_2[aa,i]M_4[aa,i]=(d[aa,i]^4)*nd[aa,i]*delta_d
M4[aa]=M4[aa]+M_4[aa,i]M_6[aa,i]=(d[aa,i]^6)*nd[aa,i]*delta_d
M6[aa]=M6[aa]+M_6[aa,i]endfor ; for i
endfor ; for aa위에서 계산된 M2,M4,M6 를 가지고 아래 eta, mu 계산에서 문제가 발생하여 알려주신대로 double을 썼을때와
쓰지 않았을때의 값을 print 하였을때;=== M246 moments ===
; eta=M4^2/(M2*M6)
; mu=((7-11*eta)-(eta^2+14*eta+1)^0.5)/(2*(eta-1))
; lamda=[(M2/M4)*(mu+3)*(mu+4)]^0.5
; No=(M2*lamda^(mu+3))/gamma(mu+3)eta=(M4^2)/(M2*M6) => 기존포맷
eta=double(M4)^2/(double(M2)*double(M6)) => 알려주신 포맷두 각각의 포맷을 print 하였을때 둘다 1값을 프린트하며
mu=((7-11*eta)-(eta^2+14*eta+1)^0.5)/(2*(eta-1))
mu 값에서는 두 값이 서로 미세하지만 차이를 보였습니다.
두 값모두 이론적으로 계산하였을때의 값과 차이를 보여 1)알려주신 방법을 사용하는 것이 맞는것인지 아님
또다른 방법이 있는지 해서 다시 질문드립니다.
2) 또한 위에서 구해진 mu를 통해 lamda, No가 계산되어지는데 아래 부분에서도 double을 써서 계산해야하는것이
맞는지 궁금하니다. 알려주신 double을 쓰면 ‘program caused….’와 같은 에러명령어는 사라지고
잘 수행되는데 값에 대한 신뢰성에 대해 판단이 잘 들지 않아서요 ㅠㅠlamda=((M2/M4)*(mu+3)*(mu+4))^0.5
No=(M2*lamda^(mu+3))/gamma(mu+3)inmyta27회원말씀해주신것처럼 cgimage 다음에
plot, e_diameter, hlk, xrange=xrange, yrange=yrange, xstyle=1, ystyle=1, color=’black’,$
position=[0.1, 0.1, 0.87, 0.87],/NOERASE 를 넣었더니 image 위에 plot이 겹쳐집니다.그런데.. 이렇게하면 x,y 축 눈금들도 모두 중첩되서 써져서요.. 이건 어떻게 해결이 안되는지 궁금합니다.
또 제가 plot을 하나만 추가가 아니라 몇개 더 추가를 하는데
noerase를 쓸때는 oplot이 적용이 안되서위와 같은 형태로 plot을 색깔과 선스타일만 바꿔서 했더니 모두 다 나타나지만.. 역시나 x, y축 눈금들도
모두 중첩되고 색깔도 하나는 빨강, 주황, 파랑 …등.. 모두 중첩되는데
혹시 이 부분을 해결할 방법을 아는지 궁금합니다.inmyta27회원읽을때 모두 실수형 숫자로 읽었습니다. 다른값들은 각 열마다 잘 읽히는데..
NaN이 있는 KDP 열만 문제가 있엇습니다.
제가 가령
D Kd_5 Kd
1 10.1 -NaN
2 10.2 -NaN
3 10.3 10.4
4 10.5 10.6
5 10.7 10.8위 데이터처럼 두 Kd_5와 kd의 regress나 scatter plot 으로 서로 비교하고싶다면
위에 말씀하신 where을 쓰면 kd는 D가 1,2,3 에 해당하는 값을… Kd_5도 D가 1,2,3에 해당되는
값들만 regress 하고 plot도 되게 되는건가요?inmyta27회원exp(-98)이 사실상 0에 가까운 수인것은 알지만
그럼에 따라 A,B 값이 같게 나와.. 혹시나 저것을 차이를 보고 싶어서 질문드렸습니다inmyta27회원조언해주신대로 수정해 보았더니 값이 잘 산출됩니다.
너무 감사드려요~~ -
글쓴이글