请为这段代码作详细注释：close all clear clc n=(0:1:127); x=cos(0.48*pi*n)+cos(0.52*pi*n); subplot(3,1,1); stem(n,x); title('x(n) (0<=n<=127)'); xlabel('n'); axis([0,127,-2.5,2.5]); w = linspace(0,2*pi,length(n)); xw = x*exp(-j*[1:length(n)]'*w); magx=abs(xw); subplot(3,1,2); plot(w,magx); title('DTFT');xlabel('w'); axis([0,2*pi,0,70]); subplot(3,1,3); x1=fft(x); magx1=abs(x1); stem(n,abs(magx1)); axis([0,128,0,10]); title('DFT'); xlabel('frequency in pi units');

这段代码的主要作用是生成一个离散时间信号，并对其进行离散时间傅里叶变换（DTFT）和离散傅里叶变换（DFT），并进行可视化展示。

具体注释如下：

close all     % 关闭所有已打开的图形窗口
clear         % 清除工作区变量
clc           % 清除命令行窗口

n=(0:1:127);  % 定义一个包含128个元素的离散时间序列，n从0到127
x=cos(0.48*pi*n)+cos(0.52*pi*n); % 生成一个由两个余弦函数组成的离散时间信号

subplot(3,1,1); % 将画布分成3行1列，并将当前位置设置为第1个子图
stem(n,x); % 绘制离散时间信号x(n)的幅度图，使用离散点表示
title('x(n) (0<=n<=127)'); % 设置子图标题
xlabel('n');   % 设置x轴标签
axis([0,127,-2.5,2.5]); % 设置坐标轴范围

w = linspace(0,2*pi,length(n)); % 生成一个从0到2π的角频率序列，长度为n的长度
xw = x*exp(-j*[1:length(n)]'*w); % 计算离散时间信号x(n)的DTFT，使用复指数表示

magx=abs(xw); % 计算DTFT的幅度
subplot(3,1,2); % 将画布分成3行1列，并将当前位置设置为第2个子图
plot(w,magx); % 绘制DTFT的幅度谱，使用连续曲线表示
title('DTFT');xlabel('w'); % 设置子图标题和x轴标签
axis([0,2*pi,0,70]); % 设置坐标轴范围

subplot(3,1,3); % 将画布分成3行1列，并将当前位置设置为第3个子图
x1=fft(x); % 计算离散时间信号x(n)的DFT
magx1=abs(x1); % 计算DFT的幅度
stem(n,abs(magx1)); % 绘制DFT的幅度谱，使用离散点表示
axis([0,128,0,10]); % 设置坐标轴范围
title('DFT'); xlabel('frequency in pi units'); % 设置子图标题和x轴标签

注：DTFT和DFT都是傅里叶变换的变体，用于将离散时间信号转换为连续时间信号或有限长度序列的频域表示。其中DTFT是将离散时间信号视为周期信号来计算的，而DFT是将离散时间信号视为有限长度序列来计算的。