24小时热门版块排行榜

返回列表

明00媛

金虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1330.1
散金: 5
帖子: 35
在线: 11小时
虫号: 4066124
注册: 2015-09-11
性别: GG
专业: 基础物理学

[求助] 寻峰的函数！！

我这里的数据是pgm的，我将其处理成多个高斯拟合的形式，现阶段只能将其最大的那个拟合出来，其他的高斯拟合我需要找到其峰值的位置！
我把前边的语句先列举上：
function [ OutArr ] = Graph( NameFile, PlotOut, Angl)
InputVar = nargin;
DAT=importdata(NameFile);

if (InputVar>2)
   Angle = Angl;
else
      Angle = 0;
end

DAT = imrotate(DAT, Angle, 'nearest');
IMG=medfilt2(DAT,[3,3]);
[md, nd] = size(IMG);
disp({'M=' ,md , 'N=', nd});

Noise = mean(IMG(1:40));

i=0:nd-1;
X2D=i;

if (InputVar >1)&(PlotOut == 'Y')
Hendl=figure;
imagesc(IMG);
end
[Max, INmax] = max(IMG);
[Max, IMmax] = max(Max);
Nmax = INmax(IMmax);
disp({'Nmax = ', Nmax; 'Mmax = ', IMmax});

Top=Nmax+20;
Bottom=Nmax-20;

Y2D(i+1)=mean(IMG(Bottom:Top,i+1));
Y2D = Y2D-Noise;
Y2D = Y2D/max(Y2D);
gaussEqn = 'a*exp(-((x-b)/c)^2)+d';
startPoints = [1, IMmax, 1, 0];
f = fit(X2D', Y2D', gaussEqn, 'Start', startPoints);
FitGaussD = f(X2D);
Cint = confint(f, 0.95);

希望有大神能帮忙写出寻峰的语句！

回复此楼

» 猜你喜欢

不骄不躁，脚踏实地

1楼 2016-06-21 18:10:33

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主明00媛的主题更新

返回列表