スプライン小区分多項式を手動で作成する

配列xとyで表される値の配列でサンプリングされた関数y（x）を考えてみましょう。私はすでに係数とXを持っている場合は：私は別のx値x0を持っている場合、私は今、スプラインスプライン小区分多項式を手動で作成する

y0 = spline(x,y,x0);

を使用して、Y（X0）を評価することができ、私はまた

pp = spline(x,y); 
y0 = ppval(pp,x0);

私の質問を書くことができますmy_coefs（size(my_coefs) = [length(y),4]）とxの場合、pp.coefs = my_coefsとそのy0 = ppval(My_pp,x0)のような区分的多項式My_ppを作成するにはどうすればよいですか？

2016-11-25 Amir Sagiv

です。スプラインは多項式ではないことにご注意ください。それは細分化a（立方晶）多項式の部分区間にしかないのですか？ - あなたが望むのは、 'my_coefs'の行に与えられた多項式の配列を評価することです。行列にベクトル '[1、x、x^2、x^3]'を掛けるだけではどうですか？ – LutzL

私が知っているように、私が書いたように、私は 'polyval'ではなく' ppval'を使っています。 'spline'が返す' struct'についての素晴らしい点は、それがすでにポイントと多項式の独自の配列を持つ区分的多項式であることです。 –

それで、あなたが構築したいのは、スカラ値のスプラインです。基本的には、空のスプラインオブジェクトを生成し、そのデータフィールドを配列で埋める必要がありますか？ – LutzL

OK、「スプラインオブジェクト」ではなく、区分的な多項式オブジェクトがあります。 my_coefsがブレークポイントmy_xによって達成されたのであれば、その後、必要なコードは、その後彼らは寸法がここにdazzelingている場合には

my_spline = mkpp(my_x,my_coefs); 
y0 = ppval(my_spline, x0);

、

my_coefsは、n *

my_x 4ですnは

y0は

012 N

ありますはN

2016-11-25 19:17:01

答えて