2014-12-17 2 views
0

gnuplot에서 동일한 데이터에 대해 2 법칙을 만들 수 있는지 알고 싶습니다. 내 그래픽에는 2 개의 지역이 있으며 두 가지 모두의 법칙이 필요합니다. 나는 xmgrace를 만드는 방법을 알고 있지만 xmgrace의 최종 예술은 좋아하지 않으며,이 그림의 엡플 래이 트를 만들고 싶습니다. 모든 시간 나는이gnuplot의 멱 법칙을 적용한 그래프의 두 영역

https://www.dropbox.com/s/hdyjm1c36evi1ph/Global.eps?dl=0

EDIT처럼 xmgrace에 내가 그 지역에서 전력이 낮은 시작을 만들 수 있습니다, 두 지역에서, 제로의 전원 법에게 그것의 시작을하려고 : 여기 모든 데이터를 게시 지금.

이것은 내 데이터입니다.

0.000000 0.000180 
0.002000 0.000180 
0.004000 0.000180 
0.006000 0.000360 
0.008000 0.000240 
0.010000 0.000360 
0.012000 0.000420 
0.014000 0.000480 
0.016000 0.000540 
0.018000 0.000720 
0.020000 0.000660 
0.022000 0.000840 
0.024000 0.000960 
0.026000 0.001080 
0.028000 0.001200 
0.030000 0.001320 
0.032000 0.001320 
0.034000 0.001500 
0.036000 0.001560 
0.038000 0.001740 
0.040000 0.001800 
0.042000 0.002040 
0.044000 0.002040 
0.046000 0.002400 
0.048000 0.002460 
0.050000 0.002340 
0.052000 0.002640 
0.054000 0.002700 
0.056000 0.002760 
0.058000 0.002940 
0.060000 0.002760 
0.062000 0.002940 
0.064000 0.003180 
0.066000 0.003000 
0.068000 0.003180 
0.070000 0.003120 
0.072000 0.003120 
0.074000 0.003060 
0.076000 0.003120 
0.078000 0.003180 
0.080000 0.003300 
0.082000 0.003120 
0.084000 0.003180 
0.086000 0.003300 
0.088000 0.003300 
0.090000 0.003300 
0.092000 0.003420 
0.094000 0.003300 
0.096000 0.003480 
0.098000 0.003420 
0.100000 0.003540 
0.102000 0.003360 
0.104000 0.003360 
0.106000 0.003540 
0.108000 0.003480 
0.110000 0.003540 
0.112000 0.003360 
0.114000 0.003540 
0.116000 0.003420 
0.118000 0.003480 
0.120000 0.003420 
0.122000 0.003600 
0.124000 0.003480 
0.126000 0.003480 
0.128000 0.003900 
0.130000 0.003540 
0.132000 0.003900 
0.134000 0.003600 
0.136000 0.003720 
0.138000 0.003780 
0.140000 0.003660 
0.142000 0.003720 
0.144000 0.003780 
0.146000 0.003660 
0.148000 0.003720 
0.150000 0.003900 
0.152000 0.003780 
0.154000 0.003840 
0.156000 0.003960 
0.158000 0.003900 
0.160000 0.003900 
0.162000 0.003840 
0.164000 0.003960 
0.166000 0.004020 
0.168000 0.003960 
0.170000 0.004140 
0.172000 0.004020 
0.174000 0.004020 
0.176000 0.004140 
0.178000 0.004080 
0.180000 0.004020 
0.182000 0.004020 
0.184000 0.004200 
0.186000 0.004200 
0.188000 0.004200 
0.190000 0.004260 
0.192000 0.004320 
0.194000 0.004080 
0.196000 0.004200 
0.198000 0.004140 

gnuplot에서 2 데이터의 데이터를 분리합니다. 각 지역을 시작하고

set term postscript eps enhanced color 
set output "powerlaw.eps" 
set ylabel "{/Symbol D}{/Symbol w}" 
set xlabel "d" 
set logscale x 
set logscale y 
f(x) = a*x**b 
g(x) = c*x**d 
h(x) = x < 0.057 ? f(x) : g(x) 
fit h(x) 'data.dat' u 1:2 via a,b,c,d 
pl "data.dat" w p tittle "", h(x) 

내가 A, B, C에 대해 가지고 그들을, D

a    = 0.0898041  +/- 0.3695  (411.4%) 
b    = 1.20131   +/- 1.299  (108.2%) 
c    = -0.626148  +/- 2.368e+19 (3.782e+21%) 
d    = 28.5183   +/- 2.372e+19 (8.316e+19%) 

편집 할 톱 : 이유를 알고 2015년 10월 1일 내가 해결책을 찾았을하지만 해달라고, 나는이 함께

set term postscript eps enhanced color 
set output "powerlaw.eps" 
set ylabel "{/Symbol D}{/Symbol w}" 
set xlabel "d" 
set logscale x 
set logscale y 
f(x) = a*x**b+e 
g(x) = c*x**d+f 
h(x) = x < 0.057 ? f(x) : g(x) 
fit h(x) 'data.dat' u 1:2 via a,b,c,d,e,f 
pl "data.dat" w p, h(x) 

의 코드를 변경 내가

a    = 0.156195   +/- 0.03023  (19.35%) 
b    = 1.40894   +/- 0.0678  (4.812%) 
c    = 0.00781003  +/- 0.00111  (14.21%) 
d    = 0.803641   +/- 0.1621  (20.17%) 
e    = 0.000178048  +/- 4.158e-05 (23.35%) 
f    = 0.00217076  +/- 0.0002812 (12.95%) 
있어

c와 d의 차이점을 확인 하시겠습니까? 그리고 e f는 거의 0입니다.

도움 주셔서 감사합니다.

+0

plot 'data.dat' w p, h(x) 불 벽 (중국에 있습니다. 당황 스럽네 ~~~). 어쩌면 gnuplot의 두 번째 축을 시도해 볼 수 있습니다. 예를 들어'plot sin (x) axis x1y1, 100 * sin (10 * x) axis x2y2' – DragonHu

+0

지금 가지고있는 스크립트를 포함 해주십시오. – Christoph

+0

각'fit' 및'plot' 명령 앞에 xrange를 설정하십시오. 나는 그것이 효과가 있을지 모르지만 시도 할 수 있습니다. – DragonHu

답변

0

나는 Great Fire Wall 뒤에 있기 때문에 여전히 데이터에 액세스 할 수 없습니다. 그러나 당신의 대본에서 분별할 때, 나는 당신이 원하는 것을 짐작할 수 있습니다. fit h(x) 'data.dat' via a,b,c,d 마지막으로

, 플롯 데이터와 시간 (X : H (x)를 사용하여 데이터에 맞게 정의, 첫째 원계 기능 그리고 h(x) = x < num ? f(x) : g(x)

:

지금은 대담하게 내 솔루션을 게시) : 나는 위대한 뒤에이기 때문에 내가 당신의 링크를 방문하지 못할

enter image description here

+0

당신이 접근 할 수있는 그림 또는 둘 다 할 수없는 그림? 당신이 볼 수있게 내가 이것을 게시 할 수있는 방법을 안다? –

+0

아니요. 질문에 그림을 게시 할 수있는 권한이 있습니까? – DragonHu

+0

아니, 적어도 10 명 이상 필요합니다. –