AM、SSB语音信号频谱图matlab实现代码

以下是一个简单的实现AM和SSB语音信号频谱图的Matlab代码：

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

% 创建时间轴
dt = 1/Fs;
t = 0:dt:(length(y)*dt)-dt;

% 绘制原始波形图
subplot(3,1,1);
plot(t, y);
xlabel('时间 (s)');
ylabel('幅度');
title('原始语音信号波形图');

% 使用调幅（AM）调制
fc = 1000; % 载波频率
Ac = 1; % 载波幅度
ym = Ac*(1+0.5*y).*cos(2*pi*fc*t); % 调制后的信号

% 绘制AM调制后的波形图
subplot(3,1,2);
plot(t, ym);
xlabel('时间 (s)');
ylabel('幅度');
title('AM调制后的语音信号波形图');

% 计算AM调制后的频谱
Ym = fft(ym);
L = length(ym);
P2 = abs(Ym/L);
P1 = P2(1:L/2+1);
P1(2:end-1) = 2*P1(2:end-1);
f = Fs*(0:(L/2))/L;

% 绘制AM调制后的频谱图
subplot(3,1,3);
plot(f,20*log10(P1));
xlabel('频率 (Hz)');
ylabel('幅度 (dB)');
title('AM调制后的语音信号频谱图');

% 使用单边带抑制载波（SSB-SC）调制
fc = 1000; % 载波频率
Ac = 1; % 载波幅度
yssb = y.*cos(2*pi*fc*t) - hilbert(y).*sin(2*pi*fc*t); % 调制后的信号

% 绘制SSB调制后的波形图
figure;
subplot(3,1,1);
plot(t, yssb);
xlabel('时间 (s)');
ylabel('幅度');
title('SSB调制后的语音信号波形图');

% 计算SSB调制后的频谱
Yssb = fft(yssb);
L = length(yssb);
P2 = abs(Yssb/L);
P1 = P2(1:L/2+1);
P1(2:end-1) = 2*P1(2:end-1);
f = Fs*(0:(L/2))/L;

% 绘制SSB调制后的频谱图
subplot(3,1,2);
plot(f,20*log10(P1));
xlabel('频率 (Hz)');
ylabel('幅度 (dB)');
title('SSB调制后的语音信号频谱图');

请注意，此代码仅提供了一个简单的实现，并且可能需要进行调整以适应您的特定需求。