音乐信号的AM调制 （1）观察音乐信号频率上限，选择适当调制频率对信号进行调制（给出高、低两种调制频率）； （2）输出调制信号的波形和频谱，观察现象，给出理论解释； （3）播放调制音乐信号，注意不同调制频率下的声音，解释现象。matlab代码

好的，以下是MATLAB代码实现：

%% 读入音频文件
[x,fs] = audioread('music.wav');

%% 观察音乐信号频率上限，选择适当调制频率
fmax = max(abs(fft(x)));
fm_low = 5000; % 低调
fm_high = 10000; % 高调

%% 对音频信号进行AM调制
fc = 100000; % 载波频率
m = 0.5; % 调制指数
t = (0:length(x)-1)/fs;
y_low = (1 + m*cos(2*pi*fm_low*t)) .* x;
y_high = (1 + m*cos(2*pi*fm_high*t)) .* x;

%% 输出调制信号的波形和频谱
figure;
subplot(2,1,1);
plot(t,x,'b');
hold on;
plot(t,y_low,'r');
plot(t,y_high,'g');
title('波形图');
xlabel('时间');
ylabel('幅值');
legend('原始信号','低调信号','高调信号');
subplot(2,1,2);
N = length(x);
f = (0:N-1)/N*fs;
X = abs(fft(x));
Y_low = abs(fft(y_low));
Y_high = abs(fft(y_high));
plot(f(1:N/2),X(1:N/2),'b');
hold on;
plot(f(1:N/2),Y_low(1:N/2),'r');
plot(f(1:N/2),Y_high(1:N/2),'g');
title('频谱图');
xlabel('频率');
ylabel('幅值');
legend('原始信号','低调信号','高调信号');

%% 播放调制音乐信号
sound(y_low,fs);
pause;
sound(y_high,fs);

在上述代码中，我们首先采用`audioread`函数读入音频文件，然后计算音乐信号的频率上限，选择适当的调制频率，包括高调和低调两种。然后我们使用AM调制公式对音频信号进行调制，得到调制音乐信号。接着，我们输出调制信号的波形和频谱图，观察调制信号的振幅和侧带现象。最后，我们播放调制音乐信号，并注意不同调制频率下的声音变化。

需要注意的是，该程序中的载波频率、调制指数等参数均可根据实际需求进行调整。