基于MATLAB R2022a提取wav音频文件的基音频率并输出图像结果

首先，你需要使用MATLAB内置函数`audioread`读取音频文件。然后，你可以使用自定义函数`pitchtrack`来提取基音频率，并将其绘制成一个图像。

下面是一个示例代码：

% 读取音频文件
[signal, fs] = audioread('example.wav');

% 提取基音频率
pitch = pitchtrack(signal, fs);

% 绘制图像
plot(pitch);
title('基音频率图像');
xlabel('时间');
ylabel('频率');

其中`pitchtrack`函数的代码如下：

function p = pitchtrack(x, fs)

    % 设置分析窗口长度
    winlen = round(0.03 * fs); % 30ms

    % 设置频率范围
    minf0 = 80;
    maxf0 = 600;

    % 计算自相关函数
    autocorr = xcorr(x, winlen, 'coeff');
    autocorr = autocorr(winlen+1:end);

    % 计算基音周期
    f0 = zeros(1, length(x));
    for n = winlen+1:length(x)
        % 计算自相关峰值
        [peaks, locs] = findpeaks(autocorr(n-winlen:n));
        if isempty(peaks)
            continue;
        end
        % 选择最大峰值
        [~, maxidx] = max(peaks);
        loc = locs(maxidx);
        % 计算基音周期
        period = loc / fs;
        if period < 1/minf0 || period > 1/maxf0
            continue;
        end
        f0(n) = 1 / period;
    end

    % 进行插值处理
    p = interp1(find(f0), f0(f0~=0), 1:length(x), 'linear', 'extrap');
end

运行代码后，会生成一个基音频率的图像，其中横轴表示时间，纵轴表示频率。你可以根据需要修改代码来调整分析窗口长度和频率范围。