例4：在sinc(t)信号中叠加噪声，分析对应频谱



fs=100 采样频率,必须大于两倍基带信号最高频率

ts=1/fs 采样时间间隔

T=2 时间窗大小

-

x=x+0.1*randn(1,N)

x=awgn(x,-10, measured )



figure

plot(t,x)

title( 时域信号图 )

xlabel( t /s )



y=fft(x,N)



figure

if mod(N,2)~=0

　　N=N-1

end

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

plot(f,abs(y(1:N/2))*2/N)

title( 叠加噪声的sinc(t)信号频谱 )

xlabel( f /Hz )



figure

y0=y(1:N/2)

P=y0.*conj(y0)/N

P=10*log10(P)

plot(f,P)

title( 功率谱 )

xlabel( f /Hz )