Settings21-1453205080.20000000000000001800.50311-1User DefinedGaussPeak_SigmoidBg_o2023b (User)0-1111101011510000::{Auto=1}175::{Auto=1}199001000::{Auto=1}344.101::{Auto=1}18221.5001000::{Auto=1}::{Auto=1}001000::{Auto=1}::{Auto=1}001000::{Auto=1}::{Auto=1}001000::{Auto=1}::{Auto=1}001000::{Auto=1}::{Auto=1}001000::{Auto=1}::{Auto=1}001000::{Auto=1}::{Auto=1}001000::{Auto=1}::{Auto=1}00095-10111001-10[Graph2]1!1"B"[176:200]-14001.0000000000000001E-9-10.00500000000000000014.9999999999999997E-304.9999999999999998E3000101-1//Gauss Peak: y=y0 + (A/(w*sqrt(pi/2)))*exp(-2*((x-xc)/w)^2) double gauss,sigmd; double ampG,ampS; double argG,argS; double tmp1,tmp2; ampG = A/(w*sqrt(pi/2)); ampS = ampG*Bg; argG = (x-xc)/(w/sqrt(2)); //Sigmoid argS = argG*2.116; //Gauss Peak gauss = ampG * exp(-argG^2); //Bg made by peak sigmd = ampS * (1 - 1/(1+exp(-argS))); if(TEST == 0) { // all y = y0 + gauss + sigmd; }else{ // y0 + sigmoid only y = TEST * (y0 + sigmd); } // partielle Ableitungen tmp1 = gauss*2/w*argG; tmp2 = A*Bg/w*exp(-argS)/(1+exp(-argS))^2; dy_y0 = 1; dy_A = (gauss+sigmd)/A; dy_xc = tmp1 - tmp2; dy_w = tmp1*argG - tmp2*(x-xc)/w; dy_Bg = sigmd/Bg;=0) ? abs(A):-abs(A); // peak area w *= 0.85; Bg = (yatxmin(x_y_curve) - yatxmax(x_y_curve)) / h; // Sigmoide amplitude // Folgende Fitter-Menü-Einstellungen sind für Doppel-Peak erforderlich: // Diese Einstellungen können im Fitter-Menü über Dialogdesign (oben, Dreieck rechts) // unter einem Namen gespeichert werden. // Einstellungen -> Datenauswahl -> Gewichtungen: keine Gewichtung // Einstellungen -> Erweitert -> // Kopien -> // Anzahl der Kopien: 1 // Einstellungen der Peaksuche -> // Suchmethode Peak: Residuum nach erster Ableitung // Peakrichtung: positiv // individuelle Kurve: ja // kumulative Kurve: ja // Eigenschaften -> Zusammenfassung des Fits -> Chi Quadr. reduziert: ja // Parameter -> // x width: teilen // sigmoid Bg: teilen // *** alte Initialisierung *** if(1==2) { // Get Fitter data smooth(x_y_curve); // Background y0 = yatxmax(x_y_curve); // Constant Bg = (yatxmin(x_y_curve) - yatxmax(x_y_curve)); // Sigmoide // Peak A = (max(y_data) - Bg/2 - y0); // Amplitude w = 0.6 * fwhm(x_y_curve, Bg/2+y0); // Width xc = xatymax(x_y_curve); // Center // Sigmoide background Bg = Bg / A; }]]>/*Enter general linear constraints here*//*Code to be executed before fitting, a good place for complicated initialization.*///Get NLFit tree from the fitting report sheet getnlr tr:=tt iw:=__REPORT$; int i, j; // Abgeleitete Parameter // Rechnung für A = Amplitude! Hier Area! double Area = tt.A * 2 * tt.w * sqrt(PI/2); double FWHM = sqrt(ln(4)) * 2 *tt.w; // Fehlerrechnung Area // Fehlertransformations-Vektor Create Tr -t 5; Tr[1] = 0; Tr[2] = 2*tt.w*sqrt(PI/2); Tr[3] = 0; Tr[4] = tt.A*2*sqrt(PI/2); Tr[5] = 0; // Hilfsvektor Create dstmp -t 5; for(i=1; i<=5; i++) { dstmp[i] = 0; }; // Fehlertransformation * Kovarianzmatrix * Fehlertransformation getresults tr:=tt2 iw:=__REPORT$; for(i=1; i<=5; i++) { for(j=1; j<=5; j++) { dstmp[i] = dstmp[i] + tt2.Covariance.Covariance1.col$(i-1).row$(j-1) * Tr[j]; // dstmp[i] = dstmp[i] + tt.Covar(j, i) * Tr[j]; }; }; double ErrA = 0; for(i=1; i<=5; i++) { ErrA = ErrA + Tr[i] * dstmp[i]; }; ErrA = sqrt(ErrA); // Fehlerrechnung FWHM double ErrF = sqrt(ln(4))*2 * tt2.Covariance.Covariance1.col$(3).row$(3) * sqrt(ln(4))*2; ErrF = sqrt(ErrF); //type ; //type "Abgeleitete Parameter fuer Gauss-Fit"; //type "mit Sigmoid-Untergrund"; //type Derived Parameters for; //type Gauss Peak Fit with Sigmoide Background; //type Peak Area: $(Area) +- $(ErrA); //type Peak FWHM: $(FWHM) +- $(ErrF); //type Wichtung: %(tt.Data1.wmethod1$); //type ; //if ( $(tt2.Covariance.Covariance1.col$(0).row$(0) ) == "--" ) { // type "Fehler koennen nicht berechnet werden"; // type "da die Kovarianzmatrix nicht zur Verfuegung steht!"; // type "-> NLFit -> Einstellungen -> Erweitert -> Kovarianzmatrix"; //}TEST=0-1-10010010095111000-111101111000100011100101101000-11000-1110-10001200-10000000-1372.3031763555108614213.02444417538187.003192589565462.76493940328681860.02178039271676197255543.657015656441188.99819758875208-10000000-10011001-11111111-1-1.23456789e-300-1.23456789e-300-1.23456789e-300-1.23456789e-300-1.23456789e-300-1.23456789e-300-1.23456789e-300-10000000-10000000-1372.3031763555108614213.02444417538187.003192589565462.76493940328681860.02178039271676197255543.657015656441188.99819758875208-10000000-1-10-1100000100000010-1010010-13*model*equ*plot*pm*rcs*cod*adj::{Auto=1||1}1010-100Name<auto>0-1<auto>Ba133<new>FitNL100<none>-1<auto>Ba133<new>FitNLCurve1-1<fittedvalue>Ba133<fittedvalue>FitNLCurve1-1<fittedvalue>Ba133<new>FitNLFindXfromY1FitNLFindYfromX1-10110020260929311522026094798688