Matlab代码编写一段梳妆频带激励，该频带由25个成分组成，从100Hz到1299Hz，频率间隔为48Hz,激励幅值为0.05，并进行整周期截断进行频谱分析绘制频谱图

以下是实现该要求的 Matlab 代码：

% 设置参数
fs = 8000; % 采样率为8kHz
f_start = 100; % 频率起点为100Hz
f_end = 1299; % 频率终点为1299Hz
f_step = 48; % 频率间隔为48Hz
n_components = 25; % 25个成分
amplitude = 0.05; % 激励幅值为0.05

% 生成梳妆频带信号
t = 0:1/fs:1-1/fs; % 生成1秒的时间序列
signal = zeros(size(t)); % 初始化信号
for n = 1:n_components
    f = f_start + (n-1)*f_step; % 计算频率
    signal = signal + amplitude*cos(2*pi*f*t); % 加入成分
end

% 整周期截断信号
N = length(signal); % 信号长度
signal_truncated = signal(1:N-mod(N,n_components));

% 绘制频谱图
N_fft = 2048; % FFT长度为2048
freq = (0:N_fft/2-1)/N_fft*fs; % 频率轴
spectrum = abs(fft(signal_truncated,N_fft)).^2/N_fft/fs; % 计算功率谱密度
spectrum = spectrum(1:N_fft/2); % 取一半频谱
plot(freq,spectrum);
xlabel('Frequency (Hz)');
ylabel('Power Spectral Density');
title('Spectrum of Comb Filter Signal');

运行上述代码，将会得到一个名为“Spectrum of Comb Filter Signal”的频谱图，显示了梳妆频带信号的功率谱密度。