MATLAB实现余弦函数的dtft

下面是MATLAB实现余弦函数（$x[n]=\cos(\omega_0 n)$）的DTFT的代码：

% 设置参数
N = 64; % DTFT计算点数
w = linspace(-pi, pi, N); % 角频率向量
w0 = pi/4; % 余弦函数的角频率

% 计算余弦函数
n = 0:N-1;
x = cos(w0*n);

% 计算DTFT
X = dtft(x, w);

% 绘制幅度响应图像
figure;
plot(w, abs(X));
xlabel('角频率 \omega');
ylabel('幅度');
title('余弦函数的DTFT幅度响应');

代码中，我们先设置了DTFT计算点数`N`和角频率向量`w`，然后计算出余弦函数的离散时间序列`x`。最后使用`dtft`函数计算`x`的DTFT，绘制出DTFT的幅度响应图像。

执行上述代码，会得到如下的图像：