Verwendung von fit() mit Bedingungen für die erste Ableitung

Question

Ist es möglich, eine Kurve durch einige Punkte anzupassen und eine Bedingung für die resultierende Anpassung anzuwenden?

Ich mag die erste Ableitung des Sitzes eine bestimmte Wert (Kraft) bei x = 0 haben. Hier

ist eine Probe des Code:

options = fitoptions; 
options.Weights = [1000,1000,1,1,1,1,1,1,1]; 
myfittype = fittype('(exp(g*x))*(a*(x^5) + b*(x^4) + c*(x^3) + d*(x^2) + e*(x) + f)','coeff',{'g','a','b','c','d','f'},'problem',{'e'}); 
[moment_fit,gof1] = fit(xdata_fit',momentarray1',myfittype,'problem',{-force}); 
[shear_fit_values,soilrctn_fit_values] = differentiate(moment_fit,points'); 

Zuvor bevor die exponentielle in der benutzerdefinierten Gleichung Einführen I verwendet, um die Ableitung der Funktion zu zwingen, einen bestimmten Wert zu sein, indem einer der Koeffizienten, wie Sie steuern kann oben sehen (wenn es ein Polynom fünften Grades ohne Exponential ist, ist der Wert von e die Ableitung der Funktion bei 0).

Ist es möglich? Ich habe die Berechnungen auf der Hand und die Bedingung durchgeführt, die auf den Koeffizienten gesetzt werden muss, ist g*f + e = force

Danke

Answer 1

Warum nicht einfach manuell die Ableitung berechnen und die Koeffizienten in Plug?

myfit_derivative =@(x) exp(myfittype.g .* x) .* (myfittype.f.*myfittype.g + myfittype.e.*(1+myfittype.g.*x)+myfittype.d.*(2.*x+myfittype.g.*x^2)+myfittype.c.*(3.*x^2+myfittype.g.*x^3)+myfittype.b.*(4.*x^3+myfittype.g.*x^4)+myfittype.a.*(5.*x^4+myfittype.g.*x^5));