matlab写一个能计算多峰半高宽的代码

以下是一个简单的MATLAB代码，用于计算多峰半高宽：

% 生成多峰高斯分布数据
x = linspace(-10,10,1000);
y = gaussmf(x,[1,-2])+gaussmf(x,[1,2])+gaussmf(x,[1,5]);

% 计算峰值位置
[pks,locs] = findpeaks(y);

% 计算半峰宽
for i = 1:length(pks)
    % 找到左右两侧的半高位置
    half_height = pks(i)/2;
    left_idx = find(y(1:locs(i))<=half_height,1,'last');
    right_idx = find(y(locs(i):end)<=half_height,1,'first')+locs(i)-1;
    
    % 计算半峰宽
    fwhm(i) = x(right_idx)-x(left_idx);
end

% 输出结果
disp(fwhm)

该代码首先生成一个多峰高斯分布数据，然后使用MATLAB内置函数`findpeaks`找到所有峰值位置。接下来，使用循环计算每个峰的半峰宽。在循环中，使用`find`函数找到左右两侧的半高位置，然后计算半峰宽。最后输出结果，即每个峰的半峰宽。