读入一段 10 秒左右的音乐，利用 M 倍抽取语句 xD=x(1:M:end)模拟对连续信 号的抽样，编程完成以下任务： a. 利用函数 ctft1 分析原信号的频谱，画出其幅度谱，大致确定原信号的最 高频率； b. 分别用 4 倍、2 倍、1 倍和 0.5 倍的最高频率对原信号进行抽样，记录抽 样频率，画出相应的时域波形，播放抽样后的音乐，描述抽样频率对声音 效果的影响。

好的，我明白了。这是一个信号处理问题，需要使用 MATLAB 进行编程。以下是我的回答：

a. 代码如下：

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

% M 倍抽样语句
M = 4;
xD = x(1:M:end);

% 利用 ctft1 分析频谱
N = length(xD);
f = (0:N-1)/N*Fs/M;
X = ctft1(xD);
plot(f, abs(X));
xlabel('Frequency (Hz)');
ylabel('Magnitude');

运行上述代码可以得到原信号的幅度谱，从图中我们可以大致确定原信号的最高频率。

b. 代码如下：

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

% 最高频率
f_max = 5000; % 假设最高频率为 5000 Hz

% 抽样频率
Fs4 = f_max * 4;
Fs2 = f_max * 2;
Fs1 = f_max * 1;
Fs05 = f_max * 0.5;

% 抽样语句
x4 = x(1:round(Fs/Fs4):end);
x2 = x(1:round(Fs/Fs2):end);
x1 = x(1:round(Fs/Fs1):end);
x05 = x(1:round(Fs/Fs05):end);

% 时域波形
t4 = (0:length(x4)-1)/Fs4;
t2 = (0:length(x2)-1)/Fs2;
t1 = (0:length(x1)-1)/Fs1;
t05 = (0:length(x05)-1)/Fs05;

% 播放音频
sound(x4, Fs4);
pause(length(x4)/Fs4);
sound(x2, Fs2);
pause(length(x2)/Fs2);
sound(x1, Fs1);
pause(length(x1)/Fs1);
sound(x05, Fs05);
pause(length(x05)/Fs05);

% 画图
subplot(4,1,1);
plot(t4, x4);
xlabel('Time (s)');
ylabel('Amplitude');
title('4x Sampling');

subplot(4,1,2);
plot(t2, x2);
xlabel('Time (s)');
ylabel('Amplitude');
title('2x Sampling');

subplot(4,1,3);
plot(t1, x1);
xlabel('Time (s)');
ylabel('Amplitude');
title('1x Sampling');

subplot(4,1,4);
plot(t05, x05);
xlabel('Time (s)');
ylabel('Amplitude');
title('0.5x Sampling');

运行上述代码可以得到抽样后的音乐波形图，并且播放抽样后的音频。从波形图上可以看出，抽样频率越低，音乐的高频部分越模糊，失去了细节和清晰度。而抽样频率越高，音乐的高频部分越清晰，但是文件大小也越大。