matlab - MATLAB中拟合高斯的最大值

Question

我需要找到我已经拟合的最大高斯，下面是我的示例代码（忽略它对高斯的拟合非常糟糕的事实，它们只是我在变量托盘中随意使用的两个备用矩阵）

x=[10;2;6;1;7;5;3;4;8;9];
z1=[0;0;0;0;0;6;7;8;9;10];
cf1=fit(x,z1,'gauss1');
plot(x,z1,'.k')
hold on
plot(cf1,'r')

你们能给我的任何帮助将不胜感激。

Answer 1

这是非常基本的数学。看看输出

>> cf1

cf1 = 
     General model Gauss1:
     cf1(x) =  a1*exp(-((x-b1)/c1)^2)
     Coefficients (with 95% confidence bounds):
       a1 =       5.187  (-0.4711, 10.85)
       b1 =       6.834  (-0.768, 14.44)
       c1 =       5.945  (-8.833, 20.72)

现在，借助有关 Gaussians 的维基百科文章，找到最大值很容易：

maximum_x = cf1.b1;
maximum_y = cf1.a1;

对于您用来拟合相同函数的任何其他工具，同样如此——分别定义最大值的系数和a1位置。b1yx