Matlab - polyvalデータの微分方程式

Question

私はカーブフィッティング法でいくつかの問題に取り組んでおり、立方体フィットの50点を見つける必要があります。質問は、与えられた細菌の増殖率を求める質問です。私の現在のコードは、私がgrowthRateまたは誘導体は、この問題のために正しいかどうかを疑問に思って、または、彼らは両方とも間違っている場合

time = [1,2,4,5,7,9]; 
bacteria = [2000,4500,7500,15000,31000,64000]; 

rcubic = polyfit(time,bacteria,3); 
newTime = linspace(1,7,50); 
vrcubic = polyval(rcubic,newTime); 

growthRate = [diff(vrcubic)./diff(newTime)]; 
derivative = diff(vrcubic) 

です。彼らは著しく異なる価値を与えたからです。また、ベクトルの長さが短くなると、時間に対するvrcubicの

乾杯

Answer 1

誘導体は確かにあなたのderivative式は、隣接データポイント間だけの違いである

growthRate = [diff(vrcubic)./diff(newTime)]; 

で与えられます。

また、それがNEWTIMEポイントが非常にまばらにサンプリングされている場合であっても正しい誘導体であるように私はいつも後者を好むだろう

ここ

growthRate = polyval(polyder(rcubic), newTime); 

を分析分化を使用することができます。