matlab 光谱相似度指数

在MATLAB中，可以使用spectrum_similarity函数来计算两个信号的光谱相似度指数。该函数需要输入两个信号的频谱，可以通过MATLAB中的fft函数来获取。具体用法如下：

% 假设有两个信号x和y，采样频率为fs
% 分别计算它们的频谱
X = abs(fft(x));
Y = abs(fft(y));

% 计算光谱相似度指数
similarity = spectrum_similarity(X, Y, fs);

其中，fs为信号的采样频率，spectrum_similarity函数的定义如下：

function similarity = spectrum_similarity(X, Y, fs)
    % 计算两个信号的光谱相似度指数
    % X, Y: 信号的频谱
    % fs: 信号的采样频率

    % 计算频率向量
    f = linspace(0, fs/2, length(X));

    % 计算相似度指数
    similarity = sum(min(X, Y).^2) / sum(X.^2 + Y.^2 - min(X, Y).^2);
end

该函数返回的相似度指数在0到1之间，数值越大表示两个信号的相似度越高。