gnuplot fit paraméter korlátozása

Triviálisan nem tudod megcsinálni.

Azzal próbálkozhatsz, hogy két paramétert használsz, az egyik értéke tetszőleges lehet és ezt adod meg a 'fit' parancs paramétereként, a másik pedig a valódi, fizikai jelentéssel bíró paraméter, ami korlátos - a kettő között pedig egy alkalmas leképezés teremt kapcsolatot.

Lásd: http://groups.google.com/group/comp.graphics.apps.gnuplot/browse_thread…

gnuplottal en sem tudnam meg csinani (!= nem lehet), de az en 5 centem:
root.cern.ch -> download, feltenni kn 20 perc (ha ferdited), aztan a
tutorials/fit konyvtarban szetnezegetni

gyakorlatilag irsz egy fuggvenyt, ami kiszamolja a fuggvenyedet


double fuggveny(double* x, double* par){
//az x az azert array, mert hat tobb dimenzios a fuggvenyed, akkor x[0]=x, x[1]=y ...


a = par[0];
b = par[1];
x = x[0];
return (a/pow(b,a))*pow(x,(a-1))*exp(-pow((x/b),a));
}

double sumfuggveny(double* x, double* par){
c1 = par[0];
c2 = par[1];
if(0 <= (1-c1-c2) && (1-c1-c2) < 1){
return c1*fuggveny(x, &par[2]) +
c2*fuggveny(x, &par[4]) +
(1-c1-c2)*fuggveny(x, &par[6]);
}else
return 0;
}

a fit meg a kovetkezokeppen megy:
vannak az osszetartozo ertekparok (esetleg error?) az X es Y vektorokban, amelyeknek N eleme van


fit(){
TF1* f1 = new TF1("f1", sumfuggveny, xmin, xmax, 7);
//fuggveny ertelemezesi tartomanya [xmin, xmax] es 7 parametere van


f1->SetParLimit(0, 0, 1);
f1->SetParLimit(1, 0, 1);
// ez adja meg a limiteket a ket amplitudora
// (melyikparameter, alsolimit,felsolimit)

//f1->SetParameter(2, ...);//setting initial parameter a fitnek
//f1->FixParameter(3, ...);//Fix parameter
//f1->FixParameter(4, ...);//Release parameter

TGarph* gr = new TGraph(N, X, Y);
//X es Y N elemu tombok ... lasd feljebb

tg->Fit(f1, "ler");
//fit log-likelihood-dal ('l'),
//"jobb" error-becslessel, ('e').
//a fuggveny altal definialt ertelmezesi tartomanyon ('r');
}

A fit a hiresneves MINUIT-ot hasznalja a MIGRAD-dal (dokumentaciot a neten konnyen talalsz)

Persze 7 parameterhez trukkozni kell (fixalni elengedni parametereket ...) a MIGDRAD gradiebnst szamol, szoval
ha a kezdo parameterek nem eleg jok, a fuggveny nem linearis, es tul sok szabad parameter van, akkor konnyen elszall ...
ajanlom a dokumentacio fittolasrol szolo reszet (vsz valami uj, de nagyon konnyen bele lehet jonni), es a TFitPanel
nevu grafikus applikaciot ...

... ja ... es hogyan kell futtatni ...

megirod a beolvaso-algoritmust, ami megkrealja az X es Y tombot (*), osszemasolod a fuggvenyeket, a beolvasoreszt meg a fittolas kodot
egy fajlba, (ahol akkor van a 'fuggveny', 'sumfuggveny', 'fit' es a beolvaso fuggveneyd ... legyen a fajl neve fit.C) beinditod a root-ot ("root" a command-line-ba) es az interpreterjebol futtatod a kododat (.x fit.C [a kodnak abban a konyvtarban kellene lennie, ahonnan futtatod a root-ot ...])
ekkor az adataidat es fitfuggvenyt is megkapod egy canvasban ... ami szinezheto, feliratozhato, miegymas ...

[*van erre 1-soros megoldas, de a konvenciok miatt maceras leirni ...]

Szoval root-site: root.cern.ch
forum : http://root.cern.ch/phpBB2/
doksi : http://root.cern.ch/drupal/content/users-guide
fittolasi peldak: http://root.cern.ch/root/html/tutorials/fit/index.html
doxygen-szeru doksi: http://root.cern.ch/drupal/content/reference-guide

amennyiben hibaid is vannak, akkor ajanlom a TGraphError-t a TGraph helyett ...

remelem szereztem egy uj rajongot a ROOT-nak ...

k

a level hosszu lett, de ha osszeszamolod kb 10-15 sor a kod amit implementalni kell (meg az adataid beolvasasat ...)

A minuit elegge szeles korben hasznalt, legalabbis fizikaban, letezik neki standalone verzioja, de annak a hasznalata egy kicsit komplikaltabb (habar ha az ember ipari meretekben fittol, akkor hasznos abba is belejonni ...): http://c-minuit.sourceforge.net/

A root meg egy kicsit agyuval verebre, de egy rakat dolgot megcsinaltak (pld a minuit kore) ami meg az is lehet hogy hasznos.

na nem csinalok tobb reklamot neki ...
Probladd feltelepiteni, es adjal neki 20 percet (tutorials/fit konyvtarat vegigjatszva ...) ha nem tetszik, akkor nem tetszik ...

Maskulonben debian repositorikban benne szok' lenni ... a tobbirol meg nem tudok ...

k.