哪位matlab高手可以帮忙写个程序啊，生成一个矩形窗函数，然后傅立叶变换求它的频谱图。

N=21;
wp=0.3*pi;
ws=0.5*pi;
wc=ws-wp;
fc=wc/pi;
n=[0:N-1];
r=(N-1)/2;
m=n-r;
hd=fc*sinc(fc*m);%产生理想低通滤波器的脉冲响应

wn=hamming(N);%产生海明窗函数序列
hn=hd.*wn';%对理想低通滤波器加海明窗
[H,w]=freqz(hn,1,512,'whole',2);
db=20*log10(abs(H)/max(abs(H)));
subplot(221);
plot(w,db);
grid;
xlabel('归一化频率w');
ylabel('相对幅度/db');
title('(a) N=21时幅度响应特性');

这是我做的21阶海明窗函数，矩形窗只需要把hamming(N)函数换成triang(N)就行了，其他窗函数也是一样。
[H,w]=freqz(hn,1,512,'whole',2);这句话是求频谱。
我本科毕业设计就是做的基于matlab的FIR滤波器。很有多论文需要在知网付费下载。
有需要的话另外提问，多给些悬赏分，留下邮箱，我发给你。不过我的毕业论文不给哈，但是有很多很多的参考材料

你是刘斌吗？

明明是三角形窗