如何将拟合曲线的结果从“cfit”转换为 Matlab 中的实际数组数据？

Question

我想拟合某个曲线。我的数据存储在“x”和“y2”中。为了拟合这条曲线，有一些预定义的函数，比如“fit”。

我的例子：

x = [20 30 40 50 60 70 80 90 100 200];
y = [4.229 2.514 1.714 1.143 0.8 0.629 0.514 0.4 0.349 0.057];
p_0 = 2*10^(-5);
y2 = 10.^(y/20)*p_0;

[xData,yData] = prepareCurveData(x,y2);
fit_type = 'pchipinterp';

[fitresult] = fit(xData,yData,ft);

结果“fitresult”现在作为“1x1 cfit”对象存储在 MATLAB 工作区中。 有什么方法可以从 fitresult 中获取一个我可以手动定义的数组长度的实际数组？ （例如 MATLAB 工作区：fitresult: 1x3841 double）

Answer 1

编辑你的脚本：

fit_type = 'pchipinterp';
[fitresult] = fit(xData,yData,fit_type);

fitresult 有两个类 Coefficients 和 Shape-preserving (pchip) interpolant，每一个都有一个分量 p 和 拟合结果（x） 拟合结果 =

 Shape-preserving (pchip) interpolant:
   fitresult(x) = piecewise polynomial computed from p
 Coefficients:
   p = coefficient structure

要提取结果，请将这些代码行添加到您的脚本中（名称由您决定）：

1. piecewise_polynomial_computed_from_p = fitresult(x)      
2. coefficient_structure = fitresult.p

这里的 coefficient_structure 是一个构造函数：

coefficient_structure = 
      form: 'pp'
    breaks: [20 30 40 50 60 70 80 90 100 200]
     coefs: [9x4 double]
    pieces: 9
     order: 4
       dim: 1

要从中提取结果，例如：

Breaks = coefficient_structure.breaks

结果是：

Breaks =
    20    30    40    50    60    70    80    90   100   200