时域离散系统的单位脉冲响应为：h(n)={1,2,3,4,5,6,6,5,4,3,2,1},画出其幅频特性曲线和相频特性曲线。使用MATLAB写出程序

时域离散系统的单位脉冲响应 \( h(n) = \{1, 2, 3, 4, 5, 6, 6, 5, 4, 3, 2, 1\} \) 描述了一个有限长序列，我们可以先将其转换成Z变换来计算其系统函数H(z)，然后利用Z变换的性质来得到幅频特性和相频特性。

幅频特性是指H(z)在单位圆上（即|z|=1）的模值，而相频特性则是H(z)在单位圆上对应的幅角导数。对于离散时间系统，幅频特性通常通过DTFT (Discrete-Time Fourier Transform) 计算，相频特性则通过求导计算。

不过，由于这里给出的是具体的数值序列，我们无法直接画出幅频特性曲线和相频特性曲线，因为这需要进行连续的频率分析，而实际应用中我们可能会使用MATLAB的内置函数如`freqz()`来进行这个过程。以下是基本的MATLAB代码示例：

% 定义离散系统的单位脉冲响应
h = [1 2 3 4 5 6 6 5 4 3 2 1];

% 计算系统函数H(z)
H = zp2tf(h); % 将序列转换为零极点形式

% 使用 MATLAB 的 freqz 函数计算幅频特性
[Hmag, w] = freqz(H, 1, 'whole'); % 'whole' 参数表示考虑所有频率

% 真实幅度谱
plot(w, abs(Hmag), 'b', 'LineWidth', 2);
xlabel('Frequency');
ylabel('Magnitude');
title('Magnitude Spectrum');

% 相位谱
phase_data = unwrap(angle(Hmag));
plot(w, phase_data, 'g', 'LineWidth', 2);
xlabel('Frequency');
ylabel('Phase (radians)');
title('Phase Spectrum');

% 如果你想要计算相频特性，通常需要对上面的相位数据求导：
% [dPhase, dW] = diff(phase_data, 1);
% plot(dW, dPhase, 'm', 'LineWidth', 2);
% title('Group Delay');

% 注意：以上代码假设你已经安装了 Control System Toolbox，并且知道如何处理零极点转换（zp2tf）。如果没有安装，需查找替代方法。

请注意，由于MATLAB中的`freqz()`函数默认是在0到π之间取样，如果你需要更全面的频域分析，可能需要调整参数或者采用其他工具。此外，实际应用中可能还需要对结果进行一些滤波和可视化优化。