t=linspace(-100,100)
t0=2a=1
T=(t-t0)/a
y = sinc(t)
plot(t,y)
xlabel('Time (sec)')ylabel('Amplitude')title('Sinc Function')
运行结果
可以这样写:clear
fs=10000
t0=0.1
t=0:0.0001:t0
m=sinc(200*t)
mk=fft(m,50000) %通过fft计算频谱
mw=2*pi/1000*abs(fftshift(mk))%频谱搬移
fw=[-25000:24999]/50000*fs
plot(fw,mw)grid
xlim([-500,500])
其中,fft的抽样点数等参数是可以调的。
希望帮到你,欢迎追问~
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)