matlab计算isnr的程序

以下是计算ISNR的Matlab程序：

function isnr = calculate_isnr(original_signal, denoised_signal)
% 计算信噪比增益比(ISNR)
% original_signal: 原始信号
% denoised_signal: 经过降噪处理后的信号
% isnr: ISNR

% 计算信号噪声能量
noise_energy = sum((original_signal - denoised_signal) .^ 2);
% 计算降噪后信号的信噪比
denoised_snr = 10 * log10(sum(original_signal .^ 2) / noise_energy);
% 计算原始信号的信噪比
original_snr = 10 * log10(sum(original_signal .^ 2) / sum(denoised_signal .^ 2));
% 计算ISNR
isnr = denoised_snr - original_snr;
end

使用方法：

% 假设original_signal和denoised_signal都是长度为N的列向量
isnr = calculate_isnr(original_signal, denoised_signal);
fprintf('ISNR: %.2f dB\n', isnr);

其中，ISNR表示信噪比增益比，单位为dB。