matlab窗函数截断及其频谱

clc
clear all;
f0=20;                  %信号频率为20Hz
Ts=0.005;               %采样周期0.005s
%截断长度为0.8s时的结果
T1=0.8;                 %截取长度为0.8s
N1=T1/Ts;               %采样点数
n1=0:N1-1;
t1=n1*Ts;
ft1=cos（2*pi*f0*t1）;
subplot（121）;plot（t1ft1）grid on;
title（‘a.截取长度为0.8s时的截取波形‘）xlabel（‘时间t/s‘）ylabel（‘幅值/ft1‘）axis（[00.9-1.11.1]）;
%截断长度为0.81s时的结果
T2=0.81;                %截取长度为0.81s
N2=T2/Ts;               %采样点数
n2=0:N2-1;
t2=n2*Ts;
ft2=cos（2*pi*f0*t2）;
subplot（122）;plot（t2ft2）grid on;
title（‘b.截取长度为0.81s时的截取波形‘）xlabel（‘时间t/s‘）ylabel（‘幅值/ft2‘）axis（[00.9-1.11.1]）;
%矩形窗函数截断
%截取长度为0.8s时的截取波形
window1_1=boxcar（N1）;   
f_boxcar1=window1_1.*（ft1）‘;
figure;subplot（221）; plot（t1f_boxcar1）grid on;
title（‘a.截取长度为0.8s时的矩形窗截断波形‘）xlabel（‘时间t/s‘）ylabel（‘幅值/ft1‘）axis（[