曲面近似のhorzcatエラー

13x22行列で表示される離散データ eta_w を曲面近似しようと試みているのですが，

エラー: horzcat

連結する配列の次元が一致しません。

エラー: b_03_working_efficiency (行 13)

sf = fit([Lb(:,1), d(1,:)],eta_w,'poly23');

% Lb(:,1)は13行のベクトル, d(1,:)は22列のベクトル

と表示されてしまいます．

eta_wの一部にNaN表示されるデータのない点があることが原因でしょうか？解決方法を教えていただきたいです．

参考に，データ点は図のようになっています．