clc; clear; close all; clear all; fs = 8000; [y, fs] = audioread("qingtian.wav"); sigPlayer = audioplayer(y(:, 1), fs); myFigure = figure('Name', '标题'); waveAxes = subplot(4, 1, 1); axis(waveAxes, [-inf inf -1 1]); waveAxes.NextPlot = 'replacechildren'; fftAxes = subplot(4, 1, 2); axis(fftAxes, [-inf inf 0 0.5]); fftAxes.NextPlot = 'replacechildren'; windLength = fs / 20; sigPlayer.TimerPeriod = windLength / fs; sigPlayer.TimerFcn = {@fleshPlot, waveAxes, fftAxes, y(:, 1), fs, windLength}; play(sigPlayer); % 添加混响效果 reverbObject = reverberator('PreDelay', 0.5, 'WetDryMix', 0.125); sound(y); while 3 > 2 [audioIn, Overrun] = step(sigPlayer); % 混响处理 reverbAudio = reverb(reverbObject, audioIn); ydata_fft = fft(reverbAudio); ydata_abs = abs(ydata_fft(1:windLength/2)); set(pic, 'ydata', audioIn); set(pic2, 'ydata', log(ydata_abs)); drawnow; end function fleshPlot(plr, ~, axwave, axfft, sigData, fs, winL) persistent ind; persistent indbase; persistent wavePlot; persistent fftPlot; if isempty(indbase) indbase = 1:winL; end ind = fix((plr.CurrentSample) / winL) * winL + indbase; segData = sigData(ind); wavePlot = plot(axwave, segData); [absFFT, fqax] = fftc2(segData, fs); fftPlot = bar(axfft, fqax, absFFT); end function [X, freqRange] = fftc2(x, fs) N = length(x); X = abs(fft(x)) / N; rawRange = 0:fix(N / 2); freqRange = (rawRange) * fs / N; LfftInd = rawRange + 1; X = X(LfftInd); X(2:end) = 2 * X(2:end); end

这段代码是用于播放音频文件并添加混响效果的。首先，代码使用`audioread`函数读取名为"qingtian.wav"的音频文件，并创建一个`audioplayer`对象用于播放音频。然后，创建一个图形窗口和两个子图用于显示波形图和频谱图。接下来，设置混响效果并使用`reverberator`函数创建一个混响对象。然后，使用`sound`函数播放原始音频。在一个循环中，通过`step`函数从`audioplayer`对象中获取音频数据，并对音频数据进行混响处理。然后，使用FFT算法计算音频数据的频谱，并更新波形图和频谱图的数据。最后，使用`drawnow`函数更新图形窗口的显示。

值得注意的是，代码中的函数`fleshPlot`用于更新波形图和频谱图的数据。函数`fftc2`用于计算FFT并返回幅度谱和频率范围。

请注意，在运行此代码之前，您需要确保已经安装了Matlab并且具有相应的音频文件。