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

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

% 定义一个包含0到9的整数的向量n
n = (0:1:9);
% 计算y(n) = cos(0.48*pi*n) + cos(0.52*pi*n)
y = cos(0.48*pi*n) + cos(0.52*pi*n);

% 定义一个包含0到99的整数的向量n1
n1 = (0:1:99);
% 构造一个长度为100的向量x，前10个元素为y，后90个元素为0
x = [y(1:1:10),zeros(1,90)];

% 在一个3行1列的图形窗口中，绘制第1个图形
subplot(3,1,1);
% 用离散点图绘制向量x的值
stem(n1,x);
% 添加标题
title('x(n) (0<=n<=99)');
% 添加x轴标签
xlabel('n');
% 设置坐标轴范围
axis([0,100,-2.5,2.5]);

% 在频域中，构造一个与x等长的向量w
w = linspace(0,2*pi,length(x));
% 将x在频域中的值计算出来，并存储在向量xw中
xw = x*exp(-j*[1:length(x)]'*w);
% 计算x在频域中的幅度
magx = abs(xw);
% 在一个3行1列的图形窗口中，绘制第2个图形
subplot(3,1,2);
% 用连续线条绘制w与magx之间的关系
plot(w,magx);
% 添加标题
title('DTFT');
% 添加x轴标签
xlabel('w');
% 设置坐标轴范围
axis([0,2*pi,0,10]);

% 在一个3行1列的图形窗口中，绘制第3个图形
subplot(3,1,3);
% 计算x的离散傅里叶变换（DFT）
x1 = fft(x);
% 计算x1的幅度
magx1 = abs(x1);
% 用离散点图绘制magx1
stem(n1,abs(magx1));
% 设置坐标轴范围
axis([0,100,0,10]);
% 添加x轴标签
xlabel('frequency in pi units');
% 添加标题
title('DFT');

此代码演示了如何在MATLAB中进行离散时间傅里叶变换（DTFT）和离散傅里叶变换（DFT）的计算和可视化。具体来说，它生成一个长度为100的向量x，其中前10个元素是通过cos函数计算得到的，其余元素为0。然后，它计算x的DTFT和DFT，并将结果绘制在一个3行1列的图形窗口中。第一个图形是离散时间信号x(n)的离散点图，第二个图形是x在频域中的幅度谱（DTFT），第三个图形是x的离散傅里叶变换（DFT）的幅度谱。