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

这里提供一个MATLAB代码示例，实现对音乐信号的抽样和分析：

% 读取音乐文件
[x, fs] = audioread('music.wav'); % x为音频数据，fs为采样率
dur = length(x) / fs; % 音频时长

% 分析原信号的频谱
N = length(x);
X = ctft1(x, N, fs);
f = linspace(0, fs, N+1);
f = f(1:end-1);
figure;
plot(f, abs(X));
xlabel('Frequency (Hz)');
ylabel('Magnitude');
title('Original Signal Spectrum');

% 设置抽样倍数
M = [4, 2, 1, 0.5];

% 对原信号进行抽样并绘制波形图
figure;
for k = 1:length(M)
    % 抽样
    xD = x(1:M(k):end);
    fsD = fs / M(k); % 抽样后的采样率
    durD = length(xD) / fsD; % 抽样后的音频时长
    % 绘制波形图
    subplot(length(M), 1, k);
    tD = linspace(0, durD, length(xD));
    plot(tD, xD);
    xlabel('Time (s)');
    ylabel('Amplitude');
    title(sprintf('Sampled Signal (M = %d)', M(k)));
    % 播放抽样后的音乐
    sound(xD, fsD);
    pause(durD); % 等待音乐播放完毕
end

在运行该代码前，请将当前目录下的名为`music.wav`的音乐文件替换为你自己的音乐文件，并确保MATLAB可以正常读取该文件。

运行该代码后，将会分析原信号的频谱并绘制幅度谱图，然后对原信号进行四种不同倍数的抽样，并分别绘制波形图和播放抽样后的音乐。你可以通过观察波形图和听取音乐，来感受抽样倍数对声音效果的影响。