多次元配列に対するPolyfitの高速化

Question

Kouta Amano 2019-9-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/478773-polyfit

评论： Kouta Amano 2019-9-4

采纳的回答： Kazuya

在 MATLAB Online 中打开

私は測定水準x, 測定値yに対し y=a0 + a1*x + a2*x^2 の形でカーブフィッティングを行いたいと考えています。

例えばxの次元は100万×3, yの次元は100万×3であり、１行１行について上記カーブフィッティングを行いたいです。

arrayfunを用いて下記のような関数を作成すると計算は可能でしたが実行速度が大変遅いです。

実行速度を改善するために、よりよい実装方法はないでしょうか。

function coefs = multiPolyfit(x,y,dim)
%%%%%%%%%%%%%%%%%%%%%%%
% x: array型。サンプルを行方向、列方向は水準。
% y: array型。サンプルを行方向、列方向は水準。
% dim: int型
%%%%%%%%%%%%%%%%%%%%%%%
s = table2struct(table(x,y,'VariableNames',{'x','y'}));
coefs = arrayfun(@(z) polyfit(z.x,z.y,dim),s,'UniformOutput',false);
end