在matelab中讲音频文件进行FM调制加上噪声进行解调的代码，显示调制的图片，显示加上噪声后解调后的图片

在MATLAB中进行音频文件的FM调制并加上噪声解调的过程可以分为几个步骤：

1. 读取音频文件。
2. 执行FM调制。
3. 添加噪声。
4. 执行FM解调。
5. 显示调制前后及解调后的信号。

以下是一个简化版的MATLAB代码示例，用于说明这个过程：

% 读取音频文件
[signal, Fs] = audioread('audiofile.wav'); % 替换为你的音频文件名

% FM调制参数
fc = 1000; % 载波频率
kf = 1;    % 频率偏移常数
t = (0:length(signal)-1)/Fs; % 时间向量

% 调制信号
carrier = cos(2*pi*fc*t);
modulated = cos(2*pi*fc*t + 2*pi*kf*signal);

% 显示调制信号的频谱
figure;
fft_modulated = fft(modulated);
f = Fs*(0:(length(fft_modulated)/2))/length(fft_modulated);
plot(f, abs(fft_modulated(1:length(fft_modulated)/2)));
title('FM Modulated Signal Spectrum');
xlabel('Frequency (Hz)');
ylabel('Magnitude');

% 添加噪声
noisy_modulated = modulated + 0.5*randn(size(modulated));

% FM解调
% 使用简单的包络检测法进行解调（实际应用中可能需要更复杂的解调器）
demodulated = abs(hilbert(noisy_modulated));

% 显示解调信号的频谱
figure;
fft_demodulated = fft(demodulated);
plot(f, abs(fft_demodulated(1:length(fft_demodulated)/2)));
title('FM Demodulated Signal Spectrum (with Noise)');
xlabel('Frequency (Hz)');
ylabel('Magnitude');

请注意，上述代码是一个简化的示例，它使用了包络检测法作为FM解调方法。在实际应用中，可能需要使用更精确的同步检测法或相位锁定环（PLL）解调器。

此外，对于解调的处理，这里没有考虑频率偏移和噪声的影响，这可能会导致解调信号失真。在实际解调时，需要更精细的处理来确保解调信号的质量。