用MATLAB实现正弦信号y=cos（1/4pit)波形及频谱，画图、对每一行代码进行标注。写出完整代码，注意给出标注

在MATLAB中，我们可以轻松地生成一个余弦信号`y = cos(1/4 * pi * t)`的波形，并利用`fft`函数计算其频率谱。以下是完整的代码示例，包括每一步的注释：

% 初始化并设定信号的采样率和时间范围
Fs = 1000; % 采样频率
t = 0:1/Fs:1-1/Fs; % 时间向量，从0到1秒

% 定义余弦函数
% y = cos(1/4 * pi * t)
y = cos(1/4 * pi * t);

% 绘制波形
% plot(t, y, 'r'); % 红色线绘制波形
% xlabel('Time (s)');
% ylabel('Amplitude');
% title('Cosine Signal with Frequency 1/4 pi Hz');
% grid on;

% 计算快速傅里叶变换 (FFT)
Y = fft(y); % 对信号进行离散傅立叶变换
Pyy = abs(Y/Fs); % 计算功率谱密度，除以采样频率以归一化

% 将频率轴从样本索引转换为实际频率 (Hz)
f = Fs*(0:(length(Y)-1))/length(Y); % 频率向量

% 绘制频谱图
plot(f, Pyy, 'b'); % 蓝色线绘制频谱
hold on; % 保持当前图形以便叠加
xlabel('Frequency (Hz)');
ylabel('Power Spectrum Density');
title('Spectrum of Cosine Signal with Peak at 1/4 pi Hz');
grid on;
legend('Signal', 'Spectrum');

% 断开坐标轴网格
hold off;

% 显示图像
figure;

这里我们首先生成了一个从0到1秒的均匀时间序列`t`，然后计算了余弦函数`y`。接着，通过`plot`函数绘制出波形。为了得到频谱，我们使用`fft`函数将信号转换到频率域，并显示频谱图。