例3：分析占空比为ta的矩形脉冲频谱



T=1 时间窗大小

ta=0.5 占空比

N=256

t=linspace(-T,T,N)

ts=2*T/N

fs=1/ts

-x=zeros(1,N)

x(N/2-round(N/2*ta):N/2+round(N/2*ta))=ones(1,2*round(N/2*ta)+1)

y=fft(x,N)



figure

plot(t,x)

title([ 时域信号图 占空比为 ,num2str(ta)])

axis([-T T-0.2 1.2])

xlabel( t /s )



figure

subplot(211)

f=linspace(-fs/2,fs/2,N)

plot(f,abs(fftshift(y))/N)

title([ 矩形脉冲双边谱　脉冲时间t=[- ,num2str(ta*T), ,　,num2str(ta*T), ] ])

xlabel( f /Hz )