用matlab绘制一段语音信号中800Hz成分短时相位随时间的变化曲线

好的，下面是用 Matlab 绘制一段语音信号中 800Hz 成分短时相位随时间的变化曲线的代码：

% 读入语音信号
[x, fs] = audioread('sample.wav');

% 定义窗口长度
window_length = 0.02;

% 定义窗口数
num_windows = floor(length(x)/(window_length*fs));

% 初始化相位数组
phi = zeros(1,length(x));

% 计算每个窗口内的相位变化
for i = 1:num_windows
    start_index = round((i-1)*window_length*fs)+1;
    end_index = round(i*window_length*fs);
    x_window = x(start_index:end_index);
    f, P = periodogram(x_window,[],[],fs,'power');
    [~,ind] = max(P);
    f0 = f(ind);
    phi(start_index:end_index) = 2*pi*f0*(0:1/fs:(window_length-1/fs));
end

% 绘制 800Hz 成分短时相位随时间的变化曲线
t = (0:length(x)-1)/fs;
plot(t,phi);
xlabel('时间 (s)');
ylabel('相位 (rad)');
title('800Hz 成分短时相位随时间的变化曲线');

在这个代码中，我们首先读入了语音信号，然后定义了窗口长度为 0.02 秒，也就是每个窗口内有 0.02*fs 个样本点，其中 fs 是采样率。接着，我们计算了窗口数，以及每个窗口内的 800Hz 成分的相位变化。具体地，我们使用了 periodogram 函数计算了每个窗口内的功率谱密度，然后找到了最大功率密度对应的频率，即为该窗口内的 800Hz 成分频率。最后，我们使用 plot 函数绘制了 800Hz 成分短时相位随时间的变化曲线。

需要注意的是，这里我们假设语音信号中只有一个主频率成分。如果存在多个主频率成分，则需要在计算相位变化时对每个主频率成分分别进行处理。

希望这个代码可以帮助你！