用MATLAB完成以下内容：利用fft计算窗函数的频谱，N=128，通过上机编程绘制|W(k)|的曲线以及w（n）的波形）w(n)=1/2*[1-cos(2pi*n/(N-1))]*RN(n)，并标出旁瓣峰值衰减

以下是MATLAB代码实现：

N = 128;
n = 0:N-1;
win = 1/2*(1-cos(2*pi*n/(N-1))).*rectwin(N)'; % 定义窗函数
W = fft(win); % 计算窗函数的频谱
W_mag = abs(W); % 取幅值
W_mag_dB = 20*log10(W_mag); % 转换为分贝

% 绘制频谱
k = 0:N-1;
plot(k, W_mag_dB);
title('Window Function Spectrum');
xlabel('Frequency (k)');
ylabel('Magnitude (dB)');

% 绘制时域波形
figure;
stem(n, win);
title('Window Function Time Domain Waveform');
xlabel('Sample (n)');
ylabel('Amplitude');

% 计算旁瓣峰值衰减
[max_peak, max_peak_index] = max(W_mag_dB);
peak_indexes = [1:max_peak_index-1, max_peak_index+1:N];
peak_values = W_mag_dB(peak_indexes);
peak_value = max(peak_values);
side_lobe_level = peak_value - max_peak;
fprintf('Side Lobe Level: %.2f dB\n', side_lobe_level);

上述代码首先定义了窗函数，然后通过FFT计算窗函数的频谱，绘制了频谱曲线和时域波形。最后通过计算旁瓣峰值衰减来评估窗函数的性能。可以通过修改窗函数的定义和N的值来计算不同的窗函数频谱。