MATLAB绘制序列g1(n)=cos(0.6πn)的DTFT幅频特性曲线

使用MATLAB绘制序列g1(n)=cos(0.6πn)的DTFT幅频特性曲线

在数字信号处理中，离散时间傅里叶变换 (DTFT) 是分析离散时间信号频率特性的重要工具。本文将介绍如何使用 MATLAB 绘制序列 g1(n) = cos(0.6πn) 的 DTFT 幅频特性曲线，并提供详细的代码和解释。

MATLAB 代码：

% 定义序列 g1(n) = cos(0.6*pi*n)
n = -100:100;
g1 = cos(0.6*pi*n);

% 计算序列的DTFT
w = -pi:0.01:pi;
G1 = freqz(g1, 1, w);

% 绘制幅频特性曲线
plot(w/pi, abs(G1))
xlabel('归一化频率 (×π rad/sample)')
ylabel('幅度')
title('序列 g1(n) 的幅频响应')

% 设置坐标轴范围
xlim([-1 1])

代码解释：

定义序列：
- n = -100:100; 定义了序列的时间范围为 -100 到 100。
- g1 = cos(0.6*pi*n); 定义了序列 g1(n) = cos(0.6πn)。
计算 DTFT：
- w = -pi:0.01:pi; 定义了频率范围为 -π 到 π，步长为 0.01π。
- G1 = freqz(g1, 1, w); 使用 freqz 函数计算序列 g1 的 DTFT，返回结果存储在变量 G1 中。
绘制幅频特性曲线：
- plot(w/pi, abs(G1)) 使用 plot 函数绘制归一化频率 (w/pi) 与 DTFT 的幅度 (abs(G1)) 的关系曲线。
- xlabel('归一化频率 (×π rad/sample)') 设置 x 轴标签为 '归一化频率 (×π rad/sample)'。
- ylabel('幅度') 设置 y 轴标签为 '幅度'。
- title('序列 g1(n) 的幅频响应') 设置图像标题为 '序列 g1(n) 的幅频响应'。
设置坐标轴范围：
- xlim([-1 1]) 设置 x 轴范围为 -1 到 1，以便更好地观察幅频特性曲线的细节。

运行以上代码，即可得到序列 g1(n) 的 DTFT 幅频特性曲线。该曲线将展示该序列在不同频率上的幅度响应。