在MATLAB中,你可以使用awgn
函数给信号添加高斯白噪声来实现给矩阵信号添加特定信噪比的噪声,并使用plot
函数可视化添加噪声前后的信号。以下是示例代码:
matlab
% 创建原始信号
Fs = 1000; % 采样率
t = 0:1/Fs:1-1/Fs; % 时间轴
signal = sin(2*pi*10*t); % 频率为10Hz的正弦信号
% 添加噪声前的信号绘图
subplot(1, 2, 1)
plot(t, signal)
title('Original Signal')
% 添加噪声
snr_db = 10; % 信噪比为10dB
noisy_signal = awgn(signal, snr_db, 'measured');
% 添加噪声后的信号绘图
subplot(1, 2, 2)
plot(t, noisy_signal)
title(['Noisy Signal (SNR = ' num2str(snr_db) 'dB)'])
sgtitle('Signal with Noise')
% 设置subplot间距
sgtitle('Signal with Noise')
% 设置subplot间距
subplotSpacing = 0.05;
subplotPosition = get(gcf, 'DefaultAxesPosition');
subplotPosition(3) = subplotPosition(3) + subplotSpacing;
set(gcf, 'DefaultAxesPosition', subplotPosition)
% 调整窗口大小以适应subplot间距
set(gcf, 'Position', get(0, 'Screensize'))
% 等待按键继续执行下一条命令
pause
运行上述代码,你将会看到一个显示原始信号和添加噪声后的信号的图表,噪声的信噪比为10dB。你可以根据需要修改信噪比值来观察不同信噪比下的效果。