clc ;clf reset;close all;clear all;warning off MATLAB:divideByZero
m = 4;T = 1;
k = m - 10:m + 10;
xt = cos((m*T/4 + 1e-5 - k*T/4)*pi/T)./(1 - 4*(m*T/4 + 1e-5 - k*T/4).^2/T^2).*sin(pi*(m*T/4 + 1e-5 - k*T/4)/T)./(pi*(m*T/4 + 1e-5 - k*T/4)/T)
xt = cos((m*T/2 + 1e-5 - k*T/2)*pi/T)./(1 - 4*(m*T/2 + 1e-5 - k*T/2).^2/T^2).*sin(pi*(m*T/2 + 1e-5 - k*T/2)/T)./(pi*(m*T/2 + 1e-5 - k*T/2)/T);
