请为我写一段matlab代码，完成本地音频的PM调制与解调

好的，以下是一个简单的MATLAB代码，用于实现本地音频的PM调制和解调：

% 读取音频文件
filename = 'test_audio.wav';
[y, Fs] = audioread(filename);

% 载波信号
fc = 5000;
t = 0:1/Fs:length(y)/Fs-1/Fs;
carrier = sin(2*pi*fc*t);

% PM调制
kf = 500;
pm = sin(2*pi*fc*t + kf*y');

% PM解调
demod = zeros(size(pm));
for i = 2:length(pm)
    demod(i) = atan2(sin(pm(i)-pm(i-1)), cos(pm(i)-pm(i-1)));
end
demod = demod/max(abs(demod)); % 归一化

% 播放原始音频和解调音频
sound(y, Fs);
pause(length(y)/Fs);
sound(demod, Fs);

% 绘制结果
subplot(3, 1, 1);
plot(y);
title('Original Audio Signal');
subplot(3, 1, 2);
plot(pm);
title('PM Modulated Signal');
subplot(3, 1, 3);
plot(demod);
title('Demodulated Audio Signal');

在这个示例中，我们首先使用`audioread`函数读取本地音频文件，并将其存储在`y`中。然后，我们使用5000Hz的载波信号调制音频信号，生成调制后的信号。在PM解调中，我们使用相移检测器来恢复原始音频信号。最后，我们使用`sound`函数播放原始音频文件和解调后的音频文件，并绘制原始音频、调制信号和解调音频的图形。