2016-09-20 7 views
0

いくつかの点に曲線をフィットさせ、結果として得られるフィットに条件を適用できますか?fit()を1次導関数の条件に使用する

Iは、x = 0で一定(力)を有するようにフィットの一次導関数したいです。ここ

コードのサンプルである:

私はあなたのように係数のいずれかを制御することにより、特定の値に関数の導関数を強制するために使用されるカスタム式に指数関数を導入する前に以前
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'); 

上記を見ることができます(指数関数のない5次多項式の場合、eの値は0での関数の微分値です)。

それは可能ですか?私は手で計算を行っていると係数に配置する必要がある条件は

g*f + e = force

答えて

1

だけ手動で微分を計算し、係数を接続しないのはなぜでありがとうございますでしょうか?

myfit_derivative [email protected](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)); 
+0

あなたが行ったことを理解していません。コードを説明してください。 @xとは何ですか? –

+0

あなたが書いたコードは、私が公式化する方法を知っている派生物です。私はその方程式をx = 0で一定の値にします。私は元の方程式を当てはめている間にこれを強制することができますか? –

+0

これはフィッティングの問題です。これらの係数の背後に物理的/数学的な推論がありますか?これ以外の点では解決策の無限大があるので – user1543042

関連する問題