当前位置：首页 > news >正文

基于MATLAB的音频信号AM调制与解调实现

news 2026/5/11 21:54:03

AM（幅度调制）是一种通过改变载波幅度传递信息的经典调制技术。

一、MATLAB实现步骤

1. 读取音频信号

% 读取WAV文件（单声道）
[filename, fs] = audioread('audio.wav');
audio = audio(:, 1); % 提取单声道
t = (0:length(audio)-1)/fs; % 时间向量

2. 参数设置与信号预处理

% 载波参数
fc = 192e3; % 载波频率（Hz）
fs_carrier = 1.92e6; % 载波采样率（需为音频采样率的整数倍）
t_carrier = (0:length(audio)-1)/fs_carrier; % 载波时间向量% 调制参数
A0 = 0.5; % 直流偏置
beta = 0.8; % 调制指数% 插值调整音频采样率
audio_interp = interp1(t, audio, t_carrier, 'linear');

3. AM调制

% 生成AM信号
carrier = cos(2*pi*fc*t_carrier);
am_signal = (A0 + beta*audio_interp) .* carrier;

4. 信道传输（可选加噪）

% 添加高斯白噪声（SNR=20dB）
noise_power = var(am_signal)/10^(20/10);
noise = sqrt(noise_power)*randn(size(am_signal));
received_signal = am_signal + noise;

5. AM解调（包络检波）

% 包络检波步骤
envelope = abs(hilbert(received_signal)); % 希尔伯特变换提取包络
demod_signal = envelope - mean(envelope); % 去除直流分量% 低通滤波（恢复原始信号）
[b, a] = butter(6, 2 * 4.8e3/fs_carrier); % 截止频率为音频最高频率
filtered_signal = filtfilt(b, a, demod_signal);

6. 结果对比与可视化

% 时域对比
figure;
subplot(3,1,1); plot(t, audio); title('原始音频');
subplot(3,1,2); plot(t_carrier, am_signal); title('AM调制信号');
subplot(3,1,3); plot(t, filtered_signal); title('解调信号');% 频谱分析
N = length(audio);
f = (-N/2:N/2-1)*(fs_carrier/N);
S_am = fftshift(fft(am_signal));
S_demod = fftshift(fft(filtered_signal));
figure;
subplot(2,1,1); plot(f, abs(S_am)); title('AM信号频谱');
subplot(2,1,2); plot(f, abs(S_demod)); title('解调信号频谱');

二、常见问题与解决方案

频谱混叠

原因：载波频率或采样率设置不当（如载波频率超过奈奎斯特极限）。

解决：通过resample函数调整音频采样率，确保载波频率满足\(f_c≥2fmax\)。
包络失真

原因：调制指数\(β>1\)。

解决：根据音频信号带宽计算最大允许\(β\)，例如4.8kHz音频的\(β_{max}=0.8\)。
解调信号延迟

原因：希尔伯特变换引入的相位延迟。

解决：使用filtfilt函数进行零相位滤波补偿。

三、扩展应用

多音频混合调制

% 混合两路音频信号
audio2 = audioread('music.wav');
mixed_signal = 0.5*audio + 0.5*audio2;
am_signal = (A0 + beta*mixed_signal) .* carrier;

抗干扰仿真

添加多径干扰或脉冲噪声，评估误码率（BER）：

% 添加多径干扰
delay = round(0.1*fs_carrier);
multipath = [zeros(1, delay), am_signal(1:end-delay)];
received_signal = am_signal + 0.3*multipath;

参考代码 matlab模拟音频信号的AM调制与解调 www.youwenfan.com/contentcnn/95841.html

四、完整代码示例

% AM调制解调完整流程
[filename, fs] = audioread('audio.wav');
audio = audio(:, 1);
t = (0:length(audio)-1)/fs;% 参数设置
fc = 192e3; fs_carrier = 1.92e6;
A0 = 0.5; beta = 0.8;% 信号处理
audio_interp = interp1(t, audio, (0:length(audio)-1)/fs_carrier, 'linear');
carrier = cos(2*pi*fc*(0:length(audio_interp)-1)/fs_carrier);
am_signal = (A0 + beta*audio_interp) .* carrier;% 解调
envelope = abs(hilbert(am_signal));
demod_signal = envelope - mean(envelope);
[b, a] = butter(6, 2 * 4.8e3/fs_carrier);
filtered_signal = filtfilt(b, a, demod_signal);% 播放对比
sound(audio, fs);
pause(length(audio)/fs + 1);
sound(filtered_signal, fs);