在matlab中1、编写FFT谱分析程序（1） x1(n)=cos(n*pi/4) （2） x1(n)=sin(n*pi/8) 要求：N=64时的程序清单，并绘制出其幅频特性曲线。

在MATLAB中，你可以使用fft函数来计算离散傅立叶变换（Discrete Fourier Transform，DFT），然后利用plot函数来绘制幅度频率特性（Magnitude Frequency Response）。以下是满足您需求的MATLAB代码示例：

% 设置信号长度
N = 64;

% 创建两个信号
x1 = cos(pi * (0:N-1) / 4); % 正弦波，频率π/4
x2 = sin(pi * (0:N-1) / 8); % 余弦波，频率π/8

% 计算每个信号的FFT
X1 = fft(x1);
X2 = fft(x2);

% 获取信号的长度（包括DC分量）
X1_length = length(X1);
X2_length = length(X2);

% 对结果取绝对值并将其归一化到0到N-1之间
X1_mag = abs(X1(1:X1_length/2+1)) ./ N; % 假设偶数长度，去掉负半部分
X2_mag = abs(X2(1:X2_length/2+1)) ./ N;

% 绘制幅频特性曲线
figure;
subplot(2,1,1)
plot(freqz(X1), 'LineWidth', 2); % 使用freqz函数获取频率轴数据
title('Signal 1 (Cosine with π/4 Frequency)');
xlabel('Frequency Index');
ylabel('Magnitude');

subplot(2,1,2)
plot(freqz(X2), 'LineWidth', 2);
title('Signal 2 (Sine with π/8 Frequency)');
xlabel('Frequency Index');
ylabel('Magnitude');

在这个代码中，`freqz`函数用于同时返回频率和对应的复数响应，我们只关注幅值部分。最后，通过`subplot`创建两个子图分别展示两个信号的幅频特性。