matlab Square函数不能正确产生50%方波问题

最近电赛培训需要对方波进行傅里叶变换求谐波频率及幅度，因为各种误差需要对单片机得到的结果再乘以一个系数来补偿，所以需要matlab来模拟生成方波看不同输入幅值时准确的谐波幅度，得到的结果却一直都不对

问题：频率为1kHz的方波，用8kHz采样，画出图像不是50%方波

错误做法：

clc;
Fs=8000;%采样频率
N=1024;   %采样点数
n=0:N-1;
t=0:N-1; 
k=1;%调整方波幅度
x=0.5*k+0.5*k*square(2*pi*1000*t/Fs,50); %注意不能是1000
subplot(2,1,1)
stem(t,x);
axis([0 50 -inf inf]);
xfft=fft(x,N);
subplot(2,1,2)
f=n*Fs/N;   %转换位频率
Amp=abs(xfft)*2/N;
stem(f,Amp);

分析：

 可以看到上图的时域和频域图像都正确，没有生成50%占空比的方波。如图所示，虚线部分就是不确定的状态，可能为上面也可能为下面，所以就导致了上述结果

解决：信号频率改为1000.1即可

clc;
Fs=8000;%采样频率
N=1024;   %采样点数
n=0:N-1;
t=0:N-1; 
k=1;%调整方波幅度
x=0.5*k+0.5*k*square(2*pi*1000.1*t/Fs,50); %注意不能是1000
subplot(2,1,1)
stem(t,x);
axis([0 50 -inf inf]);
xfft=fft(x,N);
subplot(2,1,2)
f=n*Fs/N;   %转换位频率
Amp=abs(xfft)*2/N;
stem(f,Amp);