2016-07-30 2 views
1

У меня есть вопрос относительно построения экспоненциальной подгонки в Gnuplot и наличия y-axis, установленного в logscale. Я использую Gnuplot 4.6.3. я представлю несколько способов, я попытался установить диапазон, чтобы поместиться в:

функция Я хочу, чтобы поместиться на данных е (х) = а + Ь * х

Первое:

set log y 
f(x)=a+b*x 
fit [4:8] f(x) 'CT2A_OH_R_log.dat' using (log($1)):($2) via a,b 

Resutls:

  Read 15 points 
     Skipped 15 points outside range [x=4:8] 
     No data to fit 
     error during fit 

Второе:

fit f(x) 'CT2A_OH_R_log.dat' using (log($1)):($2) via a,b 

Результаты:

Final set of parameters    Asymptotic Standard Error 
a   =16.9858     +/- 6.299   (37.08%) 
b   =-8.43215     +/- 3.502   (41.53%) 
gnuplot> plot 'CT2A_OH_R_log.dat', f(x) via a,b 

Очевидно, что подгонка смешно, но это единственный способ я знаю, что на самом деле дает мне истерику. Я знаю, что проблема, скорее всего, лежит в определении xrange при установке, но что еще я могу написать там, чтобы исправить это?

Ниже я покажу изображение моего участка (е (х) не показывают, как она изображена в отрицательной зоне.)

ответ

1

«вызывает logscale из у axis требует логарифма оси у при подгонка ;-)

fit [4:8] f(x) 'CT2A_OH_R_log.dat' u 1:(log($2)) via a,b 

ПРИМЕЧАНИЕ

Если вы хотите использовать логарифм на основе 10, вы должны использовать log($2)/log(10) или log10($2).

+0

Спасибо! Я забыл такую ​​глупую ошибку! Тем не менее, я сделал это сейчас, и я все еще получаю странную форму: gnuplot> fit [4: 8] [0,1: 60] f (x) 'CT2A_OH_R_log.dat', используя 1: (log ($ 2)) через , b gnuplot> plot 'CT2A_OH_R_log.dat', f (x) через a, b На этот раз подгонка экспоненциальна. Я хочу добавить картинку, но я не могу этого сделать в разделе комментариев. Итак, см. Комментарий ниже: – Omega

0

@tom Solid

gnuplot> fit [4:8] [0.1:60] f(x) 'CT2A_OH_R_log.dat' using 1:(log($2)) 
via a,b 
gnuplot> plot 'CT2A_OH_R_log.dat', f(x) via a,b 

На этот раз подгонка экспоненциальный как.

Итак, я вывешу картину того, что я получаю в настоящее время:

enter image description here