| 查看: 1174 | 回复: 0 | ||
[求助]
寻峰的函数！！
|
|
我这里的数据是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); 希望有大神能帮忙写出寻峰的语句！ |
回复此楼» 猜你喜欢
国自然申请面上模板最新2026版出了吗？
已经有12人回复
-
售中科院一区文章，我:8 O 5 5 1 O 5 4，含JCR-TOP
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有5人回复
-
溴的反应液脱色
已经有7人回复
-
推荐一本书
已经有12人回复
-
基金申报
已经有4人回复
-
计算机、0854电子信息（085401-058412）调剂
已经有4人回复
-
纳米粒子粒径的测量
已经有7人回复
-
常年博士招收(双一流，工科)
已经有4人回复
-
参与限项
已经有5人回复
