matlab对于连续信号xa(t)=cos(2πf1t)+5cos(2πf2t)+cos(2πf3t) ，其中f1=8.5kHz, f2=9kHz, f3=13kHz, （1）以采样频率fs=32 kHz对其进行采样，对xa(t)信号采集16点样本，分别作16点和补零到256点的FFT，并分别绘出对应的幅频特性曲线。 （2）分别以采样频率fs1=16 kHz和fs2=32 kHz对其进行采样，对xa(t)信号采集256点样本，分别作FFT，并分别绘出对应的幅频特性曲线。

好的，我会根据您的要求回答问题。对于您的第一个问题，可以按照以下步骤进行操作：

1. 生成信号xa(t)

Fs = 32000; % 采样频率
t = 0:1/Fs:15/Fs; % 时间范围
f1 = 8500; % 信号频率1
f2 = 9000; % 信号频率2
f3 = 13000; % 信号频率3
xa = cos(2*pi*f1*t) + 5*cos(2*pi*f2*t) + cos(2*pi*f3*t); % 生成信号

2. 对xa(t)信号采集16点样本，进行16点FFT

N = 16; % FFT点数
X = fft(xa(1:N)); % 进行FFT计算
f = (0:N-1)*Fs/N; % 计算频率

3. 绘制16点FFT幅频特性曲线

plot(f, abs(X));
xlabel('Frequency (Hz)');
ylabel('Magnitude');
title('16-point FFT of xa(t)');

4. 对xa(t)信号采集16点样本，进行补零到256点的FFT

N = 256; % FFT点数
X = fft(xa(1:16), N); % 进行FFT计算
f = (0:N-1)*Fs/N; % 计算频率

5. 绘制补零到256点FFT幅频特性曲线

plot(f, abs(X));
xlabel('Frequency (Hz)');
ylabel('Magnitude');
title('256-point FFT of xa(t) with zero-padding');

对于您的第二个问题，可以按照以下步骤进行操作：

1. 生成信号xa(t)

Fs1 = 16000; % 采样频率1
Fs2 = 32000; % 采样频率2
t = 0:1/Fs2:255/Fs2; % 时间范围
f1 = 8500; % 信号频率1
f2 = 9000; % 信号频率2
f3 = 13000; % 信号频率3
xa = cos(2*pi*f1*t) + 5*cos(2*pi*f2*t) + cos(2*pi*f3*t); % 生成信号

2. 对xa(t)信号采集256点样本，进行16000Hz采样频率的FFT

N = 256; % FFT点数
X1 = fft(xa, N); % 进行FFT计算
f1 = (0:N-1)*Fs1/N; % 计算频率

3. 绘制16000Hz采样频率FFT幅频特性曲线

plot(f1, abs(X1));
xlabel('Frequency (Hz)');
ylabel('Magnitude');
title('FFT of xa(t) sampled at 16 kHz');

4. 对xa(t)信号采集256点样本，进行32000Hz采样频率的FFT

X2 = fft(xa, N); % 进行FFT计算
f2 = (0:N-1)*Fs2/N; % 计算频率

5. 绘制32000Hz采样频率FFT幅频特性曲线

plot(f2, abs(X2));
xlabel('Frequency (Hz)');
ylabel('Magnitude');
title('FFT of xa(t) sampled at 32 kHz');