求大神说说这个表达式怎么弄出来图呀

fs=1000; %采样频率
w=-20:1/fs:20; %定义时间样本向量
for n=-20:1:20
Fw=(w+1-n*2).*((w-n*2)>=-1 & (w-n*2)<=0)+(1-w+n*2).*((w-n*2)>0 & (w-n*2)<=1); %定义 F(w-2n)
Pw=2*pi/3*(sin(n*pi/3))/(n*pi/3); %定义P(w)
if n==0
Pw=2*pi/3;
end
y1=1/2/pi*Pw.*Fw; %定义卷积后的结果
plot(w,y1,'b');%绘制图像
axis([-25,25,-0.2,0.5]); %设置坐标轴范围
hold on;%保持图形
end
xlabel('ω'); %定义横坐标为ω
ylabel('Fs(ω)'); %定义纵坐标为F（ω）
grid on; %图形显示网格线
title('Fs(ω)频谱图形');